Mapping quantitative trait loci for omega-3 fatty acids in Asian seabass

Candidate genes associated with fatty acid compositions in north American Atlantic salmon (Salmo salar)

Fatty acids are a requirement for normal development, however, since humans are unable to de novo produce essential fatty acids, they must be obtained from diet. Atlantic salmon is a major dietary source of nutritious and digestible fatty acids. Here, we set out to uncover the genomic basis of individual fatty acids and indices (saturated fatty acids, monounsaturated fatty acids, polyunsaturated fatty acids, n-3, and n-6) in 208 North American Atlantic salmon, to understand selection potential toward increasing relative quantities of essential fatty acids and to identify candidate genes for future research. Total n-6 (pro-inflammatory) was higher than total n-3 (anti-inflammatory) fatty acids with a ratio of 1 : 1.31 (n-3 : n-6). Heritability of fatty acids ranged from 0 to 0.99, however, most fatty acids and indices had moderate to high heritabilities (ranged from 0.20 to 0.88), implying that selection for improvement of traits could be possible. We found the same significant markers on chromosome 23 (based on false discovery rate thresholds of 2.0e-6 and suggestive significant thresholds of 2.0e-5 in Manhattan plots) in four fatty acids (γ-linoleic acid, stearidonic acid, dihimo-γ-linolenic acid, and eicosatrienoic acid), where three genes (sin3b, acbd6, and fads2) are known to be involved in lipid metabolism. These genes, fads2 in particular, would all make ideal candidates for future functional studies. In addition, there were four fatty acids with loci over the suggestive significant threshold with a variety of markers on different chromosomes (lauric acid, stearic acid, eicosatetraenoic acid (ETA), and docosadienoic acid), with associated genes that had relevant functions to fatty acids or adipose cells in general.

Age-Dependent Growth-Related QTL Variations in Pacific Abalone, Haliotis discus hannai

Pacific abalone is a high-value, commercially important marine invertebrate. It shows low growth as well as individual and yearly growth variation in aquaculture. Marker-assisted selection breeding could potentially resolve the problem of low and variable growth and increase genetic gain. Expression of quantitative trait loci (QTLs) for growth-related traits, viz., body weight, shell length, and shell width were analyzed at the first, second, and third year of age using an F1 cross population. A total of 37 chromosome-wide QTLs were identified in linkage groups 01, 02, 03, 04, 06, 07, 08, 10, 11, 12, and 13 at different ages. None of the QTLs detected at any one age were expressed in all three age groups. This result suggests that growth-related traits at different ages are influenced by different QTLs in each year. However, multiple-trait QTLs (where one QTL affects all three traits) were detected each year that are also age-specific. Eleven multiple-trait QTLs were detected at different ages: two QTLs in the first year; two QTLs in the second year; and seven QTLs in the third year. As abalone hatcheries use three-year-old abalone for breeding, QTL-linked markers that were detected at the third year of age could potentially be used in marker-assisted selection breeding programs.

Utilization of marine diatom Thalassiosira weissflogii as a feed additive in seawater-tolerant Nile tilapia (Oreochromis niloticus, Linnaeus 1758) strain

The low omega-3 content of tilapia flesh, when compared to marine fish, affects its marketability. In marine animals, the highly unsaturated fatty acids (HUFAs) can be linked to the oil produced by marine diatoms. Among the marine diatoms, the genus Thalassiosira is known to exhibit high content of HUFAs such as eicosapentaenoic acid (EPA) and docosahexaenoic acid (DHA). Thus, in this study, the use of marine diatom Thalassiosira weissflogii as a dietary additive in the seawater-tolerant Nile Tilapia strain was evaluated. One hundred ninety-two, 1.40 ± 0.05g seawater tilapia were randomly allocated into 4 treatment groups in 4 replicates. The first treatment group was fed with a control diet (D0), without the diatoms while treatments 1, 2, and 3 were each fed with diets supplemented with T. weissflogii paste at 2.55% (D1), 6% (D2), and 12% (D3), respectively for 60 days. The diets were isonitrogenous, isolipodic and the omega-3 and omega-6 requirements were satisfied. Results demonstrated that D1 had the highest percent weight gain among treatments. Although not significantly different, other parameters such as percent survival, specific growth rate (SGR), protein efficiency ratio (PER), feed conversion ratio (FCR), and feed intake had desirable results in D1. The proximate composition of seawater tilapia showed that % crude protein was highest in D0 but % crude lipid was highest in D1. The fatty acid composition of tilapia in D1 had the highest omega-3 content at 9.29 mg/g tissue and also had the highest n3:n6 at 2.19. Muscle growth-related genes (MyoD and MYG) were up-regulated while liver genes involved in long-chain polyunsaturated fatty acid synthesis (oni-fads2 and elvol5) were down-regulated in D1 as compared to D0. Feeding the diatom-supplemented diet to tilapia had no significant effects on hepatic cells and intestinal morphology. The results suggested that a 2.55% supplementation dose of T. weissflogii could promote growth and enhance the tissue content of omega-3 fatty acids of the seawater strain Oreochromis niloticus.

