Genome-wide association study (GWAS) for growth rate and age at sexual maturation in Atlantic salmon (Salmo salar)

Optimizing reproductive performance in pangasius catfish broodstock: A review of dietary and molecular strategies

Pangasius catfish, a significant player in the global whitefish market, encounters challenges in aquaculture production sustainability. Quality broodstock maintenance and seed production are impeded by growth, maturation, and fecundity issues. This review investigates the efficacy of strategic nutrient composition and molecular strategies in enhancing broodstock conditions and reproductive performance across various fish species. A notable knowledge gap for Pangasius catfish hampers aquaculture progress. The review assesses nutrient manipulation's impact on reproductive physiology, emphasizing pangasius broodstock. A systematic review analysis following PRISMA guidelines was conducted to identify research trends and hotspots quantitatively, revealing a focus on P. bocourti and fertilization techniques. Addressing this gap, the review offers insights into dietary nutrients manipulation and genetic tool utilization for improved seed production, contributing to pangasius catfish aquaculture sustainability.

Multifaceted framework for defining conservation units: An example from Atlantic salmon (Salmo salar) in Canada

Conservation units represent important components of intraspecific diversity that can aid in prioritizing and protecting at-risk populations, while also safeguarding unique diversity that can contribute to species resilience. In Canada, identification and assessments of conservation units is done by the Committee on the Status of Endangered Wildlife in Canada (COSEWIC). COSEWIC can recognize conservation units below the species level (termed "designatable units"; DUs) if the unit has attributes that make it both discrete and evolutionarily significant. There are various ways in which a DU can meet criteria of discreteness and significance, and increasing access to "big data" is providing unprecedented information that can directly inform both criteria. Specifically, the incorporation of genomic data for an increasing number of non-model species is informing more COSEWIC assessments; thus, a repeatable, robust framework is needed for integrating these data into DU characterization. Here, we develop a framework that uses a multifaceted, weight of evidence approach to incorporate multiple data types, including genetic and genomic data, to inform COSEWIC DUs. We apply this framework to delineate DUs of Atlantic salmon (Salmo salar, L.), an economically, culturally, and ecologically significant species, that is also characterized by complex hierarchical population structure. Specifically, we focus on an in-depth example of how our approach was applied to a previously data limited region of northern Canada that was defined by a single large DU. Application of our framework with newly available genetic and genomic data led to subdividing this DU into three new DUs. Although our approach was developed to meet criteria of COSEWIC, it is widely applicable given similarities in the definitions of a conservation unit.

Genome-wide association study reveals genomic loci of sex differentiation and gonadal development in Plectropomus leopardus

Introduction: Plectropomus leopardus, a commercially significant marine fish, is primarily found in the Western Pacific regions and along the coast of Southeast Asia. A thorough analysis of the molecular mechanisms involved in sex differentiation is crucial for gaining a comprehensive understanding of gonadal development and improving sex control breeding. However, the relevant fundamental studies of P. leopardus are relatively lacking. Methods: In this study, a genome-wide association study (GWAS) was conducted to investigate the genetic basis mechanism of sex differentiation and gonadal developmental traits in P. leopardus utilizing about 6,850,000 high-quality single-nucleotide polymorphisms (SNPs) derived from 168 individuals (including 126 females and 42 males) by the genome-wide efficient mixed-model association (GEMMA) algorithm. Results: The results of these single-trait GWASs showed that 46 SNP loci (-log10 p > 7) significantly associated with sex differentiation, and gonadal development traits were distributed in multiple different chromosomes, which suggested the analyzed traits were all complex traits under multi-locus control. A total of 1,838 potential candidate genes were obtained by considering a less-stringent threshold (-log10 p > 6) and ±100 kb regions surrounding the significant genomic loci. Moreover, 31 candidate genes were identified through a comprehensive analysis of significant GWAS peaks, gene ontology (GO) annotations, and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses, including taf7, ddx6, apoeb, sgk1, a2m, usf1, hsd3b7, dll4, xbp1, tet3, esr1, and gli3. These trait-associated genes have been shown to be involved in germline development, male sex differentiation, gonad morphogenesis, hormone receptor binding, oocyte development, male gonad development, steroidogenesis, estrogen-synthetic pathway, etc. Discussion: In the present study, multiple genomic loci of P. leopardus associated with sex differentiation and gonadal development traits were identified for the first time by using GWAS, providing a valuable resource for further research on the molecular genetic mechanism and sex control in P. leopardus. Our results also can contribute to understanding the genetic basis of the sex differentiation mechanism and gonadal development process in grouper fish.

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 Association Study Reveals SNPs and Candidate Genes Related to Growth and Body Shape in Bighead Carp (Hypophthalmichthys nobilis)

Growth is an economically important trait in bighead carp and other aquaculture species that affects production efficiency. Interestingly, the head of the bighead carp has a high market value in China; therefore, it is important to study the genetic bases of both growth and body shape traits. A genome-wide association study was performed based on 2b-RAD sequencing of 776 individuals to identify SNPs associated with growth and body shape traits, including body weight, body length, body height, and deheaded body length. In total, 26 significant and 19 suggestive SNPs were identified, and more than half of these significant SNPs were clustered in LG16. Two LGs (LG16 and LG21) contained QTLs associated with body weight. Fourteen SNPs of LG16 and two LG21 SNPs were found to be associated with body length. For body height, 12 significantly associated SNPs were identified in LG16. Additionally, 12 SNPs of LG16 and 3 SNPs of LG21 were found to be associated with deheaded body length. Forty-three genes were significantly or suggestively associated with body shape/growth traits based on GWAS results, 18 of which were candidate genes for all BW, BL, BH, and DBL traits. One of these genes, fndc5b, was selected for further analyses. Association analysis revealed that one SNP (g.245 C > T) in the introns of fndc5b was significantly associated with growth-related traits in growth-extreme samples. The mRNA levels of fndc5b in the brains of the lightweight group were significantly higher than those of the heavy-weight group. This study helps to reveal the genetic structure of growth and body development in fish and provides candidate genes for future molecular marker-assisted selection for fast growth and better body conformation in bighead carp.

Genomic signatures of artificial selection in fecundity of Pacific white shrimp, Penaeus vannamei

Penaeusvannamei is the most important economic shrimp in the world. Many selective breeding programs are carried out to improve its production and performance traits. Although significant differences in the reproductive ability of female P. vannamei under artificial breeding conditions have been reported, the genome-wide adaption of the reproductive ability of domesticated female P. vannamei is less investigated. In this study, whole-genome analysis was performed along with pooled DNA sequencing on two fecundity separated bulks, high fecundity bulk (HB), and low fecundity bulk (LB). Each bulk contained 30 individuals from 3 commercial populations. A sequencing depth of >30× was achieved for each bulk, leading to the identification of 625,181 and 629,748 single nucleotide polymorphisms (SNPs) in HB and LB, respectively. Fixation index (Fst) combined with p ratio allowed for the identification of 145 selective sweep regions, with a sequence length of 14.5 Mb, accounting for 0.59% of the genome. Among the 145 selective sweep regions, a total of 64,046 SNPs were identified, and further verification was performed by genotyping 50 candidate SNPs on 60 samples from the offspring of the three populations. Furthermore, 121 genes were screened from the sweep regions. GO annotation and KEGG enrichment analyses showed that partial genes were essential for fecundity regulation. This study provides important information for in-depth investigation of genomic characteristics for long-term selective breeding on the fecundity of female P. vannamei and will also be important for genome-assisted breeding of P. vannamei in the future.

Fine-mapping and association analysis of candidate genes for papilla number in sea cucumber, Apostichopus japonicus

The papilla number is one of the most economically important traits of sea cucumber in the China marketing trade. However, the genetic basis for papilla number diversity in holothurians is still scarce. In the present study, we conducted genome-wide association studies (GWAS) for the trait papilla number of sea cucumbers utilizing a set of 400,186 high-quality SNPs derived from 200 sea cucumbers. Two significant trait-associated SNPs that passed Bonferroni correction (P < 1.25E-7) were located in the intergenic region near PATS1 and the genic region of EIF4G, which were reported to play a pivotal role in cell growth and proliferation. The fine-mapping regions around the top two lead SNPs provided precise causative loci/genes related to papilla formation and cellular activity, including PPP2R3C, GBP1, and BCAS3. Potential SNPs with P < 1E-4 were acquired for the following GO and KEGG enrichment analysis. Moreover, the two lead SNPs were verified in another population of sea cucumber, and the expressive detection of three potential candidate genes PATS1, PPP2R3C, and EIF4G that near or cover the two lead SNPs was conducted in papilla tissue of TG (Top papilla number group) and BG (Bottom papilla number group) by qRT-PCR. We found the significantly higher expression profile of PATS1 (3.34-fold), PPP2R3C (4.90-fold), and EIF4G (4.23-fold) in TG, implying their potential function in papilla polymorphism. The present results provide valuable information to decipher the phenotype differences of the papilla trait and will provide a scientific basis for selective breeding in sea cucumbers.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-022-00139-w.

