GWAS analysis of QTL for enteric septicemia of catfish and their involved genes suggest evolutionary conservation of a molecular mechanism of disease resistance

GWAS with principal component analysis identify QTLs associated with main peanut flavor-related traits

Peanut flavor is a complex and important trait affected by raw material and processing technology owing to its significant impact on consumer preference. In this research, principal component analysis (PCA) on 33 representative traits associated with flavor revealed that total sugars, sucrose, and total tocopherols provided more information related to peanut flavor. Genome-wide association studies (GWAS) using 102 U.S. peanut mini-core accessions were performed to study associations between 12,526 single nucleotide polymorphic (SNP) markers and the three traits. A total of 7 and 22 significant quantitative trait loci (QTLs) were identified to be significantly associated with total sugars and sucrose, respectively. Among these QTLs, four and eight candidate genes for the two traits were mined. In addition, two and five stable QTLs were identified for total sugars and sucrose in both years separately. No significant QTLs were detected for total tocopherols. The results from this research provide useful knowledge about the genetic control of peanut flavor, which will aid in clarifying the molecular mechanisms of flavor research in peanuts.

Mapping of a major QTL for increased robustness and detection of genome assembly errors in Asian seabass (Lates calcarifer)

BACKGROUND

For Asian seabass (Lates calcarifer, Bloch 1790) cultured at sea cages various aquatic pathogens, complex environmental and stress factors are considered as leading causes of disease, causing tens of millions of dollars of annual economic losses. Over the years, we conducted farm-based challenges by exposing Asian seabass juveniles to complex natural environmental conditions. In one of these challenges, we collected a total of 1,250 fish classified as either 'sensitive' or 'robust' individuals during the 28-day observation period.

RESULTS

We constructed a high-resolution linkage map with 3,089 SNPs for Asian seabass using the double digest Restriction-site Associated DNA (ddRAD) technology and a performed a search for Quantitative Trait Loci (QTL) associated with robustness. The search detected a major genome-wide significant QTL for increased robustness in pathogen-infected marine environment on linkage group 11 (ASB_LG11; 88.9 cM to 93.6 cM) with phenotypic variation explained of 81.0%. The QTL was positioned within a > 800 kb genomic region located at the tip of chromosome ASB_LG11 with two Single Nucleotide Polymorphism markers, R1-38468 and R1-61252, located near to the two ends of the QTL. When the R1-61252 marker was validated experimentally in a different mass cross population, it showed a statistically significant association with increased robustness. The majority of thirty-six potential candidate genes located within the QTL have known functions related to innate immunity, stress response or disease. By utilizing this ddRAD-based map, we detected five mis-assemblies corresponding to four chromosomes, namely ASB_LG8, ASB_LG9, ASB_LG15 and ASB_LG20, in the current Asian seabass reference genome assembly.

CONCLUSION

According to our knowledge, the QTL associated with increased robustness is the first such finding from a tropical fish species. Depending on further validation in other stocks and populations, it might be potentially useful for selecting robust Asian seabass lines in selection programs.

Whole-genome sequencing identifies interferon-induced protein IFI6/IFI27-like as a strong candidate gene for VNN resistance in European sea bass

BACKGROUND

Viral nervous necrosis (VNN) is a major disease that affects European sea bass, and understanding the biological mechanisms that underlie VNN resistance is important for the welfare of farmed fish and sustainability of production systems. The aim of this study was to identify genomic regions and genes that are associated with VNN resistance in sea bass.

RESULTS

We generated a dataset of 838,451 single nucleotide polymorphisms (SNPs) identified from whole-genome sequencing (WGS) in the parental generation of two commercial populations (A: 2371 individuals and B: 3428 individuals) of European sea bass with phenotypic records for binary survival in a VNN challenge. For each population, three cohorts were submitted to a red-spotted grouper nervous necrosis virus (RGNNV) challenge by immersion and genotyped on a 57K SNP chip. After imputation of WGS SNPs from their parents, quantitative trait loci (QTL) were mapped using a Bayesian sparse linear mixed model (BSLMM). We found several QTL regions that were specific to one of the populations on different linkage groups (LG), and one 127-kb QTL region on LG12 that was shared by both populations and included the genes ZDHHC14, which encodes a palmitoyltransferase, and IFI6/IFI27-like, which encodes an interferon-alpha induced protein. The most significant SNP in this QTL region was only 1.9 kb downstream of the coding sequence of the IFI6/IFI27-like gene. An unrelated population of four large families was used to validate the effect of the QTL. Survival rates of susceptible genotypes were 40.6% and 45.4% in populations A and B, respectively, while that of the resistant genotype was 66.2% in population B and 78% in population A.

CONCLUSIONS

We have identified a genomic region that carries a major QTL for resistance to VNN and includes the ZDHHC14 and IFI6/IFI27-like genes. The potential involvement of the interferon pathway, a well-known anti-viral defense mechanism in several organisms (chicken, human, or fish), in survival to VNN infection is of particular interest. Our results can lead to major improvements for sea bass breeding programs through marker-assisted genomic selection to obtain more resistant fish.

Reference genomes of channel catfish and blue catfish reveal multiple pericentric chromosome inversions

BACKGROUND

Channel catfish and blue catfish are the most important aquacultured species in the USA. The species do not readily intermate naturally but F₁ hybrids can be produced through artificial spawning. F₁ hybrids produced by mating channel catfish female with blue catfish male exhibit heterosis and provide an ideal system to study reproductive isolation and hybrid vigor. The purpose of the study was to generate high-quality chromosome level reference genome sequences and to determine their genomic similarities and differences.

RESULTS

We present high-quality reference genome sequences for both channel catfish and blue catfish, containing only 67 and 139 total gaps, respectively. We also report three pericentric chromosome inversions between the two genomes, as evidenced by long reads across the inversion junctions from distinct individuals, genetic linkage mapping, and PCR amplicons across the inversion junctions. Recombination rates within the inversional segments, detected as double crossovers, are extremely low among backcross progenies (progenies of channel catfish female × F₁ hybrid male), suggesting that the pericentric inversions interrupt postzygotic recombination or survival of recombinants. Identification of channel catfish- and blue catfish-specific genes, along with expansions of immunoglobulin genes and centromeric Xba elements, provides insights into genomic hallmarks of these species.

CONCLUSIONS

