Chromosome-level genome sequence of the Genetically Improved Farmed Tilapia (GIFT, Oreochromis niloticus) highlights regions of introgression with O. mossambicus.

BACKGROUND: The Nile tilapia (Oreochromis niloticus) is the third most important freshwater fish for aquaculture. Its success is directly linked to continuous breeding efforts focusing on production traits such as growth rate and weight. Among those elite strains, the Genetically Improved Farmed Tilapia (GIFT) programme initiated by WorldFish is now distributed worldwide. To accelerate the development of the GIFT strain through genomic selection, a high-quality reference genome is necessary. RESULTS: Using a combination of short (10X Genomics) and long read (PacBio HiFi, PacBio CLR) sequencing and a genetic map for the GIFT strain, we generated a chromosome level genome assembly for the GIFT. Using genomes of two closely related species (O. mossambicus, O. aureus), we characterised the extent of introgression between these species and O. niloticus that has occurred during the breeding process. Over 11 Mb of O. mossambicus genomic material could be identified within the GIFT genome, including genes associated with immunity but also with traits of interest such as growth rate. CONCLUSION: Because of the breeding history of elite strains, current reference genomes might not be the most suitable to support further studies into the GIFT strain. We generated a chromosome level assembly of the GIFT strain, characterising its mixed origins, and the potential contributions of introgressed regions to selected traits.

Exploring genetic resistance to infectious salmon anaemia virus in Atlantic salmon by genome-wide association and RNA sequencing.

BACKGROUND: Infectious Salmonid Anaemia Virus (ISAV) causes a notifiable disease that poses a large threat for Atlantic salmon (Salmo salar) aquaculture worldwide. There is no fully effective treatment or vaccine, and therefore selective breeding to increase resistance to ISAV is a promising avenue for disease prevention. Genomic selection and potentially genome editing can be applied to enhance host resistance, and these approaches benefit from improved knowledge of the genetic and functional basis of the target trait. The aim of this study was to characterise the genetic architecture of resistance to ISAV in a commercial Atlantic salmon population and study its underlying functional genomic basis using RNA Sequencing. RESULTS: A total of 2833 Atlantic salmon parr belonging to 194 families were exposed to ISAV in a cohabitation challenge in which cumulative mortality reached 63% over 55 days. A total of 1353 animals were genotyped using a 55 K SNP array, and the estimate of heritability for the trait of binary survival was 0.13-0.33 (pedigree-genomic). A genome-wide association analysis confirmed that resistance to ISAV was a polygenic trait, albeit a genomic region in chromosome Ssa13 was significantly associated with resistance and explained 3% of the genetic variance. RNA sequencing of the heart of 16 infected (7 and 14 days post infection) and 8 control fish highlighted 4927 and 2437 differentially expressed genes at 7 and 14 days post infection respectively. The complement and coagulation pathway was down-regulated in infected fish, while several metabolic pathways were up-regulated. The interferon pathway showed little evidence of up-regulation at 7 days post infection but was mildly activated at 14 days, suggesting a potential crosstalk between host and virus. Comparison of the transcriptomic response of fish with high and low breeding values for resistance highlighted TRIM25 as being up-regulated in resistant fish. CONCLUSIONS: ISAV resistance shows moderate heritability with a polygenic architecture, but a significant QTL was detected on chromosome 13. A mild up-regulation of the interferon pathway characterises the response to the virus in heart samples from this population of Atlantic salmon, and candidate genes showing differential expression between samples with high and low breeding values for resistance were identified.

Development and testing of a combined species SNP array for the European seabass (Dicentrarchus labrax) and gilthead seabream (Sparus aurata).