Detecting Local Adaptation between North and South European Atlantic Salmon Populations

Pollution and other anthropogenic effects have driven a decrease in Atlantic salmon (Salmo salar) in the Iberian Peninsula. The restocking effort carried out in the 1980s, with salmon from northern latitudes with the aim of mitigating the decline of native populations, failed, probably due to the deficiency in adaptation of foreign salmon from northern Europe to the warm waters of the Iberian Peninsula. This result would imply that the Iberian populations of Atlantic salmon have experienced local adaptation in their past evolutionary history, as has been described for other populations of this species and other salmonids. Local adaptation can occur by divergent selections between environments, favoring the fixation of alleles that increase the fitness of a population in the environment it inhabits relative to other alleles favored in another population. In this work, we compared the genomes of different populations from the Iberian Peninsula (Atlantic and Cantabric basins) and Scotland in order to provide tentative evidence of candidate SNPs responsible for the adaptive differences between populations, which may explain the failures of restocking carried out during the 1980s. For this purpose, the samples were genotyped with a 220,000 high-density SNP array (Affymetrix) specific to Atlantic salmon. Our results revealed potential evidence of local adaptation for North Spanish and Scottish populations. As expected, most differences concerned the comparison of the Iberian Peninsula with Scotland, although there were also differences between Atlantic and Cantabric populations. A high proportion of the genes identified are related to development and cellular metabolism, DNA transcription and anatomical structure. A particular SNP was identified within the NADP-dependent malic enzyme-2 (_mMEP-2*), previously reported by independent studies as a candidate for local adaptation in salmon from the Iberian Peninsula. Interestingly, the corresponding SNP within the _mMEP-2* region was consistent with a genomic pattern of divergent selection.

A High-Quality Genome Assembly of Striped Catfish (Pangasianodon hypophthalmus) Based on Highly Accurate Long-Read HiFi Sequencing Data

The HiFi sequencing technology yields highly accurate long-read data with accuracies greater than 99.9% that can be used to improve results for complex applications such as genome assembly. Our study presents a high-quality chromosome-scale genome assembly of striped catfish (Pangasianodon hypophthalmus), a commercially important species cultured mainly in Vietnam, integrating HiFi reads and Hi-C data. A 788.4 Mb genome containing 381 scaffolds with an N50 length of 21.8 Mb has been obtained from HiFi reads. These scaffolds have been further ordered and clustered into 30 chromosome groups, ranging from 1.4 to 57.6 Mb, based on Hi-C data. The present updated assembly has a contig N50 of 14.7 Mb, representing a 245-fold and 4.2-fold improvement over the previous Illumina and Illumina-Nanopore-Hi-C based version, respectively. In addition, the proportion of repeat elements and BUSCO genes identified in our genome is remarkably higher than in the two previously released striped catfish genomes. These results highlight the power of using HiFi reads to assemble the highly repetitive regions and to improve the quality of genome assembly. The updated, high-quality genome assembled in this work will provide a valuable genomic resource for future population genetics, conservation biology and selective breeding studies of striped catfish.

Genome-wide association study reveals the genetic basis of growth trait in yellow catfish with sexual size dimorphism

Sexual size dimorphism has been widely observed in a large number of animals including fish species. Genome-wide association study (GWAS) is a powerful tool to dissect the genetic basis of complex traits, whereas the sex-differences in the genomics of animal complex traits have been ignored in the GWAS analysis. Yellow catfish (Pelteobagrus fulvidraco) is an important aquaculture fish in China with significant sexual size dimorphism. In this study, GWAS was conducted to identify candidate SNPs and genes related to body length (BL) and body weight (BW) in 125 female yellow catfish from a breeding population. In total, one BL-related SNP and three BW-related SNPs were identified to be significantly associated with the traits. Besides, one of these SNPs (Chr15:19195072) was shared in both the BW and BL traits in female yellow catfish, which was further validated in 185 male individuals and located on the exon of stat5b gene. Transgenic yellow catfish and zebrafish that expressed yellow catfish stat5b showed increased growth rate and reduction of sexual size dimorphism. These results not only reveal the genetic basis of growth trait and sexual size dimorphism in fish species, but also provide useful information for the marker-assisted breeding in yellow catfish.

Projecting global mariculture production and adaptation pathways under climate change

The sustainability of global seafood supply to meet increasing demand is facing several challenges, including increasing consumption levels due to a growing human population, fisheries resources over-exploitation and climate change. Whilst growth in seafood production from capture fisheries is limited, global mariculture production is expanding. However, climate change poses risks to the potential seafood production from mariculture. Here, we apply a global mariculture production model that accounts for changing ocean conditions, suitable marine area for farming, fishmeal and fish oil production, farmed species dietary demand, farmed fish price and global seafood demand to project mariculture production under two climate and socio-economic scenarios. We include 85 farmed marine fish and mollusc species, representing about 70% of all mariculture production in 2015. Results show positive global mariculture production changes by the mid and end of the 21^st century relative to the 2000s under the SSP1-2.6 scenario with an increase of 17%±5 and 33%±6, respectively. However, under the SSP5-8.5 scenario, an increase of 8%±5 is projected, with production peaking by mid-century and declining by 16%±5 towards the end of the 21^st century. More than 25% of mariculture-producing nations are projected to lose 40%-90% of their current mariculture production potential under SSP5-8.5 by mid-century. Projected impacts are mainly due to the direct ocean warming effects on farmed species and suitable marine areas, and the indirect impacts of changing availability of forage fishes supplies to produce aquafeed. Fishmeal replacement with alternative protein can lower climate impacts on a subset of finfish production. However, such adaptation measures do not apply to regions dominated by non-feed-based farming (i.e. molluscs) and regions losing substantial marine areas suitable for mariculture. Our study highlights the importance of strong mitigation efforts and the need for different climate adaptation options tailored to the diversity of mariculture systems, to support climate-resilient mariculture development.

Genomic Selection and Genome-wide Association Study for Feed-Efficiency Traits in a Farmed Nile Tilapia (Oreochromis niloticus) Population

Nile tilapia is a key aquaculture species with one of the highest production volumes globally. Genetic improvement of feed efficiency via selective breeding is an important goal, and genomic selection may expedite this process. The aims of this study were to 1) dissect the genetic architecture of feed-efficiency traits in a Nile tilapia breeding population, 2) map the genomic regions associated with these traits and identify candidate genes, 3) evaluate the accuracy of breeding value prediction using genomic data, and 4) assess the impact of the genetic marker density on genomic prediction accuracies. Using an experimental video recording trial, feed conversion ratio (FCR), body weight gain (BWG), residual feed intake (RFI) and feed intake (FI) traits were recorded in 40 full-sibling families from the GIFT (Genetically Improved Farmed Tilapia) Nile tilapia breeding population. Fish were genotyped with a ThermoFisher Axiom 65 K Nile tilapia SNP array. Significant heritabilities, ranging from 0.12 to 0.22, were estimated for all the assessed traits using the genomic relationship matrix. A negative but favourable genetic correlation was found between BWG and the feed-efficiency related traits; -0.60 and -0.63 for FCR and RFI, respectively. While the genome-wide association analyses suggested a polygenic genetic architecture for all the measured traits, there were significant QTL identified for BWG and FI on chromosomes seven and five respectively. Candidate genes previously found to be associated with feed-efficiency traits were located in these QTL regions, including ntrk3a, ghrh and eif4e3. The accuracy of breeding value prediction using the genomic data was up to 34% higher than using pedigree records. A SNP density of approximately 5,000 SNPs was sufficient to achieve similar prediction accuracy as the full genotype data set. Our results highlight the potential of genomic selection to improve feed efficiency traits in Nile tilapia breeding programmes.

Increased accuracy of genomic predictions for growth under chronic thermal stress in rainbow trout by prioritizing variants from GWAS using imputed sequence data