We generated high-quality reference genome sequences for both blue catfish and channel catfish and identified major chromosomal inversions on chromosomes 6, 11, and 24. These perimetric inversions were validated by additional sequencing analysis, genetic linkage mapping, and PCR analysis across the inversion junctions. The reference genome sequences, as well as the contrasted chromosomal architecture should provide guidance for the interspecific breeding programs.

Applying genetic technologies to combat infectious diseases in aquaculture

Disease and parasitism cause major welfare, environmental and economic concerns for global aquaculture. In this review, we examine the status and potential of technologies that exploit genetic variation in host resistance to tackle this problem. We argue that there is an urgent need to improve understanding of the genetic mechanisms involved, leading to the development of tools that can be applied to boost host resistance and reduce the disease burden. We draw on two pressing global disease problems as case studies-sea lice infestations in salmonids and white spot syndrome in shrimp. We review how the latest genetic technologies can be capitalised upon to determine the mechanisms underlying inter- and intra-species variation in pathogen/parasite resistance, and how the derived knowledge could be applied to boost disease resistance using selective breeding, gene editing and/or with targeted feed treatments and vaccines. Gene editing brings novel opportunities, but also implementation and dissemination challenges, and necessitates new protocols to integrate the technology into aquaculture breeding programmes. There is also an ongoing need to minimise risks of disease agents evolving to overcome genetic improvements to host resistance, and insights from epidemiological and evolutionary models of pathogen infestation in wild and cultured host populations are explored. Ethical issues around the different approaches for achieving genetic resistance are discussed. Application of genetic technologies and approaches has potential to improve fundamental knowledge of mechanisms affecting genetic resistance and provide effective pathways for implementation that could lead to more resistant aquaculture stocks, transforming global aquaculture.

Genome-wide association study of host resistance to the ectoparasite Ichthyophthirius multifiliis in the Amazon fish Colossoma macropomum

BACKGROUND

Tambaqui, Colossoma macropomum, is the most important native fish species farmed in South America, particularly in Brazil, where its production is limited in the southern and southeastern regions due to disease outbreaks caused by the parasite Ichthyophthirius multifiliis. Therefore, genome level analysis to understand the genetic architecture of the host resistance against I. multifiliis is fundamental to improve this trait in tambaqui. The objective of the present study was to map QTL (quantitative trait loci) associated with resistance to I. multifiliis in tambaqui by GWAS (genome-wide association study).

METHODS AND RESULTS

Individuals belonging to seven families, which were previously submitted to an experimental challenge to assess the natural resistance to the parasite I. multifiliis, were used for genomic analysis. A total of 7717 SNPs were identified in this population by ddRAD (double digest restriction site associated DNA). GWAS revealed four SNPs significantly associated in the LGs (linkage groups) 2, 9, 11 and 20 for the traits time of death and parasite load. The SNPs explained a low proportion of the variance to I. multifiliis resistance for time of death and parasite load (about 0.622% and 0.375%, respectively). The SNPs were close to 11 genes related to the immune system: abcf3, znf830, ccr9, gli3, ackr4, tbata, ndr2, tgfbr3, nhej1, znf644b, and cldn10a.

CONCLUSIONS

In conclusion, the resistance to I. multifiliis is probably under polygenic control in tambaqui, in which different QTLs of low variance can be involved in the immune responses against this ectoparasite.

Genetic and epigenetic regulation of growth, reproduction, disease resistance and stress responses in aquaculture

Major progress has been made with genomic and genetic studies in aquaculture in the last decade. However, research on epigenetic regulation of aquaculture traits is still at an early stage. It is apparent that most, if not all, aquaculture traits are regulated at both genetic and epigenetic levels. This paper reviews recent progress in understanding of genetic and epigenetic regulation of important aquaculture traits such as growth, reproduction, disease resistance, and stress responses. Although it is challenging to make generalized statements, DNA methylation is mostly correlated with down-regulation of gene expression, especially when at promoters and enhancers. As such, methylation of growth factors and their receptors is negatively correlated with growth; hypomethylation of genes important for stress tolerance is correlated with increased stress tolerance; hypomethylation of genes important for male or female sex differentiation leads to sex differentiation into males or females, respectively. It is apparent that environmental regulation of aquaculture traits is mediated at the level of epigenetic regulation, and such environment-induced epigenetic changes appeared to be intergenerationally inherited, but evidences for transgenerational inheritance are still limited.

Chromosome-level assembly and annotation of the blue catfish Ictalurus furcatus, an aquaculture species for hybrid catfish reproduction, epigenetics, and heterosis studies

Background

The blue catfish is of great value in aquaculture and recreational fisheries. The F₁ hybrids of female channel catfish (Ictalurus punctatus) × male blue catfish (Ictalurusfurcatus) have been the primary driver of US catfish production in recent years because of superior growth, survival, and carcass yield. The channel–blue hybrid also provides an excellent model to investigate molecular mechanisms of environment-dependent heterosis. However, transcriptome and methylome studies suffered from low alignment rates to the channel catfish genome due to divergence, and the genome resources for blue catfish are not publicly available.

Results

The blue catfish genome assembly is 841.86 Mbp in length with excellent continuity (8.6 Mbp contig N50, 28.2 Mbp scaffold N50) and completeness (98.6% Eukaryota and 97.0% Actinopterygii BUSCO). A total of 30,971 protein-coding genes were predicted, of which 21,781 were supported by RNA sequencing evidence. Phylogenomic analyses revealed that it diverged from channel catfish approximately 9 million years ago with 15.7 million fixed nucleotide differences. The within-species single-nucleotide polymorphism (SNP) density is 0.32% between the most aquaculturally important blue catfish strains (D&B and Rio Grande). Gene family analysis discovered significant expansion of immune-related families in the blue catfish lineage, which may contribute to disease resistance in blue catfish.

Conclusions

We reported the first high-quality, chromosome-level assembly of the blue catfish genome, which provides the necessary genomic tool kit for transcriptome and methylome analysis, SNP discovery and marker-assisted selection, gene editing and genome engineering, and reproductive enhancement of the blue catfish and hybrid catfish.

Allelically and Differentially Expressed Genes After Infection of Edwardsiella ictaluri in Channel Catfish as Determined by Bulk Segregant RNA-Seq

