1
|
Tralamazza SM, Gluck-Thaler E, Feurtey A, Croll D. Copy number variation introduced by a massive mobile element facilitates global thermal adaptation in a fungal wheat pathogen. Nat Commun 2024; 15:5728. [PMID: 38977688 PMCID: PMC11231334 DOI: 10.1038/s41467-024-49913-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/17/2023] [Accepted: 06/25/2024] [Indexed: 07/10/2024] Open
Abstract
Copy number variation (CNV) can drive rapid evolution in changing environments. In microbial pathogens, such adaptation is a key factor underpinning epidemics and colonization of new niches. However, the genomic determinants of such adaptation remain poorly understood. Here, we systematically investigate CNVs in a large genome sequencing dataset spanning a worldwide collection of 1104 genomes from the major wheat pathogen Zymoseptoria tritici. We found overall strong purifying selection acting on most CNVs. Genomic defense mechanisms likely accelerated gene loss over episodes of continental colonization. Local adaptation along climatic gradients was likely facilitated by CNVs affecting secondary metabolite production and gene loss in general. One of the strongest loci for climatic adaptation is a highly conserved gene of the NAD-dependent Sirtuin family. The Sirtuin CNV locus localizes to an ~68-kb Starship mobile element unique to the species carrying genes highly expressed during plant infection. The element has likely lost the ability to transpose, demonstrating how the ongoing domestication of cargo-carrying selfish elements can contribute to selectable variation within populations. Our work highlights how standing variation in gene copy numbers at the global scale can be a major factor driving climatic and metabolic adaptation in microbial species.
Collapse
Affiliation(s)
- Sabina Moser Tralamazza
- Laboratory of Evolutionary Genetics, Institute of Biology, University of Neuchâtel, CH-2000, Neuchâtel, Switzerland
| | - Emile Gluck-Thaler
- Laboratory of Evolutionary Genetics, Institute of Biology, University of Neuchâtel, CH-2000, Neuchâtel, Switzerland
- Department of Plant Pathology, University of Wisconsin-Madison, Madison, WI, USA
| | - Alice Feurtey
- Laboratory of Evolutionary Genetics, Institute of Biology, University of Neuchâtel, CH-2000, Neuchâtel, Switzerland
- Plant Pathology, D-USYS, ETH Zurich, CH-8092, Zurich, Switzerland
| | - Daniel Croll
- Laboratory of Evolutionary Genetics, Institute of Biology, University of Neuchâtel, CH-2000, Neuchâtel, Switzerland.
| |
Collapse
|
2
|
Melton AE, Novak SJ, Buerki S. Utilizing a comparative approach to assess genome evolution during diploidization in Artemisia tridentata, a keystone species of western North America. AMERICAN JOURNAL OF BOTANY 2024; 111:e16353. [PMID: 38826031 DOI: 10.1002/ajb2.16353] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/07/2023] [Revised: 04/03/2024] [Accepted: 04/03/2024] [Indexed: 06/04/2024]
Abstract
PREMISE Polyploidization is often followed by diploidization. Diploidization is generally studied using synthetic polyploid lines and/or crop plants, but rarely using extant diploids or nonmodel plants such as Artemisia tridentata. This threatened western North American keystone species has a large genome compared to congeneric Artemisia species; dominated by diploid and tetraploid cytotypes, with multiple origins of tetraploids with genome size reduction. METHODS The genome of an A. tridentata sample was resequenced to study genome evolution and compared to that of A. annua, a diploid congener. Three diploid genomes of A. tridentata were compared to test for multiple diploidization events. RESULTS The A. tridentata genome had many chromosomal rearrangements relative to that of A. annua, while large-scale synteny of A. tridentata chromosome 3 and A. annua chromosome 4 was conserved. The three A. tridentata genomes had similar sizes (4.19-4.2 Gbp), heterozygosity (2.24-2.25%), and sequence (98.73-99.15% similarity) across scaffolds, and in k-mer analyses, similar patterns of diploid heterozygous k-mers (AB = 41%, 47%, and 47%), triploid heterozygous k-mers (AAB = 18-21%), and tetraploid k-mers (AABB = 13-17%). Biallelic SNPs were evenly distributed across scaffolds for all individuals. Comparisons of transposable element (TE) content revealed differential enrichment of TE clades. CONCLUSIONS Our findings suggest population-level TE differentiation after a shared polyploidization-to-diploidization event(s) and exemplify the complex processes of genome evolution. This research approached provides new resources for exploration of abiotic stress response, especially the roles of TEs in response pathways.
Collapse
Affiliation(s)
- Anthony E Melton
- Department of Biological Sciences, Boise State University, Boise, 83725, ID, USA
| | - Stephen J Novak
- Department of Biological Sciences, Boise State University, Boise, 83725, ID, USA
| | - Sven Buerki
- Department of Biological Sciences, Boise State University, Boise, 83725, ID, USA
| |
Collapse
|
3
|
Sjodin BMF, Schmidt DA, Galbreath KE, Russello MA. Putative climate adaptation in American pikas (Ochotona princeps) is associated with copy number variation across environmental gradients. Sci Rep 2024; 14:8568. [PMID: 38609461 PMCID: PMC11014952 DOI: 10.1038/s41598-024-59157-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/17/2023] [Accepted: 04/08/2024] [Indexed: 04/14/2024] Open
Abstract
Improved understanding of the genetic basis of adaptation to climate change is necessary for maintaining global biodiversity moving forward. Studies to date have largely focused on sequence variation, yet there is growing evidence that suggests that changes in genome structure may be an even more significant source of adaptive potential. The American pika (Ochotona princeps) is an alpine specialist that shows some evidence of adaptation to climate along elevational gradients, but previous work has been limited to single nucleotide polymorphism based analyses within a fraction of the species range. Here, we investigated the role of copy number variation underlying patterns of local adaptation in the American pika using genome-wide data previously collected across the entire species range. We identified 37-193 putative copy number variants (CNVs) associated with environmental variation (temperature, precipitation, solar radiation) within each of the six major American pika lineages, with patterns of divergence largely following elevational and latitudinal gradients. Genes associated (n = 158) with independent annotations across lineages, variables, and/or CNVs had functions related to mitochondrial structure/function, immune response, hypoxia, olfaction, and DNA repair. Some of these genes have been previously linked to putative high elevation and/or climate adaptation in other species, suggesting they may serve as important targets in future studies.
Collapse
Affiliation(s)
- Bryson M F Sjodin
- Department of Biology, The University of British Columbia, 3247 University Way, Kelowna, BC, V1V 1V7, Canada
| | - Danielle A Schmidt
- Department of Biology, The University of British Columbia, 3247 University Way, Kelowna, BC, V1V 1V7, Canada
| | - Kurt E Galbreath
- Department of Biology, Northern Michigan University, 1401 Presque Isle Ave, Marquette, MI, 49855, USA
| | - Michael A Russello
- Department of Biology, The University of British Columbia, 3247 University Way, Kelowna, BC, V1V 1V7, Canada.
| |
Collapse
|
4
|
Lecomte L, Árnyasi M, Ferchaud A, Kent M, Lien S, Stenløkk K, Sylvestre F, Bernatchez L, Mérot C. Investigating structural variant, indel and single nucleotide polymorphism differentiation between locally adapted Atlantic salmon populations. Evol Appl 2024; 17:e13653. [PMID: 38495945 PMCID: PMC10940791 DOI: 10.1111/eva.13653] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/23/2023] [Revised: 12/14/2023] [Accepted: 01/13/2024] [Indexed: 03/19/2024] Open
Abstract
Genomic structural variants (SVs) are now recognized as an integral component of intraspecific polymorphism and are known to contribute to evolutionary processes in various organisms. However, they are inherently difficult to detect and genotype from readily available short-read sequencing data, and therefore remain poorly documented in wild populations. Salmonid species displaying strong interpopulation variability in both life history traits and habitat characteristics, such as Atlantic salmon (Salmo salar), offer a prime context for studying adaptive polymorphism, but the contribution of SVs to fine-scale local adaptation has yet to be explored. Here, we performed a comparative analysis of SVs, single nucleotide polymorphisms (SNPs) and small indels (<50 bp) segregating in the Romaine and Puyjalon salmon, two putatively locally adapted populations inhabiting neighboring rivers (Québec, Canada) and showing pronounced variation in life history traits, namely growth, fecundity, and age at maturity and smoltification. We first catalogued polymorphism using a hybrid SV characterization approach pairing both short- (16X) and long-read sequencing (20X) for variant discovery with graph-based genotyping of SVs across 60 salmon genomes, along with characterization of SNPs and small indels from short reads. We thus identified 115,907 SVs, 8,777,832 SNPs and 1,089,321 short indels, with SVs covering 4.8 times more base pairs than SNPs. All three variant types revealed a highly congruent population structure and similar patterns of F ST and density variation along the genome. Finally, we performed outlier detection and redundancy analysis (RDA) to identify variants of interest in the putative local adaptation of Romaine and Puyjalon salmon. Genes located near these variants were enriched for biological processes related to nervous system function, suggesting that observed variation in traits such as age at smoltification could arise from differences in neural development. This study therefore demonstrates the feasibility of large-scale SV characterization and highlights its relevance for salmonid population genomics.
Collapse
Affiliation(s)
- Laurie Lecomte
- Institut de Biologie Intégrative et des Systèmes (IBIS)Université LavalQuébecCanada
- Département de BiologieUniversité LavalQuébecCanada
| | - Mariann Árnyasi
- Department of Animal and Aquacultural Sciences (IHA), Faculty of Life Sciences (BIOVIT), Centre for Integrative Genetics (CIGENE)Norwegian University of Life Sciences (NMBU)ÅsNorway
| | - Anne‐Laure Ferchaud
- Institut de Biologie Intégrative et des Systèmes (IBIS)Université LavalQuébecCanada
- Département de BiologieUniversité LavalQuébecCanada
- Present address:
Parks Canada, Office of the Chief Ecosystem ScientistQuébecQCCanada
| | - Matthew Kent
- Department of Animal and Aquacultural Sciences (IHA), Faculty of Life Sciences (BIOVIT), Centre for Integrative Genetics (CIGENE)Norwegian University of Life Sciences (NMBU)ÅsNorway
| | - Sigbjørn Lien
- Department of Animal and Aquacultural Sciences (IHA), Faculty of Life Sciences (BIOVIT), Centre for Integrative Genetics (CIGENE)Norwegian University of Life Sciences (NMBU)ÅsNorway
| | - Kristina Stenløkk
- Department of Animal and Aquacultural Sciences (IHA), Faculty of Life Sciences (BIOVIT), Centre for Integrative Genetics (CIGENE)Norwegian University of Life Sciences (NMBU)ÅsNorway
| | - Florent Sylvestre
- Institut de Biologie Intégrative et des Systèmes (IBIS)Université LavalQuébecCanada
- Département de BiologieUniversité LavalQuébecCanada
| | - Louis Bernatchez
- Institut de Biologie Intégrative et des Systèmes (IBIS)Université LavalQuébecCanada
- Département de BiologieUniversité LavalQuébecCanada
| | - Claire Mérot
- Institut de Biologie Intégrative et des Systèmes (IBIS)Université LavalQuébecCanada
- Département de BiologieUniversité LavalQuébecCanada
- Present address:
UMR 6553 Ecobio, OSUR, CNRSUniversité de RennesRennesFrance
| |
Collapse
|
5
|
Layton KKS, Brieuc MSO, Castilho R, Diaz-Arce N, Estévez-Barcia D, Fonseca VG, Fuentes-Pardo AP, Jeffery NW, Jiménez-Mena B, Junge C, Kaufmann J, Leinonen T, Maes SM, McGinnity P, Reed TE, Reisser CMO, Silva G, Vasemägi A, Bradbury IR. Predicting the future of our oceans-Evaluating genomic forecasting approaches in marine species. GLOBAL CHANGE BIOLOGY 2024; 30:e17236. [PMID: 38519845 DOI: 10.1111/gcb.17236] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/18/2023] [Revised: 02/21/2024] [Accepted: 02/25/2024] [Indexed: 03/25/2024]
Abstract
Climate change is restructuring biodiversity on multiple scales and there is a pressing need to understand the downstream ecological and genomic consequences of this change. Recent advancements in the field of eco-evolutionary genomics have sought to include evolutionary processes in forecasting species' responses to climate change (e.g., genomic offset), but to date, much of this work has focused on terrestrial species. Coastal and offshore species, and the fisheries they support, may be even more vulnerable to climate change than their terrestrial counterparts, warranting a critical appraisal of these approaches in marine systems. First, we synthesize knowledge about the genomic basis of adaptation in marine species, and then we discuss the few examples where genomic forecasting has been applied in marine systems. Next, we identify the key challenges in validating genomic offset estimates in marine species, and we advocate for the inclusion of historical sampling data and hindcasting in the validation phase. Lastly, we describe a workflow to guide marine managers in incorporating these predictions into the decision-making process.
Collapse
Affiliation(s)
- K K S Layton
- Department of Biology, University of Toronto Mississauga, Mississauga, Ontario, Canada
- School of Biological Sciences, University of Aberdeen, Aberdeen, UK
| | | | - R Castilho
- University of the Algarve, Faro, Portugal
- Centre for Marine Sciences, University of the Algarve, Faro, Portugal
- Pattern Institute, Faro, Portugal
| | - N Diaz-Arce
- AZTI Marine Research, Basque Research and Technology Alliance (BRTA), Sukarrieta, Spain
| | - D Estévez-Barcia
- Department of Fish and Shellfish, Greenland Institute of Natural Resources, Nuuk, Greenland
| | - V G Fonseca
- Centre for Environment, Fisheries and Aquaculture Science, Weymouth, UK
| | - A P Fuentes-Pardo
- Department of Immunology, Genetics and Pathology, SciLifeLab Data Centre, Uppsala University, Uppsala, Sweden
| | - N W Jeffery
- Bedford Institute of Oceanography, Fisheries and Oceans Canada, Dartmouth, Nova Scotia, Canada
| | - B Jiménez-Mena
- Section for Marine Living Resources, National Institute of Aquatic Resources, Technical University of Denmark, Silkeborg, Denmark
| | - C Junge
- Institute of Marine Research, Tromso, Norway
| | | | - T Leinonen
- Natural Resources Institute Finland, Helsinki, Finland
| | - S M Maes
- Flanders Research Institute for Agriculture, Fisheries and Food, Ostend, Belgium
| | - P McGinnity
- School of Biological, Earth & Environmental Sciences, University College Cork, Cork, Ireland
| | - T E Reed
- School of Biological, Earth & Environmental Sciences, University College Cork, Cork, Ireland
| | - C M O Reisser
- MARBEC, University of Montpellier, CNRS, IFREMER, IRD, Montpellier, France
| | - G Silva
- MARE-Marine and Environmental Sciences Centre/ARNET-Aquatic Research Network, ISPA-Instituto Universitário, Lisbon, Portugal
| | - A Vasemägi
- Swedish University of Agricultural Sciences, Drottningholm, Sweden
- Estonian University of Life Sciences, Tartu, Estonia
| | - I R Bradbury
- Northwest Atlantic Fisheries Centre, Fisheries and Oceans Canada, St. John's, Newfoundland and Labrador, Canada
| |
Collapse
|
6
|
Smeds L, Huson LSA, Ellegren H. Structural genomic variation in the inbred Scandinavian wolf population contributes to the realized genetic load but is positively affected by immigration. Evol Appl 2024; 17:e13652. [PMID: 38333557 PMCID: PMC10848878 DOI: 10.1111/eva.13652] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/09/2023] [Revised: 01/08/2024] [Accepted: 01/16/2024] [Indexed: 02/10/2024] Open
Abstract
When populations decrease in size and may become isolated, genomic erosion by loss of diversity from genetic drift and accumulation of deleterious mutations is likely an inevitable consequence. In such cases, immigration (genetic rescue) is necessary to restore levels of genetic diversity and counteract inbreeding depression. Recent work in conservation genomics has studied these processes focusing on the genetic diversity of single nucleotide polymorphisms. In contrast, our knowledge about structural genomic variation (insertions, deletions, duplications and inversions) in endangered species is limited. We analysed whole-genome, short-read sequences from 212 wolves from the inbred Scandinavian population and from neighbouring populations in Finland and Russia, and detected >35,000 structural variants (SVs) after stringent quality and genotype frequency filtering; >26,000 high-confidence variants remained after manual curation. The majority of variants were shorter than 1 kb, with a distinct peak in the length distribution of deletions at 190 bp, corresponding to insertion events of SINE/tRNA-Lys elements. The site frequency spectrum of SVs in protein-coding regions was significantly shifted towards rare alleles compared to putatively neutral variants, consistent with purifying selection. The realized genetic load of SVs in protein-coding regions increased with inbreeding levels in the Scandinavian population, but immigration provided a genetic rescue effect by lowering the load and reintroducing ancestral alleles at loci fixed for derived SVs. Our study shows that structural variation comprises a common type of in part deleterious mutations in endangered species and that establishing gene flow is necessary to mitigate the negative consequences of loss of diversity.
Collapse
Affiliation(s)
- Linnéa Smeds
- Department of Ecology and Genetics, Evolutionary BiologyUppsala UniversityUppsalaSweden
| | - Lars S. A. Huson
- Department of Ecology and Genetics, Evolutionary BiologyUppsala UniversityUppsalaSweden
| | - Hans Ellegren
- Department of Ecology and Genetics, Evolutionary BiologyUppsala UniversityUppsalaSweden
| |
Collapse
|
7
|
Li Y, Yao J, Sang H, Wang Q, Su L, Zhao X, Xia Z, Wang F, Wang K, Lou D, Wang G, Waterhouse RM, Wang H, Luo S, Sun C. Pan-genome analysis highlights the role of structural variation in the evolution and environmental adaptation of Asian honeybees. Mol Ecol Resour 2024; 24:e13905. [PMID: 37996991 DOI: 10.1111/1755-0998.13905] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/07/2023] [Revised: 10/20/2023] [Accepted: 11/09/2023] [Indexed: 11/25/2023]
Abstract
The Asian honeybee, Apis cerana, is an ecologically and economically important pollinator. Mapping its genetic variation is key to understanding population-level health, histories and potential capacities to respond to environmental changes. However, most efforts to date were focused on single nucleotide polymorphisms (SNPs) based on a single reference genome, thereby ignoring larger scale genomic variation. We employed long-read sequencing technologies to generate a chromosome-scale reference genome for the ancestral group of A. cerana. Integrating this with 525 resequencing data sets, we constructed the first pan-genome of A. cerana, encompassing almost the entire gene content. We found that 31.32% of genes in the pan-genome were variably present across populations, providing a broad gene pool for environmental adaptation. We identified and characterized structural variations (SVs) and found that they were not closely linked with SNP distributions; however, the formation of SVs was closely associated with transposable elements. Furthermore, phylogenetic analysis using SVs revealed a novel A. cerana ecological group not recoverable from the SNP data. Performing environmental association analysis identified a total of 44 SVs likely to be associated with environmental adaptation. Verification and analysis of one of these, a 330 bp deletion in the Atpalpha gene, indicated that this SV may promote the cold adaptation of A. cerana by altering gene expression. Taken together, our study demonstrates the feasibility and utility of applying pan-genome approaches to map and explore genetic feature variations of honeybee populations, and in particular to examine the role of SVs in the evolution and environmental adaptation of A. cerana.
Collapse
Affiliation(s)
- Yancan Li
- State Key Laboratory of Resource Insects, Institute of Apicultural Research, Chinese Academy of Agricultural Sciences, Beijing, China
- Western Research Institute, Chinese Academy of Agricultural Sciences, Changji, China
| | - Jun Yao
- State Key Laboratory of Resource Insects, Institute of Apicultural Research, Chinese Academy of Agricultural Sciences, Beijing, China
| | - Huiling Sang
- State Key Laboratory of Resource Insects, Institute of Apicultural Research, Chinese Academy of Agricultural Sciences, Beijing, China
| | - Quangui Wang
- State Key Laboratory of Resource Insects, Institute of Apicultural Research, Chinese Academy of Agricultural Sciences, Beijing, China
| | - Long Su
- Institute of Plant Protection, Shandong Academy of Agricultural Sciences, Jinan, China
| | - Xiaomeng Zhao
- College of Animal Science, Shanxi Agricultural University, Shanxi, China
| | - Zhenyu Xia
- State Key Laboratory of Resource Insects, Institute of Apicultural Research, Chinese Academy of Agricultural Sciences, Beijing, China
| | - Feiran Wang
- State Key Laboratory of Resource Insects, Institute of Apicultural Research, Chinese Academy of Agricultural Sciences, Beijing, China
- Western Research Institute, Chinese Academy of Agricultural Sciences, Changji, China
| | - Kai Wang
- State Key Laboratory of Resource Insects, Institute of Apicultural Research, Chinese Academy of Agricultural Sciences, Beijing, China
| | - Delong Lou
- Shandong Provincial Animal Husbandry Station, Jinan, China
| | - Guizhi Wang
- Department of Animal Science, Shandong Agricultural University, Taian, China
| | - Robert M Waterhouse
- Department of Ecology and Evolution, University of Lausanne, and SIB Swiss Institute of Bioinformatics, Lausanne, Switzerland
| | - Huihua Wang
- State Key Laboratory of Resource Insects, Institute of Apicultural Research, Chinese Academy of Agricultural Sciences, Beijing, China
| | - Shudong Luo
- State Key Laboratory of Resource Insects, Institute of Apicultural Research, Chinese Academy of Agricultural Sciences, Beijing, China
- Western Research Institute, Chinese Academy of Agricultural Sciences, Changji, China
| | - Cheng Sun
- College of Life Sciences, Capital Normal University, Beijing, China
| |
Collapse
|
8
|
Dallaire X, Bouchard R, Hénault P, Ulmo-Diaz G, Normandeau E, Mérot C, Bernatchez L, Moore JS. Widespread Deviant Patterns of Heterozygosity in Whole-Genome Sequencing Due to Autopolyploidy, Repeated Elements, and Duplication. Genome Biol Evol 2023; 15:evad229. [PMID: 38085037 PMCID: PMC10752349 DOI: 10.1093/gbe/evad229] [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] [Accepted: 11/30/2023] [Indexed: 12/28/2023] Open
Abstract
Most population genomic tools rely on accurate single nucleotide polymorphism (SNP) calling and filtering to meet their underlying assumptions. However, genomic complexity, resulting from structural variants, paralogous sequences, and repetitive elements, presents significant challenges in assembling contiguous reference genomes. Consequently, short-read resequencing studies can encounter mismapping issues, leading to SNPs that deviate from Mendelian expected patterns of heterozygosity and allelic ratio. In this study, we employed the ngsParalog software to identify such deviant SNPs in whole-genome sequencing (WGS) data with low (1.5×) to intermediate (4.8×) coverage for four species: Arctic Char (Salvelinus alpinus), Lake Whitefish (Coregonus clupeaformis), Atlantic Salmon (Salmo salar), and the American Eel (Anguilla rostrata). The analyses revealed that deviant SNPs accounted for 22% to 62% of all SNPs in salmonid datasets and approximately 11% in the American Eel dataset. These deviant SNPs were particularly concentrated within repetitive elements and genomic regions that had recently undergone rediploidization in salmonids. Additionally, narrow peaks of elevated coverage were ubiquitous along all four reference genomes, encompassed most deviant SNPs, and could be partially associated with transposons and tandem repeats. Including these deviant SNPs in genomic analyses led to highly distorted site frequency spectra, underestimated pairwise FST values, and overestimated nucleotide diversity. Considering the widespread occurrence of deviant SNPs arising from a variety of sources, their important impact in estimating population parameters, and the availability of effective tools to identify them, we propose that excluding deviant SNPs from WGS datasets is required to improve genomic inferences for a wide range of taxa and sequencing depths.
Collapse
Affiliation(s)
- Xavier Dallaire
- Institut de biologie intégrative et des systèmes, Université Laval, Québec, Canada
- Centre d'Études Nordiques, Université Laval, Québec, Canada
| | - Raphael Bouchard
- Institut de biologie intégrative et des systèmes, Université Laval, Québec, Canada
- Ressources Aquatique Québec, Université de Rimouski, Rimouski, Canada
| | - Philippe Hénault
- Institut de biologie intégrative et des systèmes, Université Laval, Québec, Canada
- Ressources Aquatique Québec, Université de Rimouski, Rimouski, Canada
| | - Gabriela Ulmo-Diaz
- Institut de biologie intégrative et des systèmes, Université Laval, Québec, Canada
- Ressources Aquatique Québec, Université de Rimouski, Rimouski, Canada
| | - Eric Normandeau
- Institut de biologie intégrative et des systèmes, Université Laval, Québec, Canada
- Ressources Aquatique Québec, Université de Rimouski, Rimouski, Canada
- Plateforme de bio-informatique de l’IBIS, Université Laval, Québec, Canada
| | - Claire Mérot
- CNRS, UMR 6553 ECOBIO, Université de Rennes, Rennes, France
| | - Louis Bernatchez
- Institut de biologie intégrative et des systèmes, Université Laval, Québec, Canada
- Ressources Aquatique Québec, Université de Rimouski, Rimouski, Canada
| | - Jean-Sébastien Moore
- Institut de biologie intégrative et des systèmes, Université Laval, Québec, Canada
- Centre d'Études Nordiques, Université Laval, Québec, Canada
- Ressources Aquatique Québec, Université de Rimouski, Rimouski, Canada
| |
Collapse
|
9
|
Wold JR, Guhlin JG, Dearden PK, Santure AW, Steeves TE. The promise and challenges of characterizing genome-wide structural variants: A case study in a critically endangered parrot. Mol Ecol Resour 2023. [PMID: 36916824 DOI: 10.1111/1755-0998.13783] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/15/2022] [Revised: 02/24/2023] [Accepted: 03/09/2023] [Indexed: 03/15/2023]
Abstract
There is growing interest in the role of structural variants (SVs) as drivers of local adaptation and speciation. From a biodiversity genomics perspective, the characterization of genome-wide SVs provides an exciting opportunity to complement single nucleotide polymorphisms (SNPs). However, little is known about the impacts of SV discovery and genotyping strategies on the characterization of genome-wide SV diversity within and among populations. Here, we explore a near whole-species resequence data set, and long-read sequence data for a subset of highly represented individuals in the critically endangered kākāpō (Strigops habroptilus). We demonstrate that even when using a highly contiguous reference genome, different discovery and genotyping strategies can significantly impact the type, size and location of SVs characterized genome-wide. Further, we found that the mean number of SVs in each of two kākāpō lineages differed both within and across generations. These combined results suggest that genome-wide characterization of SVs remains challenging at the population-scale. We are optimistic that increased accessibility to long-read sequencing and advancements in bioinformatic approaches including multireference approaches like genome graphs will alleviate at least some of the challenges associated with resolving SV characteristics below the species level. In the meantime, we address caveats, highlight considerations, and provide recommendations for the characterization of genome-wide SVs in biodiversity genomic research.
Collapse
Affiliation(s)
- Jana R Wold
- University of Canterbury, Christchurch, New Zealand
| | - Joseph G Guhlin
- Genomics Aotearoa and Biochemistry Department, University of Otago, Dunedin, New Zealand
| | - Peter K Dearden
- Genomics Aotearoa and Biochemistry Department, University of Otago, Dunedin, New Zealand
| | | | | |
Collapse
|
10
|
Non-synonymous variation and protein structure of candidate genes associated with selection in farm and wild populations of turbot (Scophthalmus maximus). Sci Rep 2023; 13:3019. [PMID: 36810752 PMCID: PMC9944912 DOI: 10.1038/s41598-023-29826-z] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/02/2022] [Accepted: 02/10/2023] [Indexed: 02/24/2023] Open
Abstract
Non-synonymous variation (NSV) of protein coding genes represents raw material for selection to improve adaptation to the diverse environmental scenarios in wild and livestock populations. Many aquatic species face variations in temperature, salinity and biological factors throughout their distribution range that is reflected by the presence of allelic clines or local adaptation. The turbot (Scophthalmus maximus) is a flatfish of great commercial value with a flourishing aquaculture which has promoted the development of genomic resources. In this study, we developed the first atlas of NSVs in the turbot genome by resequencing 10 individuals from Northeast Atlantic Ocean. More than 50,000 NSVs where detected in the ~ 21,500 coding genes of the turbot genome, and we selected 18 NSVs to be genotyped using a single Mass ARRAY multiplex on 13 wild populations and three turbot farms. We detected signals of divergent selection on several genes related to growth, circadian rhythms, osmoregulation and oxygen binding in the different scenarios evaluated. Furthermore, we explored the impact of NSVs identified on the 3D structure and functional relationship of the correspondent proteins. In summary, our study provides a strategy to identify NSVs in species with consistently annotated and assembled genomes to ascertain their role in adaptation.
Collapse
|
11
|
Le Moan A, Panova M, De Jode A, Ortega-Martinez O, Duvetorp M, Faria R, Butlin R, Johannesson K. An allozyme polymorphism is associated with a large chromosomal inversion in the marine snail Littorina fabalis. Evol Appl 2023; 16:279-292. [PMID: 36793696 PMCID: PMC9923470 DOI: 10.1111/eva.13427] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/09/2022] [Revised: 05/04/2022] [Accepted: 05/26/2022] [Indexed: 11/26/2022] Open
Abstract
Understanding the genetic targets of natural selection is one of the most challenging goals of population genetics. Some of the earliest candidate genes were identified from associations between allozyme allele frequencies and environmental variation. One such example is the clinal polymorphism in the arginine kinase (Ak) gene in the marine snail Littorina fabalis. While other enzyme loci do not show differences in allozyme frequencies among populations, the Ak alleles are near differential fixation across repeated wave exposure gradients in Europe. Here, we use this case to illustrate how a new sequencing toolbox can be employed to characterize the genomic architecture associated with historical candidate genes. We found that the Ak alleles differ by nine nonsynonymous substitutions, which perfectly explain the different migration patterns of the allozymes during electrophoresis. Moreover, by exploring the genomic context of the Ak gene, we found that the three main Ak alleles are located on different arrangements of a putative chromosomal inversion that reaches near fixation at the opposing ends of two transects covering a wave exposure gradient. This shows Ak is part of a large (3/4 of the chromosome) genomic block of differentiation, in which Ak is unlikely to be the only target of divergent selection. Nevertheless, the nonsynonymous substitutions among Ak alleles and the complete association of one allele with one inversion arrangement suggest that the Ak gene is a strong candidate to contribute to the adaptive significance of the inversion.
Collapse
Affiliation(s)
- Alan Le Moan
- Tjärnö Marine Laboratory, Department of Marine Sciences University of Gothenburg Strömstad Sweden
| | - Marina Panova
- Tjärnö Marine Laboratory, Department of Marine Sciences University of Gothenburg Strömstad Sweden
| | - Aurélien De Jode
- Tjärnö Marine Laboratory, Department of Marine Sciences University of Gothenburg Strömstad Sweden
| | - Olga Ortega-Martinez
- Tjärnö Marine Laboratory, Department of Marine Sciences University of Gothenburg Strömstad Sweden
| | - Mårten Duvetorp
- Tjärnö Marine Laboratory, Department of Marine Sciences University of Gothenburg Strömstad Sweden
| | - Rui Faria
- InBIO Laboratório Associado, CIBIO, Centro de Investigação em Biodiversidade e Recursos Genéticos Universidade do Porto Vairão Portugal.,BIOPOLIS Program in Genomics, Biodiversity and Land Planning, CIBIO Campus de Vairão Vairão Portugal
| | - Roger Butlin
- Tjärnö Marine Laboratory, Department of Marine Sciences University of Gothenburg Strömstad Sweden.,Ecology and Evolutionary Biology, School of Biosciences University of Sheffield Sheffield UK
| | - Kerstin Johannesson
- Tjärnö Marine Laboratory, Department of Marine Sciences University of Gothenburg Strömstad Sweden
| |
Collapse
|
12
|
Venney CJ, Cayuela H, Rougeux C, Laporte M, Mérot C, Normandeau E, Leitwein M, Dorant Y, Præbel K, Kenchington E, Clément M, Sirois P, Bernatchez L. Genome-wide DNA methylation predicts environmentally driven life history variation in a marine fish. Evolution 2023; 77:186-198. [PMID: 36622671 DOI: 10.1093/evolut/qpac028] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/18/2022] [Revised: 09/21/2022] [Accepted: 11/16/2022] [Indexed: 01/10/2023]
Abstract
Epigenetic modifications are thought to be one of the molecular mechanisms involved in plastic adaptive responses to environmental variation. However, studies reporting associations between genome-wide epigenetic changes and habitat-specific variations in life history traits (e.g., lifespan, reproduction) are still scarce, likely due to the recent application of methylome resequencing methods to non-model species. In this study, we examined associations between whole genome DNA methylation and environmentally driven life history variation in 2 lineages of a marine fish, the capelin (Mallotus villosus), from North America and Europe. In both lineages, capelin harbor 2 contrasting life history tactics (demersal vs. beach-spawning). Performing whole genome and methylome sequencing, we showed that life history tactics are associated with epigenetic changes in both lineages, though the effect was stronger in European capelin. Genetic differentiation between the capelin harboring different life history tactics was negligible, but we found genome-wide methylation changes in both lineages. We identified 9,125 European and 199 North American differentially methylated regions (DMRs) due to life history. Gene ontology (GO) enrichment analysis for both lineages revealed an excess of terms related to neural function. Our results suggest that environmental variation causes important epigenetic changes that are associated with contrasting life history tactics in lineages with divergent genetic backgrounds, with variable importance of genetic variation in driving epigenetic variation. Our study emphasizes the potential role of genome-wide epigenetic variation in adaptation to environmental variation.
Collapse
Affiliation(s)
- Clare J Venney
- Institut de Biologie Intégrative et des Systèmes (IBIS), Université Laval, Quebec City, QC, Canada
| | - Hugo Cayuela
- Institut de Biologie Intégrative et des Systèmes (IBIS), Université Laval, Quebec City, QC, Canada.,Department of Ecology and Evolution, University of Lausanne, Lausanne, Switzerland.,University of Lyon, Université Claude Bernard Lyon 1, CNRS, UMR 5558, Laboratoire de Biométrie et Biologie Evolutive, Villeurbanne, France
| | - Clément Rougeux
- Department of Biological Sciences, University of Calgary, Calgary, AB, Canada
| | - Martin Laporte
- Institut de Biologie Intégrative et des Systèmes (IBIS), Université Laval, Quebec City, QC, Canada
| | - Claire Mérot
- Institut de Biologie Intégrative et des Systèmes (IBIS), Université Laval, Quebec City, QC, Canada
| | - Eric Normandeau
- Institut de Biologie Intégrative et des Systèmes (IBIS), Université Laval, Quebec City, QC, Canada
| | - Maëva Leitwein
- Institut de Biologie Intégrative et des Systèmes (IBIS), Université Laval, Quebec City, QC, Canada
| | - Yann Dorant
- Institut de Biologie Intégrative et des Systèmes (IBIS), Université Laval, Quebec City, QC, Canada
| | - Kim Præbel
- Norwegian College of Fishery Science, Faculty of Biosciences, Fisheries and Economics, UiT The Arctic University of Norway, Tromsø, Norway
| | - Ellen Kenchington
- Department of Fisheries and Oceans, Bedford Institute of Oceanography, Dartmouth, NS, Canada
| | - Marie Clément
- Centre for Fisheries Ecosystems Research, Fisheries and Marine Institute of Memorial University of Newfoundland, St. John's, NL, Canada.,Labrador Institute, Memorial University of Newfoundland, Happy Valley-Goose Bay, NL, Canada
| | - Pascal Sirois
- Département des sciences fondamentales, Université du Québec à Chicoutimi, Chicoutimi, QC, Canada
| | - Louis Bernatchez
- Institut de Biologie Intégrative et des Systèmes (IBIS), Université Laval, Quebec City, QC, Canada
| |
Collapse
|
13
|
Canales‐Aguirre CB, Larson WA, McKinney GJ, Claure CE, Rocha JD, Ceballos SG, Cádiz MI, Yáñez JM, Gomez‐Uchida D. Neutral and adaptive loci reveal fine‐scale population structure in Eleginops maclovinus from north Patagonia. Ecol Evol 2022; 12:e9343. [PMID: 36225825 PMCID: PMC9530513 DOI: 10.1002/ece3.9343] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/19/2022] [Revised: 09/05/2022] [Accepted: 09/06/2022] [Indexed: 11/10/2022] Open
Abstract
Patagonia is an understudied area, especially when it comes to population genomic studies with relevance to fishery management. However, the dynamic and heterogeneous landscape in this area can harbor an important but cryptic genetic population structure. Once such information is revealed, it can be integrated into the management of infrequently investigated species. Eleginops maclovinus is a protandrous hermaphrodite species with economic importance for local communities that are currently managed as a single genetic unit. In this study, we sampled five locations distributed across a salinity cline from Northern Patagonia to investigate the genetic population structure of E. maclovinus. We used restriction site‐associated DNA (RAD) sequencing and outlier tests to obtain neutral and adaptive loci, using FST and GEA approaches. We identified a spatial pattern of structuration with gene flow and spatial selection by environmental association. Neutral and adaptive loci showed two and three genetic groups, respectively. The effective population sizes estimated ranged from 572 (Chepu) to 14,454 (Chaitén) and were influenced more by locality than by salinity cline. We found loci putatively associated with salinity suggesting that salinity may act as a selective driver in E. maclovinus populations. These results suggest a complex interaction between genetic drift, gene flow, and natural selection in this area. Our findings also suggest several evolutionary significant units in this area, and the information should be integrated into the management of this species. We discussed the significance of these results for fishery management and suggest future directions to improve our understanding of how E. maclovinus has adapted to the dynamic waters of Northern Patagonia.
Collapse
Affiliation(s)
| | - Wesley A. Larson
- National Oceanographic and Atmospheric Administration, National Marine Fisheries Service, Alaska Fisheries Science CenterAuke Bay LaboratoriesJuneauAlaskaUSA
| | | | - C. Eliza Claure
- Centro i~mar, Universidad de Los LagosPuerto MonttChile,Núcleo Milenio INVASALConcepciónChile
| | - J. Dellis Rocha
- Centro i~mar, Universidad de Los LagosPuerto MonttChile,Núcleo Milenio INVASALConcepciónChile
| | - Santiago G. Ceballos
- Centro Austral de Investigaciones Científicas (CADIC‐CONICET)UshuaiaTierra del FuegoArgentina,Universidad Nacional de Tierra del Fuego (ICPA‐UNTDF)UshuaiaArgentina
| | - María I. Cádiz
- Núcleo Milenio INVASALConcepciónChile,Department of BiologyAarhus UniversityAarhus CDenmark
| | - José M. Yáñez
- Núcleo Milenio INVASALConcepciónChile,Facultad de Ciencias Veterinarias y PecuariasUniversidad de ChileLa PintanaSantiagoChile
| | - Daniel Gomez‐Uchida
- Núcleo Milenio INVASALConcepciónChile,Genomics in Ecology, Evolution & Conservation Lab (GEECLAB), Departamento de Zoología. Facultad de Ciencias Naturales y OceanográficasUniversidad de ConcepciónConcepciónChile
| |
Collapse
|
14
|
Wold J, Koepfli KP, Galla SJ, Eccles D, Hogg CJ, Le Lec MF, Guhlin J, Santure AW, Steeves TE. Expanding the conservation genomics toolbox: Incorporating structural variants to enhance genomic studies for species of conservation concern. Mol Ecol 2021; 30:5949-5965. [PMID: 34424587 PMCID: PMC9290615 DOI: 10.1111/mec.16141] [Citation(s) in RCA: 21] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/20/2021] [Revised: 07/28/2021] [Accepted: 08/18/2021] [Indexed: 12/28/2022]
Abstract
Structural variants (SVs) are large rearrangements (>50 bp) within the genome that impact gene function and the content and structure of chromosomes. As a result, SVs are a significant source of functional genomic variation, that is, variation at genomic regions underpinning phenotype differences, that can have large effects on individual and population fitness. While there are increasing opportunities to investigate functional genomic variation in threatened species via single nucleotide polymorphism (SNP) data sets, SVs remain understudied despite their potential influence on fitness traits of conservation interest. In this future-focused Opinion, we contend that characterizing SVs offers the conservation genomics community an exciting opportunity to complement SNP-based approaches to enhance species recovery. We also leverage the existing literature-predominantly in human health, agriculture and ecoevolutionary biology-to identify approaches for readily characterizing SVs and consider how integrating these into the conservation genomics toolbox may transform the way we manage some of the world's most threatened species.
Collapse
Affiliation(s)
- Jana Wold
- School of Biological Sciences, University of Canterbury, Christchurch, New Zealand
| | - Klaus-Peter Koepfli
- Smithsonian-Mason School of Conservation, Front Royal, Virginia, USA.,Centre for Species Survival, Smithsonian Conservation Biology Institute, National Zoological Park, Washington, District of Columbia, USA.,Computer Technologies Laboratory, ITMO University, Saint Petersburg, Russia
| | - Stephanie J Galla
- School of Biological Sciences, University of Canterbury, Christchurch, New Zealand.,Department of Biological Sciences, Boise State University, Boise, Idaho, USA
| | - David Eccles
- Malaghan Institute of Medical Research, Wellington, New Zealand
| | - Carolyn J Hogg
- School of Life and Environmental Sciences, The University of Sydney, Sydney, NSW, Australia
| | - Marissa F Le Lec
- Department of Biochemistry, University of Otago, Dunedin, Otago, New Zealand
| | - Joseph Guhlin
- Department of Biochemistry, University of Otago, Dunedin, Otago, New Zealand.,Genomics Aotearoa, Dunedin, Otago, New Zealand
| | - Anna W Santure
- School of Biological Sciences, The University of Auckland, Auckland, New Zealand
| | - Tammy E Steeves
- School of Biological Sciences, University of Canterbury, Christchurch, New Zealand
| |
Collapse
|
15
|
Cayuela H, Dorant Y, Forester BR, Jeffries DL, Mccaffery RM, Eby LA, Hossack BR, Gippet JMW, Pilliod DS, Chris Funk W. Genomic signatures of thermal adaptation are associated with clinal shifts of life history in a broadly distributed frog. J Anim Ecol 2021; 91:1222-1238. [PMID: 34048026 PMCID: PMC9292533 DOI: 10.1111/1365-2656.13545] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/30/2021] [Accepted: 05/17/2021] [Indexed: 12/14/2022]
Abstract
Temperature is a critical driver of ectotherm life‐history strategies, whereby a warmer environment is associated with increased growth, reduced longevity and accelerated senescence. Increasing evidence indicates that thermal adaptation may underlie such life‐history shifts in wild populations. Single nucleotide polymorphisms (SNPs) and copy number variants (CNVs) can help uncover the molecular mechanisms of temperature‐driven variation in growth, longevity and senescence. However, our understanding of these mechanisms is still limited, which reduces our ability to predict the response of non‐model ectotherms to global temperature change. In this study, we examined the potential role of thermal adaptation in clinal shifts of life‐history traits (i.e. life span, senescence rate and recruitment) in the Columbia spotted frog Rana luteiventris along a broad temperature gradient in the western United States. We took advantage of extensive capture–recapture datasets of 20,033 marked individuals from eight populations surveyed annually for 14–18 years to examine how mean annual temperature and precipitation influenced demographic parameters (i.e. adult survival, life span, senescence rate, recruitment and population growth). After showing that temperature was the main climatic predictor influencing demography, we used RAD‐seq data (50,829 SNPs and 6,599 putative CNVs) generated for 352 individuals from 31 breeding sites to identify the genomic signatures of thermal adaptation. Our results showed that temperature was negatively associated with annual adult survival and reproductive life span and positively associated with senescence rate. By contrast, recruitment increased with temperature, promoting the long‐term viability of most populations. These temperature‐dependent demographic changes were associated with strong genomic signatures of thermal adaptation. We identified 148 SNP candidates associated with temperature including three SNPs located within protein‐coding genes regulating resistance to cold and hypoxia, immunity and reproduction in ranids. We also identified 39 CNV candidates (including within 38 transposable elements) for which normalized read depth was associated with temperature. Our study indicates that both SNPs and structural variants are associated with temperature and could eventually be found to play a functional role in clinal shifts in senescence rate and life‐history strategies in R. luteiventris. These results highlight the potential role of different sources of molecular variation in the response of ectotherms to environmental temperature variation in the context of global warming.
Collapse
Affiliation(s)
- Hugo Cayuela
- Department of Ecology and Evolution, University of Lausanne, Lausanne, Switzerland
| | - Yann Dorant
- Institut de Biologie Intégrative et des Systèmes (IBIS), Université Laval, Québec, QC, Canada
| | - Brenna R Forester
- Department of Biology, Graduate Degree Program in Ecology, Colorado State University, Fort Collins, CO, USA
| | - Dan L Jeffries
- Department of Ecology and Evolution, University of Lausanne, Lausanne, Switzerland
| | - Rebecca M Mccaffery
- US Geological Survey, Forest and Rangeland Ecosystem Science Center, Port Angeles, WA, USA
| | - Lisa A Eby
- Wildlife Biology Program, W. A. Franke College of Forestry and Conservation, University of Montana, Missoula, MT, USA
| | - Blake R Hossack
- US Geological Survey, Northern Rocky Mountain Science Center, Missoula, MT, USA
| | - Jérôme M W Gippet
- Department of Ecology and Evolution, University of Lausanne, Lausanne, Switzerland
| | - David S Pilliod
- US Geological Survey, Forest and Rangeland Ecosystem Science Center, Boise, ID, USA
| | - W Chris Funk
- Department of Biology, Graduate Degree Program in Ecology, Colorado State University, Fort Collins, CO, USA
| |
Collapse
|