Through imputation of genotypes, genome‐wide association study (GWAS) and genomic prediction (GP) using whole‐genome sequencing (WGS) data are cost‐efficient and feasible in aquaculture breeding schemes. The objective was to dissect the genetic architecture of growth traits under chronic heat stress in rainbow trout (Oncorhynchus mykiss) and to assess the accuracy of GP based on imputed WGS and different preselected single nucleotide polymorphism (SNP) arrays. A total of 192 and 764 fish challenged to a heat stress experiment for 62 days were genotyped using a customized 1 K and 26 K SNP panels, respectively, and then, genotype imputation was performed from a low‐density chip to WGS using 102 parents (36 males and 66 females) as the reference population. Imputed WGS data were used to perform GWAS and test GP accuracy under different preselected SNP scenarios. Heritability was estimated for body weight (BW), body length (BL) and average daily gain (ADG). Estimates using imputed WGS data ranged from 0.33 ± 0.05 to 0.55 ± 0.05 for growth traits under chronic heat stress. GWAS revealed that the top five cumulatively SNPs explained a maximum of 0.94%, 0.86% and 0.51% of genetic variance for BW, BL and ADG, respectively. Some important functional candidate genes associated with growth‐related traits were found among the most important SNPs, including signal transducer and activator of transcription 5B and 3 (STAT5B and STAT3, respectively) and cytokine‐inducible SH2‐containing protein (CISH). WGS data resulted in a slight increase in prediction accuracy compared with pedigree‐based method, whereas preselected SNPs based on the top GWAS hits improved prediction accuracies, with values ranging from 1.2 to 13.3%. Our results support the evidence of the polygenic nature of growth traits when measured under heat stress. The accuracies of GP can be improved using preselected variants from GWAS, and the use of WGS marginally increases prediction accuracy.

Toward Genome-Based Selection in Asian Seabass: What Can We Learn From Other Food Fishes and Farm Animals?

Due to the steadily increasing need for seafood and the plateauing output of fisheries, more fish need to be produced by aquaculture production. In parallel with the improvement of farming methods, elite food fish lines with superior traits for production must be generated by selection programs that utilize cutting-edge tools of genomics. The purpose of this review is to provide a historical overview and status report of a selection program performed on a catadromous predator, the Asian seabass (Lates calcarifer, Bloch 1790) that can change its sex during its lifetime. We describe the practices of wet lab, farm and lab in detail by focusing onto the foundations and achievements of the program. In addition to the approaches used for selection, our review also provides an inventory of genetic/genomic platforms and technologies developed to (i) provide current and future support for the selection process; and (ii) improve our understanding of the biology of the species. Approaches used for the improvement of terrestrial farm animals are used as examples and references, as those processes are far ahead of the ones used in aquaculture and thus they might help those working on fish to select the best possible options and avoid potential pitfalls.

Identification of Growth-Associated Genes by Genome-Wide Association Study and Their Potential Application in the Breeding of Pacific White Shrimp (Litopenaeus vannamei)

The Pacific white shrimp (Litopenaeus vannamei) is the most widely cultured shrimp in the world. A great attention has been paid to improve its body weight (BW) at harvest through genetic selection for decades. Genome-wide association study (GWAS) is a tool to dissect the genetic basis of the traits. In this study, a GWAS approach was conducted to find genes related to BW through genotyping 94,113 single nucleotide polymorphisms (SNPs) in 200 individuals from a breeding population. Four BW-related SNPs located in LG19 and LG39 were identified. Through further candidate gene association analysis, the SNPs in two candidate genes, deoxycytidylate deaminase and non-receptor protein tyrosine kinase, were found to be related with the body weight of the shrimp. Marker-assisted best linear unbiased prediction (MA-BLUP) based on the SNPs in these two genes was used to estimate the breeding values, and the result showed that the highest prediction accuracy of MA-BLUP was increased by 9.4% than traditional BLUP. These results will provide useful information for the marker-assisted breeding in L. vannamei.

First genomic prediction and genome‐wide association for complex growth‐related traits in Rock Bream (Oplegnathus fasciatus)

Rock Bream (Oplegnathus fasciatus) is an important aquaculture species for offshore cage aquaculture and fish stocking of marine ranching in East Asia. Genomic selection has the potential to expedite genetic gain for the key target traits of a breeding program, but has not yet been evaluated in Oplegnathus. The purposes of the present study were to explore the performance of genomic selection to improve breeding value accuracy through real data analyses using six statistical models and to carry out genome‐wide association studies (GWAS) to dissect the genetic architecture of economically vital growth‐related traits (body weight, total length, and body depth) in the O. fasciatus population. After quality control, genotypes for 16,162 SNPs were acquired for 455 fish. Heritability was estimated to be moderate for the three traits (0.38 for BW, 0.33 for TL, and 0.24 for BD), and results of GWAS indicated that the underlying genetic architecture was polygenic. Six statistic models (GBLUP, BayesA, BayesB, BayesC, Bayesian Ridge‐Regression, and Bayesian LASSO) showed similar performance for the predictability of genomic estimated breeding value (GEBV). The low SNP density (around 1 K selected SNP based on GWAS) is sufficient for accurate prediction on the breeding value for the three growth‐related traits in the current studied population, which will provide a good compromise between genotyping costs and predictability in such standard breeding populations advanced. These consequences illustrate that the employment of genomic selection in O. fasciatus breeding could provide advantages for the selection of breeding candidates to facilitate complex economic growth traits.

Sexually Dimorphic Growth Stimulation in a Strain of Growth Hormone Transgenic Coho Salmon (Oncorhynchus kisutch)

Growth hormone (GH) transgenic fish often exhibit remarkable transformations in growth rate and other phenotypes relative to wild-type. The 5750A transgenic coho salmon strain exhibits strong sexually dimorphic growth, with females possessing growth stimulation at a level typical of that seen for both sexes in other strains harbouring the same gene construct (e.g. M77), while males display a modest level of growth stimulation. GH mRNA levels were significantly higher in females than in males of the 5750A strain but equivalent in the M77 strain, indicating sex and transgene insertion locus altered transgene expression. We found that acute estradiol treatments did not influence GH expression in either strain (5750A and M77) or the transgene promoter (metallothionein-B), suggesting that estradiol level was not a significant factor influencing transgene activity. The feminization of XX and XY fish of the 5750A and M77 strains generated all-female groups and resulted in equalized growth of the two genetic sexes, suggesting that the presence of the Y chromosome was not directly capable of influencing the GH transgene-mediated growth in a physiological female conditions. These data suggest that the difference in growth rate seen between the sexes in the 5750A strain arises from non-estradiol-mediated sex influences on gene regulation at the transgene locus. This study shows how genetic factors and transgene insertion sites can influence transgene expression with significant consequent effects on phenotype.

Genomic predictions for fillet yield and firmness in rainbow trout using reduced-density SNP panels

Background

One of the most important goals for the rainbow trout aquaculture industry is to improve fillet yield and fillet quality. Previously, we showed that a 50 K transcribed-SNP chip can be used to detect quantitative trait loci (QTL) associated with fillet yield and fillet firmness. In this study, data from 1568 fish genotyped for the 50 K transcribed-SNP chip and ~ 774 fish phenotyped for fillet yield and fillet firmness were used in a single-step genomic BLUP (ssGBLUP) model to compute the genomic estimated breeding values (GEBV). In addition, pedigree-based best linear unbiased prediction (PBLUP) was used to calculate traditional, family-based estimated breeding values (EBV).

Results

The genomic predictions outperformed the traditional EBV by 35% for fillet yield and 42% for fillet firmness. The predictive ability for fillet yield and fillet firmness was 0.19–0.20 with PBLUP, and 0.27 with ssGBLUP. Additionally, reducing SNP panel densities indicated that using 500–800 SNPs in genomic predictions still provides predictive abilities higher than PBLUP.

Conclusion

These results suggest that genomic evaluation is a feasible strategy to identify and select fish with superior genetic merit within rainbow trout families, even with low-density SNP panels.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07404-9.

Multi-trait GWAS using imputed high-density genotypes from whole-genome sequencing identifies genes associated with body traits in Nile tilapia

Background

Body traits are generally controlled by several genes in vertebrates (i.e. polygenes), which in turn make them difficult to identify through association mapping. Increasing the power of association studies by combining approaches such as genotype imputation and multi-trait analysis improves the ability to detect quantitative trait loci associated with polygenic traits, such as body traits.

Results

A multi-trait genome-wide association study (mtGWAS) was performed to identify quantitative trait loci (QTL) and genes associated with body traits in Nile tilapia (Oreochromis niloticus) using genotypes imputed to whole-genome sequences (WGS). To increase the statistical power of mtGWAS for the detection of genetic associations, summary statistics from single-trait genome-wide association studies (stGWAS) for eight different body traits recorded in 1309 animals were used. The mtGWAS increased the statistical power from the original sample size from 13 to 44%, depending on the trait analyzed. The better resolution of the WGS data, combined with the increased power of the mtGWAS approach, allowed the detection of significant markers which were not previously found in the stGWAS. Some of the lead single nucleotide polymorphisms (SNPs) were found within important functional candidate genes previously associated with growth-related traits in other terrestrial species. For instance, we identified SNP within the α1,6-fucosyltransferase (FUT8), solute carrier family 4 member 2 (SLC4A2), A disintegrin and metalloproteinase with thrombospondin motifs 9 (ADAMTS9) and heart development protein with EGF like domains 1 (HEG1) genes, which have been associated with average daily gain in sheep, osteopetrosis in cattle, chest size in goats, and growth and meat quality in sheep, respectively.

Conclusions

