Assessing the genetic diversity of farmed and wild Rufiji tilapia (Oreochromis urolepis urolepis) populations using ddRAD sequencing

Genetic diversity patterns in farmed rainbow trout (Oncorhynchus mykiss) populations using genome-wide SNP and haplotype data

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.

Genetic diversity and genome-scale population structure of wild Indian major carp, Labeo catla (Hamilton, 1822), revealed by genotyping-by-sequencing

Labeo catla (catla) is the second most commercially important and widely cultured Indian major carp (IMC). It is indigenous to the Indo-Gangetic riverine system of India and the rivers of Bangladesh, Nepal, Myanmar, and Pakistan. Despite the availability of substantial genomic resources in this important species, detailed information on the genome-scale population structure using SNP markers is yet to be reported. In the present study, the identification of genome-wide single nucleotide polymorphisms (SNPs) and population genomics of catla was undertaken by re-sequencing six catla populations of riverine origin from distinct geographical regions. DNA isolated from 100 samples was used to perform genotyping-by-sequencing (GBS). A published catla genome with 95% genome coverage was used as the reference for mapping reads using BWA software. From a total of 472 million paired-end (150 × 2 bp) raw reads generated in this study, we identified 10,485 high-quality polymorphic SNPs using the STACKS pipeline. Expected heterozygosity (He) across the populations ranged from 0.162 to 0.20, whereas observed heterozygosity (Ho) ranged between 0.053 and 0.06. The nucleotide diversity (π) was the lowest (0.168) in the Ganga population. The within-population variation was found to be higher (95.32%) than the among-population (4.68%) variation. However, genetic differentiation was observed to be low to moderate, with F_st values ranging from 0.020 to 0.084, and the highest between Brahmani and Krishna populations. Bayesian and multivariate techniques were used to further evaluate the population structure and supposed ancestry in the studied populations using the structure and discriminant analysis of principal components (DAPC), respectively. Both analyses revealed the existence of two separate genomic clusters. The maximum number of private alleles was observed in the Ganga population. The findings of this study will contribute to a deeper understanding of the population structure and genetic diversity of wild populations of catla for future research in fish population genomics.

Genome-Wide SNP Discovery and Population Genetic Analysis of Mesocentrotus nudus in China Seas

Mesocentrotus nudus is an important commercially aquatic species because of its high edible and medicinal values. However, wild stocks have dramatically decreased in recent decades. Understanding the population structure and genetic diversity can provide vital information for genetic conservation and improvement. In the present study, the genotyping-by-sequencing (GBS) approach was adopted to identify the genome-wide single-nucleotide polymorphisms (SNPs) from a collection of 80 individuals consisting of five geographical populations (16 individuals from each population), covering the natural habitats of M. nudus in China seas. An average of 0.96-Gb clean reads per sample were sequenced, and a total of 51,738 biallelic SNPs were identified. Based on these SNPs, diversity index analysis showed that all populations have a similar pattern with positive F_is (0.136) and low Ne (724.3). Low genetic differentiation and high genetic connectivity among five geographical populations were detected by pairwise F_st, principal component analysis (PCA), admixture, and phylogenetic analysis. Besides, two YWL individuals originating from an isolated ancestor may imply that there is a genetically differentiated population in the adjacent sea. Overall, the results showed that GBS is an effective method to detect genome-wide SNPs for M. nudus and suggested that the protective measures and the investigation with larger spatial scale and sample size for M. nudus should be carried out in the future.

A genomic‐based vision on the genetic diversity and key performance traits in selectively bred Arctic charr ( Salvelinus alpinus )