Genetic architecture and genomic selection of fatty acid composition predicted by Raman spectroscopy in rainbow trout

Background

In response to major challenges regarding the supply and sustainability of marine ingredients in aquafeeds, the aquaculture industry has made a large-scale shift toward plant-based substitutions for fish oil and fish meal. But, this also led to lower levels of healthful n−3 long-chain polyunsaturated fatty acids (PUFAs)—especially eicosapentaenoic (EPA) and docosahexaenoic (DHA) acids—in flesh. One potential solution is to select fish with better abilities to retain or synthesise PUFAs, to increase the efficiency of aquaculture and promote the production of healthier fish products. To this end, we aimed i) to estimate the genetic variability in fatty acid (FA) composition in visceral fat quantified by Raman spectroscopy, with respect to both individual FAs and groups under a feeding regime with limited n-3 PUFAs; ii) to study the genetic and phenotypic correlations between FAs and processing yields- and fat-related traits; iii) to detect QTLs associated with FA composition and identify candidate genes; and iv) to assess the efficiency of genomic selection compared to pedigree-based BLUP selection.

Results

Proportions of the various FAs in fish were indirectly estimated using Raman scattering spectroscopy. Fish were genotyped using the 57 K SNP Axiom™ Trout Genotyping Array. Following quality control, the final analysis contained 29,652 SNPs from 1382 fish. Heritability estimates for traits ranged from 0.03 ± 0.03 (n-3 PUFAs) to 0.24 ± 0.05 (n-6 PUFAs), confirming the potential for genomic selection. n-3 PUFAs are positively correlated to a decrease in fat deposition in the fillet and in the viscera but negatively correlated to body weight. This highlights the potential interest to combine selection on FA and against fat deposition to improve nutritional merit of aquaculture products. Several QTLs were identified for FA composition, containing multiple candidate genes with indirect links to FA metabolism. In particular, one region on Omy1 was associated with n-6 PUFAs, monounsaturated FAs, linoleic acid, and EPA, while a region on Omy7 had effects on n-6 PUFAs, EPA, and linoleic acid. When we compared the effectiveness of breeding programmes based on genomic selection (using a reference population of 1000 individuals related to selection candidates) or on pedigree-based selection, we found that the former yielded increases in selection accuracy of 12 to 120% depending on the FA trait.

Conclusion

This study reveals the polygenic genetic architecture for FA composition in rainbow trout and confirms that genomic selection has potential to improve EPA and DHA proportions in aquaculture species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08062-7.

The evolutionary ecology of fatty-acid variation: Implications for consumer adaptation and diversification

The nutritional diversity of resources can affect the adaptive evolution of consumer metabolism and consumer diversification. The omega-3 long-chain polyunsaturated fatty acids eicosapentaenoic acid (EPA; 20:5n-3) and docosahexaenoic acid (DHA; 22:6n-3) have a high potential to affect consumer fitness, through their widespread effects on reproduction, growth and survival. However, few studies consider the evolution of fatty acid metabolism within an ecological context. In this review, we first document the extensive diversity in both primary producer and consumer fatty acid distributions amongst major ecosystems, between habitats and amongst species within habitats. We highlight some of the key nutritional contrasts that can shape behavioural and/or metabolic adaptation in consumers, discussing how consumers can evolve in response to the spatial, seasonal and community-level variation of resource quality. We propose a hierarchical trait-based approach for studying the evolution of consumers' metabolic networks and review the evolutionary genetic mechanisms underpinning consumer adaptation to EPA and DHA distributions. In doing so, we consider how the metabolic traits of consumers are hierarchically structured, from cell membrane function to maternal investment, and have strongly environment-dependent expression. Finally, we conclude with an outlook on how studying the metabolic adaptation of consumers within the context of nutritional landscapes can open up new opportunities for understanding evolutionary diversification.