Identification of genetic markers associated with resistance against enteric septicemia of catfish (ESC) is of great interest for genetic enhancement programs of catfish. In the present study, bulk segregant RNA-Seq analysis was applied to determine differentially expressed genes and alleles after ESC infection. Here we report three genomic regions on LG1, LG12, and LG26, containing significant single-nucleotide polymorphisms (SNPs). These genomic regions aligned well with quantitative trait loci (QTL) previously identified. Within the QTL regions, eleven genes were found to be differentially regulated between phenotypic bulks. Importantly, the QTL on linkage group 1 (LG1) were found to be expressed in the liver, whereas the QTL on LG12 and LG26 were expressed in the intestine, suggesting multiple mechanisms of ESC resistance. It is apparent that apolipoproteins may be important for ESC resistance as the QTL on LG1 included the 14-kDa apolipoprotein genes that are both allelically expressed and differentially expressed between the resistant and susceptible bulks. Traf2 and NCK-interacting protein kinase (TNIK) were found in the QTL on LG12, and it was downregulated in resistant fish, suggesting the importance of NCK downregulation in ESC resistance, as previously reported. In addition, we observed divergent gene expression patterns between the liver and intestine after infection. Immune/inflammatory-related processes were overrepresented from liver DEGs, while those DEGs identified from intestine were enriched for proteolysis and wounding processes. Taken together, the BSR-Seq analysis presented here advanced the knowledge of ESC resistance, providing information of not only positions of QTL but also genes and their differential expression between resistant and susceptible fish, making it one step closer to the identification of the causal genes for ESC resistance.

An association of CEP78, MEF2C, VPS13A and ARRDC3 genes with survivability to heat stress in an F2 chicken population

Heat stress is a serious problem in the poultry industry. An effective tool for improving heat tolerance can be genomic selection based on single nucleotide polymorphisms. This study was performed to identify genomic regions controlling survivability to heat stress in a population of F₂ chickens that accidentally experienced acute heat stress, using Illumina 60K Chicken SNP Bead Chip. After quality control in markers, 47,730 SNPs remained for genome-wide association study (GWAS). The GWAS results indicated that markers Gga_rs16111480 (p = 8.503e-08), GGaluGA354375 (p = 5.99e-07) and Gga_rs14748694 (p = 7.085e-07) located on Z chromosome showed significant association with heat stress tolerance trait. The Gga_rs16111480 marker was located inside the CEP78 gene. The marker GGaluGA354375 was located inside the LOC101752071 gene and next to the MEF2C gene. The Gga_rs14748694 marker was adjacent to LOC101752071 and MEF2C genes. Moreover, the SNP maker of Gga_rs16111480 was located on 243 kb downstream of the VPS13A gene, and the GGaluGA354375 and Gga_rs14748694 SNPs were located on 947 kb and 888 kb downstream of the ARRDC3 gene, respectively. The results of this study suggest that apart from the gene LOC101752071, which its function was unknown, each of the two MEF2C and CEP78 genes were found to be closely related to heat stress resistance in bird.

Evolution and expression analysis of STAT family members in blunt snout bream (Megalobrama amblycephala)

The Janus kinase/signal transducer and activator of transcription (JAK/STAT) pathway is involved in regulating the body's immunity, cell proliferation, differentiation, and apoptosis. Members of the STAT family have been extensively studied in different mammalian species. However, there are few studies on the STAT family genes in farmed economic fish. In this study, eight STAT genes including STAT1a, STAT1b, STAT2, STAT3, STAT4, STAT5a, STAT5b and STAT6, in blunt snout bream (Megalobrama amblycephala), an economically important fish in China, were identified and characterized. Analyses of gene location, phylogeny and conserved synteny were conducted to infer the evolutionary origin of these STAT family genes. Furthermore, the evolutionary origin model of STATs was constructed based on the 2R hypothesis and teleost genome duplication (TGD) hypothesis, which clarified the evolutionary origin of the eight STATs in blunt snout bream. Besides, expression of the eight STATs was detected in 10 tissues of healthy blunt snout bream, which showed different expression patterns, and all had the highest level in the blood. In addition, expression of the STATs was significantly induced in the spleen, liver, and kidney after infection of Aeromonas hydrophila, suggesting that they play an important role in protecting the host from pathogens. In general, the evolution of cytokine-related genes parallels that of the immune system, which has likely been a main evolutionary driver. Therefore, the evolutionary model of STAT genes, constructed in the non-model organism pioneeringly, may provide some enlightenment for the evolution of the fish STAT family genes and their involvement in the immune function.

Whole-genome sequencing of brown-marbled grouper (Epinephelus fuscoguttatus) provides insights into adaptive evolution and growth differences

The brown-marbled grouper (Epinephelus fuscoguttatus) is an important species of fish in the coral reef ecosystem and marine aquaculture industry. In this study, a high-quality chromosome-level genome of brown-marbled grouper was assembled using Oxford Nanopore technology and Hi-C technology. The GC content and heterozygosity were approximately 42% and 0.35%, respectively. A total of 230 contigs with a total length of 1047 Mb and contig N50 of 13.8 Mb were assembled, and 228 contigs (99.13%) were anchored into 24 chromosomes. A total of 24,005 protein-coding genes were predicted, among which 23,862 (99.4%) predicted genes were annotated. Phylogenetic analysis showed that brown-marbled grouper and humpback grouper were clustered into one clade that separated approximately 11-23 million years ago. Collinearity analyses showed that there was no obvious duplication of large fragments between chromosomes in the brown-marbled grouper. Genomes of the humpback grouper and giant grouper showed a high collinearity with that of the brown-marbled grouper. A total of 305 expanded gene families were detected in the brown-marbled grouper genome, which is mainly involved in disease resistance. In addition, a genetic linkage map with 3061.88 cM was constructed. Based on the physical and genetic map, one growth-related quantitative trait loci was detected in 32,332,447 bp of chromosome 20, and meox1 and etv4 were considered candidate growth-related genes. This study provides pivotal genetic resources for further evolutionary analyses and artificial breeding of groupers.

Genome-Wide Association Study Reveals PC4 as the Candidate Gene for Thermal Tolerance in Bay Scallop (Argopecten irradians irradians)