SNP arrays are powerful tools for high-resolution studies of the genetic basis of complex traits, facilitating both selective breeding and population genomic research. The European seabass (Dicentrarchus labrax) and the gilthead seabream (Sparus aurata) are the two most important fish species for Mediterranean aquaculture. While selective breeding programmes increasingly underpin stock supply for this industry, genomic selection is not yet widespread. Genomic selection has major potential to expedite genetic gain, particularly for traits practically impossible to measure on selection candidates, such as disease resistance and fillet characteristics. The aim of our study was to design a combined-species 60 K SNP array for European seabass and gilthead seabream, and to test its performance on farmed and wild populations from numerous locations throughout the species range. To achieve this, high coverage Illumina whole-genome sequencing of pooled samples was performed for 24 populations of European seabass and 27 populations of gilthead seabream. This resulted in a database of ~20 million SNPs per species, which were then filtered to identify high-quality variants and create the final set for the development of the 'MedFish' SNP array. The array was then tested by genotyping a subset of the discovery populations, highlighting a high conversion rate to functioning polymorphic assays on the array (92% in seabass; 89% in seabream) and repeatability (99.4-99.7%). The platform interrogates ~30 K markers in each species, includes features such as SNPs previously shown to be associated with performance traits, and is enriched for SNPs predicted to have high functional effects on proteins. The array was demonstrated to be effective at detecting population structure across a wide range of fish populations from diverse geographical origins, and to examine the extent of haplotype sharing among Mediterranean farmed fish populations. In conclusion, the new MedFish array enables efficient and accurate high-throughput genotyping for genome-wide distributed SNPs for each fish species, and will facilitate stock management, population genomics approaches, and acceleration of selective breeding through genomic selection.

Quantitative trait loci and genes associated with salmonid alphavirus load in Atlantic salmon: implications for pancreas disease resistance and tolerance.

Salmonid alphavirus infection results in pancreas disease causing severe economic losses for Atlantic salmon aquaculture. Knowledge about genes and pathways contributing to resistance is limited. A 54 K SNP panel was used to genotype 10 full-sibling families each consisting of ~ 110 offspring challenged with salmonid alphavirus subtype 3. Relative heart viral load was assessed at 4- and 10-weeks post-infection using quantitative PCR. A moderate genomic heritability of viral load at 4 weeks (0.15-0.21) and a high positive correlation with survival (0.91-0.98) were detected. Positions of QTL detected on chromosome 3 matched those for survival detected by other studies. The SNP of highest significance occurred in the 3' untranslated region of gig1, a fish-specific antiviral effector. Locus B of immunoglobulin heavy chain mapped to an area containing multiple SNPs with genome-wide association. Heart mRNA-seq comparing parr from families with high- versus low-genomic breeding value, and matching sample genotypes for SNPs, identified two eQTL for salmonid alphavirus load. Immune genes associated with trans-eQTL were numerous and spread throughout the genome. QTL regions contained several genes with known or predicted immune functions, some differentially expressed. The putative functional genes and variants identified could help improve marker-based selection for pancreas disease resistance.

Potential of genomic selection for improvement of resistance to ostreid herpesvirus in Pacific oyster (Crassostrea gigas).

In genomic selection (GS), genome-wide SNP markers are used to generate genomic estimated breeding values for selection candidates. The application of GS in shellfish looks promising and has the potential to help in dealing with one of the main issues currently affecting Pacific oyster production worldwide, which is the 'summer mortality syndrome'. This causes periodic mass mortality in farms worldwide and has mainly been attributed to a specific variant of the ostreid herpesvirus (OsHV-1). In the current study, we evaluated the potential of genomic selection for host resistance to OsHV-1 in Pacific oysters, and compared it with pedigree-based approaches. An OsHV-1 disease challenge was performed using an immersion-based virus exposure treatment for oysters for 7 days. A total of 768 samples were genotyped using the medium-density SNP array for oysters. A GWAS was performed for the survival trait using a GBLUP approach in blupf90 software. Heritability ranged from 0.25 ± 0.05 to 0.37 ± 0.05 (mean ± SE) based on pedigree and genomic information respectively. Genomic prediction was more accurate than pedigree prediction, and SNP density reduction had little impact on prediction accuracy until marker densities dropped below approximately 500 SNPs. This demonstrates the potential for GS in Pacific oyster breeding programmes, and importantly, demonstrates that a low number of SNPs might suffice to obtain accurate genomic estimated breeding values, thus potentially making the implementation of GS more cost effective.

Assessment of genetic diversity and population structure in cultured Australian Pacific oysters.

The recent development of Pacific oyster (Crassostrea gigas) SNP genotyping arrays has allowed detailed characterisation of genetic diversity and population structure within and between oyster populations. It also raises the potential of harnessing genomic selection for genetic improvement in oyster breeding programmes. The aim of this study was to characterise a breeding population of Australian oysters through genotyping and analysis of 18 027 SNPs, followed by comparison with genotypes of oyster sampled from Europe and Asia. This revealed that the Australian populations had similar population diversity (HE ) to oysters from New Zealand, the British Isles, France and Japan. Population divergence was assessed using PCA of genetic distance and revealed that Australian oysters were distinct from all other populations tested. Australian Pacific oysters originate from planned introductions sourced from three Japanese populations. Approximately 95% of these introductions were from geographically, and potentially genetically, distinct populations from the Nagasaki oysters assessed in this study. Finally, in preparation for the application of genomic selection in oyster breeding programmes, the strength of LD was evaluated and subsets of loci were tested for their ability to accurately infer relationships. Weak LD was observed on average; however, SNP subsets were shown to accurately reconstitute a genomic relationship matrix constructed using all loci. This suggests that low-density SNP panels may have utility in the Australian population tested, and the findings represent an important first step towards the design and implementation of genomic approaches for applied breeding in Pacific oysters.