Genome-Wide Marker Analysis for Traits of Economic Importance in Asian Seabass Lates calcarifer

To date, it is not known whether animal breeding values in Asian seabass (Lates calcarifer) can be estimated using single nucleotide polymorphisms (SNPs) generated from new high-throughput genotyping by sequencing platforms. The principal aim of the present study was to assess the genomic prediction accuracy for growth traits, survival, cannibalism, and disease resistance against Streptococcus iniae in this species L. calcarifer. Additionally, this study attempted to identify markers associated with the five traits studied as well as to understand if the genotype data can be used to estimate genetic parameters for these complex traits. The genomic best linear unbiased prediction (gBLUP) method was used to analyze 11,084 SNPs and showed that the prediction accuracies for growth traits (weight and length) were high (0.67–0.75). By contrast, these estimates for survival were low (0.25). Multi-locus mixed model analyses identified four SNPs significantly associated with body weight (p < 5 × 10−8 or −log10 p ≥ 5). There were, however, no significant associations detected for other traits. Similarly, the SNP heritability was moderate, while the estimates for other traits were approximated to zero and not significant. Genetic correlations between body weight and standard length were close to unity. Collectively, the results obtained from this study suggest that genotyping by sequencing platforms can provide informative DNA markers to conduct genome-wide association analysis, estimation of genetic parameters, and evaluation of genomic prediction accuracy for complex traits in Asian seabass.

A High-Density Genetic Linkage Map and QTL Mapping for Sex and Growth-Related Traits of Large-Scale Loach (Paramisgurnus dabryanus)

Large-scale loach (Paramisgurnus dabryanus) is a commercially important species in East Asia; however, the cultured population that exhibited degradation of germplasm resource cannot meet the market needs, and the genome resources for P. dabryanus are still lacking. In this study, the first high-density genetic map of P. dabryanus was constructed using 15,830 SNP markers based on high-throughput sequencing with an improved SLAF-seq strategy. The quantitative trait locus (QTL) mapping for sex, growth, and morphology traits was performed for the first time. The genetic map spanned 4,657.64 cM in length with an average inter-marker distance of 0.30 cM. QTL mapping and association analysis identified eight QTLs of growth traits, nine QTLs of morphology traits, and five QTLs of sex-related traits, respectively. Interestingly, the most significant QTLs for almost all the traits were concentrated on the same linkage group LG11. Seven candidate markers and 12 potentially key genes, which were associated with sex determination and growth, were identified within the overlapped QTL regions on LG11. Further, the first genome survey analysis of P. dabryanus was performed which represents the first step toward fully decoding the P. dabryanus genome. The genome scaffolds were anchored to the high-density linkage map, spanning 960.27 Mb of P. dabryanus reference genome. The collinearity analysis revealed a high level of collinearity between the genetic map and the reference genome of P. dabryanus. Moreover, a certain degree of homology was observed between large-scale loach and zebrafish using comparative genomic analysis. The constructed high-density genetic map was an important basis for QTL fine mapping, genome assembly, and genome comparison. The present study will provide a valuable resource for future marker-assisted breeding, and further genetic and genomic researches in P. dabryanus.

Mapping and Validating QTL for Fatty Acid Compositions and Growth Traits in Asian Seabass

Asian seabass is an important food fish species. While improving growth, increasing the nutritional value is important, omega-3 fatty acids are indispensable to human health. Identifying and validating DNA markers associated with traits is the first step towards marker-assisted selection (MAS). We quantified 13 different fatty acids and three growth traits in 213 F₂ Asian seabass from a family at the age 270 days post hatch, and screened QTL for these traits. The content of total fatty acids in 100 g flesh was 2.57 ± 0.80 g, while the proportions of docosahexaenoic acid (DHA) and eicosapentaenoic acid (EPA) were 16.96 ± 2.20% and 5.42 ± 0.90%, respectively. A linkage map with 2424 SNPs was constructed and used for QTL mapping. For fatty acid compositions, 14 significant QTL were identified on three linkage groups (LG5, LG11 and LG14), with phenotypic variance explained (PVE) from 12.8 to 24.6%. Thirty-nine suggestive QTL were detected on 16 LGs. Two significant QTL for EPA were identified on LG5 and LG14, with PVE of 15.2% and 15.1%, respectively. No significant QTL was identified for DHA. For growth traits, six significant and 13 suggestive QTL were identified on two and seven LGs, respectively. Only a few significant QTL for fatty acids overlapped with previously mapped QTL for these traits, suggesting that most QTL detected in a family are family-specific and could only be used in MAS in the family per se. To facilitate population-wide molecular breeding, more powerful methods (e.g. GWAS) should be used to identify SNPs for genomic selection.