The high-resolution mtGWAS presented here allowed the identification of significant SNPs, linked to strong functional candidate genes, associated with body traits in Nile tilapia. These results provide further insights about the genetic variants and genes underlying body trait variation in cichlid fish with high accuracy and strong statistical support.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07341-z.

Discordance for genotypic sex in phenotypic female Atlantic salmon (Salmo salar) is related to a reduced sdY copy number

The master sex determinant in rainbow trout (Oncorhynchus mykiss), sexually dimorphic on the Y chromosome (sdY), is strongly but not perfectly associated with male phenotype in several other species from the family Salmonidae. Currently, the cause and implications of discordance for sdY-predicted genotypic sex and phenotypic sex in these species is unclear. Using an established multiplex PCR test for exons 2 and 3 of sdY, we demonstrated that sdY-predicted genotypic sex was discordant with histologically evidenced phenotypic sex in 4% of 176 Tasmanian Atlantic salmon. All discordant individuals were phenotypic females presenting a male genotype. Using real-time qPCR assays that we developed and validated for exons 2, 3 and 4 of sdY, all genotype-phenotype discordant females were confirmed to possess sdY, albeit at a reduced number of copies when compared to phenotypic males. The real-time qPCR assays also demonstrated reduced levels of sdY in 30% of phenotypic females that the established multiplex PCR-based test indicated to be devoid of sdY. These findings suggest sdY may be reduced in copy number or mosaicked in the genomic DNA of sdY-positive phenotypic female Atlantic salmon and highlight the importance of understanding the effects of reduced sdY copies on the development of phenotypic sex.

Genome Wide Analysis for Growth at Two Growth Stages in A New Fast-Growing Common Carp Strain (Cyprinus carpio L.)

In order to identify candidate genes or loci associated with growth performance of the newly established common carp strain, Xinlong, we conducted a genome-wide association analysis using 2b-RAD technology on 123 individuals. We constructed two sets of libraries associated with growth-related parameters (weight, length, width and depth) measured at two different grow-out stages. Among the 413,059 SNPs identified using SOAP SNP calling, 147,131 were tested for GWAS after quality filtering. Finally, 39 overlapping SNPs, assigned to four genomic locations, were associated with growth traits in two stages. These loci were assigned to functional classes related to immune response, response to stress, neurogenesis, cholesterol metabolism and development, and proliferation and differentiation of cells. By overlapping results of Plink and EMMAX analyses, we identified three genes: TOX, PLK2 and CD163 (both methods P < 0.05). Our study results could be used for marker-assisted selection to further improve the growth of the Xinlong strain, and illustrate that largely different sets of genes drive the growth of carp in the early and late grow-out stages.

Genome-wide association mapping and accuracy of predictions for amoebic gill disease in Atlantic salmon (Salmo salar)

Amoebic gill disease (AGD) is a parasitic disease caused by the amoeba Paramoeba perurans, which colonizes the gill tissues and causes distress for the host. AGD can cause high morbidity and mortalities in salmonid and non-salmonid fish species. To understand the genetic basis of AGD and improve health status of farmed A. salmon, a population of ~ 6,100 individuals belonging to 150 full-sib families was monitored for development of AGD in the sea of Ireland. The population was followed for two rounds of AGD infections, and fish were gill scored to identify severity of disease in first (N = 3,663) and the second (N = 3,511) infection with freshwater treatment after the first gill-scoring. A subset of this gill-scored population (N = 1,141) from 119 full-sib families were genotyped with 57,184 SNPs using custom-made Affymetrix SNP-chip. GWAS analyses were performed which resulted in five significantly associated SNP variants distributed over chromosome 1, 2 and 5. Three candidate genes; c4, tnxb and slc44a4 were found within QTL region of chromosome 2. The tnxb and c4 genes are known to be a part of innate immune system, and may play a role in resistance to AGD. The gain in prediction accuracy obtained by involving genomic information was 9–17% higher than using traditional pedigree information.

Identification of Candidate Growth-Related SNPs and Genes Using GWAS in Brown-Marbled Grouper (Epinephelus fuscoguttatus)

Brown-marbled grouper, Epinephelus fuscoguttatus, is not only an important commercial fish species, but also an important crossbreeding parent in grouper industry. Improvement of growth traits of this species contributes to the development of grouper breeding. Currently, the development of molecular marker associated with growth of brown-marbled grouper is rare. Thus, we performed the first genome-wide association study (GWAS) for five growth traits in 172 brown-marbled groupers with 43,688 SNPs detected by ddRAD-seq. We identified a total of 5 significant and 18 suggestive QTLs located in multiple chromosomes associated with growth traits. In the 20 kb window of the significant SNPs and suggestive SNPs, 5 and 14 potential candidate genes affecting growth were detected, respectively. Five potential candidate genes near the significantly associated SNPs were selected for expression analysis. Among of which, bmp2k, wasf1, and acyp2 involved in bone development, maintenance of mitochondrion structure, and metabolism were differentially expressed. Interestingly, the SNP 23:29601315 located in the intron of bmp2k was significantly associated with body weight, body length, body height, and body thickness and suggestively associated with total length. We verified the locus using another new group including 123 individuals. The results showed that individuals with CC genotype have better growth traits comparing other individuals. Our findings not only contribute to understanding the molecular mechanism of growth regulation, but also promote the advance of marker-assisted selection in brown-marbled grouper.

Genome-wide identification of loci associated with growth in rainbow trout

Background

Growth is a major economic production trait in aquaculture. Improvements in growth performance will reduce time and cost for fish to reach market size. However, genes underlying growth have not been fully explored in rainbow trout.

Results

A previously developed 50 K gene-transcribed SNP chip, containing ~ 21 K SNPs showing allelic imbalances potentially associated with important aquaculture production traits including body weight, muscle yield, was used for genotyping a total of 789 fish with available phenotypic data for bodyweight gain. Genotyped fish were obtained from two consecutive generations produced in the NCCCWA growth-selection breeding program. Weighted single-step GBLUP (WssGBLUP) was used to perform a genome-wide association (GWA) analysis to identify quantitative trait loci (QTL) associated with bodyweight gain. Using genomic sliding windows of 50 adjacent SNPs, 247 SNPs associated with bodyweight gain were identified. SNP-harboring genes were involved in cell growth, cell proliferation, cell cycle, lipid metabolism, proteolytic activities, chromatin modification, and developmental processes. Chromosome 14 harbored the highest number of SNPs (n = 50). An SNP window explaining the highest additive genetic variance for bodyweight gain (~ 6.4%) included a nonsynonymous SNP in a gene encoding inositol polyphosphate 5-phosphatase OCRL-1. Additionally, based on a single-marker GWA analysis, 33 SNPs were identified in association with bodyweight gain. The highest SNP explaining variation in bodyweight gain was identified in a gene coding for thrombospondin-1 (THBS1) (R² = 0.09).

Conclusion

The majority of SNP-harboring genes, including OCRL-1 and THBS1, were involved in developmental processes. Our results suggest that development-related genes are important determinants for growth and could be prioritized and used for genomic selection in breeding programs.

Beyond large-effect loci: large-scale GWAS reveals a mixed large-effect and polygenic architecture for age at maturity of Atlantic salmon

BACKGROUND

Understanding genetic architecture is essential for determining how traits will change in response to evolutionary processes such as selection, genetic drift and/or gene flow. In Atlantic salmon, age at maturity is an important life history trait that affects factors such as survival, reproductive success, and growth. Furthermore, age at maturity can seriously impact aquaculture production. Therefore, characterizing the genetic architecture that underlies variation in age at maturity is of key interest.

RESULTS

Here, we refine our understanding of the genetic architecture for age at maturity of male Atlantic salmon using a genome-wide association study of 11,166 males from a single aquaculture strain, using imputed genotypes at 512,397 single nucleotide polymorphisms (SNPs). All individuals were genotyped with a 50K SNP array and imputed to higher density using parents genotyped with a 930K SNP array and pedigree information. We found significant association signals on 28 of 29 chromosomes (P-values: 8.7 × 10^-133-9.8 × 10^-8), including two very strong signals spanning the six6 and vgll3 gene regions on chromosomes 9 and 25, respectively. Furthermore, we identified 116 independent signals that tagged 120 candidate genes with varying effect sizes. Five of the candidate genes found here were previously associated with age at maturity in other vertebrates, including humans.

DISCUSSION

These results reveal a mixed architecture of large-effect loci and a polygenic component that consists of multiple smaller-effect loci, suggesting a more complex genetic architecture of Atlantic salmon age at maturity than previously thought. This more complex architecture will have implications for selection on this key trait in aquaculture and for management of wild salmon populations.

Diversity in reproductive seasonality in the three-spined stickleback, Gasterosteus aculeatus