Atlantic salmon (Salmo salar L.) genetics in the 21st century: taking leaps forward in aquaculture and biological understanding.

Atlantic salmon (Salmo salar L.) is among the most iconic and economically important fish species and was the first member of Salmonidae to have a high-quality reference genome assembly published. Advances in genomics have become increasingly central to the genetic improvement of farmed Atlantic salmon as well as conservation of wild salmon stocks. The salmon genome has also been pivotal in shaping our understanding of the evolutionary and functional consequences arising from an ancestral whole-genome duplication event characterising all Salmonidae members. Here, we provide a review of the current status of Atlantic salmon genetics and genomics, focussed on progress made from genome-wide research aimed at improving aquaculture production and enhancing understanding of salmonid ecology, physiology and evolution. We present our views on the future direction of salmon genomics, including the role of emerging technologies (e.g. genome editing) in elucidating genetic features that underpin functional variation in traits of commercial and evolutionary importance.

Mapping and validation of a major QTL affecting resistance to pancreas disease (salmonid alphavirus) in Atlantic salmon (Salmo salar).

Pancreas disease (PD), caused by a salmonid alphavirus (SAV), has a large negative economic and animal welfare impact on Atlantic salmon aquaculture. Evidence for genetic variation in host resistance to this disease has been reported, suggesting that selective breeding may potentially form an important component of disease control. The aim of this study was to explore the genetic architecture of resistance to PD, using survival data collected from two unrelated populations of Atlantic salmon; one challenged with SAV as fry in freshwater (POP 1) and one challenged with SAV as post-smolts in sea water (POP 2). Analyses of the binary survival data revealed a moderate-to-high heritability for host resistance to PD in both populations (fry POP 1 h(2)~0.5; post-smolt POP 2 h(2)~0.4). Subsets of both populations were genotyped for single nucleotide polymorphism markers, and six putative resistance quantitative trait loci (QTL) were identified. One of these QTL was mapped to the same location on chromosome 3 in both populations, reaching chromosome-wide significance in both the sire- and dam-based analyses in POP 1, and genome-wide significance in a combined analysis in POP 2. This independently verified QTL explains a significant proportion of host genetic variation in resistance to PD in both populations, suggesting a common underlying mechanism for genetic resistance across lifecycle stages. Markers associated with this QTL are being incorporated into selective breeding programs to improve PD resistance.

Single nucleotide polymorphisms in the insulin-like growth factor 1 (IGF1) gene are associated with growth-related traits in farmed Atlantic salmon.

Understanding the genetic basis of variation in traits related to growth and fillet quality in Atlantic salmon is of importance to the aquaculture industry. Several growth-related QTL have been identified via the application of genetic markers. The IGF1 gene is considered a highly conserved and crucial growth-regulating gene in salmonid species. However, the association between polymorphisms in the IGF1 gene and growth-related traits in Atlantic salmon is unknown. Therefore, in this study, regions of the Atlantic salmon IGF1 gene were sequenced, aligned and compared across individuals. Three SNPs were identified in the putative promoter (SNP1, g.5763G>T; GenBank no. AGKD01012745), intron 1 (SNP2, g.7292C>T; GenBank no. AGKD01012745) and intron 3 (SNP3, g.4671A>C; GenBank no. AGKD01133398) regions respectively. These SNPs were genotyped in a population of 4800 commercial Atlantic salmon with data on several weight and fillet traits measured at harvest (at approximately 3 years of age). In a mixed model, association analysis of individual SNPs, SNP1 and SNP3 were both significantly associated with several weight traits (P < 0.05). The estimated additive effect on overall harvest weight was approximately 35 and 110 g for SNPs 1 and 3 respectively. A haplotype analysis confirmed the association between genetic variation in the IGF1 gene with overall body weight (P < 0.05) and fillet component traits (P < 0.05). Our findings suggest the identified nucleotide polymorphisms of the IGF1 gene may either affect farmed Atlantic salmon growth directly or be in population-wide linkage disequilibrium with causal variation, highlighting their possible utility as candidates for marker-assisted selection in the aquaculture industry.