Genome-wide association analysis of nutrient traits in the oyster Crassostrea gigas: genetic effect and interaction network

Background

Oyster is rich in glycogen and free amino acids and is called “the milk of sea”. To understand the main genetic effects of these traits and the genetic networks underlying their correlation, we have conducted the whole genome resequencing with 427 oysters collected from the world-wide scale.

Results

After association analysis, 168 clustered significant single nucleotide polymorphism (SNP) loci were identified for glycogen content and 17 SNPs were verified with 288 oyster individuals in another wide populations. These were the most important candidate loci for oyster breeding. Among 24 genes in the 100-kb regions of the leading SNP loci, cytochrome P450 17A1 (CYP17A1) contained a non-synonymous SNP and displayed higher expressions in high glycogen content individuals. This might enhance the gluconeogenesis process by the transcriptionally regulating the expression of phosphoenolpyruvate carboxykinase (PEPCK) and glucose 6-phosphatase (G6Pase). Also, for amino acids content, 417 clustered significant SNPs were identified. After genetic network analysis, three node SNP regions were identified to be associated with glycogen, protein, and Asp content, which might explain their significant correlation.

Conclusion

Overall, this study provides insights into the genetic correlation among complex traits, which will facilitate future oyster functional studies and breeding through molecular design.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5971-z) contains supplementary material, which is available to authorized users.

GWAS Analysis Indicated Importance of NF-κB Signaling Pathway in Host Resistance Against Motile Aeromonas Septicemia Disease in Catfish

Motile Aeromonas septicemia (MAS) disease caused by a bacterial pathogen, Aeromonas hydrophila, is an emerging but severe disease of catfish. Genetic enhancement of disease resistance is considered to be effective to control the disease. To provide an insight into the genomic basis of MAS disease resistance, in this study, we conducted a genome-wide association study (GWAS) to identify quantitative trait loci (QTL). A total of 1820 interspecific backcross catfish of 7 families were challenged with A. hydrophila, and 382 phenotypic extremes were selected for genotyping with the catfish 690 K SNP arrays. Three QTL on linkage group (LG) 2, 26 and 29 were identified to be significantly associated with MAS resistance. Within these regions, a total of 24 genes had known functions in immunity, 10 of which were involved in NF-κB signaling pathway, suggesting the importance of NF-κB signaling pathway in MAS resistance. In addition, three suggestively significant QTL were identified on LG 11, 17, and 20. The limited numbers of QTL involved in MAS resistance suggests that marker-assisted selection may be a viable approach for catfish breeding.

Genomic, Transcriptomic, and Epigenomic Features Differentiate Genes That Are Relevant for Muscular Polyunsaturated Fatty Acids in the Common Carp

Polyunsaturated fatty acids (PUFAs) are a set of important nutrients that mainly include arachidonic acid (ARA4), docosahexaenoic acid (DHA), eicosapentaenoic acid (EPA), and α-linolenic acid (ALA). Recently, fish-derived PUFAs have been associated with cardiovascular health, fetal development, and improvement of brain functions. Studies have shown that fish muscular tissues are rich in PUFAs, which are influenced by various factors, including genetic variations, regulatory profiles, and methylation status of desaturase genes during fatty acid desaturation and elongation processes. However, the genetic mechanism and the pathways involved in fatty acid metabolism in fishes remain unclear. The overall aim of this study was to assess differences in gene expression responses among fishes with different fatty acid levels. To achieve this goal, we conducted genome-wide association analysis (GWAS) using a 250K SNP array in a population of 203 samples of common carp (Cyprinus carpio) and identified nine SNPs and 15 genes associated with muscular PUFA content. Then, RNA-Seq and whole genome bisulfite sequencing (WGBS) of different groups with high and low EPA, DHA, ARA4, and ALA contents in muscle, liver and brain tissues were conducted, resulting in 6,750 differentially expressed genes and 5,631 genes with differentially methylated promoters. Gene ontology and KEGG pathway enrichment analyses of RNA-Seq and WGBS results identified enriched pathways for fatty acid metabolism, which included the adipocytokine signaling pathway, ARA4 and linoleic acid metabolism pathway, and insulin signaling pathway. Integrated analysis indicated significant correlations between gene expression and methylation status among groups with high and low PUFA contents in muscular tissues. Taken together, these multi-level results uncovered candidate genes and pathways that are associated with fatty acid metabolism and paved the way for further genomic selection and carp breeding for PUFA traits.

