1
|
Doublet M, Degalez F, Lagarrigue S, Lagoutte L, Gueret E, Allais S, Lecerf F. Variant calling and genotyping accuracy of ddRAD-seq: Comparison with 20X WGS in layers. PLoS One 2024; 19:e0298565. [PMID: 39058708 PMCID: PMC11280156 DOI: 10.1371/journal.pone.0298565] [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: 01/26/2024] [Accepted: 05/23/2024] [Indexed: 07/28/2024] Open
Abstract
Whole Genome Sequencing (WGS) remains a costly or unsuitable method for routine genotyping of laying hens. Until now, breeding companies have been using or developing SNP chips. Nevertheless, alternatives methods based on sequencing have been developed. Among these, reduced representation sequencing approaches can offer sequencing quality and cost-effectiveness by reducing the genomic regions covered by sequencing. The aim of this study was to evaluate the ability of double digested Restriction site Associated DNA sequencing (ddRAD-seq) to identify and genotype SNPs in laying hens, by comparison with a presumed reliable WGS approach. Firstly, the sensitivity and precision of variant calling and the genotyping reliability of ddRADseq were determined. Next, the SNP Call Rate (CRSNP) and mean depth of sequencing per SNP (DPSNP) were compared between both methods. Finally, the effect of multiple combinations of thresholds for these parameters on genotyping reliability and amount of remaining SNPs in ddRAD-seq was studied. In raw form, the ddRAD-seq identified 349,497 SNPs evenly distributed on the genome with a CRSNP of 0.55, a DPSNP of 11X and a mean genotyping reliability rate per SNP of 80%. Considering genomic regions covered by expected enzymatic fragments (EFs), the sensitivity of the ddRAD-seq was estimated at 32.4% and its precision at 96.4%. The low CRSNP and DPSNP values were explained by the detection of SNPs outside the EFs theoretically generated by the ddRAD-seq protocol. Indeed, SNPs outside the EFs had significantly lower CRSNP (0.25) and DPSNP (1X) values than SNPs within the EFs (0.7 and 17X, resp.). The study demonstrated the relationship between CRSNP, DPSNP, genotyping reliability and the number of SNPs retained, to provide a decision-support tool for defining filtration thresholds. Severe quality control over ddRAD-seq data allowed to retain a minimum of 40% of the SNPs with a CcR of 98%. Then, ddRAD-seq was defined as a suitable method for variant calling and genotyping in layers.
Collapse
Affiliation(s)
| | | | | | | | - Elise Gueret
- MGX-Montpellier GenomiX, Univ. Montpellier, CNRS, INSERM, Montpellier, France
| | | | | |
Collapse
|
2
|
Paul K, Restoux G, Phocas F. Genome-wide detection of positive and balancing signatures of selection shared by four domesticated rainbow trout populations (Oncorhynchus mykiss). Genet Sel Evol 2024; 56:13. [PMID: 38389056 PMCID: PMC10882880 DOI: 10.1186/s12711-024-00884-9] [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: 03/09/2023] [Accepted: 02/12/2024] [Indexed: 02/24/2024] Open
Abstract
BACKGROUND Evolutionary processes leave footprints along the genome over time. Highly homozygous regions may correspond to positive selection of favorable alleles, while maintenance of heterozygous regions may be due to balancing selection phenomena. We analyzed data from 176 fish from four disconnected domestic rainbow trout populations that were genotyped using a high-density Axiom Trout genotyping 665K single nucleotide polymorphism array, including 20 from the US and 156 from three French lines. Using methods based on runs of homozygosity and extended haplotype homozygosity, we detected signatures of selection in these four populations. RESULTS Nine genomic regions that included 253 genes were identified as being under positive selection in all four populations Most were located on chromosome 2 but also on chromosomes 12, 15, 16, and 20. In addition, four heterozygous regions that contain 29 genes that are putatively under balancing selection were also shared by the four populations. These were located on chromosomes 10, 13, and 19. Regardless of the homozygous or heterozygous nature of the regions, in each region, we detected several genes that are highly conserved among vertebrates due to their critical roles in cellular and nuclear organization, embryonic development, or immunity. We identified new candidate genes involved in rainbow trout fitness, as well as 17 genes that were previously identified to be under positive selection, 10 of which in other fishes (auts2, atp1b3, zp4, znf135, igf-1α, brd2, col9a2, mrap2, pbx1, and emilin-3). CONCLUSIONS Using material from disconnected populations of different origins allowed us to draw a genome-wide map of signatures of positive selection that are shared between these rainbow trout populations, and to identify several regions that are putatively under balancing selection. These results provide a valuable resource for future investigations of the dynamics of genetic diversity and genome evolution during domestication.
Collapse
Affiliation(s)
- Katy Paul
- Université Paris-Saclay, INRAE, AgroParisTech, GABI, 78350, Jouy-en-Josas, France
| | - Gwendal Restoux
- Université Paris-Saclay, INRAE, AgroParisTech, GABI, 78350, Jouy-en-Josas, France
| | - Florence Phocas
- Université Paris-Saclay, INRAE, AgroParisTech, GABI, 78350, Jouy-en-Josas, France.
| |
Collapse
|
3
|
Longo A, Kurta K, Vanhala T, Jeuthe H, de Koning DJ, Palaiokostas C. Genetic diversity patterns in farmed rainbow trout (Oncorhynchus mykiss) populations using genome-wide SNP and haplotype data. Anim Genet 2024; 55:87-98. [PMID: 37994156 DOI: 10.1111/age.13378] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/10/2023] [Revised: 10/25/2023] [Accepted: 11/01/2023] [Indexed: 11/24/2023]
Abstract
Rainbow trout is one of the most popular aquaculture species worldwide, with a long history of domestication. However, limited information exists about the genetic diversity of farmed rainbow trout populations globally, with most available reports relying on low-throughput genotyping technologies. Notably, no information exists about the genetic diversity status of farmed rainbow trout in Sweden. Double-digest restriction-site-associated DNA sequencing was performed on more than 500 broodfish from two leading producers in Sweden and from the country's national breeding program. Following the detection of single nucleotide polymorphisms (SNPs), genetic diversity was studied by using either individual SNPs (n = 8680; one SNP retained per 300 bp sequence reads) or through SNP haplotypes (n = 20 558; all SNPs retained in 300 bp sequence reads). Similar amounts of genetic diversity were found amongst the three populations when individual SNPs were used. Furthermore, principal component analysis and discriminant analysis of principal components suggested two genetic clusters with the two industry populations grouped together. Genetic differentiation based on the FST fixation index was ~0.01 between the industry populations and ~0.05 when those were compared with the breeding program. Preliminary estimates of effective population size (Ne ) and inbreeding (based on runs of homozygosity; FROH ) were similar amongst the three populations (Ne ≈ 50-80; median FROH ≈ 0.11). Finally, the haplotype-based analysis suggested that animals from the breeding program had higher shared coancestry levels than those from the other two populations. Overall, our study provides novel insights into the genetic diversity and structure of Sweden's three main farmed rainbow trout populations, which could guide their future management.
Collapse
Affiliation(s)
- Alessio Longo
- Department of Animal Breeding and Genetics, Swedish University of Agricultural Sciences, Uppsala, Sweden
| | - Khrystyna Kurta
- Department of Animal Breeding and Genetics, Swedish University of Agricultural Sciences, Uppsala, Sweden
| | - Tytti Vanhala
- Department of Animal Breeding and Genetics, Swedish University of Agricultural Sciences, Uppsala, Sweden
| | - Henrik Jeuthe
- Department of Animal Breeding and Genetics, Swedish University of Agricultural Sciences, Uppsala, Sweden
- Aquaculture Center North, Kälarne, Sweden
| | - Dirk-Jan de Koning
- Department of Animal Breeding and Genetics, Swedish University of Agricultural Sciences, Uppsala, Sweden
| | - Christos Palaiokostas
- Department of Animal Breeding and Genetics, Swedish University of Agricultural Sciences, Uppsala, Sweden
| |
Collapse
|
4
|
Criscione A, Chessari G, Cesarani A, Ablondi M, Asti V, Bigi D, Bordonaro S, Ciampolini R, Cipolat-Gotet C, Congiu M, De Palo P, Landi V, Macciotta NPP, Matassino D, Portolano B, Riggio S, Sabbioni A, Sardina MT, Senczuk G, Tumino S, Vasini M, Ciani E, Mastrangelo S. Analysis of ddRAD-seq data provides new insights into the genomic structure and patterns of diversity in Italian donkey populations. J Anim Sci 2024; 102:skae165. [PMID: 38874306 PMCID: PMC11214105 DOI: 10.1093/jas/skae165] [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: 03/07/2024] [Accepted: 06/13/2024] [Indexed: 06/15/2024] Open
Abstract
With more than 150 recognized breeds, donkeys assume relevant economic importance, especially in developing countries. Even if the estimated number of heads worldwide is 53M, this species received less attention than other livestock species. Italy has traditionally been considered one of the cradles of European donkey breeding, and despite a considerable loss of biodiversity, today still counts nine autochthonous populations. A total of 220 animals belonging to nine different populations were genotyped using the double-digest restriction site associated DNA (ddRAD) sequencing to investigate the pattern of diversity using a multi-technique approach. A total of 418,602,730 reads were generated and successfully demultiplexed to obtain a medium-density SNP genotypes panel with about 27K markers. The diversity indices showed moderate levels of variability. The genetic distances and relationships, largely agree with the breeding history of the donkey populations under investigation. The results highlighted the separation of populations based on their genetic origin or geographical proximity between breeding areas, showed low to moderate levels of admixture, and indicated a clear genetic difference in some cases. For some breeds, the results also validate the success of proper management conservation plans. Identified runs of homozygosity islands, mapped within genomic regions related to immune response and local adaptation, are consistent with the characteristics of the species known for its rusticity and adaptability. This study is the first exhaustive genome-wide analysis of the diversity of Italian donkey populations. The results emphasized the high informativeness of genome-wide markers retrieved through the ddRAD approach. The findings take on great significance in designing and implementing conservation strategies. Standardized genotype arrays for donkey species would make it possible to combine worldwide datasets to provide further insights into the evolution of the genomic structure and origin of this important genetic resource.
Collapse
Affiliation(s)
- Andrea Criscione
- Dipartimento di Agricoltura, Alimentazione e Ambiente, University of Catania, Catania, Italy
| | - Giorgio Chessari
- Dipartimento di Agricoltura, Alimentazione e Ambiente, University of Catania, Catania, Italy
- Department of Animal Sciences, Georg-August-University Göttingen, Göttingen, Germany
| | - Alberto Cesarani
- Dipartimento di Agraria, University of Sassari, Sassari, Italy
- Department of Animal and Dairy Science, University of Georgia, Athens, USA
| | - Michela Ablondi
- Dipartimento di Scienze Medico-Veterinarie, University of Parma, Parma, Italy
| | - Vittoria Asti
- Dipartimento di Scienze Medico-Veterinarie, University of Parma, Parma, Italy
| | - Daniele Bigi
- Dipartimento di Scienze e Tecnologie Agro-Alimentari, University of Bologna, Bologna, Italy
| | - Salvatore Bordonaro
- Dipartimento di Agricoltura, Alimentazione e Ambiente, University of Catania, Catania, Italy
| | | | | | - Michele Congiu
- Dipartimento di Agraria, University of Sassari, Sassari, Italy
| | - Pasquale De Palo
- Dipartimento di Medicina Veterinaria, University of Bari Aldo Moro, Valenzano, Italy
| | - Vincenzo Landi
- Dipartimento di Medicina Veterinaria, University of Bari Aldo Moro, Valenzano, Italy
| | | | - Donato Matassino
- Consorzio per la Sperimentazione, Divulgazione e Applicazione di Biotecniche Innovative, Benevento, Italy
| | - Baldassare Portolano
- Dipartimento di Scienze Agrarie, Alimentari e Forestali, University of Palermo, Palermo, Italy
| | - Silvia Riggio
- Dipartimento di Scienze Agrarie, Alimentari e Forestali, University of Palermo, Palermo, Italy
| | - Alberto Sabbioni
- Dipartimento di Scienze Medico-Veterinarie, University of Parma, Parma, Italy
| | - Maria Teresa Sardina
- Dipartimento di Scienze Agrarie, Alimentari e Forestali, University of Palermo, Palermo, Italy
| | - Gabriele Senczuk
- Dipartimento di Agricoltura, Ambiente e Alimenti, University of Molise, Campobasso, Italy
| | - Serena Tumino
- Dipartimento di Agricoltura, Alimentazione e Ambiente, University of Catania, Catania, Italy
| | - Matteo Vasini
- Associazione Nazionale Allevatori delle Razze Equine ed Asinine Italiane, ANAREAI, Roma, Italy
| | - Elena Ciani
- Dipartimento di Bioscienze, Biotecnologie e Biofarmaceutica, University of Bari Aldo Moro, Bari, Italy
| | - Salvatore Mastrangelo
- Dipartimento di Scienze Agrarie, Alimentari e Forestali, University of Palermo, Palermo, Italy
| |
Collapse
|
5
|
Vera M, Aparicio E, Heras S, Abras A, Casanova A, Roldán MI, García-Marin JL. Regional environmental and climatic concerns on preserving native gene pools of a least concern species: Brown trout lineages in Mediterranean streams. THE SCIENCE OF THE TOTAL ENVIRONMENT 2023; 862:160739. [PMID: 36502686 DOI: 10.1016/j.scitotenv.2022.160739] [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: 08/31/2022] [Revised: 11/10/2022] [Accepted: 12/03/2022] [Indexed: 06/17/2023]
Abstract
The European brown trout, Salmo trutta, is a cold-adapted fish reported as a Least Concern species in the IUCN Red List. This species colonized new territories from southern refuges during the last glacial melting, but during the 20th century suffered from anthropic impacts on its habitats. The long-time survival of the species relies on the genetic diversity within and among populations. Brown trout is among the genetically most diverse vertebrate species; however, native populations in Mediterranean rivers have dramatically suffered of introgressive hybridization from extensive releases of evolutionary distant non-native Atlantic stocks. In addition, in Mediterranean rivers climate change will result in unsuitable conditions for the species during the 21st century. Using brown trout populations at the headstreams of a Pyrenean river as a model, this paper revised how hatchery releases have affected the native gene pools and how environmental and climatic variables controlled the amount of local introgression at intra-basin level. Introgressive hybridization was detected in all studied sites. Ten times larger divergence was observed among populations at tributaries than among populations along the main stem. A highly impacted population distributed in a long transect in the main stem suggested that hatchery fish move towards the main stem wherever released. From already highly impacted populations and despite the cessation of hatchery releases, warmer temperatures and lower precipitation expected from climate change will extend the introgressive hybridization along the basin, contributing to the extinction of the native gene pools. Based on available morphological distinction of native, hatchery and hybrid brown trout, we advocate the involvement of regional social groups (e.g. riverside dwellers, anglers, conservationists, hikers) in citizen science programs to detect the spread of non-native phenotypes along the rivers. These are cheap and fast methods to collaborate with fishery managers in the preservation and recovery of the regional native populations.
Collapse
Affiliation(s)
- Manuel Vera
- Departamento de Zoología, Xenética e Antropología Física, Campus Lugo, Universidad de Santiago de Compostela, 27002 Lugo, Spain.
| | - Enric Aparicio
- GRECO, Institute of Aquatic Ecology, University of Girona, 17071 Girona, Spain
| | - Sandra Heras
- Laboratori d'Ictiologia Genètica, Campus Montilivi, Universitat de Girona, 17003 Girona, Spain.
| | - Alba Abras
- Laboratori d'Ictiologia Genètica, Campus Montilivi, Universitat de Girona, 17003 Girona, Spain.
| | - Adrián Casanova
- Laboratori d'Ictiologia Genètica, Campus Montilivi, Universitat de Girona, 17003 Girona, Spain
| | - Maria-Inés Roldán
- Laboratori d'Ictiologia Genètica, Campus Montilivi, Universitat de Girona, 17003 Girona, Spain.
| | - Jose-Luis García-Marin
- Laboratori d'Ictiologia Genètica, Campus Montilivi, Universitat de Girona, 17003 Girona, Spain.
| |
Collapse
|
6
|
Valette T, Leitwein M, Lascaux JM, Desmarais E, Berrebi P, Guinand B. Redundancy analysis, genome-wide association studies and the pigmentation of brown trout (Salmo trutta L.). JOURNAL OF FISH BIOLOGY 2023; 102:96-118. [PMID: 36218076 DOI: 10.1111/jfb.15243] [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: 09/26/2021] [Accepted: 10/04/2022] [Indexed: 06/16/2023]
Abstract
The association of molecular variants with phenotypic variation is a main issue in biology, often tackled with genome-wide association studies (GWAS). GWAS are challenging, with increasing, but still limited, use in evolutionary biology. We used redundancy analysis (RDA) as a complimentary ordination approach to single- and multitrait GWAS to explore the molecular basis of pigmentation variation in brown trout (Salmo trutta) belonging to wild populations impacted by hatchery fish. Based on 75,684 single nucleotide polymorphic (SNP) markers, RDA, single- and multitrait GWAS allowed the extraction of 337 independent colour patterning loci (CPLs) associated with trout pigmentation traits, such as the number of red and black spots on flanks. Collectively, these CPLs (i) mapped onto 35 out of 40 brown trout linkage groups indicating a polygenic genomic architecture of pigmentation, (ii) were found to be associated with 218 candidate genes, including 197 genes formerly mentioned in the literature associated to skin pigmentation, skin patterning, differentiation or structure notably in a close relative, the rainbow trout (Onchorhynchus mykiss), and (iii) related to functions relevant to pigmentation variation (e.g., calcium- and ion-binding, cell adhesion). Annotated CPLs include genes with well-known pigmentation effects (e.g., PMEL, SLC45A2, SOX10), but also markers associated with genes formerly found expressed in rainbow or brown trout skins. RDA was also shown to be useful to investigate management issues, especially the dynamics of trout pigmentation submitted to several generations of hatchery introgression.
Collapse
|
7
|
Genetic and Phenotypic Characteristics of the Salmo trutta Complex in Italy. APPLIED SCIENCES-BASEL 2022. [DOI: 10.3390/app12073219] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/05/2023]
Abstract
Salmonid fish have become ecological and research models of study in the field of conservation genetics and genomics. Over the last decade, brown trout have received a high level of interest in research and publications. The term Salmo trutta complex is used to indicate the large number of geographic forms present in the species Salmo trutta. In Europe, the S. trutta complex consists (based on mitochondrial DNA control region analysis) of seven major evolutionary lineages: Atlantic (AT), Mediterranean (ME), Adriatic (AD), Danubian (DA), Marmoratus (MA), Duero (DU) and Tigris (TI). In several nations, the difficulty of identifying some lineages derives from their wide phenotypic and geographic plasticity and the presence of mixed lineages (due to introgressive hybridization with domestic AT populations). In Italy, the S. trutta complex populations living in the Tyrrhenian area and on the main islands (Sicily, Sardinia and Corsica) showed high genetic diversity. Currently, on the Italian Red List, the protected (near threatened) populations are the AD and ME lineages. Recent studies based on traditional (mitochondrial and nuclear markers) and NGS (next-generation sequencing) analyses have clarified some genetic differences between the populations of the Tyrrhenian region, Sicily, Sardinia and Corsica. Native populations in Sardinia belong to the AD lineage, while those living in Corsica are mainly characterized by the AD, MA and ME haplotypes. In Sicily, in the area of the Iblei mountains, an AT lineage (North African) exists. According to some authors, the term Salmo macrostigma should only be used for populations in North Africa. The use of genotyping methods based on mtDNA and nuclear markers and the latest generation sequencing techniques can improve the study of populations and evolutionary lineages in areas where there are overlaps and hybridization phenomena.
Collapse
|