The annual timing of reproduction is a key life history trait with a large effect on fitness. Populations often vary in the timing and duration of reproduction to adapt to different seasonality of ecological and environmental variables between habitats. However, little is known about the molecular genetic mechanisms underlying interpopulation variation in reproductive seasonality. Here, we demonstrate that the three-spined stickleback (Gasterosteus aculeatus) is a good model for molecular genetic analysis of variations in reproductive seasonality. We first compiled data on reproductive seasons of diverse ecotypes, covering marine-anadromous, lake and stream ecotypes, of three-spined stickleback inhabiting a wide range of latitudes. Our analysis showed that both ecotype and latitude significantly contribute to variation in reproductive seasons. Stream ecotypes tend to start breeding earlier and end later than other ecotypes. Populations from lower latitudes tend to start breeding earlier than those from higher latitudes in all three ecotypes. Additionally, stream ecotypes tend to have extended breeding seasons at lower latitudes than at higher latitudes, leading to nearly year-round reproduction in the most southern stream populations. A review of recent progress in our understanding of the physiological mechanisms underlying seasonal reproduction in the three-spined stickleback indicates that photoperiod is an important external cue that stimulates and/or suppresses reproduction in this species. Taking advantage of genomic tools available for this species, the three-spined stickleback will be a good model to investigate what kinds of genes and mutations underlie variations in the physiological signalling pathways that regulate reproduction in response to photoperiod.

High-Throughput Single Nucleotide Polymorphism (SNP) Discovery and Validation Through Whole-Genome Resequencing in Nile Tilapia (Oreochromis niloticus)

Nile tilapia (Oreochromis niloticus) is the second most important farmed fish in the world and a sustainable source of protein for human consumption. Several genetic improvement programs have been established for this species in the world. Currently, the estimation of genetic merit of breeders is typically based on genealogical and phenotypic information. Genome-wide information can be exploited to efficiently incorporate traits that are difficult to measure into the breeding goal. Thus, single nucleotide polymorphisms (SNPs) are required to investigate phenotype-genotype associations and determine the genomic basis of economically important traits. We performed de novo SNP discovery in three different populations of farmed Nile tilapia. A total of 29.9 million non-redundant SNPs were identified through Illumina (HiSeq 2500) whole-genome resequencing of 326 individual samples. After applying several filtering steps, including removing SNP based on genotype and site quality, presence of Mendelian errors, and non-unique position in the genome, a total of 50,000 high-quality SNPs were selected for the development of a custom Illumina BeadChip SNP panel. These SNPs were highly informative in the three populations analyzed showing between 43,869 (94%) and 46,139 (99%) SNPs in Hardy-Weinberg Equilibrium; 37,843 (76%) and 45,171(90%) SNPs with a minor allele frequency (MAF) higher than 0.05; and 43,450 (87%) and 46,570 (93%) SNPs with a MAF higher than 0.01. The 50K SNP panel developed in the current work will be useful for the dissection of economically relevant traits, enhancing breeding programs through genomic selection, as well as supporting genetic studies in farmed populations of Nile tilapia using dense genome-wide information.

Genetic segregation for male body coloration and female mate preference in the guppy

OBJECTIVE

The purpose of this study was to segregate the genetic lines responsible for the orange area of coloration in males and the response to orange coloration exhibited by females in the guppy (Poecilia reticulata) through artificial selection. This study is part of a project that uses QTL-seq to search for candidate genes involved in male orange coloration and female response to male coloration. We created two lines: high-selected lines of males having large areas of orange spots and of females with high response to male orange coloration; and low-selected lines of males having small areas of orange spots and of females with low response to male orange coloration.

RESULTS

The male orange area and the female response became significantly different between high- and low-selected lines after three generations of artificial selection. This indicates that the differences in the frequencies of alleles at loci affecting the orange area and the female response between the lines increased over the generations through selection.

Genome-Wide Association Study of Growth and Body-Shape-Related Traits in Large Yellow Croaker (Larimichthys crocea) Using ddRAD Sequencing

Large yellow croaker (Larimichthys crocea) is an economically important marine fish species of China. Due to overfishing and marine pollution, the wild stocks of this croaker have collapsed in the past decades. Meanwhile, the cultured croaker is facing the difficulties of reduced genetic diversity and low growth rate. To explore the molecular markers related to the growth traits of croaker and providing the related SNPs for the marker-assisted selection, we used double-digest restriction-site associated DNA (ddRAD) sequencing to dissect the genetic bases of growth traits in a cultured population and identify the SNPs that associated with important growth traits by GWAS. A total of 220 individuals were genotyped by ddRAD sequencing. After quality control, 27,227 SNPs were identified in 220 samples and used for GWAS analysis. We identified 13 genome-wide significant associated SNPs of growth traits on 8 chromosomes, and the beta P of these SNPs ranged from 0.01 to 0.86. Through the definition of candidate regions and gene annotation, candidate genes related to growth were identified, including important regulators such as fgf18, fgf1, nr3c1, cyp8b1, fabp2, cyp2r1, ppara, and ccm2l. We also identified SNPs and candidate genes that significantly associated with body shape, including bmp7, col1a1, col11a2, and col18a1, which are also economically important traits for large yellow croaker aquaculture. The results provided insights into the genetic basis of growth and body shape in large yellow croaker population and would provide reliable genetic markers for molecular marker-assisted selection in the future. Meanwhile, the result established a basis for our subsequent fine mapping and related gene study.

Genome-Wide Patterns of Population Structure and Linkage Disequilibrium in Farmed Nile Tilapia (Oreochromis niloticus)

Nile tilapia (Oreochromis niloticus) is one of the most produced farmed fish in the world and represents an important source of protein for human consumption. Farmed Nile tilapia populations are increasingly based on genetically improved stocks, which have been established from admixed populations. To date, there is scarce information about the population genomics of farmed Nile tilapia, assessed by dense single nucleotide polymorphism (SNP) panels. The patterns of linkage disequilibrium (LD) may affect the success of genome-wide association studies (GWAS) and genomic selection (GS), and also provide key information about demographic history of farmed Nile tilapia populations. The objectives of this study were to provide further knowledge about the population structure and LD patterns, as well as, estimate the effective population size (N e ) for three farmed Nile tilapia populations, one from Brazil (POP A) and two from Costa Rica (POP B and POP C). A total of 55 individuals from each population, were genotyped using a 50K SNP panel selected from a whole-genome sequencing (WGS) experiment. The first two principal components explained about 20% of the total variation and clearly differentiated between the three populations. Population genetic structure analysis showed evidence of admixture, especially for POP C. The contemporary N e estimated, based on LD values, ranged from 78 to 159. No differences were observed in the LD decay among populations, with a rapid decrease of r 2 with increasing inter-marker distance. Average r 2 between adjacent SNP pairs ranged from 0.19 to 0.03 for both POP A and C, and 0.20 to 0.03 f or POP B. Based on the number of independent chromosome segments in the Nile tilapia genome, at least 9.4, 7.6, and 4.6K SNPs for POP A, POP B, and POP C respectively, are required for the implementation of GS in the present farmed Nile tilapia populations.

Association Mapping Based on a Common-Garden Migration Experiment Reveals Candidate Genes for Migration Tendency in Brown Trout

A better understanding of the environmental and genetic contribution to migratory behavior and the evolution of traits linked to migration is crucial for fish conservation and fisheries management. Up to date, a few genes with unequivocal influence on the adoption of alternative migration strategies have been identified in salmonids. Here, we used a common garden set-up to measure individual migration distances of generally highly polymorphic brown trout Salmo trutta from two populations. Fish from the assumedly resident population showed clearly shorter migration distances than the fish from the assumed migratory population at the ages of 2 and 3 years. By using two alternative analytical pipelines with 22186 and 18264 SNPs obtained through RAD-sequencing, we searched for associations between individual migration distance, and both called genotypes and genotype probabilities. None of the SNPs showed statistically significant individual effects on migration after correction for multiple testing. By choosing a less stringent threshold, defined as an overlap of the top 0.1% SNPs identified by the analytical pipelines, GAPIT and Angsd, we identified eight candidate genes that are potentially linked to individual migration distance. While our results demonstrate large individual and population level differences in migration distances, the detected genetic associations were weak suggesting that migration traits likely have multigenic control.

Variation in migration pattern, broodstock origin, and family productivity of coho salmon hatchery populations in British Columbia, Canada, derived from parentage-based tagging