The FTO Gene Is Associated with Growth and Omega-3/-6 Ratio in Asian Seabass

Polymorphisms in the FTO gene are associated with obesity and body mass index in humans and livestock. Little information of whether FTO plays an important role in aquaculture fish species is available. We cloned and characterized the FTO gene in an economically important food fish species: Asian seabass (Lates calcarifer). The full-length cDNA of the gene is 3679 bp, containing an ORF of 1935 bp encoding 644 amino acids, a 216 bp 5' UTR and a 1538 bp 3' UTR. The gene consisted of nine exons and eight introns and was 117,679 bp in length. Phylogenetic analysis revealed that the gene in Asian seabass was closely related to those of Japanese flounder and Nile tilapia. Analysis of its expressions using qRT-PCR showed that it was expressed ubiquitously, but was higher in the liver, stomach and intestine. Comparative analysis of the genomic sequences of part of intron 1 of the gene among 10 unrelated individuals identified two SNPs. Analysis of associations between SNPs and traits (i.e. growth, oil content, omega-3 and -6 contents) in an F₂ family demonstrated that the two SNPs were significantly associated with growth, oil content, omega-3 content and omega-3/-6 ratio. Altogether, our data suggest that the gene or/and its linked genes play an important role in growth and fatty acid synthesis, and that the SNPs associated with traits may be used as markers for selecting quicker growth and higher omega-3/-6 ratio at the fingerling stage.

Mapping QTL for Omega-3 Content in Hybrid Saline Tilapia

Tilapia is one of most important foodfish species. The low omega-3 to omega-6 fatty acid ratio in freshwater tilapia meat is disadvantageous for human health. Increasing omega-3 content is an important task in breeding to increase the nutritional value of tilapia. However, conventional breeding to increase omega-3 content is difficult and slow. To accelerate the increase of omega-3 through marker-assisted selection (MAS), we conducted QTL mapping for fatty acid contents and profiles in a F₂ family of saline tilapia generated by crossing red tilapia and Mozambique tilapia. The total omega-3 content in F₂ hybrid tilapia was 2.5 ± 1.0 mg/g, higher than that (2.00 mg/g) in freshwater tilapia. Genotyping by sequencing (GBS) technology was used to discover and genotype SNP markers, and microsatellites were also genotyped. We constructed a linkage map with 784 markers (151 microsatellites and 633 SNPs). The linkage map was 2076.7 cM long and consisted of 22 linkage groups. Significant and suggestive QTL for total lipid content were mapped on six linkage groups (LG3, -4, -6, -8, -13, and -15) and explained 5.8-8.3% of the phenotypic variance. QTL for omega-3 fatty acids were located on four LGs (LG11, -18, -19, and -20) and explained 5.0 to 7.5% of the phenotypic variance. Our data suggest that the total lipid and omega-3 fatty acid content were determined by multiple genes in tilapia. The markers flanking the QTL for omega-3 fatty acids can be used in MAS to accelerate the genetic improvements of these traits in salt-tolerant tilapia.

Construction of the First High-Density Genetic Linkage Map and Analysis of Quantitative Trait Loci for Growth-Related Traits in Sinonovacula constricta

The razor clam (Sinonovacula constricta) is an important aquaculture species, for which a high-density genetic linkage map would play an important role in marker-assisted selection (MAS). In this study, we constructed a high-density genetic map and detected quantitative trait loci (QTLs) for Sinonovacula constricta with an F₁ cross population by using the specific locus amplified fragment sequencing (SLAF-seq) method. A total of 315,553 SLAF markers out of 467.71 Mreads were developed. The final linkage map was composed of 7516 SLAFs (156.60-fold in the parents and 20.80-fold in each F₁ population on average). The total distance of the linkage map was 2383.85 cM, covering 19 linkage groups with an average inter-marker distance of 0.32 cM. The proportion of gaps less than 5.0 cM was on average 96.90%. A total of 16 suggestive QTLs for five growth-related traits (five QTLs for shell height, six QTLs for shell length, three QTLs for shell width, one QTL for total body weight, and one QTL for soft body weight) were identified. These QTLs were distributed on five linkage groups, and the regions showed overlapping on LG9 and LG13. In conclusion, the high-density genetic map and QTLs for S. constricta provide a valuable genetic resource and a basis for MAS.