Routine implementation of genomic information for guiding selection decisions is not yet common in the majority of aquaculture species. Reduced representation sequencing approaches offer a cost‐effective solution for obtaining genome‐wide information in species with a limited availability of genomic resources. In the current study, we implemented double‐digest restriction site‐associated DNA sequencing (ddRAD‐seq) on an Arctic charr strain with the longest known history of selection (approximately 40 years) aiming to improve selection decisions. In total, 1730 animals reared at four different farms in Sweden and spanning from year classes 2013–2017 were genotyped using ddRAD‐seq. Approximately 5000 single nucleotide polymorphisms (SNPs) were identified, genetic diversity‐related metrics were estimated, and genome‐wide association studies (GWAS) for body length at different time points and age of sexual maturation were conducted. Low genetic differentiation amongst animals from the different farms was observed based on both the results from pairwise F_st values and principal component analysis (PCA). The existence of associations was investigated between the mean genome‐wide heterozygosity of each full‐sib family (year class 2017) and the corresponding inbreeding coefficient or survival to the eyed stage. A moderate correlation (−0.33) was estimated between the mean observed heterozygosity of each full‐sib family and the corresponding inbreeding coefficient, while no linear association was obtained with the survival to the eyed stage. GWAS did not detect loci with major effect for any of the studied traits. However, genomic regions explaining more than 1% of the additive genetic variance for either studied traits were suggested across 14 different chromosomes. Overall, key insights valuable for future selection decisions of Arctic charr have been obtained, suggesting ddRAD as an attractive genotyping platform for obtaining genome‐wide information in a cost‐effective manner.

Comparing divergence landscapes from reduced-representation and whole genome resequencing in the yellow-rumped warbler (Setophaga coronata) species complex

Researchers seeking to generate genomic data for non-model organisms are faced with a number of trade-offs when deciding which method to use. The selection of reduced representation approaches versus whole genome resequencing will ultimately affect the marker density, sequencing depth, and the number of individuals that can multiplexed. These factors can affect researchers' ability to accurately characterize certain genomic features, such as landscapes of divergence-how F_ST varies across the genomes. To provide insight into the effect of sequencing method on the estimation of divergence landscapes, we applied an identical bioinformatic pipeline to three generations of sequencing data (GBS, ddRAD, and WGS) produced for the same system, the yellow-rumped warbler species complex. We compare divergence landscapes generated using each method for the myrtle warbler (Setophaga coronata coronata) and the Audubon's warbler (S. c. auduboni), and for Audubon's warblers with deeply divergent mtDNA resulting from mitochondrial introgression. We found that most high-F_ST peaks were not detected in the ddRAD data set, and that while both GBS and WGS were able to identify the presence of large peaks, WGS was superior at a finer scale. Comparing Audubon's warblers with divergent mitochondrial haplotypes, only WGS allowed us to identify small (10-20 kb) regions of elevated differentiation, one of which contained the nuclear-encoded mitochondrial gene NDUFAF3. We calculated the cost per base pair for each method and found it was comparable between GBS and WGS, but significantly higher for ddRAD. These comparisons highlight the advantages of WGS over reduced representation methods when characterizing landscapes of divergence.

Genotyping Strategies Using ddRAD Sequencing in Farmed Arctic Charr (Salvelinus alpinus)

Incorporation of genomic technologies into fish breeding programs is a modern reality, promising substantial advances regarding the accuracy of selection, monitoring the genetic diversity and pedigree record verification. Single nucleotide polymorphism (SNP) arrays are the most commonly used genomic tool, but the investments required make them unsustainable for emerging species, such as Arctic charr (Salvelinus alpinus), where production volume is low. The requirement to genotype a large number of animals for breeding practices necessitates cost effective genotyping approaches. In the current study, we used double digest restriction site-associated DNA (ddRAD) sequencing of either high or low coverage to genotype Arctic charr from the Swedish national breeding program and performed analytical procedures to assess their utility in a range of tasks. SNPs were identified and used for deciphering the genetic structure of the studied population, estimating genomic relationships and implementing an association study for growth-related traits. Missing information and underestimation of heterozygosity in the low coverage set were limiting factors in genetic diversity and genomic relationship analyses, where high coverage performed notably better. On the other hand, the high coverage dataset proved to be valuable when it comes to identifying loci that are associated with phenotypic traits of interest. In general, both genotyping strategies offer sustainable alternatives to hybridization-based genotyping platforms and show potential for applications in aquaculture selective breeding.