In salmonid parentage-based tagging (PBT) applications, entire hatchery broodstocks are genotyped, and subsequently, progeny can be nonlethally sampled and assigned back to their parents using parentage analysis, thus identifying their hatchery of origin and brood year (i.e., age). Inter- and intrapopulation variability in migration patterns, life history traits, and fishery contributions can be determined from PBT analysis of samples derived from both fisheries and escapements (portion of a salmon population that does not get caught in fisheries and returns to its natal river to spawn). In the current study of southern British Columbia coho salmon (Oncorhynchus kisutch) populations, PBT analysis provided novel information on intrapopulation heterogeneity among males in the total number of progeny identified in fisheries and escapements, the proportion of progeny sampled from fisheries versus escapement, the proportion of two-year-old progeny (jacks) produced, and the within-season return time of progeny. Fishery recoveries of coho salmon revealed heterogeneity in migration patterns among and within populations, with recoveries from north and central coast fisheries distinguishing "northern migrating" from "resident" populations. In northern migrating populations, the mean distance between fishery captures of sibs (brothers and sisters) was significantly less than the mean distance between nonsibs, indicating the possible presence of intrapopulation genetic heterogeneity for migration pattern. Variation among populations in productivity and within populations in fish catchability indicated that population selection and broodstock management can be implemented to optimize harvest benefits from hatcheries. Application of PBT provided valuable information for assessment and management of hatchery-origin coho salmon in British Columbia.

Genome-Wide Association Study and Cost-Efficient Genomic Predictions for Growth and Fillet Yield in Nile Tilapia (Oreochromis niloticus)

Fillet yield (FY) and harvest weight (HW) are economically important traits in Nile tilapia production. Genetic improvement of these traits, especially for FY, are lacking, due to the absence of efficient methods to measure the traits without sacrificing fish and the use of information from relatives to selection. However, genomic information could be used by genomic selection to improve traits that are difficult to measure directly in selection candidates, as in the case of FY. The objectives of this study were: (i) to perform genome-wide association studies (GWAS) to dissect the genetic architecture of FY and HW, (ii) to evaluate the accuracy of genotype imputation and (iii) to assess the accuracy of genomic selection using true and imputed low-density (LD) single nucleotide polymorphism (SNP) panels to determine a cost-effective strategy for practical implementation of genomic information in tilapia breeding programs. The data set consisted of 5,866 phenotyped animals and 1,238 genotyped animals (108 parents and 1,130 offspring) using a 50K SNP panel. The GWAS were performed using all genotyped and phenotyped animals. The genotyped imputation was performed from LD panels (LD0.5K, LD1K and LD3K) to high-density panel (HD), using information from parents and 20% of offspring in the reference set and the remaining 80% in the validation set. In addition, we tested the accuracy of genomic selection using true and imputed genotypes comparing the accuracy obtained from pedigree-based best linear unbiased prediction (PBLUP) and genomic predictions. The results from GWAS supports evidence of the polygenic nature of FY and HW. The accuracy of imputation ranged from 0.90 to 0.98 for LD0.5K and LD3K, respectively. The accuracy of genomic prediction outperformed the estimated breeding value from PBLUP. The use of imputation for genomic selection resulted in an increased relative accuracy independent of the trait and LD panel analyzed. The present results suggest that genotype imputation could be a cost-effective strategy for genomic selection in Nile tilapia breeding programs.

Scanning of Genetic Variants and Genetic Mapping of Phenotypic Traits in Gilthead Sea Bream Through ddRAD Sequencing

Gilthead sea bream (Sparus aurata) is a teleost of considerable economic importance in Southern European aquaculture. The aquaculture industry shows a growing interest in the application of genetic methods that can locate phenotype-genotype associations with high economic impact. Through selective breeding, the aquaculture industry can exploit this information to maximize the financial yield. Here, we present a Genome Wide Association Study (GWAS) of 112 samples belonging to seven different sea bream families collected from a Greek commercial aquaculture company. Through double digest Random Amplified DNA (ddRAD) Sequencing, we generated a per-sample genetic profile consisting of 2,258 high-quality Single Nucleotide Polymorphisms (SNPs). These profiles were tested for association with four phenotypes of major financial importance: Fat, Weight, Tag Weight, and the Length to Width ratio. We applied two methods of association analysis. The first is the typical single-SNP to phenotype test, and the second is a feature selection (FS) method through two novel algorithms that are employed for the first time in aquaculture genomics and produce groups with multiple SNPs associated to a phenotype. In total, we identified 9 single SNPs and 6 groups of SNPs associated with weight-related phenotypes (Weight and Tag Weight), 2 groups associated with Fat, and 16 groups associated with the Length to Width ratio. Six identified loci (Chr4:23265532, Chr6:12617755, Chr:8:11613979, Chr13:1098152, Chr15:3260819, and Chr22:14483563) were present in genes associated with growth in other teleosts or even mammals, such as semaphorin-3A and neurotrophin-3. These loci are strong candidates for future studies that will help us unveil the genetic mechanisms underlying growth and improve the sea bream aquaculture productivity by providing genomic anchors for selection programs.

Genome Scan for Genomic Regions and Genes Associated with Growth Trait in Pacific White Shrimp Litopeneaus vannamei

The Pacific white shrimp Litopeneaus vannmei (L. vannmei) is a predominant aquaculture shrimp species worldwide, and it is considered as the aquaculture species with the highest single output value. Advances in selective breeding have accelerated the development of L. vannmei aquaculture. Recently, the genome-wide association studies (GWAS) have been applied in aquaculture animals and markers associated with economic traits were identified. In this study, we focused on the growth trait of L. vannamei and performed GWAS to identify SNPs or genes associated with growth. Genomic regions in linkage group 7, 27, 33, and 38 were identified to be associated with body weight and body length of the shrimp. Further, candidate gene association analysis was performed in two independent populations and the result demonstrated that the SNPs in the genes protein kinase C delta type and ras-related protein Rap-2a were significantly associated with the growth trait of L. vannamei. This study showed that GWAS analysis is an efficient approach for screening trait-related markers or genes. The genomic regions and genes identified in this study are essential for further fine mapping of growth-related genes. The identified markers will provide useful information for marker-assisted selection in L. vannamei.

A Novel Candidate Gene Associated With Body Weight in the Pacific White Shrimp Litopenaeus vannamei

Improvements of growth traits are always the focus in selective breeding programs for the Pacific white shrimp Litopenaeus vannamei (L. vannamei). Identification of growth-related genes or markers can contribute to the application of modern breeding technologies, and thus accelerate the genetic improvement of growth traits. The aim of this study was to identify the genes and molecular markers associated with the growth traits of L. vannamei. A population of 200 individuals was genotyped using 2b-RAD techniques for genome-wide linkage disequilibrium (LD) analysis and genome-wide association study (GWAS). The results showed that the LD decayed fast in the studied population, which suggest that it is feasible to fine map the growth-related genes with GWAS in L. vannamei. One gene designated as LvSRC, encoding the class C scavenger receptor (SRC), was identified as a growth-related candidate gene by GWAS. Further targeted sequencing of the candidate gene in another population of 322 shrimps revealed that several non-synonymous mutations within LvSRC were significantly associated with the body weight (P < 0.01), and the most significant marker (SRC_24) located in the candidate gene could explain 13% of phenotypic variance. The current results provide not only molecular markers for genetic improvement in L. vannamei, but also new insights for understanding the growth regulation mechanism in penaeid shrimp.

Whole Genome Linkage Disequilibrium and Effective Population Size in a Coho Salmon (Oncorhynchus kisutch) Breeding Population Using a High-Density SNP Array

The estimation of linkage disequilibrium between molecular markers within a population is critical when establishing the minimum number of markers required for association studies, genomic selection, and inferring historical events influencing different populations. This work aimed to evaluate the extent and decay of linkage disequilibrium in a coho salmon breeding population using a high-density SNP array. Linkage disequilibrium was estimated between a total of 93,502 SNPs found in 64 individuals (33 dams and 31 sires) from the breeding population. The markers encompass all 30 coho salmon chromosomes and comprise 1,684.62 Mb of the genome. The average density of markers per chromosome ranged from 48.31 to 66 per 1 Mb. The minor allele frequency averaged 0.26 (with a range from 0.22 to 0.27). The overall average linkage disequilibrium among SNPs pairs measured as r² was 0.10. The Average r² value decreased with increasing physical distance, with values ranging from 0.21 to 0.07 at a distance lower than 1 kb and up to 10 Mb, respectively. An r² threshold of 0.2 was reached at distance of approximately 40 Kb. Chromosomes Okis05, Okis15 and Okis28 showed high levels of linkage disequilibrium (>0.20 at distances lower than 1 Mb). Average r² values were lower than 0.15 for all chromosomes at distances greater than 4 Mb. An effective population size of 43 was estimated for the population 10 generations ago, and 325, for 139 generations ago. Based on the effective number of chromosome segments, we suggest that at least 74,000 SNPs would be necessary for an association mapping study and genomic predictions. Therefore, the SNP panel used allowed us to capture high-resolution information in the farmed coho salmon population. Furthermore, based on the contemporary N_e, a new mate allocation strategy is suggested to increase the effective population size.

The influence of vgll3 genotypes on sea age at maturity is altered in farmed mowi strain Atlantic salmon

BACKGROUND

In Atlantic salmon in the wild, age at maturity is strongly influenced by the vgll3 locus. Under farming conditions, light, temperature and feeding regimes are known significantly advance or delay age at maturity. However, the potential influence of the vgll3 locus on the maturation of salmon reared under farming conditions has been rarely investigated, especially in females.