The increasing sea temperature caused by global warming has resulted in severe mortalities in maricultural scallops. Therefore, improving thermal tolerance has become an active research area in the scallop farming industry. Bay scallop (Argopecten irradians irradians) was introduced into China in 1982 and has developed into a vast aquaculture industry in northern China. To date, genetic studies on thermal tolerance in bay scallops are limited, and no systematic screening of thermal tolerance-related loci or genes has been conducted in this species. In the present study, we conducted a genome-wide association study (GWAS) for thermal tolerance using the Arrhenius break temperature (ABT) indicators of 435 bay scallops and 38,011 single nucleotide polymorphism (SNP) markers. The GWAS identified 1,906 significant thermal tolerance-associated SNPs located in 16 chromosomes of bay scallop. Gene ontology and Kyoto Encyclopedia of Genes and Genomes pathway analyses showed that 638 genes were enriched in 42 GO terms, while 549 annotated genes were enriched in aggregation pathways. Additionally, the SNP (15-5091-20379557-1) with the lowest P value was located in the transcriptional coactivator p15 (PC4) gene, which is involved in regulating DNA damage repair and stabilizing genome functions. Further analysis in another population identified two new thermal tolerance-associated SNPs in the first coding sequence of PC4 in bay scallops (AiPC4). Moreover, AiPC4 expression levels were significantly correlated (r = 0.675–0.962; P < 0.05) with the ABT values of the examined bay scallops. Our data suggest that AiPC4 might be a positive regulator of thermal tolerance and a potential candidate gene for molecular breeding in bay scallop aiming at thermal tolerance improvement.

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

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

Genome-Wide Association Analysis Reveals the Genetic Architecture of Parasite (Cryptocaryon irritans) Resistance in Large Yellow Croaker (Larimichthys crocea)

Large yellow croaker is an important marine culture species in China. Recently, the large yellow croaker industry is threatened by various disease problems, especially for the white spot disease, which is caused by parasite Cryptocaryon irritans. In the current study, we conducted a genome-wide association study (GWAS) for C. irritans resistance in two large yellow croaker populations (n = 264 and n = 480, respectively). We identified 15 QTL with explained genetic variance ranging from 1 to 8% in the two populations. One QTL on chromosome 23 was shared by the two populations, and three QTL had been reported in the previous study. We identified a lot of biological pathways associated with C. irritans resistance, such as hormone transport, response to bacterium, apoptotic process, acute inflammatory response to antigenic stimulus, and NF-kappa B signaling pathway. The genes casp8 and traf6 involved in regulatory network for apoptosis and inflammation were identified to be candidate genes for C. irritans resistance. Our results showed the complex polygenic architecture of resistance of large yellow croaker against C. irritans. These results would be helpful for the researches of the molecular mechanism of C. irritans resistance and genome-assisted breeding of large yellow croaker.

A high-density SNP-based genetic map and several economic traits-related loci in Pelteobagrus vachelli

BACKGROUND

A high-density genetic linkage map is essential for QTL fine mapping, comparative genome analysis, identification of candidate genes and marker-assisted selection in aquaculture species. Pelteobagrus vachelli is a very popular commercial species in Asia. However, some specific characters hindered achievement of the traditional selective breeding based on phenotypes, such as lack of large-scale genomic resource and short of markers tightly associated with growth, sex determination and hypoxia tolerance related traits.

RESULTS

By making use of 5059 ddRAD markers in P. vachelli, a high-resolution genetic linkage map was successfully constructed. The map' length was 4047.01 cM by using an interval of 0.11 cm, which is an average marker standard. Comparative genome mapping revealed that a high proportion (83.2%) of markers with a one-to-one correspondence were observed between P. vachelli and P. fulvidraco. Based on the genetic map, 8 significant genome-wide QTLs for 4 weight, 1 body proportion, 2 sex determination, and 1 hypoxia tolerance related traits were detected on 4 LGs. Some SNPs from these significant genome-wide QTLs were observably associated with these phenotypic traits in other individuals by Kompetitive Allele Specific PCR. In addition, two candidate genes for weight, Sipa1 and HSD11B2, were differentially expressed between fast-, medium- and slow-growing P. vachelli. Sema7a, associated with hypoxia tolerance, was induced after hypoxia exposure and reoxygenation.

CONCLUSIONS

We mapped a set of suggestive and significant QTLs as well as candidate genes for 12 growth, 1 sex determination and 1 hypoxia tolerance related traits based on a high-density genetic linkage map by making use of SNP markers for P. fulvidraco. Our results have offered a valuable method about the much more efficient production of all-male, fast growth and hypoxia tolerance P. vachelli for the aquaculture industry.

Effect of dietary isoleucine on skin mucus barrier and epithelial physical barrier functions of hybrid bagrid catfish Pelteobagrus vachelli × Leiocassis longirostris

The study investigated the effects of dietary isoleucine (Ile) on skin mucus barrier and epithelial physical barrier functions of hybrid bagrid catfish Pelteobagrus vachelli × Leiocassis longirostris. A total of 630 fish (33.11 ± 0.09 g) were fed semi-purified isonitrogenous diets containing 5.0 (control), 7.5, 10.0, 12.5, 15.0, 17.5, and 20.0 g Ile kg ^-1 diet for 8 weeks. The results indicated that dietary Ile increased (P < 0.05) in skin (1) mucus protein content and antimicrobial activity against three gram-negative bacteria (Aeromonas hydrophila, Escherichia coli, and Yersinia ruckeri) and two gram-positive bacteria (Streptococcus agalactiae and Staphylococcus aureus), (2) mucus lysofew information is available about the influencezyme (LZM), acid phosphatase (ACP), and alkaline phosphatase (AKP) activities, and complement 3 and 4 (C3 and C4) and immunoglobulin M (IgM) contents, (3) intelectin 1 (intl1), intelectin 2 (intl2), c-type-lysozyme (c-LZM), g-type-lysozyme (g-LZM), and β-defensin mRNA levels. Dietary Ile decreased (P < 0.05) reactive oxygen species (ROS), malondialdehyde (MDA), and protein carbonyl (PC) contents, and up-regulated (P < 0.05) CuZnSOD, GST, GPX1a, muc5ac, muc5b, zonula occludens-1 (ZO-1), zonula occludens-2 (ZO-2), occludin, and claudin 3 mRNA levels in skin. These results indicated that Ile improved skin mucus barrier function via increasing mucus protein, C3 and C4, and IgM contents and antibacterial factors activities, and promoted epithelial physical barrier function via decreasing skin antioxidant damage and improving tight junction structure in hybrid bagrid catfish.

