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

Genomic Signatures of Domestication in European Seabass (Dicentrarchus labrax L.) Reveal a Potential Role for Epigenetic Regulation in Adaptation to Captivity

Genome scans provide a comprehensive method to explore genome-wide variation associated with traits under study. However, linking individual genes to broader functional groupings and pathways is often challenging, yet crucial for understanding the evolutionary mechanisms underlying these traits. This task is particularly relevant for multi-trait processes such as domestication, which are influenced by complex interactions between numerous genetic and non-genetic factors, including epigenetic regulation. As various traits within the broader spectrum of domestication are selected in concert over time, this process offers an opportunity to identify broader functional overlaps and understand the integrated genetic architecture underlying these traits. In this study, we analyzed approximately 600,000 SNPs from a Pool-Seq experiment comparing eight natural-origin and 12 farmed populations of European seabass in the Mediterranean Sea region. We implemented two genome scan approaches and focused on genomic regions supported by both methods, resulting in the identification of 96 candidate genes, including nine CpG islands, which highligt potential epigenetic influences. Many of these genes and CpG islands are in linkage groups previously associated with domestication-related traits. The most significantly overrepresented molecular function was "oxidoreductase activity". Furthermore, a dense network of interactions was identified, connecting 22 of the candidate genes. Within this network, the most significantly enriched pathways and central genes were involved in "chromatin organization", highlighting another potential epigenetic mechanism. Altogether, our findings underscore the utility of interactome-assisted pathway analysis in elucidating the genomic architecture of polygenic traits and suggest that epigenetic regulation may play a crucial role in the domestication of European seabass.

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.

Probabilistic reaction norm reveals family-related variation in the association between size, condition, and sexual maturation onset in Atlantic salmon (Salmo salar)

This study investigated the relationship between the size, condition, year class, family, and sexual maturity of Atlantic salmon (Salmo salar) using data collected in an aquaculture selective breeding programme. Males that were sexually mature at 2 years of age (maiden spawn) have, on average, greater fork length and condition factor (K) at 1 year of age than their immature counterparts. For every 10-mm increase in fork length or 0.1 increase in K at 1 year of age, the odds of sexual maturity at 2 years of age increased by 1.48 or 1.22 times, respectively. Females that were sexually mature at 3 years of age (maiden spawn) have, on average, greater fork length and K at 2 years of age than their immature counterparts. For every 10-mm increase in fork length or 0.1 increase in K at 2 years of age, the odds of sexual maturity at 3 years of age increased by 1.06 or 1.44 times, respectively. The family explained 34.93% of the variation in sexual maturity among 2-year-old males that was not attributable to the average effects of fork length and K at 1 year of age and year class. The proportion of variation in sexual maturity among 3-year-old females explained by the family could not be investigated. These findings suggest that the onset of sexual maturation in Atlantic salmon is conditional on performance (with respect to energy availability) surpassing a threshold, the magnitude of which can vary between families and is determined by a genetic component. This could support the application of genetic selection to promote or inhibit the onset of sexual maturation in farmed stocks.

FishSNP: a high quality cross-species SNP database of fishes