Characterization of OAR1 and OAR18 QTL associated with muscle depth in British commercial terminal sire sheep.

This study aimed at verifying previously identified QTL affecting growth and carcass traits on ovine chromosome 18 (OAR18) in Texel sheep (n = 1844), and on OAR1 in Charollais (n = 851) and Suffolk (n = 998) sheep. The QTL were investigated using regression and variance component mapping (VCA) of body weight, muscle and fat depth measurements. In addition, the mode of inheritance of the Texel OAR18 QTL was explored, using data from 4376 Texel sheep, fitting VCA models testing for additive and imprinting effects. We also simulated a 480-sheep population with different QTL imprinting models and various available levels of marker information to understand the behaviour of the VCA results under different assumed genetic models. In summary, the previously identified QTL were successfully verified using both interval mapping and VCA in the three breeds. We propose a polar overdominance mode of inheritance for the OAR18 QTL in Texel sheep, and we present methods to dissect the QTL mode of inheritance, using the Texel OAR18 QTL as an example.

Segregation of infectious pancreatic necrosis resistance QTL in the early life cycle of Atlantic Salmon (Salmo salar).

In a previous study, three significant quantitative trait loci (QTL) associated with resistance to Infectious Pancreatic Necrosis (IPN) disease were identified by analysing challenge data from one sub-population of Landcatch Atlantic salmon (Salmo salar) smolt. While these QTL were shown to affect the resistance in seawater, their effect in freshwater was unknown. This study investigates the effect of these QTL on IPN resistance in salmon fry in freshwater. Twenty families with intermediate levels of IPN mortality were analysed from a freshwater challenge trial undertaken on a different sup-population of LNS salmon to that studied previously. Only the QTL from linkage group 21 (LG21) appeared to have a significant and large effect on resistance in freshwater; the same QTL was found to have the largest effect in seawater in the previous study. Variance component analysis showed a high heritability for the QTL: 0.45±0.07 on the liability scale and 0.25±0.05 on the observed scale. In a family where both parents were segregating for the QTL, there was a 0% vs. 100% mortality in homozygous offspring for resistant and susceptible QTL alleles. The finding that the same QTL has major effect in both freshwater and seawater has important practical implications, as this will allow the improvement of resistance in both phases through marker assisted selection by targeting this QTL. Moreover, the segregation of the LG21 QTL in a different sub-population gives further evidence of its association with IPN-resistance.

The susceptibility of Atlantic salmon fry to freshwater infectious pancreatic necrosis is largely explained by a major QTL.

Infectious pancreatic necrosis (IPN) is a viral disease with a significant negative impact on the global aquaculture of Atlantic salmon. IPN outbreaks can occur during specific windows of both the freshwater and seawater stages of the salmon life cycle. Previous research has shown that a proportion of the variation seen in resistance to IPN is because of host genetics, and we have shown that major quantitative trait loci (QTL) affect IPN resistance at the seawater stage of production. In the current study, we completed a large freshwater IPN challenge experiment to allow us to undertake a thorough investigation of the genetic basis of resistance to IPN in salmon fry, with a focus on previously identified QTL regions. The heritability of freshwater IPN resistance was estimated to be 0.26 on the observed scale and 0.55 on the underlying scale. Our results suggest that a single QTL on linkage group 21 explains almost all the genetic variation in IPN mortality under our experimental conditions. A striking contrast in mortality is seen between fry classified as homozygous susceptible versus homozygous resistant, with QTL-resistant fish showing virtually complete resistance to IPN mortality. The findings highlight the importance of the major QTL in the genetic regulation of IPN resistance across distinct physiological lifecycle stages, environmental conditions and viral isolates. These results have clear scientific and practical implications for the control of IPN.

Detection of QTL affecting harvest traits in a commercial Atlantic salmon population.

Genetic variation in performance and quality traits measured at harvest has previously been demonstrated in Atlantic salmon aquaculture populations. To map major loci underlying this variation, we utilized data from 10 families from a commercial breeding programme. Significant QTL were detected affecting harvest weight and length traits on linkage group 1, and affecting waste weight on linkage group 5. In total, 11 of the 29 linkage groups examined showed at least suggestive evidence for a QTL. These data suggest that major loci affecting economically important harvest characteristics are segregating in commercial salmon populations.