Combining Multiple Approaches and Models to Dissect the Genetic Architecture of Resistance to Infections in Fish

Infectious diseases represent a major threat for the sustainable development of fish farming. Efficient vaccines are not available against all diseases, and growing antibiotics resistance limits the use of antimicrobial drugs in aquaculture. It is therefore important to understand the basis of fish natural resistance to infections to help genetic selection and to develop new approaches against infectious diseases. However, the identification of the main mechanisms determining the resistance or susceptibility of a host to a pathogenic microbe is challenging, integrating the complexity of the variation of host genetics, the variability of pathogens, and their capacity of fast evolution and adaptation. Multiple approaches have been used for this purpose: (i) genetic approaches, QTL (quantitative trait loci) mapping or GWAS (genome-wide association study) analysis, to dissect the genetic architecture of disease resistance, and (ii) transcriptomics and functional assays to link the genetic constitution of a fish to the molecular mechanisms involved in its interactions with pathogens. To date, many studies in a wide range of fish species have investigated the genetic determinism of resistance to many diseases using QTL mapping or GWAS analyses. A few of these studies pointed mainly toward adaptive mechanisms of resistance/susceptibility to infections; others pointed toward innate or intrinsic mechanisms. However, in the majority of studies, underlying mechanisms remain unknown. By comparing gene expression profiles between resistant and susceptible genetic backgrounds, transcriptomics studies have contributed to build a framework of gene pathways determining fish responsiveness to a number of pathogens. Adding functional assays to expression and genetic approaches has led to a better understanding of resistance mechanisms in some cases. The development of knock-out approaches will complement these analyses and help to validate putative candidate genes critical for resistance to infections. In this review, we highlight fish isogenic lines as a unique biological material to unravel the complexity of host response to different pathogens. In the future, combining multiple approaches will lead to a better understanding of the dynamics of interaction between the pathogen and the host immune response, and contribute to the identification of potential targets of selection for improved resistance.

Proteomic profiling of yellow catfish (Pelteobagrus fulvidraco) skin mucus identifies differentially-expressed proteins in response to Edwardsiella ictaluri infection

Fish mucus acts as a physiological and immunological barrier for maintaining normal fish physiology and conferring defense against pathogens infection. Here we report proteomic profiling of skin mucus of yellow catfish before and after E. ictaluri infection by Label-free LC-MS/MS approach. A total of 918 non-redundant proteins were identified from 54443 spectra referring to yellow catfish genome database. Further annotation via GO and KEGG database revealed complex protein composition of yellow catfish mucus. Besides structural proteins in mucus, a lot of immune-related proteins were retrieved, such as lectins, complement components, antibacterial peptides and immunoglobins. 133 differentially-expressed proteins (DEPs), including 76 up-regulated and 57 down-regulated proteins, were identified, most of which were enriched into 17 pathways centering on "immune system" category with 33 proteins involved. Consistently, significant proliferation of mucus-secreting goblet cells and CYPA-expressing cells were observed along outside of yellow catfish skin after E. ictaluri infection, indicating an enhanced immune response to E. ictaluri infection in yellow catfish skin mucus. The proteomic data provide systematic protein information to comprehensively understand the biological function of yellow catfish skin mucus in response to bacterial infection.

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

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

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.

Polyadenylation sites and their characteristics in the genome of channel catfish (Ictalurus punctatus) as revealed by using RNA-Seq data

Polyadenylation plays important roles in gene expression regulation in eukaryotes, which typically involves cleavage and poly(A) tail addition at the polyadenylation site (PAS) of the pre-mature mRNA. Many eukaryotic genes contain more than one PASs, termed as alternative polyadenylation (APA). As a crucial post-transcriptional regulation, polyadenylation affects various aspects of RNA metabolism such as mRNA stability, translocation, and translation. However, polyadenylation has been rarely studied in teleosts. Here we conducted polyadenylation analysis in channel catfish, a commercially important aquaculture species around the world. Using RNA-Seq data, we identified 20,320 PASs which were classified into 14,500 clusters by merging adjacent PASs. Most of the PASs were found in 3' UTRs, followed by intron regions based on the annotation of channel catfish reference genome. No apparent difference in PAS distribution was observed between the sense and antisense strand of the channel catfish genome. The sequence analysis of nucleotide composition and motif around PASs yielded a highly similar profile among various organisms, suggesting the conservation and importance of polyadenylation in evolution. Using APA genes with more than two PASs, gene ontology enrichment revealed genes particularly involved in RNA binding. Reactome pathway analysis showed the enrichment of the innate immune system, especially neutrophil degranulation.

Constructing a High-Density Genetic Linkage Map for Large Yellow Croaker (Larimichthys crocea) and Mapping Resistance Trait Against Ciliate Parasite Cryptocaryon irritans

The large yellow croaker (Larimichthys crocea) is the most economically important marine cage-farming fish in China in the past decade. However, the sustainable development of large yellow croaker aquaculture has been severely hampered by several diseases, of which, the white spot disease caused by ciliate protozoan parasite Cryptocaryon irritans ranks the most damaging disease in large yellow croaker cage farms. To better understand the genetic basis of parasite infection and disease resistance to C. irritans, it is vital to map the traits and localize the underlying candidate genes in L. crocea genome. Here, we constructed a high-density genetic linkage map using double-digest restriction-site associated DNA (ddRAD)-based high-throughput SNP genotyping data of a F1 mapping family, which had been challenged with C. irritans for resistant trait measure. A total of 5261 SNPs was grouped and oriented into 24 linkage groups (LGs), representing 24 chromosomes of L. crocea. The total genetic map length was 1885.67 cM with an average inter-locus distance of 0.36 cM. Quantitative trait loci (QTL) mapping identified seven significant QTLs in four LGs linked to C. irritans disease resistance. Candidate genes underlying disease resistance were identified from the reference genome, including ifnar1, ifngr2, ikbke, and CD112. Comparative genomic analysis between large yellow croaker and the four closely related species revealed high evolutionary conservation of chromosomes, though inter-chromosomal rearrangements do exist. Especially, the croaker genome structure was closer to the medaka genome than stickleback, indicating that the croaker genome might retain the teleost ancestral genome structure. The high-density genetic linkage map provides an important tool and resource for fine mapping, comparative genome analysis, and molecular selective breeding of large yellow croaker.