Genome-Wide Association Study Identifies Loci Associated with Resistance to Viral Nervous Necrosis Disease in Asian Seabass

Viral nervous necrosis disease (VNN), caused by nervous necrosis virus (NNV), is one major threat to mariculture. Identifying loci and understanding the mechanisms associated with resistance to VNN are important in selective breeding programs. We performed a genome-wide association study (GWAS) using genotyping-by-sequencing (GBS) to study the genomic architecture of resistance to NNV infection in Asian seabass. We genotyped 986 individuals from 43 families produced by 15 founders with 44498 bi-allelic genetic variants using GBS. The GWAS identified three genome-wide significant loci on chromosomes 16, 19, and 20, respectively, and six suggestive loci on chromosomes 1, 8, 14, 15, 21, and 24, respectively, associated with resistance to NNV infection measured as binary and quantitative traits. Using the 500 most significant markers in combination with a training population of 800 samples could reach a genomic prediction accuracy of 0.7. Candidate genes significantly associated with resistance to NNV, including lysine-specific demethylase 2A, beta-defensin 1, and cystatin-B, which play important roles in immune responses against virus infection, were identified. Almost all the candidate genes were differentially expressed in different tissues against NNV infection. The significant genetic variants can be used in genomic selection and help understand the mechanism of resistance to VNN. Future studies should use populations of large effective size and whole genome resequencing to identify more useful genetic variants.

Aquaculture genomics, genetics and breeding in the United States: current status, challenges, and priorities for future research

Advancing the production efficiency and profitability of aquaculture is dependent upon the ability to utilize a diverse array of genetic resources. The ultimate goals of aquaculture genomics, genetics and breeding research are to enhance aquaculture production efficiency, sustainability, product quality, and profitability in support of the commercial sector and for the benefit of consumers. In order to achieve these goals, it is important to understand the genomic structure and organization of aquaculture species, and their genomic and phenomic variations, as well as the genetic basis of traits and their interrelationships. In addition, it is also important to understand the mechanisms of regulation and evolutionary conservation at the levels of genome, transcriptome, proteome, epigenome, and systems biology. With genomic information and information between the genomes and phenomes, technologies for marker/causal mutation-assisted selection, genome selection, and genome editing can be developed for applications in aquaculture. A set of genomic tools and resources must be made available including reference genome sequences and their annotations (including coding and non-coding regulatory elements), genome-wide polymorphic markers, efficient genotyping platforms, high-density and high-resolution linkage maps, and transcriptome resources including non-coding transcripts. Genomic and genetic control of important performance and production traits, such as disease resistance, feed conversion efficiency, growth rate, processing yield, behaviour, reproductive characteristics, and tolerance to environmental stressors like low dissolved oxygen, high or low water temperature and salinity, must be understood. QTL need to be identified, validated across strains, lines and populations, and their mechanisms of control understood. Causal gene(s) need to be identified. Genetic and epigenetic regulation of important aquaculture traits need to be determined, and technologies for marker-assisted selection, causal gene/mutation-assisted selection, genome selection, and genome editing using CRISPR and other technologies must be developed, demonstrated with applicability, and application to aquaculture industries.Major progress has been made in aquaculture genomics for dozens of fish and shellfish species including the development of genetic linkage maps, physical maps, microarrays, single nucleotide polymorphism (SNP) arrays, transcriptome databases and various stages of genome reference sequences. This paper provides a general review of the current status, challenges and future research needs of aquaculture genomics, genetics, and breeding, with a focus on major aquaculture species in the United States: catfish, rainbow trout, Atlantic salmon, tilapia, striped bass, oysters, and shrimp. While the overall research priorities and the practical goals are similar across various aquaculture species, the current status in each species should dictate the next priority areas within the species. This paper is an output of the USDA Workshop for Aquaculture Genomics, Genetics, and Breeding held in late March 2016 in Auburn, Alabama, with participants from all parts of the United States.

Characterization of a novel disease resistance gene rtp3 and its association with VNN disease resistance in Asian seabass