RESULTS

Here, we reared domesticated salmon (mowi strain) with different vgll3 genotypes under standard farming conditions until they matured at either one, two or more than two sea winters. Interestingly, and in contrast to previous findings in the wild, we were not able to identify a link between vgll3 and age at maturity in females when reared under farming conditions. For males however, we found that the probability of delaying maturation from one to two sea winters was significantly lower in fish homozygous for the early allele compared to homozygous fish for the late allele, while the probability for heterozygous fish was intermediate. These data also contrast to previous findings in the wild where the early allele has been reported as dominant. However, we found that the probability of males delaying maturation from two to three sea winters was regulated in the same manner as the wild.

CONCLUSIONS

Collectively, our data suggest that increased growth rates in mowi salmon, caused by high feed intake and artificial light and temperature regimes together with other possible genetic/epigenetic components, may significantly influence the impact that the vgll3 locus has on age at maturity, especially in females. In turn, our results show that the vgll3 locus can only to a large extent be used in selective breeding to control age at maturation in mowi males. In summary, we here show that in contrast to the situation in wild salmon, under farming conditions vgll3 does not seem to influence age at maturity in mowi females whereas in mowi males, maturing as one or two sea winters it alters the early allele effect from dominant to intermediate.

Design and characterization of an 87k SNP genotyping array for Arctic charr (Salvelinus alpinus)

We have generated a high-density, high-throughput genotyping array for characterizing genome-wide variation in Arctic charr (Salvelinus alpinus). Novel single nucleotide polymorphisms (SNPs) were identified in charr from the Fraser, Nauyuk and Tree River aquaculture strains, which originated from northern Canada and fish from Iceland using high coverage sequencing, reduced representation sequencing and RNA-seq datasets. The array was designed to capture genome-wide variation from a diverse suite of Arctic charr populations. Cross validation of SNPs from various sources and comparison with previously published Arctic charr SNP data provided a set of candidate SNPs that generalize across populations. Further candidate SNPs were identified based on minor allele frequency, association with RNA transcripts, even spacing across intergenic regions and association with the sex determining (sdY) gene. The performance of the 86,503 SNP array was assessed by genotyping Fraser, Nauyuk and Tree River strain individuals, as well as wild Icelandic Arctic charr. Overall, 63,060 of the SNPs were polymorphic within at least one group and 36.8% were unique to one of the four groups, suggesting that the array design allows for characterization of both within and across population genetic diversity. The concordance between sdY markers and known phenotypic sex indicated that the array can accurately determine the sex of individuals based on genotype alone. The Salp87k genotyping array provides researchers and breeders the opportunity to analyze genetic variation in Arctic charr at a more detailed level than previously possible.

Polygenic and sex specific architecture for two maturation traits in farmed Atlantic salmon

BACKGROUND

A key developmental transformation in the life of all vertebrates is the transition to sexual maturity, whereby individuals are capable of reproducing for the first time. In the farming of Atlantic salmon, early maturation prior to harvest size has serious negative production impacts.

RESULTS

We report genome wide association studies (GWAS) using fish measured for sexual maturation in freshwater or the marine environment. Genotypic data from a custom 50 K single nucleotide polymorphism (SNP) array was used to identify 13 significantly associated SNP for freshwater maturation with the most strongly associated on chromosomes 10 and 11. A higher number of associations (48) were detected for marine maturation, and the two peak loci were found to be the same for both traits. The number and broad distribution of GWAS hits confirmed a highly polygenetic nature, and GWAS performed separately within males and females revealed sex specific genetic behaviour for loci co-located with positional candidate genes phosphatidylinositol-binding clathrin assembly protein-like (picalm) and membrane-associated guanylate kinase, WW and PDZ domain-containing protein 2 (magi2).

CONCLUSIONS

The results extend earlier work and have implications for future applied breeding strategies to delay maturation in this important aquaculture species.

Genome-wide association analysis for body weight identifies candidate genes related to development and metabolism in rainbow trout (Oncorhynchus mykiss)

Growth is one of the most important traits from both a physiological and economic perspective in aquaculture species. Thus, identifying the genomic regions and genes underpinning genetic variation for this trait is of particular interest in several fish species, including rainbow trout. In this work, we perform a genome-wide association study (GWAS) to identify the genomic regions associated with body weight at tagging (BWT) and at 18 months (BW18M) using a dense SNP panel (57 k) and 4596 genotyped rainbow trout from 105 full-sib families belonging to a Chilean breeding population. Analysis was performed by means of single-step GBLUP approach. Genetic variance explained by 20 adjacent SNP windows across the whole genome is reported. To further explore candidate genes, we focused on windows that explained the highest proportion of genetic variance in the top 10 chromosomes for each trait. The main window from the top 10 chromosomes was explored by BLAST using the first and last SNP position of each window to determine the target nucleotide sequence. As expected, the percentage of genetic variance explained by windows was relatively low, due to the polygenic nature of body weight. The most important genomic region for BWT and BW18M were located on chromosomes 15 and 24 and they explained 2.14% and 3.02% of the genetic variance for each trait, respectively. Candidate genes including several growth factors, genes involved in development of skeletal muscle and bone tissue and nutrient metabolism were identified within the associated regions for both traits BWT and BW18M. These results indicate that body weight is polygenic in nature in rainbow trout, with the most important loci explaining as much as 3% of the genetic variance for the trait. The genes identified here represent good candidates for further functional validation to uncover biological mechanisms underlying variation for growth in rainbow trout.

Comparing genomic signatures of domestication in two Atlantic salmon (Salmo salar L.) populations with different geographical origins

Selective breeding and genetic improvement have left detectable signatures on the genomes of domestic species. The elucidation of such signatures is fundamental for detecting genomic regions of biological relevance to domestication and improving management practices. In aquaculture, domestication was carried out independently in different locations worldwide, which provides opportunities to study the parallel effects of domestication on the genome of individuals that have been selected for similar traits. In this study, we aimed to detect potential genomic signatures of domestication in two independent pairs of wild/domesticated Atlantic salmon populations of Canadian and Scottish origins, respectively. Putative genomic regions under divergent selection were investigated using a 200K SNP array by combining three different statistical methods based either on allele frequencies (LFMM, Bayescan) or haplotype differentiation (Rsb). We identified 337 and 270 SNPs potentially under divergent selection in wild and hatchery populations of Canadian and Scottish origins, respectively. We observed little overlap between results obtained from different statistical methods, highlighting the need to test complementary approaches for detecting a broad range of genomic footprints of selection. The vast majority of the outliers detected were population-specific but we found four candidate genes that were shared between the populations. We propose that these candidate genes may play a role in the parallel process of domestication. Overall, our results suggest that genetic drift may have override the effect of artificial selection and/or point toward a different genetic basis underlying the expression of similar traits in different domesticated strains. Finally, it is likely that domestication may predominantly target polygenic traits (e.g., growth) such that its genomic impact might be more difficult to detect with methods assuming selective sweeps.

Population Genomic Structure and Genome-Wide Linkage Disequilibrium in Farmed Atlantic Salmon (Salmo salar L.) Using Dense SNP Genotypes

Chilean Farmed Atlantic salmon (Salmo salar) populations were established with individuals of both European and North American origins. These populations are expected to be highly genetically differentiated due to evolutionary history and poor gene flow between ancestral populations from different continents. The extent and decay of linkage disequilibrium (LD) among single nucleotide polymorphism (SNP) impacts the implementation of genome-wide association studies and genomic selection and provides relevant information about demographic processes of fish populations. We assessed the population structure and characterized the extent and decay of LD in three Chilean commercial populations of Atlantic salmon with North American (NAM), Scottish (SCO), and Norwegian (NOR) origin. A total of 123 animals were genotyped using a 159 K SNP Axiom® myDesignTM Genotyping Array. A total of 32 K SNP markers, representing the common SNPs along the three populations after quality control were used. The principal component analysis explained 78.9% of the genetic diversity between populations, clearly discriminating between populations of North American and European origin, and also between European populations. NAM had the lowest effective population size, followed by SCO and NOR. Large differences in the LD decay were observed between populations of North American and European origin. An r 2 threshold of 0.2 was estimated for marker pairs separated by 7,800, 64, and 50 kb in the NAM, SCO, and NOR populations, respectively. In this study we show that this SNP panel can be used to detect association between markers and traits of interests and also to capture high-resolution information for genome-enabled predictions. Also, we suggest the feasibility to achieve similar prediction accuracies using a smaller SNP data set for the NAM population, compared with samples with European origin which would need a higher density SNP array.

Development and Validation of 58K SNP-Array and High-Density Linkage Map in Nile Tilapia (O. niloticus)