Genomic Selection in Aquaculture: Application, Limitations and Opportunities With Special Reference to Marine Shrimp and Pearl Oysters

Within aquaculture industries, selection based on genomic information (genomic selection) has the profound potential to change genetic improvement programs and production systems. Genomic selection exploits the use of realized genomic relationships among individuals and information from genome-wide markers in close linkage disequilibrium with genes of biological and economic importance. We discuss the technical advances, practical requirements, and commercial applications that have made genomic selection feasible in a range of aquaculture industries, with a particular focus on molluscs (pearl oysters, Pinctada maxima) and marine shrimp (Litopenaeus vannamei and Penaeus monodon). The use of low-cost genome sequencing has enabled cost-effective genotyping on a large scale and is of particular value for species without a reference genome or access to commercial genotyping arrays. We highlight the pitfalls and offer the solutions to the genotyping by sequencing approach and the building of appropriate genetic resources to undertake genomic selection from first-hand experience. We describe the potential to capture large-scale commercial phenotypes based on image analysis and artificial intelligence through machine learning, as inputs for calculation of genomic breeding values. The application of genomic selection over traditional aquatic breeding programs offers significant advantages through being able to accurately predict complex polygenic traits including disease resistance; increasing rates of genetic gain; minimizing inbreeding; and negating potential limiting effects of genotype by environment interactions. Further practical advantages of genomic selection through the use of large-scale communal mating and rearing systems are highlighted, as well as presenting rate-limiting steps that impact on attaining maximum benefits from adopting genomic selection. Genomic selection is now at the tipping point where commercial applications can be readily adopted and offer significant short- and long-term solutions to sustainable and profitable aquaculture industries.

Increased Alternative Splicing as a Host Response to Edwardsiella ictaluri Infection in Catfish

Alternative splicing is the process of generating multiple transcripts from a single pre-mRNA used by eukaryotes to regulate gene expression and increase proteomic complexity. Although alternative splicing profiles have been well studied in mammalian species, they have not been well studied in aquatic species, especially after biotic stresses. In the present study, genomic information and RNA-Seq datasets were utilized to characterize alternative splicing profiles and their induced changes after bacterial infection with Edwardsiella ictaluri in channel catfish (Ictalurus punctatus). A total of 27,476 alternative splicing events, derived from 9694 genes, were identified in channel catfish. Exon skipping was the most abundant while mutually exclusive exon was the least abundant type of alternative splicing. Alternative splicing was greatly induced by E. ictaluri infection with 21.9% increase in alternative splicing events. Interestingly, genes involved in RNA binding and RNA splicing themselves were significantly enriched in differentially alternatively spliced genes after infection. Sequence analyses of splice variants of a representative alternatively spliced gene, splicing factor srsf2, revealed that certain spliced transcripts may undergo nonsense-mediated decay (NMD), suggesting functional significance of the induced alternative splicing. Although statistical analysis was not possible with such large datasets, results from quantitative real-time PCR from representative differential alternative splicing events provided general validation of the bacterial infection-induced alternative splicing. This is the first comprehensive study of alternative splicing and its changes in response to bacterial infection in fish species, providing insights into the molecular mechanisms of host responses to biotic stresses.

Mapping and Sequencing of a Significant Quantitative Trait Locus Affecting Resistance to Koi Herpesvirus in Common Carp

Cyprinids are the most highly produced group of fishes globally, with common carp being one of the most valuable species of the group. Koi herpesvirus (KHV) infections can result in high levels of mortality, causing major economic losses, and is listed as a notifiable disease by the World Organization for Animal Health. Selective breeding for host resistance has the potential to reduce morbidity and losses due to KHV. Therefore, improving knowledge about host resistance and methods of incorporating genomic data into breeding for resistance may contribute to a decrease in economic losses in carp farming. In the current study, a population of 1,425 carp juveniles, originating from a factorial cross between 40 sires and 20 dams was challenged with KHV. Mortalities and survivors were recorded and sampled for genotyping by sequencing using Restriction Site-Associated DNA sequencing (RADseq). Genome-wide association analyses were performed to investigate the genetic architecture of resistance to KHV. A genome-wide significant QTL affecting resistance to KHV was identified on linkage group 44, explaining approximately 7% of the additive genetic variance. Pooled whole genome resequencing of a subset of resistant (n = 60) and susceptible animals (n = 60) was performed to characterize QTL regions, including identification of putative candidate genes and functional annotation of associated polymorphisms. The TRIM25 gene was identified as a promising positional and functional candidate within the QTL region of LG 44, and a putative premature stop mutation in this gene was discovered.

RADseq and mate choice assays reveal unidirectional gene flow among three lamprey ecotypes despite weak assortative mating: Insights into the formation and stability of multiple ecotypes in sympatry

Adaptive divergence with gene flow often results in complex patterns of variation within taxa exhibiting substantial ecological differences among populations. One example where this may have occurred is the parallel evolution of freshwater-resident nonparasitic lampreys from anadromous-parasitic ancestors. Previous studies have focused on transitions between these two phenotypic extremes, but here, we considered more complex evolutionary scenarios where an intermediate freshwater form that remains parasitic is found sympatrically with the other two ecotypes. Using population genomic analysis (restriction-associated DNA sequencing), we found that a freshwater-parasitic ecotype was highly distinct from an anadromous-parasitic form (Q_lake-P  = 96.8%, F_st  = 0.154), but that a freshwater-nonparasitic form was almost completely admixed in Loch Lomond, Scotland. Demographic reconstructions indicated that both freshwater populations likely derived from a common freshwater ancestor. However, while the nonparasitic ecotype has experienced high levels of introgression from the anadromous-parasitic ecotype (Q_anad-P  = 37.7%), there is no evidence of introgression into the freshwater-parasitic ecotype. Paradoxically, mate choice experiments predicted high potential for gene flow: Males from all ecotypes were stimulated to spawn with freshwater-parasitic females, which released gametes in response to all ecotypes. Differentially fixed single nucleotide polymorphisms identified genes associated with growth and development, which could possibly influence the timing of metamorphosis, resulting in significant ecological differences between forms. This suggests that multiple lamprey ecotypes can persist in sympatry following shifts in adaptive peaks, due to environmental change during their repeated colonization of post-glacial regions, followed by periods of extensive gene flow among such diverging populations.