Detection and confirmation of a major QTL affecting resistance to infectious pancreatic necrosis (IPN) in Atlantic salmon (Salmo salar).

Infectious pancreatic necrosis (IPN) is a viral disease currently presenting a major problem to the aquaculture of Atlantic salmon (Salmon salar), during both the freshwater and seawater stages of production. Genetic variation in resistance to IPN has previously been demonstrated and the purpose of this study was to determine whether this variation includes loci of major effect. The initial QTL detection methodology utilized the limited recombination seen in male salmon to detect QTL in ten large full-sib families, using a genome-wide scan of two to three markers per linkage group. QTL were then positioned by adding additional markers to the significant linkage groups in a female-based analysis. The most significant QTL was mapped to LG 21, and further confirmation of the LG 21 QTL is provided in an analysis of the QTL flanking markers in an additional nine full-sib families from the same population. The size of QTL effect is such that the QTL flanking markers can be immediately applied in marker-assisted selection programmes to improve the resistance of salmon populations to IPN, thus reducing mortality due to the disease.

The cholecystokinin type A receptor g.179A>G polymorphism affects feeding rate.

A polymorphism within the 5' untranslated region of the cholecystokinin type A receptor (CCKAR) gene has been shown to affect feed intake and growth in commercial pig lines. To further investigate the phenotype of animals carrying alternative alleles at this polymorphism, we genotyped animals from a distinct segregating commercial line and an experimental cross F(2) population, both with electronically recorded feeding pattern data. The data indicate that the daily feed intake increasing effect of the DQ496228:g.179G allele is mediated through a faster rate of feed intake, without evidence for an effect on other feeding behaviour traits.

A polymorphism in the 5'-untranslated region of the porcine cholecystokinin type a receptor gene affects feed intake and growth.

The location and utilization of quantitative trait loci (QTL) and candidate genes with significant effects on economically important traits are becoming increasingly important in livestock breeding programs. The porcine cholecystokinin type A receptor (CCKAR) is a candidate gene for performance traits, due to its known role in the physiological control of feed intake, satiety, and obesity. We investigated the association of CCKAR polymorphisms with feeding, growth, and efficiency traits in an F2 population derived from a cross between Meishan and Large White founder animals and in lines of Large White pigs that had been divergently selected on the basis of lean growth efficiency traits. In the F2 population, CCKAR genotype was significantly associated with daily feed intake and average daily gain. The effects of the polymorphisms were then assessed in a larger-scale analysis of segregating commercial lines. A newly discovered single-nucleotide polymorphism (SNP) within the 5'-untranslated region (5'-UTR) had highly significant effects on feed intake, average daily gain, and days to 110 kg, which were not seen for a previously reported SNP within the CCKAR gene. Furthermore, we provide evidence that the novel SNP disrupts the binding of the YY1 transcription factor, which raises the possibility that it is the causal variant. The 5'-UTR SNP could be utilized as a molecular genetic test for increased feed intake, faster lean growth, and reduced days to market weight in segregating commercial lines.

A melanocortin-4 receptor (MC4R) polymorphism is associated with performance traits in divergently selected Large White pig populations.

The melanocortin-4 receptor (MC4R) has a vital role in the control of energy balance and the genetic basis of obesity. A polymorphism, which results in the replacement of aspartic acid with asparagine at position 298 of the porcine MC4R gene, within the seventh transmembrane domain, has previously been described. In the current study, allele frequencies for this Asp298Asn polymorphism were investigated in lines of Large White pigs which had been divergently selected for seven generations based on lean food conversion (LFC), lean growth with ad libitum feeding (LGA), lean growth with restricted feeding (LGS) and daily feed intake (DFI). The association of the Asp298Asn polymorphism with performance traits in these lines was assessed. The frequency of Asp298 was higher (P < 0.001) in the LFC high line (0.48) than the low line (0.00), while the frequency of Asn298 was higher (P < 0.01) in the LGA high line (0.22) than the low line (0.04). When analysed across all lines, the Asp298Asn polymorphism was significantly associated with ultrasonic backfat depth, average daily gain and daily feed intake (P < 0.05). Asp298 homozygous animals had mean values of 13.3 mm, 733 g and 1933 g for backfat, average daily gain and daily feed intake respectively, compared with 14.7 mm, 805 g and 2098 g for Asn298 homozygotes. Therefore, the data support a role for the MC4R Asp298Asn polymorphism in the genetic basis of economically important traits in Large White pigs.