Asian seabass, an important food fish in Southeast Asia, has suffered from nervous necrosis virus (NNV) infection, resulting in massive mortality of Asian seabass larvae and enormous economic losses. Identification and characterization of disease resistance genes is important. Previous transcriptome analysis of Asians seabass epithelial cells after NNV infection revealed a highly inducible gene, receptor-transporting protein 3 (rtp3), indicating it could play an important role in Asian seabass - NNV interaction. To characterize this gene, we determined its expression pattern and subcellular localization. The rtp3 was highly induced in most examined tissues and organs of Asian seabass after NNV infection, and protein Rtp3 was localized in cytoplasm. Further association study in multiple families revealed that a microsatellite marker, (GT)ntt(GT)n, in the 3' UTR of rtp3 was significantly associated with VNN disease resistance in Asian seabass. Our results imply that rtp3 may be a novel disease resistance gene in Asian seabass. This data could improve our understanding of molecular interaction between Asian seabass and NNV, and has the potential to be applied in marker-assisted selection for disease resistance breeding in Asian seabass.

Genome-Wide Association Study for Muscle Fat Content and Abdominal Fat Traits in Common Carp (Cyprinus carpio)

Muscle fat content is an important phenotypic trait in fish, as it affects the nutritional, technical and sensory qualities of flesh. To identify loci and candidate genes associated with muscle fat content and abdominal fat traits, we performed a genome-wide association study (GWAS) using the common carp 250 K SNP assay in a common carp F₂ resource population. A total of 18 loci surpassing the genome-wide suggestive significance level were detected for 4 traits: fat content in dorsal muscle (MFdo), fat content in abdominal muscle (MFab), abdominal fat weight (AbFW), and AbFW as a percentage of eviscerated weight (AbFP). Among them, one SNP (carp089419) affecting both AbFW and AbFP reached the genome-wide significance level. Ten of those loci were harbored in or near known genes. Furthermore, relative expressions of 5 genes related to MFdo were compared using dorsal muscle samples with high and low phenotypic values. The results showed that 4 genes were differentially expressed between the high and low phenotypic groups. These genes are, therefore, prospective candidate genes for muscle fat content: ankyrin repeat domain 10a (ankrd10a), tetratricopeptide repeat, ankyrin repeat and coiled-coil containing 2 (tanc2), and four jointed box 1 (fjx1) and choline kinase alpha (chka). These results offer valuable insights into the complex genetic basis of fat metabolism and deposition.

BAC-pool sequencing and analysis confirms growth-associated QTLs in the Asian seabass genome

The Asian seabass is an important marine food fish that has been cultured for several decades in Asia Pacific. However, the lack of a high quality reference genome has hampered efforts to improve its selective breeding. A 3D BAC pool set generated in this study was screened using 22 SSR markers located on linkage group 2 which contains a growth-related QTL region. Seventy-two clones corresponding to 22 FPC contigs were sequenced by Illumina MiSeq technology. We co-assembled the MiSeq-derived scaffolds from each FPC contig with error-corrected PacBio reads, resulting in 187 sequences covering 9.7 Mb. Eleven genes annotated within this region were found to be potentially associated with growth and their tissue-specific expression was investigated. Correlation analysis demonstrated that SNPs in ctsb, skp1 and ppp2ca can be potentially used as markers for selecting fast-growing fingerlings. Conserved syntenies between seabass LG2 and five other teleosts were identified. This study i) provided a 10 Mb targeted genome assembly; ii) demonstrated NGS of BAC pools as a potential approach for mining candidates underlying QTLs of this species; iii) detected eleven genes potentially responsible for growth in the QTL region; and iv) identified useful SNP markers for selective breeding programs of Asian seabass.

Fine mapping QTL for resistance to VNN disease using a high-density linkage map in Asian seabass

Asian seabass has suffered from viral nervous necrosis (VNN) disease. Our previous study has mapped quantitative trait loci (QTL) for resistance to VNN disease. To fine map these QTL and identify causative genes, we identified 6425 single nucleotide polymorphisms (SNPs) from 85 dead and 94 surviving individuals. Combined with 155 microsatellites, we constructed a genetic map consisting of 24 linkage groups (LGs) containing 3000 markers, with an average interval of 1.27 cM. We mapped one significant and three suggestive QTL with phenotypic variation explained (PVE) of 8.3 to 11.0%, two significant and two suggestive QTL with PVE of 7.8 to 10.9%, for resistance in three LGs and survival time in four LGs, respectively. Further analysis one QTL with the largest effect identified protocadherin alpha-C 2-like (Pcdhac2) as the possible candidate gene. Association study in 43 families with 1127 individuals revealed a 6 bp insertion-deletion was significantly associated with disease resistance. qRT-PCR showed the expression of Pcdhac2 was significantly induced in the brain, muscle and skin after nervous necrosis virus (NNV) infection. Our results could facilitate marker-assisted selection (MAS) for resistance to NNV in Asian seabass and set up the basis for functional analysis of the potential causative gene for resistance.