Despite being the second most important aquaculture species in the world accounting for 7.4% of global production in 2015, tilapia aquaculture has lacked genomic tools like SNP-arrays and high-density linkage maps to improve selection accuracy and accelerate genetic progress. In this paper, we describe the development of a genotyping array containing more than 58,000 SNPs for Nile tilapia (Oreochromis niloticus). SNPs were identified from whole genome resequencing of 32 individuals from the commercial population of the Genomar strain, and were selected for the SNP-array based on polymorphic information content and physical distribution across the genome using the Orenil1.1 genome assembly as reference sequence. SNP-performance was evaluated by genotyping 4991 individuals, including 689 offspring belonging to 41 full-sib families, which revealed high-quality genotype data for 43,588 SNPs. A preliminary genetic linkage map was constructed using Lepmap2 which in turn was integrated with information from the O_niloticus_UMD1 genome assembly to produce an integrated physical and genetic linkage map comprising 40,186 SNPs distributed across 22 linkage groups (LGs). Around one-third of the LGs showed a different recombination rate between sexes, with the female being greater than the male map by a factor of 1.2 (1632.9 to 1359.6 cM, respectively), with most LGs displaying a sigmoid recombination profile. Finally, the sex-determining locus was mapped to position 40.53 cM on LG23, in the vicinity of the anti-Müllerian hormone (amh) gene. These new resources has the potential to greatly influence and improve the genetic gain when applying genomic selection and surpass the difficulties of efficient selection for invasively measured traits in Nile tilapia.

Genetic Mapping of Head Size Related Traits in Common Carp (Cyprinus carpio)

Head size is important economic trait for many aquaculture fish which is directly linked to their carcass yield. The genetic basis of head size trait remains unclear in many widely cultured fish species. Common carp (Cyprinus carpio) is one of the most widely studied fish due to its importance on both economic and environmental aspects. In this study, we performed genome-wide association study using 433 Yellow River carp individuals from multiple families to identify loci and genes potentially associated with head size related traits including head length (HL), head length/body length ratio (HBR), eye diameter (ED), and eye cross (EC). QTL mapping was utilized to filter the effects of population stratification and improve power for the candidates identification in the largest surveyed family with a published genetic linkage map. Twelve SNPs showed significant for head size traits in GWAS and 18 QTLs were identified in QTL mapping. Our study combining both GWAS and QTL mapping could compensate the deficiency from each other and advance our understanding of head size traits in common carp. To acquire a better understanding of the correlation between head size and body growth, we also performed comparisons between QTLs of head size traits and growth-related traits. Candidate genes underlying head size traits were identified surrounding the significant SNPs, including parvalbumin, srpk2, fsrp5, igf1, igf3, grb10, igf1r, notch2, sfrp2. Many of these genes have been identified with potential functions on bone formation and growth. Igf1 was a putative gene associated with both head size and body growth in Yellow River carp. The teleost-specific igf3 was a candidate head size related gene, related to both HL and HBR. Our study also indicated the importance of Igf signaling pathway for both growth and head size determination in common carp, which could be potentially used in future selective breeding in common carp as well as other species.

Genome-Scale Association Study of Abnormal Scale Pattern in Yellow River Carp Identified Previously Known Causative Gene in European Mirror Carp

Common carp (Cyprinus carpio) is one of the most widely studied fish species due to its great economic value and strong environmental adaptability. Scattered scale, a typical phenotype of the mirror carp that is derived from Europe, has never been observed in the Yellow River carp previously. We recently identified approximately one fourth of the F1 progenies displaying scattered scale in a full-sib Yellow River carp family in our breeding program, despite both parents that showed wild type with normal scale patterns. This family provides us unique materials to investigate the genetic basis underlying the abnormal scale mutant in Yellow River carp population. Genome-wide association study (GWAS) and association mapping were performed based on genome-wide single nucleotide polymorphisms (SNP) genotyped with common carp 250 K SNP genotyping array in 82 samples of the Yellow River carp family. We identified a 1.4 Mb genome region that was significantly associated with abnormal scattered scale patterns. We further identified a deletion mutation in fibroblast growth factor receptor 1 a1 (fgfr1a1) gene within this genome region. Amplification and sequencing analysis of this gene revealed a 311-bp deletion in intron 10 and exon 11, which proved that fgfr1a1 could be the causal gene responsible for abnormal scattered scale in the Yellow River carp family. Since similar fragment mutation with 306-bp and 310-bp deletions had been previously reported as causal mutation of scattered scale patterns in the mirror carp, we speculate that either the deletion mutation was introduced from Europe-derived mirror carp or the deletion independently occurred in the mutation hotspot in fgfr1a1 gene. The results provided insights into the genetic basis of scale pattern mutant in Yellow River carp population, which would help us to eliminate the recessive allele of the abnormal scale patterns in Yellow River carp population by molecular marker-assisted breeding.

Genomic Selection for Growth Traits in Pacific Oyster (Crassostrea gigas): Potential of Low-Density Marker Panels for Breeding Value Prediction

Pacific oysters are a key aquaculture species globally, and genetic improvement via selective breeding is a major target. Genomic selection has the potential to expedite genetic gain for key target traits of a breeding program, but has not yet been evaluated in oyster. The recent development of SNP arrays for Pacific oyster (Crassostrea gigas) raises the opportunity to test genomic selection strategies for polygenic traits. In this study, a population of 820 oysters (comprising 23 full-sibling families) were genotyped using a medium density SNP array (23 K informative SNPs), and the genetic architecture of growth-related traits [shell height (SH), shell length (SL), and wet weight (WW)] was evaluated. Heritability was estimated to be moderate for the three traits (0.26 ± 0.06 for SH, 0.23 ± 0.06 for SL and 0.35 ± 0.05 for WW), and results of a GWAS indicated that the underlying genetic architecture was polygenic. Genomic prediction approaches were used to estimate breeding values for growth, and compared to pedigree based approaches. The accuracy of the genomic prediction models (GBLUP) outperformed the traditional pedigree approach (PBLUP) by ∼25% for SL and WW, and ∼30% for SH. Further, reduction in SNP marker density had little impact on prediction accuracy, even when density was reduced to a few hundred SNPs. These results suggest that the use of genomic selection in oyster breeding could offer benefits for the selection of breeding candidates to improve complex economic traits at relatively modest cost.

Genetic (co)variation in skin pigmentation patterns and growth in rainbow trout

From a physiological-behavioral perspective, it has been shown that fish with a higher density of black eumelanin spots are more dominant, less sensitive to stress, have higher feed intake, better feed efficiency and therefore are larger in size. Thus, we hypothesized that genetic (co)variation between skin pigmentation patterns and growth exists and it is advantageous in rainbow trout. The objective of this study was to determine the genetic relationships between skin pigmentation patterns and BW in a breeding population of rainbow trout. We performed a genetic analysis of pigmentation traits including dorsal color (DC), lateral band (LB) intensity, amount of spotting above (SA) and below (SB) the lateral line, and BW at harvest (HW). Variance components were estimated using a multi-trait linear animal model fitted by restricted maximum likelihood. Estimated heritabilities were 0.08±0.02, 0.17±0.03, 0.44±0.04, 0.17±0.04 and 0.23±0.04 for DC, LB, SA, SB and HW, respectively. Genetic correlations between HW and skin color traits were 0.42±0.13, 0.32±0.14 and 0.25±0.11 for LB, SA and SB, respectively. These results indicate positive, but low to moderate genetic relationships between the amount of spotting and BW in rainbow trout. Thus, higher levels of spotting are genetically associated with better growth performance in this population.

Genomic Tools and Selective Breeding in Molluscs

The production of most farmed molluscs, including mussels, oysters, scallops, abalone, and clams, is heavily dependent on natural seed from the plankton. Closing the lifecycle of species in hatcheries can secure independence from wild stocks and enables long-term genetic improvement of broodstock through selective breeding. Genomic techniques have the potential to revolutionize hatchery-based selective breeding by improving our understanding of the characteristics of mollusc genetics that can pose a challenge for intensive aquaculture and by providing a new suite of tools for genetic improvement. Here we review characteristics of the life history and genetics of molluscs including high fecundity, self-fertilization, high genetic diversity, genetic load, high incidence of deleterious mutations and segregation distortion, and critically assess their impact on the design and effectiveness of selective breeding strategies. A survey of the results of current breeding programs in the literature show that selective breeding with inbreeding control is likely the best strategy for genetic improvement of most molluscs, and on average growth rate can be improved by 10% per generation and disease resistance by 15% per generation across the major farmed species by implementing individual or family-based selection. Rapid advances in sequencing technology have resulted in a wealth of genomic resources for key species with the potential to greatly improve hatchery-based selective breeding of molluscs. In this review, we catalog the range of genomic resources currently available for molluscs of aquaculture interest and discuss the bottlenecks, including lack of high-quality reference genomes and the relatively high cost of genotyping, as well as opportunities for applying genomics-based selection.