A Review of Molecular Responses of Catfish to Bacterial Diseases and Abiotic Stresses

Catfish is one of the major aquaculture species in the United States. However, the catfish industry is threatened by several bacterial diseases such as enteric septicemia of catfish (ESC), columnaris disease and Aeromonas disease, as well as by abiotic stresses such as high temperature and low oxygen. Research has been conducted for several decades to understand the host responses to these diseases and abiotic stresses. With the development of sequencing technologies, and the application of genome-wide association studies in aquaculture species, significant progress has been made. This review article summarizes recent progress in understanding the molecular responses of catfish after bacterial infection and stress challenges, and in understanding of genomic and genetic basis for disease resistance and stress tolerance.

Genome-wide association analysis of intra-specific QTL associated with the resistance for enteric septicemia of catfish

Disease resistance is one of the most important traits for aquaculture industry. For catfish industry, enteric septicemia of catfish (ESC), caused by the bacterial pathogen Edwardsiella ictaluri, is the most severe disease, causing enormous economic losses every year. In this study, we used three channel catfish families with 900 individuals (300 fish per family) and the 690K catfish SNP array, and conducted a genome-wide association study to detect the quantitative trait loci (QTL) associated with ESC resistance. Three significant QTL, with two of located on LG1 and one on LG26, and three suggestive QTL located on LG1, LG3, and LG21, respectively, were identified to be associated with ESC resistance. With a well-assembled- and -annotated reference genome sequence, genes around the involved QTL regions were identified. Among these genes, 37 genes had known functions in immunity, which may be involved in ESC resistance. Notably, nlrc3 and nlrp12 identified here were also found in QTL regions of ESC resistance in the channel catfish × blue catfish interspecific hybrid system, suggesting this QTL was operating within both intra-specific channel catfish populations and interspecific hybrid backcross populations. Many of the genes of the Class I MHC pathway, for mediated antigen processing and presentation, were found in the QTL regions. The positional correlation found in this study and the expressional correlation found in previous studies indicated that Class I MHC pathway was significantly associated with ESC resistance. This study validated one QTL previously identified using the second and fourth generation of the interspecific hybrid backcross progenies, and identified five additional QTL among channel catfish families. Taken together, it appears that there are only a few major QTL for ESC disease resistance, making marker-assisted selection an effective approach for genetic improvements of ESC resistance.

GWAS analysis using interspecific backcross progenies reveals superior blue catfish alleles responsible for strong resistance against enteric septicemia of catfish

Infectious diseases pose significant threats to the catfish industry. Enteric septicemia of catfish (ESC) caused by Edwardsiella ictaluri is the most devastating disease for catfish aquaculture, causing huge economic losses annually. Channel catfish and blue catfish exhibit great contrast in resistance against ESC, with channel catfish being highly susceptible and blue catfish being highly resistant. As such, the interspecific backcross progenies provide an ideal system for the identification of quantitative trait locus (QTL). We previously reported one significant QTL on linkage group (LG) 1 using the third-generation backcrosses, but the number of founders used to make the second- and third-generation backcross progenies was very small. Although the third-generation backcross progenies provided a greater power for fine mapping than the first-generation backcrosses, some major QTL for disease resistance may have been missing due to the small numbers of founders used to produce the higher generation backcrosses. In this study, we performed a genome-wide association study using first-generation backcrosses with the catfish 690 K SNP arrays to identify additional ESC disease resistance QTL, especially those at the species level. Two genomic regions on LG1 and LG23 were determined to be significantly associated with ESC resistance as revealed by a mixed linear model and family-based association test. Examination of the resistance alleles indicated their origin from blue catfish, indicating that at least two major disease resistance loci exist among blue catfish populations. Upon further validation, markers linked with major ESC disease resistance QTL should be useful for marker-assisted introgression, allowing development of highly ESC resistant breeds of catfish.

Basal polarization of the immune responses to Streptococcus agalactiae susceptible and resistant tilapia (Oreochromis niloticus)

One of the highest priority areas for improvement is the development of effective strategies for decreasing disease mortality levels in aquaculture production, a better understanding of the components of the fish immune system and their functions in the context of pathogen invasion is needed. Tilapia is the most common fish in South China, and Streptococcus agalactiae has become the most serious disease problem for tilapia industry in China. Here, we profiled gene expression differences between tilapia differing in their susceptibility to S. agalactiae both basally (before infection) and at three early timepoints post-infection (5 h, 50 h, and 7 d). Between group comparisons revealed 5756 unique genes differentially expressed greater than 2-fold at one or more timepoints. And the resistant fish showed much more strong ability in pathogen recognition, antigen presentation, immune activation, while the susceptible fish showed fast activation of apoptosis. Taken together, the immune profiles expand our knowledge for molecular mechanisms for disease resistance, as well as provide solid molecular resources for further identification of the candidate markers for disease-resistant selection and evaluation of disease prevention and treatment options for tilapia industry.

JAK and STAT members in channel catfish: Identification, phylogenetic analysis and expression profiling after Edwardsiella ictaluri infection

The Janus kinase/signal transducers and activators of transcription (JAK/STAT) signaling pathway is one of the main pleiotropic cascades used to transmit information from extracellular receptors to the nucleus, which results in DNA transcription and expression of genes involved in immunity, proliferation, differentiation, migration, apoptosis, and cell survival. Members of JAK family and STAT family have been extensively studied in different mammalian species because of their important roles in innate and adaptive immune responses. However, they have not been systematically studied among teleost fish species. In this study, five JAK family members and eight STAT family members were identified and characterized from channel catfish. Phylogenetic analysis was conducted to properly annotate these genes. Syntenic analysis was also conducted to establish orthology, and confirm the results from phylogenetic analysis. Compared to mammals, more members of the JAK and STAT family were identified in channel catfish genome. Expression of JAK and STAT family members was detected in healthy catfish tissues, but was induced in gill, liver, and intestine after bacterial challenge. Notably, the significant upregulation of STAT1b gene in catfish liver, gill and intestine after Edwardsiella ictaluri infection supported the notion that high STAT1 expression are involved in defense against pathogens. Collectively, the increased expression of JAK and STAT members in tested tissues suggested their crucial function in defending the host against pathogen invasion.