Charactering the ZFAND3 gene mapped in the sex-determining locus in hybrid tilapia (Oreochromis spp.)

Zinc finger AN1-type domain 3 (ZFAND3) is essential for spermatogenesis in mice. However, its function in teleosts remains unclear. In this study, we characterized the ZFAND3 gene (termed as OsZFAND3) in an important food fish, tilapia. The OsZFAND3 cDNA sequence is 1,050 bp in length, containing an ORF of 615 bp, which encodes a putative peptide of 204 amino acid residues. Quantitative real-time PCR revealed that the OsZFAND3 transcripts were exclusively expressed in the testis and ovary. In situ hybridization showed that the high expression of OsZFAND3 transcripts was predominantly localized in the spermatocyte and spermatid. These results suggest that OsZFAND3 is involved in male germ cell maturation. Three single nucleotide polymorphisms (SNPs) were detected in the introns of OsZFAND3. The OsZFAND3 gene was mapped in the sex-determining locus on linkage group 1 (LG1). The three SNPs in the OsZFAND3 gene were strictly associated with sex phenotype, suggesting that the OsZFAND3 gene is tightly linked to the sex-determining locus. Our study provides new insights into the functions of the OsZFAND3 gene in tilapia and a foundation for further detailed analysis of the OsZFAND3 gene in sex determination and differentiation.

Mapping QTL for Resistance Against Viral Nervous Necrosis Disease in Asian Seabass

Viral nervous necrosis disease (VNN), caused by nervous necrosis virus (NNV), leads to mass mortality in mariculture. However, phenotypic selection for resistance against VNN is very difficult. To facilitate marker-assisted selection (MAS) for resistance against VNN and understanding of the genetic architecture underlying the resistance against this disease, we mapped quantitative trait loci (QTL) for resistance against VNN in Asian seabass. We challenged fingerlings at 37 days post-hatching (dph), from a single back-cross family, with NNV at a concentration of 9 × 10(6) TCID50/ml for 2 h. Daily mortalities were recorded and collected. A panel of 330 mortalities and 190 surviving fingerlings was genotyped using 149 microsatellites with 145 successfully mapped markers covering 24 linkage groups (LGs). Analysis of QTL for both resistance against VNN and survival time was conducted using interval mapping. Five significant QTL located in four LGs and eight suggestive QTL in seven LGs were identified for resistance. Another five significant QTL in three LGs and five suggestive QTL in three LGs were detected for survival time. One significant QTL, spanning 3 cM in LG20, was identified for both resistance and survival time. These QTL explained 2.2-4.1% of the phenotypic variance for resistance and 2.2-3.3% of the phenotypic variance for survival time, respectively. Our results suggest that VNN resistance in Asian seabass is controlled by many loci with small effects. Our data provide information for fine mapping of QTL and identification of candidate genes for a better understanding of the mechanism of disease resistance.

Molecular characterization and mapping of Fgf21 gene in a foodfish species asian seabass

Fgf21 is a newly discovered fibroblast growth factor. It is typically induced by fasting and plays important roles in the regulation of glucose and lipid metabolisms and energy balance in mammals, whereas potential functions of this gene in teleosts are still unknown. We identified the Fgf21 gene and studied its functions in Asian seabass (Lates calcarifer). The cDNA of the Fgf21 encoded a protein with 206 amino acids. Analysis of DNA and amino acid sequences of Fgf21 genes revealed that the sequences and structure of the Fgf21 genes were highly conserved in vertebrates. Real-time PCR revealed that Fgf21 was exclusively expressed in the intestine and kidney, which was different from the expression profiles of mammals. Fgf21 was down-regulated under fasting, whereas it was significantly increased during the LPS challenge. Exogenous recombinant FGF21 significantly suppressed the appetite of Asian seabass. Our data suggest that Fgf21 plays a role in energy regulation and acute phase response in Asian seabass, and may have different functions in fish and mammals. In addition, we identified one SNP in Fgf21. By using this SNP, the gene was mapped on the linkage group 23, where a suggestive QTL for growth was mapped previously. Association mapping identified significant associations between Fgf21 genotypes at the SNP and growth traits. These results not only provide important information of the functions of Fgf21, but also suggest that the SNP in this gene can be used as a marker in selecting fast-growing individuals of Asian seabass.