The progress of aquaculture heavily depends on the efficient utilization of diverse genetic resources to enhance production efficiency and maximize profitability. Single nucleotide polymorphisms (SNPs) have been widely used in the study of aquaculture genomics, genetics, and breeding research since they are the most prevalent molecular markers on the genome. Currently, a large number of SNP markers from cultured fish species are scattered in individual studies, making querying complicated and data reuse problematic. We compiled relevant SNP data from literature and public databases to create a fish SNP database, FishSNP ( http://bioinfo.ihb.ac.cn/fishsnp ), and also used a unified analysis pipeline to process raw data that the author of the literature did not perform SNP calling on to obtain SNPs with high reliability. This database presently contains 45,690,243 (45 million) nonredundant SNP data for 13 fish species, with 30,288,958 (30 million) of those being high-quality SNPs. The main function of FishSNP is to search, browse, annotate and download SNPs, which provide researchers various and comprehensive associated information.

Investigating the role of genetic variation in vgll3 and six6 in the domestication of gilthead seabream (Sparus aurata Linnaeus) and European seabass (Dicentrarchus labrax Linnaeus)

Gene function conservation is crucial in molecular ecology, especially for key traits like growth and maturation in teleost fish. The vgll3 and six6 genes are known to influence age-at-maturity in Atlantic salmon, but their impact on other fish species is poorly understood. Here, we investigated the association of vgll3 and six6 in the domestication of gilthead seabream and European seabass, both undergoing selective breeding for growth-related traits in the Mediterranean. We analysed two different sets of samples using two different genotyping approaches. The first dataset comprised farmed and wild populations from Greece, genotyped for SNPs within the two genes ('gene-level genotyping'). The second dataset examined 300-600 k SNPs located in the chromosomes of the two genes, derived from a meta-analysis of a Pool-Seq experiment involving farmed and wild populations distributed widely across the Mediterranean ('chromosome-level genotyping'). The gene-level analysis revealed a statistically significant allele frequency differences between farmed and wild populations on both genes in each species. This finding was partially supported by the chromosome-level analysis, identifying highly differentiated regions may be involved in the domestication process at varying distances from the candidate genes. Noteworthy genomic features were found, such as a CpG island in gilthead seabream and novel candidate genes in European seabass, warranting further investigation. These findings support a putative role of vgll3 and six6 in the maturation and growth of gilthead seabream and European seabass, emphasizing the need for further research on their conserved function.

Atlantic salmon (Salmo salar) age at maturity is strongly affected by temperature, population and age-at-maturity genotype

Age at maturity is a key life history trait involving a trade-off between survival risk and reproductive investment, and is an important factor for population structures. In ectotherms, a warming environment may have a dramatic influence on development and life history, but this influence may differ between populations. While an increasing number of studies have examined population-dependent reactions with temperature, few have investigated this in the context of maturation timing. Atlantic salmon, a species of high conservation relevance, is a good study species for this topic as it displays considerable variation in age at maturity, of which a large proportion has been associated with a genomic region including the strong candidate gene vgll3. Until now, the effect of this gene in the context of different environments and populations has not been studied. Using a large-scale common-garden experiment, we find strong effects of temperature, population-of-origin, and vgll3 genotype on maturation in 2-year-old male Atlantic salmon (Salmo salar). With a temperature difference of 1.8°C, maturation probability was 4.8 times higher in the warm treatment than the cold treatment. This temperature effect was population-specific and was higher in the southern (60.48°N) compared to the northern (65.01°N) population. The early maturation vgll3*E allele was associated with a significantly higher maturation probability, but there was no vgll3 interaction with temperature or population. Both body condition and body mass associated with maturation. The body mass association was only present in the warm treatment. Our findings demonstrate that (i) populations can vary in their response to temperature change in terms of age at maturity, (ii) high intrinsic growth could be associated with higher thermal sensitivity for life history variation and (iii) vgll3 effects on age at maturity might be similar between populations and different thermal environments.

Leveraging transcriptome and epigenome landscapes to infer regulatory networks during the onset of sexual maturation

Background

Despite sexual development being ubiquitous to vertebrates, the molecular mechanisms underpinning this fundamental transition remain largely undocumented in many organisms. We designed a time course experiment that successfully sampled the period when Atlantic salmon commence their trajectory towards sexual maturation.

Results

Through deep RNA sequencing, we discovered key genes and pathways associated with maturation in the pituitary-ovarian axis. Analyzing DNA methylomes revealed a bias towards hypermethylation in ovary that implicated maturation-related genes. Co-analysis of DNA methylome and gene expression changes revealed chromatin remodeling genes and key transcription factors were both significantly hypermethylated and upregulated in the ovary during the onset of maturation. We also observed changes in chromatin state landscapes that were strongly correlated with fundamental remodeling of gene expression in liver. Finally, a multiomic integrated analysis revealed regulatory networks and identified hub genes including TRIM25 gene (encoding the estrogen-responsive finger protein) as a putative key regulator in the pituitary that underwent a 60-fold change in connectivity during the transition to maturation.

Conclusion

The study successfully documented transcriptome and epigenome changes that involved key genes and pathways acting in the pituitary – ovarian axis. Using a Systems Biology approach, we identified hub genes and their associated networks deemed crucial for onset of maturation. The results provide a comprehensive view of the spatiotemporal changes involved in a complex trait and opens the door to future efforts aiming to manipulate puberty in an economically important aquaculture species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08514-8.

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.

A major quantitative trait locus affecting resistance to Tilapia lake virus in farmed Nile tilapia (Oreochromis niloticus)

Enhancing host resistance to infectious disease has received increasing attention in recent years as a major goal of farm animal breeding programs. Combining field data with genomic tools can provide opportunities to understand the genetic architecture of disease resistance, leading to new opportunities for disease control. In the current study, a genome-wide association study was performed to assess resistance to the Tilapia lake virus (TiLV), one of the biggest threats affecting Nile tilapia (Oreochromis niloticus); a key aquaculture species globally. A pond outbreak of TiLV in a pedigreed population of the GIFT strain was observed, with 950 fish classified as either survivor or mortality, and genotyped using a 65 K SNP array. A significant QTL of large effect was identified on chromosome Oni22. The average mortality rate of tilapia homozygous for the resistance allele at the most significant SNP (P value = 4.51E-10) was 11%, compared to 43% for tilapia homozygous for the susceptibility allele. Several candidate genes related to host response to viral infection were identified within this QTL, including lgals17, vps52, and trim29. These results provide a rare example of a major QTL affecting a trait of major importance to a farmed animal. Genetic markers from the QTL region have potential in marker-assisted selection to improve host resistance, providing a genetic solution to an infectious disease where few other control or mitigation options currently exist.

Polygenic and major-locus contributions to sexual maturation timing in Atlantic salmon

Sexual maturation timing is a life-history trait central to the balance between mortality and reproduction. Maturation may be triggered when an underlying compound trait, called liability, exceeds a threshold. In many different species and especially fishes, this liability is approximated by growth and body condition. However, environmental vs. genetic contributions either directly or via growth and body condition to maturation timing remain unclear. Uncertainty exists also because the maturation process can reverse this causality and itself affect growth and body condition. In addition, disentangling the contributions of polygenic and major loci can be important. In many fishes, males mature before females, enabling the study of associations between male maturation and maturation-unbiased female liability traits. Using 40 Atlantic salmon families, longitudinal common-garden experimentation, and quantitative genetic analyses, we disentangled environmental from polygenic and major locus (vgll3) effects on male maturation, and sex-specific growth and condition. We detected polygenic heritabilities for maturation, growth, and body condition, and vgll3 effects on maturation and body condition but not on growth. Longitudinal patterns for sex-specific phenotypic liability, and for genetic variances and correlations between sexes suggested that early growth and condition indeed positively affected maturation initiation. However, towards spawning time, causality appeared reversed for males whereby maturation affected growth negatively and condition positively via both the environmental and genetic effects. Altogether, the results indicate that growth and condition are useful traits to study liability for maturation initiation, but only until maturation alters their expression, and that vgll3 contributes to maturation initiation via condition.

Detection of selection signatures in farmed coho salmon (Oncorhynchus kisutch) using dense genome-wide information

Animal domestication and artificial selection give rise to gradual changes at the genomic level in populations. Subsequent footprints of selection, known as selection signatures or selective sweeps, have been traced in the genomes of many animal livestock species by exploiting variation in linkage disequilibrium patterns and/or reduction of genetic diversity. Domestication of most aquatic species is recent in comparison with land animals, and salmonids are one of the most important fish species in aquaculture. Coho salmon (Oncorhynchus kisutch), cultivated primarily in Chile, has been subjected to breeding programs to improve growth, disease resistance traits, and flesh color. This study aimed to identify selection signatures that may be involved in adaptation to culture conditions and traits of productive interest. To do so, individuals of two domestic populations cultured in Chile were genotyped with 200 thousand SNPs, and analyses were conducted using iHS, XP-EHH and CLR. Several signatures of selection on different chromosomal regions were detected across both populations. Some of the identified regions under selection contained genes such anapc2, alad, chp2 and myn, which have been previously associated with body weight in Atlantic salmon, or sec24d and robo1, which have been associated with resistance to Piscirickettsia salmonis in coho salmon. Findings in our study can contribute to an integrated genome-wide map of selection signatures, to help identify the genetic mechanisms of phenotypic diversity in coho salmon.

Heterogeneous genetic basis of age at maturity in salmonid fishes

Understanding the genetic basis of repeated evolution of the same phenotype across taxa is a fundamental aim in evolutionary biology and has applications in conservation and management. However, the extent to which interspecific life-history trait polymorphisms share evolutionary pathways remains underexplored. Here, we address this gap by studying the genetic basis of a key life-history trait, age at maturity, in four species of Pacific salmonids (genus Oncorhynchus) that exhibit intra- and interspecific variation in this trait-Chinook Salmon, Coho Salmon, Sockeye Salmon, and Steelhead Trout. We tested for associations in all four species between age at maturity and two genome regions, six6 and vgll3, that are strongly associated with the same trait in Atlantic Salmon (Salmo salar). We also conducted a genome-wide association analysis in Steelhead to assess whether additional regions were associated with this trait. We found the genetic basis of age at maturity to be heterogeneous across salmonid species. Significant associations between six6 and age at maturity were observed in two of the four species, Sockeye and Steelhead, with the association in Steelhead being particularly strong in both sexes (p = 4.46 × 10^-9 after adjusting for genomic inflation). However, no significant associations were detected between age at maturity and the vgll3 genome region in any of the species, despite its strong association with the same trait in Atlantic Salmon. We discuss possible explanations for the heterogeneous nature of the genetic architecture of this key life-history trait, as well as the implications of our findings for conservation and management.

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.

Transcription Profiles of Age-at-Maturity-Associated Genes Suggest Cell Fate Commitment Regulation as a Key Factor in the Atlantic Salmon Maturation Process

Despite recent taxonomic diversification in studies linking genotype with phenotype, follow-up studies aimed at understanding the molecular processes of such genotype-phenotype associations remain rare. The age at which an individual reaches sexual maturity is an important fitness trait in many wild species. However, the molecular mechanisms regulating maturation timing processes remain obscure. A recent genome-wide association study in Atlantic salmon (Salmo salar) identified large-effect age-at-maturity-associated chromosomal regions including genes vgll3, akap11 and six6, which have roles in adipogenesis, spermatogenesis and the hypothalamic-pituitary-gonadal (HPG) axis, respectively. Here, we determine expression patterns of these genes during salmon development and their potential molecular partners and pathways. Using Nanostring transcription profiling technology, we show development- and tissue-specific mRNA expression patterns for vgll3, akap11 and six6. Correlated expression levels of vgll3 and akap11, which have adjacent chromosomal location, suggests they may have shared regulation. Further, vgll3 correlating with arhgap6 and yap1, and akap11 with lats1 and yap1 suggests that Vgll3 and Akap11 take part in actin cytoskeleton regulation. Tissue-specific expression results indicate that vgll3 and akap11 paralogs have sex-dependent expression patterns in gonads. Moreover, six6 correlating with slc38a6 and rtn1, and Hippo signaling genes suggests that Six6 could have a broader role in the HPG neuroendrocrine and cell fate commitment regulation, respectively. We conclude that Vgll3, Akap11 and Six6 may influence Atlantic salmon maturation timing via affecting adipogenesis and gametogenesis by regulating cell fate commitment and the HPG axis. These results may help to unravel general molecular mechanisms behind maturation.