Applied Hologenomics: Feasibility and Potential in Aquaculture

Aquaculture will play an essential role in feeding a growing human population, but several biological challenges impede sustainable growth of production. Emerging evidence across all areas of life has revealed the importance of the intimate biological interactions between animals and their associated gut microbiota. Based on challenges in aquaculture, we leverage current knowledge in molecular biology and host microbiota interactions to propose an applied holo-omic framework that integrates molecular data including genomes, transcriptomes, epigenomes, proteomes, and metabolomes for analyzing fish and their gut microbiota as interconnected and coregulated systems. With an eye towards aquaculture, we discuss the feasibility and potential of our holo-omic framework to improve growth, health, and sustainability in any area of food production, including livestock and agriculture.

Identification of novel genes significantly affecting growth in catfish through GWAS analysis

Growth is the most important economic trait in aquaculture. Improvements in growth-related traits can enhance production, reduce costs and time to produce market-size fish. Catfish is the major aquaculture species in the United States, accounting for 65% of the US finfish production. However, the genes underlying growth traits in catfish were not well studied. Currently, the majority of the US catfish industry uses hybrid catfish derived from channel catfish female mated with blue catfish male. Interestingly, channel catfish and blue catfish exhibit differences in growth-related traits, and therefore the backcross progenies provide an efficient system for QTL analysis. In this study, we conducted a genome-wide association study for catfish body weight using the 250 K SNP array with 556 backcross progenies generated from backcross of male F1 hybrid (female channel catfish × male blue catfish) with female channel catfish. A genomic region of approximately 1 Mb on linkage group 5 was found to be significantly associated with body weight. In addition, four suggestively associated QTL regions were identified on linkage groups 1, 2, 23 and 24. Most candidate genes in the associated regions are known to be involved in muscle growth and bone development, some of which were reported to be associated with obesity in humans and pigs, suggesting that the functions of these genes may be evolutionarily conserved in controlling growth. Additional fine mapping or functional studies should allow identification of the causal genes for fast growth in catfish, and elucidation of molecular mechanisms of regulation of growth in fish.

A Genome-Wide Association Study Reveals That Genes with Functions for Bone Development Are Associated with Body Conformation in Catfish

Body conformation is of great scientific and commercial interest for aquaculture fish species because it affects biological adaptation of the organism to environments, and is of economic importance to the aquaculture industry considering its direct effect on fillet yield. Catfish is the primary aquaculture species in the USA. Two major species used in the aquaculture industry, channel catfish and blue catfish, differ in body shape and therefore the backcross progenies serve as a good model for quantitative trait locus (QTL) analysis. Here, a genome-wide association study (GWAS) with hybrid catfish was conducted to identify the QTL for body conformation, including deheaded body length (DBL), body length (BL), body depth (BD), and body breadth (BB), which were all standardized by cubic root of body weight. Overall, the results indicate that the traits are polygenic. For DBL, linkage group (LG) 2 and LG 24 contain significant QTL, and LG 13 and LG 26 contain suggestively associated QTL (-log₁₀(P value) > 4.5). Compared with DBL, additional SNPs were identified to be associated with body length on LG 2, LG 7, and LG 18. Although no significant QTL for body depth was found, three suggestively associated QTLs were identified on LG 5, LG 13, and LG 14. No SNP for body breadth reached the threshold for suggestive association. Genes close to the associated SNPs were determined, many of which are known to be involved in bone development. This work therefore provides the basis for future identification of causal genes for the control of body conformation.

The NCK and ABI adaptor genes in catfish and their involvement in ESC disease response

Adaptor proteins non-catalytic region of tyrosine kinase (NCK) and Abelson interactor (ABI) are crucial for disease response. NCK1 was identified to be a candidate gene for enteric septicemia of catfish (ESC) disease resistance, and was speculated to play similar roles during ESC and enteropathogenic Escherichia coli (EPEC) pathogenicity. ABI1 was reported as a positional candidate gene for bacterial cold water disease (BCWD) resistance in rainbow trout. In this study, three NCK genes and six ABI genes were identified in the channel catfish (Ictalurus punctatus) genome and blue catfish (I. furcatus) transcriptome, and annotated by domain structures, phylogenetic and syntenic analyses. Their expression patterns were examined in the intestine and liver of catfish after challenge with Edwardsiella ictaluri. In the intestine, NCK1, ABI2a, ABI2b, ABI3a were differentially expressed after E. ictaluri infection. In the liver, NCK2a, NCK2b, ABI1b, ABI2a, ABI2b were significantly upregulated in ESC susceptible fish. In general, the NCK and ABI genes, with exception of ABI3a gene and NCK1 gene, were expressed at higher levels in susceptible fish after infection than in control fish, but were expressed at lower levels in resistant fish than in the control fish. Taken together, these results support the notion that NCK and ABI genes are involved in disease processes facilitating pathogenesis of the E. ictaluri bacteria.

A Genome-Wide Association Study Identifies the Genomic Region Associated with Shell Color in Yesso Scallop, Patinopecten yessoensis

The shell color polymorphism widely exists in economic shellfish, which not only results in a better visual perception but also shows great value as an economic trait for breeding. Small numbers of reddish-orange shell Yesso scallops, Patinopecten yessoensis, were found in cultured populations compared to the brown majority. In this study, a genome-wide association study was conducted to understand the genetic basis of shell color. Sixty-six 2b-RAD libraries with equal numbers of reddish-orange and brown shell individuals were constructed and sequenced using the Illumina HiSeq 2000 platform. A total of 322,332,684 high-quality reads were obtained, and the average sequencing depth was 18.4×. One genomic region on chromosome 11 that included 239 single-nucleotide polymorphisms (SNPs) was identified as significantly associated with shell color. After verification by high-resolution melting in another population, two SNPs were selected as specific loci for reddish-orange shell color. These two SNPs could be used to improve the selective breeding progress of true-breeding strains with complete reddish-orange scallops. In addition, within the significantly associated genomic region, candidate genes were identified using marker sequences to search the draft genome of Yesso scallop. Three genes (LDLR, FRIS, and FRIY) with known functions in carotenoid metabolism were identified. Further study using high-performance liquid chromatography proved that the relative level of carotenoids in the reddish-orange shells was 40 times higher than that in the brown shells. These results suggested that the accumulation of carotenoids contributes to the formation of reddish-orange shells.

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.