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

Transcriptome analysis reveals the immune responses of leopard coral grouper to nervous necrosis virus infection

The leopard coral grouper (Plectropomus leopardus) is a key marine aquaculture species that is vulnerable to nervous necrosis virus infection during its larval stage. This infection often causes viral nervous necrosis disease, resulting in significant losses in aquaculture. The brain tissue serves as a target organ for viral invasion in fish, with vacuolization being the most prominent characteristic observed in susceptible individuals. In this study, a transcriptome analysis was conducted using brain tissues from naturally resistant and susceptible P. leopardus to investigate the defensive mechanisms against nervous necrosis virus infection. Histological observations revealed irregular vacuolization in the brains of fish infected with the nervous necrosis virus. Based on the differentially expressed genes, our study revealed a significant upregulation pathway of neuroactive ligand-receptor interaction, glutamatergic synapse, and GABAergic synapse in the resistant group, while COVID-19 infection pathway and NF-kappa B signaling pathway were primarily downregulated in the susceptible group. These findings suggest that neuroprotective and repair mechanisms may enhance resistance to viral infection while suppressing the immune response, thereby preventing hyperactive damage caused by inflammation. The expression profiles of key differentially expressed genes were validated using qRT-PCR. These genes include those related to the nervous system (grin2b, id4, mpz, sema6bb, dab1a, and nmda2d), cytokine interaction (ccr7), and antigen processing and presentation (hsp90aa1). Fluorescence in situ hybridization (FISH) demonstrated virus co-localization at both the tissue and cellular levels. These findings provide insights into the immune mechanism of leopard coral grouper, aiding in the prevention and treatment of nervous necrosis virus infection.

Genome-Wide Association Study Reveals Growth-Related SNPs and Candidate Genes in Largemouth Bass (Micropterus salmoides) Adapted to Hypertonic Environments

Sustainable development of the largemouth bass industry is hindered by limited freshwater resources. Consequently, the expansion of farming space by brackish and saline water aquaculture has become imperative. Largemouth bass is an economically important freshwater fish species. However, there is presently a lack of germplasm resources with the capacity to adapt to hypertonic environments and maintain rapid growth. A genome-wide association study is a technique used for the detection of genetic variants associated with specific phenotypic traits. In this study, we firstly applied this technique to explore the potential single-nucleotide polymorphism (SNP) locus and candidate genes associated with rapid growth and adaptation to the hypertonic environment of largemouth bass. A total of 10 potential growth-related SNPs were obtained on chromosome 16, and SNP16:4120214 was a significant SNP peak. Based on these SNPs, 23 candidate genes were annotated in the genome, including Nkcc1, Mapkap1, Hmgcs1, Slc27a6, and Shroom3. Shroom3 expression was significantly higher in individuals enriched for the most growth-advantageous genotypes. Shroom3 upregulation is beneficial for fish growth in hyperosmotic environments. This study provides insight into the genetic basis of rapid growth in hypertonic environments and foundational information for the future breeding of salt-tolerant largemouth bass.

Highly Efficient Methods with a Generalized Linear Mixed Model for the Quantitative Trait Locus Mapping of Resistance to Columnaris Disease in Rainbow Trout (Oncorhynchus mykiss)

Linear mixed models (LMMs) are commonly used in genome-wide association studies (GWASs) to evaluate population structures and relatedness. However, LMMs have been shown to be ineffective in controlling false positive errors for the analysis of resistance to Columnaris disease in Rainbow Trout. To solve this problem, we conducted a series of studies using generalized linear mixed-model association software such as GMMAT (v1.4.0) (generalized linear mixed-model association tests), SAIGE (v1.4.0) (Scalable and Accurate Implementation of Generalized mixed model), and Optim-GRAMMAR for scanning a total of 25,853 SNPs. Seven different SNPs (single-nucleotide polymorphisms) associated with the trait of resistance to Columnaris were detected by Optim-GRAMMAR, four SNPs were detected by GMMAT, and three SNPs were detected by SAIGE, and all of these SNPs can explain 8.87% of the genetic variance of the trait of resistance to Columnaris disease. The heritability of the trait of resistance to Columnaris re-evaluated by GMMAT was calibrated and was found to amount to a total of 0.71 other than 0.35, which was seriously underestimated in previous research. The identification of LOC110520307, LOC110520314, and LOC110520317 associated with the resistance to Columnaris disease will provide us more genes to improve the genetic breeding by molecular markers. Finally, we continued the haplotype and gene-based analysis and successfully identified some haplotypes and a gene (TTf-2) associated with resistance to Columnaris disease.

First insight of the genome-wide association study and genomic prediction into enteritis disease (Vibrio harveyi) resistance trait in the lined seahorse (Hippocampus erectus)

Enteritis caused by Vibrio is a highly die-off disease that severely impeded substantial production in seahorse aquaculture. In the present study, challenged with LD50 of concentration of Vibrio harveyi, a total of 161 of susceptible and 166 of resistant individuals were allocated into binary survival phenotypes, thus, to firstly investigate the genetic architecture by genome-wide association study (GWAS) analysis, as well as to evaluate the feasibility of genomic selection (GS) in enteritis disease resistance trait of the lined seahorse Hippocampus erectus. Results indicated that the heritability for resistance to Vibrio harveyi was estimated to be 0.10. And a set of 10 significant/suggestive SNPs in a multiple chromosomes localization were identified, explaining 7.76% to 13.28% of genetic variance. Associated 82 of candidate genes were clustered into signal transduction, cell proliferation, response of external stress, bacteria defence, and anti-inflammatory processes. Moreover, the potential performance of genomic selection (GS) in application in selective breeding for enteritis disease resistance seahorses was assessed by genomic prediction (GP). In general, the predictive accuracy of the genomic estimated breeding value (GEBV) of BayesC exceeded the rrBLUP, BayesA, RKHS, and SVM models while with no significant difference. And the GWAS-informative SNPs was significantly superior in efficience than random selected markers by comparison of predictive performance on different selection strategies of SNPs. Overall, the genetic basis of enteritis disease resistance trait in the lined seahorse is a polygenic genetic architecture. SNPs associated with the important genes of cathepsin L1-like previously reported with respect to disease resistance are consider as potential molecular markers of genetic breeding. Furthermore, GS approach is an appropriate, effective, and less-cost application in breeding enteritis disease-resistant seahorses.

Genome-Wide Association Study of Resistance to Largemouth Bass Ranavirus (LMBV) in Micropterus salmoides

The disease caused by Largemouth bass ranavirus (LMBV) is one of the most severe viral diseases in largemouth bass (Micropterus salmoides). It is crucial to evaluate the genetic resistance of largemouth bass to LMBV and develop markers for disease-resistance breeding. In this study, 100 individuals (45 resistant and 55 susceptible) were sequenced and evaluated for resistance to LMBV and a total of 2,579,770 variant sites (SNPs-single-nucleotide polymorphisms (SNPs) and insertions-deletions (InDels)) were identified. A total of 2348 SNPs-InDels and 1018 putative candidate genes associated with LMBV resistance were identified by genome-wide association analyses (GWAS). Furthermore, GO and KEGG analyses revealed that the 10 candidate genes (MHC II, p38 MAPK, AMPK, SGK1, FOXO3, FOXO6, S1PR1, IL7R, RBL2, and GADD45) were related to intestinal immune network for IgA production pathway and FoxO signaling pathway. The acquisition of candidate genes related to resistance will help to explore the molecular mechanism of resistance to LMBV in largemouth bass. The potential polymorphic markers identified in this study are important molecular markers for disease resistance breeding in largemouth bass.

Identification of candidate SNPs and genes associated with resistance to nervous necrosis virus in leopard coral grouper (Plectropomus leopardus) using GWAS

The leopard coral grouper (Plectropomus leopardus), which has become increasingly popular in consumption due to its bright body color and great nutritional, holds a high economic and breeding potential. However, in recent years, the P.leopardus aquaculture industry has been impeded by the nervous necrosis virus (NNV) outbreak, leading to widespread mortality among fry and juvenile grouper. However, the genetic basis of resistance to NNV in P. leopardus remains to be investigated. In the present study, we conducted a genome-wide association analysis (GWAS) on 100 resistant and 100 susceptible samples to discover variants and potential genes linked with NNV resistance. For this study, 157,926 high-quality single nucleotide polymorphisms (SNPs) based on whole genome resequencing were discovered, and eighteen SNPs loci linked to disease resistance were discovered. We annotated six relevant candidate genes, including sik2, herc2, pip5k1c, npr1, mybpc3, and arhgap9, which showed important roles in lipid metabolism, oxidative stress, and neuronal survival. In the brain tissues of resistant and susceptible groups, candidate genes against NNV infection showed significant differential expression. The results indicate that regulating neuronal survival or pathways involved in lipid metabolism may result in increased resistance to NNV. Understanding the molecular mechanisms that lead to NNV resistance will be beneficial for the growth of the P. leopardus breeding sector. Additionally, the identified SNPs could be employed as biomarkers of disease resistance in P. leopardus, which will facilitate the selective breeding of grouper.

LysM-containing proteins function in the resistance of Litopenaeus vannamei against Vibrio parahaemolyticus infection

Lysin motif (LysM) is a functional domain that can bind to peptidoglycans, chitin and their derivatives. The LysM-containing proteins participate in multiple biological processes, such as the hydrolysis of bacterial cell walls and the perception of PAMPs in plants and high animals. In the present study, two genes encoding LysM-containing proteins, designated as LvLysM1 and LvLysM2, were identified in the Pacific white shrimp, Litopenaeus vannamei, and their functions during Vibrio infection were analyzed. The open-reading frame (ORF) of LvLysM1 was 795 bp, only encoding a LysM domain at the N-terminal region. The ORF of LvLysM2 was 834 bp, encoding a LysM domain at the central region and a transmembrane region at the C-terminal region. Both LvLysM1 and LvLysM2 were widely transcribed in all tested shrimp tissues. Enzyme-linked immunosorbent assay (ELISA) showed that the recombinant protein of LvLysM2 could bind to different bacterial polysaccharides, while LvLysM1 showed no direct binding activity. The transcripts of LvLysMs in gills increased significantly after infection with Vibrio parahaemolyticus. When LvLysM1 or LvLysM2 was knocked down by dsRNA, the mortality of shrimp was significantly increased after infection with Vibrio parahaemolyticus. Interestingly, some SNPs existed in these two genes were apparently correlated with the VpAHPND resistance of shrimp. These results suggested that LvLysM1 and LvLysM2 might contribute to the disease resistance of shrimp. The data provide new knowledge about the function of LysM-containing proteins in shrimp and potential genetic markers for disease resistance breeding.

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.

Combination of GWAS and FST-based approaches identified loci associated with economic traits in sugarcane

Sugarcane is a globally important plant for both sugar and biofuel production. Although conventional breeding has played an important role in increasing the productivity of sugarcane, it takes a long time to achieve breeding goals such as high yield and resistant to diseases. Molecular breeding, including marker-assisted breeding and genomic selection, can accelerate genetic improvement by selecting elites at the seedling stage with DNA markers. However, only a few DNA markers associated with important traits were identified in sugarcane. The purpose of this study was to identify DNA markers associated with sugar content, stalk diameter, and sugarcane top borer resistance. The sugarcane samples with trait records were genotyped using the restriction site-associated DNA sequencing (RADseq) technology. Using F_ST analysis and genome-wide association study (GWAS), a total of 9, 23 and 9 DNA variants (single nucleotide polymorphisms (SNPs)/insertions and deletions (indels)) were associated with sugar content, stalk diameter, and sugarcane top borer resistance, respectively. The identified genetic variants were on different chromosomes, suggesting that these traits are complex and determined by multiple genetic factors. These DNA markers identified by both approaches have the potential to be used in selecting elite clones at the seeding stage in our sugarcane breeding program to accelerate genetic improvement. Certainly, it is essential to verify the reliability of the identified DNA markers associated with traits before they are used in molecular breeding in other populations.

Genetic diversity and genome-scale population structure of wild Indian major carp, Labeo catla (Hamilton, 1822), revealed by genotyping-by-sequencing

Labeo catla (catla) is the second most commercially important and widely cultured Indian major carp (IMC). It is indigenous to the Indo-Gangetic riverine system of India and the rivers of Bangladesh, Nepal, Myanmar, and Pakistan. Despite the availability of substantial genomic resources in this important species, detailed information on the genome-scale population structure using SNP markers is yet to be reported. In the present study, the identification of genome-wide single nucleotide polymorphisms (SNPs) and population genomics of catla was undertaken by re-sequencing six catla populations of riverine origin from distinct geographical regions. DNA isolated from 100 samples was used to perform genotyping-by-sequencing (GBS). A published catla genome with 95% genome coverage was used as the reference for mapping reads using BWA software. From a total of 472 million paired-end (150 × 2 bp) raw reads generated in this study, we identified 10,485 high-quality polymorphic SNPs using the STACKS pipeline. Expected heterozygosity (He) across the populations ranged from 0.162 to 0.20, whereas observed heterozygosity (Ho) ranged between 0.053 and 0.06. The nucleotide diversity (π) was the lowest (0.168) in the Ganga population. The within-population variation was found to be higher (95.32%) than the among-population (4.68%) variation. However, genetic differentiation was observed to be low to moderate, with F_st values ranging from 0.020 to 0.084, and the highest between Brahmani and Krishna populations. Bayesian and multivariate techniques were used to further evaluate the population structure and supposed ancestry in the studied populations using the structure and discriminant analysis of principal components (DAPC), respectively. Both analyses revealed the existence of two separate genomic clusters. The maximum number of private alleles was observed in the Ganga population. The findings of this study will contribute to a deeper understanding of the population structure and genetic diversity of wild populations of catla for future research in fish population genomics.

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.

Viral nervous necrosis resistance in gilthead sea bream (Sparus aurata) at the larval stage: heritability and accuracy of genomic prediction with different training and testing settings

BACKGROUND

The gilthead sea bream (Sparus aurata) has long been considered resistant to viral nervous necrosis (VNN), until recently, when significant mortalities caused by a reassortant nervous necrosis virus (NNV) strain were reported. Selective breeding to enhance resistance against NNV might be a preventive action. In this study, 972 sea bream larvae were subjected to a NNV challenge test and the symptomatology was recorded. All the experimental fish and their parents were genotyped using a genome-wide single nucleotide polymorphism (SNP) array consisting of over 26,000 markers.

RESULTS

Estimates of pedigree-based and genomic heritabilities of VNN symptomatology were consistent with each other (0.21, highest posterior density interval at 95% (HPD95%): 0.1-0.4; 0.19, HPD95%: 0.1-0.3, respectively). The genome-wide association study suggested one genomic region, i.e., in linkage group (LG) 23 that might be involved in sea bream VNN resistance, although it was far from the genome-wide significance threshold. The accuracies (r) of the predicted estimated breeding values (EBV) provided by three Bayesian genomic regression models (Bayes B, Bayes C, and Ridge Regression) were consistent and on average were equal to 0.90 when assessed in a set of cross-validation (CV) procedures. When genomic relationships between training and testing sets were minimized, accuracy decreased greatly (r = 0.53 for a validation based on genomic clustering, r = 0.12 for a validation based on a leave-one-family-out approach focused on the parents of the challenged fish). Classification of the phenotype using the genomic predictions of the phenotype or using the genomic predictions of the pedigree-based, all data included, EBV as classifiers was moderately accurate (area under the ROC curve 0.60 and 0.66, respectively).

CONCLUSIONS

The estimate of the heritability for VNN symptomatology indicates that it is feasible to implement selective breeding programs for increased resistance to VNN of sea bream larvae/juveniles. Exploiting genomic information offers the opportunity of developing prediction tools for VNN resistance, and genomic models can be trained on EBV using all data or phenotypes, with minimal differences in classification performance of the trait phenotype. In a long-term view, the weakening of the genomic ties between animals in the training and test sets leads to decreased genomic prediction accuracies, thus periodical update of the reference population with new data is mandatory.

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.

Silencing Asian Seabass gab3 Inhibits Nervous Necrosis Virus Replication

The nervous necrosis virus (NNV) causes the viral nervous necrosis (VNN) disease in aquatic animals and has been a major threat in aquaculture. Thus, it is essential for the development of a prevention method to minimize economic losses caused by NNV such as the identification of NNV resistance genes and application of these genes in molecular breeding to increase disease resistance. gab3 is an important NNV resistance gene in Asian seabass. However, the mechanism of gab3 in NNV resistance has not been elucidated. In this study, knockdown of gab3 in NNV-infected Asian seabass cells resulted in a significant decrease in viral RNA and virus titers. Knockout of gab3 in zebrafish led to an increased survival rate and resistant time after NNV infection. Cellular localization of the GAB3 and NNV by immunofluorescence staining showed that the GAB3 was translocated from the nucleus to the cytoplasm, and finally reached the cell membrane of SB cells after 48 h post NNV infection. Our study suggests that gab3 plays an important role in NNV replication and silencing gab3 can inhibit virus replication.

Genome-wide association study of VHSV-resistance trait in Paralichthys olivaceus

In flounder aquaculture, selective breeding plays a vital role in the development of disease-resistant traits and animals with high growth rates. Moreover, superior animals are required to achieve high profits. Unlike growth-related traits, disease-resistant experiments need to be conducted in a controlled environment, as the improper measurement of traits often leads to low genetic correlation and incorrect estimation of breeding values. In this study, viral hemorrhagic septicemia virus (VHSV) resistance was studied using a genome-wide association study (GWAS), and the genetic parameters were estimated. Genotyping was performed using a high-quality 70 K single nucleotide polymorphism (SNP) Affymetrix® Axiom® myDesign™ Genotyping Array of olive flounder. A heritability of ∼0.18 for resistance to VHSV was estimated using genomic information of the fish. According to the GWAS, significant SNPs were detected in chromosomes 21, 24, and contig AGQT02032065.1. Three SNPs showed significance at the genome-wide level (p < 1 × 10^-6), while others showed significance above the suggestive cutoff (p < 1 × 10^-4). The 3% phenotypic variation was explained by the highest significant SNP, named AX-419319631. Of the important genes for disease resistance, SNPs were associated with plcg1, epha4, clstn2, pik3cb, hes6, meis3, prx6, cep164, siae, and kirrel3b. Most of the genes associated with these SNPs have been previously reported with respect to viral entry, propagation, and immune mechanisms. Therefore, our study provides helpful information regarding VHSV resistance in olive flounder, which can be used for breeding applications.

Genotyping, the Usefulness of Imputation to Increase SNP Density, and Imputation Methods and Tools

Imputation has become a standard practice in modern genetic research to increase genome coverage and improve accuracy of genomic selection and genome-wide association study as a large number of samples can be genotyped at lower density (and lower cost) and, imputed up to denser marker panels or to sequence level, using information from a limited reference population. Most genotype imputation algorithms use information from relatives and population linkage disequilibrium. A number of software for imputation have been developed originally for human genetics and, more recently, for animal and plant genetics considering pedigree information and very sparse SNP arrays or genotyping-by-sequencing data. In comparison to human populations, the population structures in farmed species and their limited effective sizes allow to accurately impute high-density genotypes or sequences from very low-density SNP panels and a limited set of reference individuals. Whatever the imputation method, the imputation accuracy, measured by the correct imputation rate or the correlation between true and imputed genotypes, increased with the increasing relatedness of the individual to be imputed with its denser genotyped ancestors and as its own genotype density increased. Increasing the imputation accuracy pushes up the genomic selection accuracy whatever the genomic evaluation method. Given the marker densities, the most important factors affecting imputation accuracy are clearly the size of the reference population and the relationship between individuals in the reference and target populations.

Genomic Predictions of Phenotypes and Pseudo-Phenotypes for Viral Nervous Necrosis Resistance, Cortisol Concentration, Antibody Titer and Body Weight in European Sea Bass

Simple Summary

Selective breeding programs based on genomic data are still not a common practice in aquaculture, although genomic selection has been widely demonstrated to be advantageous when trait phenotyping is a difficult task. In this study, we investigated the accuracy of predicting the phenotype and the estimated breeding value (EBV) of three Bayesian models and a Random Forest algorithm exploiting the information of a genome-wide SNP panel for European sea bass. The genomic predictions were developed for mortality caused by viral nervous necrosis, post-stress cortisol concentration, antibody titer against nervous necrosis virus and body weight. Selective breeding based on genomic data is a possible option for improving these traits while overcoming difficulties related to individual phenotyping of the investigated traits. Our results evidenced that the EBV used as a pseudo-phenotype enhances the predictive performances of genomic models, and that EBV can be predicted with satisfactory accuracy. The genomic prediction of the EBV for mortality might also be used to classify the phenotype for the same trait.

Abstract

In European sea bass (Dicentrarchus labrax L.), the viral nervous necrosis mortality (MORT), post-stress cortisol concentration (HC), antibody titer (AT) against nervous necrosis virus and body weight (BW) show significant heritability, which makes selective breeding a possible option for their improvement. An experimental population (N = 650) generated by a commercial broodstock was phenotyped for the aforementioned traits and genotyped with a genome-wide SNP panel (16,075 markers). We compared the predictive accuracies of three Bayesian models (Bayes B, Bayes C and Bayesian Ridge Regression) and a machine-learning method (Random Forest). The prediction accuracy of the EBV for MORT was approximately 0.90, whereas the prediction accuracies of the EBV and the phenotype were 0.86 and 0.21 for HC, 0.79 and 0.26 for AT and 0.71 and 0.38 for BW. The genomic prediction of the EBV for MORT used to classify the phenotype for the same trait showed moderate classification performance. Genome-wide association studies confirmed the polygenic nature of MORT and demonstrated a complex genetic structure for HC and AT. Genomic predictions of the EBV for MORT could potentially be used to classify the phenotype of the same trait, though further investigations on a larger experimental population are needed.

Effects of rrm1 on NNV Resistance Revealed by RNA-seq and Gene Editing

Viral nervous necrosis (VNN) disease caused by the nervous necrosis virus (NNV) is a major disease, leading to a huge economic loss in aquaculture. Previous GWAS and QTL mapping have identified a major QTL for NNV resistance in linkage group 20 in Asian seabass. However, no causative gene for NNV resistance has been identified. In this study, RNA-seq from brains of Asian seabass fingerlings challenged with NNV at four time points (5, 10, 15 and 20 days post-challenge) identified 1228, 245, 189 and 134 DEGs, respectively. Eight DEGs, including rrm1, were located in the major QTL for NNV resistance. An association study in 445 survived and 608 dead fingerlings after NNV challenge revealed that the SNP in rrm1 were significantly associated with NNV resistance. Therefore, rrm1 was selected for functional analysis, as a candidate gene for NNV resistance. The expression of rrm1 was significantly increased in the gill, liver, spleen and muscle, and was suppressed in the brain, gut and skin after NNV challenge. The rrm1 protein was localized in the nuclear membrane. Over-expression of rrm1 significantly decreased viral RNA and titer in NNV-infected Asian seabass cells, whereas knock-down of rrm1 significantly increased viral RNA and titer in NNV-infected Asian seabass cells. The rrm1 knockout heterozygous zebrafish was more susceptible to NNV infection. Our study suggests that rrm1 is one of the causative genes for NNV resistance and the SNP in the gene may be applied for accelerating genetic improvement for NNV resistance.

Selection of functional EPHB2 genotypes from ENU mutated grass carp treated with GCRV

Background

N-ethyl-N-nitrosourea (ENU) mutagenesis is a useful method for the genetic engineering of plants, and the production of functional mutants in animal models including mice and zebrafish. Grass carp reovirus (GCRV) is a haemorrhagic disease of grass carp which has caused noteworthy losses in fingerlings over the last few years. To overcome this problem, we used ENU mutant grass carp in an attempt to identify functional resistance genes for future hereditary rearing projects in grass carp.

Results

This study used ENU-mutated grass carp to identify genetic markers associated with resistance to the haemorrhagic disease caused by GCRV. Bulked segregant analysis (BSA) was performed on two homozygous gynogenetic ENU grass carp groups who were susceptible or resistant to GCRV. This analysis identified 466,162 SNPs and 197,644 InDels within the genomes of these mixed pools with a total of 170 genes annotated in the associated region, including 49 genes with non-synonymous mutations at SNP sites and 25 genes with frame shift mutations at InDel sites. Of these 170 mutated genes, 5 randomly selected immune-related genes were shown to be more strongly expressed in the resistant group as compared to the susceptible animals. In addition, we found that one immune-related gene, EPHB2, presented with two heterozygous SNP mutations which altered the animal’s responded to GCRV disease. These SNPs were found in the intron region of EPHB2 at positions 5859 (5859^G > A) and 5968 (5968^G > A) and were significantly (p = 0.002, 0.003) associated with resistance to GCRV. These SNP sites were also shown to correlate with the GCRV-resistant phenotype in these ENU grass carp. We also evaluated the mortality of the different ENU fish genotypes in response to GCRV and the SNPs in EPHB2. The outcomes of these evaluations will be useful in future selections of GCRV-resistant genes for genetic breeding in grass carp.

Conclusion

Our results provide a proof of concept for the application of BSA-sequence analysis in detecting genes responsible for specific functional genotypes and may help to develop better methods for marker-assisted selection, especially for disease resistance in response to GCRV.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07858-x.

Effects of Ocean Acidification on Transcriptomes in Asian Seabass Juveniles

Ocean acidification is changing the fate of marine organisms. It is essential to predict the biological responses and evolutionary processes driven by ocean acidification, to maintain the equilibrium of the marine ecosystem and to facilitate aquaculture. However, how marine organisms, particularly the marine fish species, respond to ocean acidification, is still poorly understood. Consequences of ocean acidification on finfish aquaculture are largely not well known. We studied the effects of ocean acidification for 7 days on growth, behaviour and gene expression profiles in the brain, gill and kidney of Asian seabass juveniles. Results showed that growth and behaviour were not affected by short-term ocean acidification. We found tissue-specific differentially expressed genes (DEGs) involving many molecular processes, such as organ development, growth, muscle development, ion homeostasis and neurogenesis and development, as well as behaviours. Most of the DEGs, which were functionally enriched in ion homeostasis, were related to calcium transport, followed by sodium/potassium channels. We found that genes associated with neurogenesis and development were significantly enriched, implying that ocean acidification has also adversely affected the neural regulatory mechanism. Our results indicate that although the short-term ocean acidification does not cause obvious phenotypic and behavioural changes, it causes substantial changes of gene expressions in all three analysed tissues. All these changes of gene expressions may eventually affect physiological fitness. The DEGs identified here should be further investigated to discover DNA markers associated with adaptability to ocean acidification to improve fish's capability to adapt to ocean acidification.

Challenges and Solutions to Viral Diseases of Finfish in Marine Aquaculture

Aquaculture is the fastest food-producing sector in the world, accounting for one-third of global food production. As is the case with all intensive farming systems, increase in infectious diseases has adversely impacted the growth of marine fish farming worldwide. Viral diseases cause high economic losses in marine aquaculture. We provide an overview of the major challenges limiting the control and prevention of viral diseases in marine fish farming, as well as highlight potential solutions. The major challenges include increase in the number of emerging viral diseases, wild reservoirs, migratory species, anthropogenic activities, limitations in diagnostic tools and expertise, transportation of virus contaminated ballast water, and international trade. The proposed solutions to these problems include developing biosecurity policies at global and national levels, implementation of biosecurity measures, vaccine development, use of antiviral drugs and probiotics to combat viral infections, selective breeding of disease-resistant fish, use of improved diagnostic tools, disease surveillance, as well as promoting the use of good husbandry and management practices. A multifaceted approach combining several control strategies would provide more effective long-lasting solutions to reduction in viral infections in marine aquaculture than using a single disease control approach like vaccination alone.

Identification of candidate SNPs and genes associated with anti-RGNNV using GWAS in the red-spotted grouper, Epinephelus akaara

The red-spotted grouper, Epinephelus akaara, has been cultured widely in China, and in several countries of Southeast Asia, due to its important economic value. However, in recent years the outbreak of disease caused by red-spotted grouper nervous necrosis virus (RGNNV) has caused mass mortality in the stage of the grouper lifecycle from fry to juvenile, resulting in considerable economic loss in commercial aquaculture. However, the molecular mechanism underlying anti-RGNNV infection in red-spotted grouper has never been fully understood. To identify the anti-RGNNV related markers and candidate genes, we performed a genome-wide association study (GWAS) on a natural population of 100 individuals for a full-genome screen of the red-spotted grouper. In this research, 36,311 single, high quality nucleotide polymorphisms (SNPs) were developed. Two significantly associated SNPs and three suggestively associated SNPs were identified at the genome level. From these identified SNPs, five candidate genes were annotated: EPHA7, Osbpl2, GPC5, CDH4 and Pou3f1. These genes are involved in nervous system development, retinal formation, and lipid metabolism regulation. In combination with studies on the characteristics of NNV infection, it was speculated that in the fry stage of the grouper lifecycle, the immune system is not fully developed. Therefore, improved resistance to RGNNV may come through regulating nervous system development or lipid metabolism related pathways. In addition, the genotypes of SNPs associated with disease-resistant traits were analyzed. The markers and genes obtained in this study may facilitate a marker-assisted selection for red-spotted grouper aiming at disease resistance to RGNNV.

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.

Identification and characterization of cystatin B from black rockfish, Sebastes schlegelii, indicating its potent immunological importance

Cystatins represent a large superfamily of proteins involved in the competitive reversible inhibition of C1 class cysteine proteases. Plant-derived papain proteases and cysteine cathepsins are the major cysteine proteases that interact with cystatins. The cystatin superfamily can be further classified into three groups: stefins, cystatins, and kininogens. Among these, cystatin B is categorized under stefins. Cystatin B lacks a signal sequence, disulfide bonds, and carbohydrate groups. However, it contains the conserved cystatin family signature, including a single cystatin-like domain, cysteine protease inhibitory signature concealing pentapeptide (QXVXG) consensus sequence, and two conserved neighboring glycine (⁸GG⁹) residues at the N-terminal. In the current study, a member of cystatin B was identified from Korean black rockfish (Sebastes schlegeli) using a cDNA database and designated as RfCytB. The full-length cDNA of RfCytB was 573 bp long, with a coding region of 294 bp. The 5'-untranslated region (UTR) comprised 55 bp, and the 263-bp-long 3'-UTR included a polyadenylation signal sequence and a poly-A tail. The coding sequence encodes a polypeptide comprising 97 amino acids, with a predicted molecular weight of 11 kDa and theoretical isoelectric point of 6.3. RfCytB shared homology features with similar molecules from other teleost and vertebrate species, and was clustered with Cystatin family 1 in our phylogenetic reconstruction. RfCytB was ubiquitously expressed in all tissue types of healthy animals, with the highest levels of expression observed in gill and spleen. Temporal expression of RfCytB displayed significant up-regulation upon infection with Aeromonas salmonicida. Recombinantly expressed RfCytB showed a concentration-dependent inhibitory activity towards papain, with a high thermal stability. Transient expression of RfCytB in LPS activated murine macrophages, thereby inducing the expression of genes related to pro-inflammatory conditions, such as iNOS and TNF α. These results provide evidence for its protease inhibitory and immunity relevant roles in hosts.

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.

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

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

Validation of growth-related quantitative trait loci markers in different Exopalaemon carinicauda families for marker-assisted selection

We detected growth-related QTL and associated markers from the backcross population of Exopalaemon carinicauda in the previous study. Based on our previous study, the 47 SNP markers associated with candidate growth trait QTL were selected to analyze the association between these markers and body weight (BW), body length and abdominal segment length traits in four different populations including wild population, a full-sib family, a half-sib family and a backcross population for evaluating their potential application of marker-assisted selection in E. carinicauda. The general linear model (GLM) and mixed linear model were applied and the associations between SNP loci and three growth-related traits verified. The results showed that the Marker79268 and Marker100644 were significantly associated with the BW trait in more than three populations by the GLM method. The Marker100644 was significantly associated with BW in the full-sib family, half-sib family and backcross populations by the GLM and mixed linear model methods. Our findings will provide useful SNP markers to go forward to improve growth performance through more refined marker-assisted selection in E. carinicauda.

Quantitative Trait Loci (QTL) Mapping

Quantitative trait loci (QTL) are genetic regions that influence phenotypic variation of a complex trait, often through genetic interactions with each other and the environment. These are commonly identified through a statistical genetic analysis known as QTL mapping. Here, I present a step-by-step, practical approach to QTL mapping along with a sample data file. I focus on methods commonly used and discoveries that have been made in fishes, and utilize a multiple QTL mapping (MQM) approach in the free software package R/qtl.

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.

First Genome-wide Association Analysis for Growth Traits in the Largest Coral Reef-Dwelling Bony Fishes, the Giant Grouper (Epinephelus lanceolatus)

The giant grouper, Epinephelus lanceolatus, is the largest coral reef-dwelling bony fish species. However, despite extremely fast growth performance and the considerable economic importance in this species, its genetic regulation of growth remains unknown. Here, we performed the first genome-wide association study (GWAS) for five growth traits in 289 giant groupers using 42,323 single nucleotide polymorphisms (SNPs) obtained by genotyping-by-sequencing (GBS). We identified a total of 36 growth-related SNPs, of which 11 SNPs reached a genome-wide significance level. The phenotypic variance explained by these SNPs varied from 7.09% for body height to 18.42% for body length. Moreover, 22 quantitative trait loci (QTLs) for growth traits, including nine significant QTLs and 13 suggestive QTLs, were found on multiple chromosomes. Interestingly, the QTL (LG17: 6934451) was shared between body weight and body height, while two significant QTLs (LG7: 22596399 and LG15: 11877836) for body length were consistent with the associated regions of total length at the genome-wide suggestive level. Eight potential candidate genes close to the associated SNPs were selected for expression analysis, of which four genes (phosphatidylinositol transfer protein cytoplasmic 1, protein tyrosine phosphatase receptor type E, alpha/beta hydrolase domain-containing protein 17C, and vascular endothelial growth factor A-A) were differentially expressed and involved in metabolism, development, response stress, etc. This study improves our understanding of the complex genetic architecture of growth in the giant grouper. The results contribute to the selective breeding of grouper species and the conservation of coral reef fishes.

Functional differentiation of three phosphatidylinositol 3-kinase catalytic subunit alpha (PIK3CA) in response to Vibrio anguillarum infection in turbot (Scophthalmus maximus)

PIK3CA has been extensively investigated from its molecular mechanism perspective and association with its mutations in different types of cancers. However, little has been reported regarding the pathological significance of PIK3CA expression in teleost. Here, in our present study, three PIK3CA genes termed SmPIK3CAa, SmPIK3CAb and SmPIK3CA-like were firstly identified in the genome of turbot S. maximus. Although these three genes located in different chromosomes, all of them share the same five domains. Phylogenetic and synteny analysis indicated that SmPIK3CAa, SmPIK3CAb and SmPIK3CA-like were three paralogs that may originate from duplication of the same ancestral PIK3CA gene. Subcellular localization analysis confirmed the cytoplasm distribution of these three paralogs. All three SmPIK3CA were ubiquitously expressed in examined tissues in turbot, with the higher expression levels in immune-related tissues such as blood, spleen, kidney, gills and intestines. Upon Vibrio anguillarum challenge, SmPIK3CAa and SmPIK3CA-like transcripts were significantly induced in spleen, intestine and blood despite of differential expression levels and responsive time points. Additionally, individuals in resistant group showed significantly higher expression level of both two genes than in the susceptible group. Moreover, four SNPs (102, 2530, 3027 and 3060) and one haplotype (Hap2) located in exon region of SmPIK3CA-like were identified and confirmed to be associated with V. anguillarum resistance in turbot by association analysis in different populations. Taken together, these results suggested that functional differentiation occurred in three SmPIK3CA paralogs with Vibrio anguillarum resistance and SmPIK3CAa and SmPIK3CA-like probable play potential roles in innate immune response to pathogenic invasions in turbot.

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

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

QTL Mapping for Red Blotches in Malaysia Red Tilapia (Oreochromis spp.)

Body color is an interesting economic trait in fish. Red tilapia with red blotches may decrease its commercial values. Conventional selection of pure red color lines is a time-consuming and labor-intensive process. To accelerate selection of pure lines through marker-assisted selection, in this study, double-digest restriction site-associated DNA sequencing (ddRAD-seq) technology was applied to genotype a full-sib mapping family of Malaysia red tilapia (Oreochromis spp.) (N = 192). Genome-wide significant quantitative trait locus (QTL)-controlling red blotches were mapped onto two chromosomes (chrLG5 and chrLG15) explaining 9.7% and 8.2% of phenotypic variances by a genome-wide association study (GWAS) and linkage-based QTL mapping. Six SNPs from the chromosome chrLG5 (four), chrLG15 (one), and unplaced supercontig GL831288-1 (one) were significantly associated to the red blotch trait in GWAS analysis. We developed nine microsatellite markers and validated significant correlations between genotypes and blotch data (p < 0.05). Our study laid a foundation for exploring a genetic mechanism of body colors and carrying out genetic improvement for color quality in tilapia.

Constructing High-Density Genetic Maps and Developing Sexing Markers in Northern Snakehead (Channa argus)

High-density genetic maps are essential for mapping QTL, improving genome assembly, comparative genomics, and studying sex chromosome evolution. The northern snakehead (Channa argus) is an economically important foodfish species with significant sexual dimorphism, where the males grow much faster and bigger than the females. However, to date, the sex determination pattern is still not clear, limiting identification of sex chromosomes, even sex determination genes and development of monosex populations that are valuable for both sex evolution of vertebrates and aquaculture practices. Here, a sex-averaged map and two sex-specific genetic maps were constructed with 2974, 2323, and 2338 SNPs, respectively. Little difference was observed in the pattern of sex-specific recombination between female- and male-specific genetic maps. Genome scan identified a major locus for sex determination at LG16. Females and males are, respectively, homogametic and heterogametic, suggesting an XY sex determination system for this species. By resequencing genomes, InDels in the sex-associated QTL region were discovered and used for developing sex-specific PCR assays for fast sexing of snakehead. These high-density genetic maps provide useful resources for future genomic studies in snakehead and its related species. The PCR assays for sexing are of importance in developing all male populations for aquaculture.

Whole-genome resequencing from bulked-segregant analysis reveals gene set based association analyses for the Vibrio anguillarum resistance of turbot (Scophthalmus maximus)

Many achievements have been made to develop quantitative trait loci (QTLs) and gene-associated single nucleotide polymorphisms (SNPs) to facilitate practical marker-assisted selection (MAS) in aquatic animals. However, the systematic studies of SNPs associated with extreme threshold traits were poor in populations lacking of parental genomic information. Coupling next generation sequencing with bulked segregant analysis (BSA) should allow identification of numerous associated SNPs with extreme phenotypes. In the present study, using combination of SNP frequency difference and Euclidean distance, we conducted linkage analysis of SNPs located in genes involved in immune responses, and identified markers associated with Vibrio anguillarum resistance in turbot (Scophthalmus maximus). A total of 221 SNPs was found as candidate SNPs between resistant and susceptible individuals. Among these SNPs, 35 loci located in immune related genes were genotyped in verification population and 7 of them showed significant association with V. anguillarum resistance in both alleles and genotypes (P < 0.05). Among these 7 genes, PIK3CA-like, CYLD, VCAM1, RhoB and RhoGEF are involved in PI3K/Akt/mTOR pathway and NF-κB pathway, which influence the efficiency of bacteria entering the host and inflammation. SNP-SNP interaction analysis was performed by generalized multifactor dimensionality reduction (GMDR). The combination of SNP loci in RhoB, PIK3CA-like and ADCY3 showed a significant effect on V. anguillarum resistance with the verification rate in the sequencing population up to 70.8%. Taken all, our findings demonstrated the feasibility of BSA-seq approach in identifying genes responsible for the extreme phenotypes and will aid in performing MAS in turbot.

Predicting Growth Traits with Genomic Selection Methods in Zhikong Scallop (Chlamys farreri)

Selective breeding is a common and effective approach for genetic improvement of aquaculture stocks with parental selection as the key factor. Genomic selection (GS) has been proposed as a promising tool to facilitate selective breeding. Here, we evaluated the predictability of four GS methods in Zhikong scallop (Chlamys farreri) through real dataset analyses of four economical traits (e.g., shell length, shell height, shell width, and whole weight). Our analysis revealed that different GS models exhibited variable performance in prediction accuracy depending on genetic and statistical factors, but non-parametric method, including reproducing kernel Hilbert spaces regression (RKHS) and sparse neural networks (SNN), generally outperformed parametric linear method, such as genomic best linear unbiased prediction (GBLUP) and BayesB. Furthermore, we demonstrated that the predictability relied mainly on the heritability regardless of GS methods. The size of training population and marker density also had considerable effects on the predictive performance. In practice, increasing the training population size could better improve the genomic prediction than raising the marker density. This study is the first to apply non-linear model and neural networks for GS in scallop and should be valuable to help develop strategies for aquaculture breeding programs.

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

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

Genomic Selection Using BayesCπ and GBLUP for Resistance Against Edwardsiella tarda in Japanese Flounder (Paralichthys olivaceus)

The Japanese flounder is one of the most widely farmed economic flatfish species throughout eastern Asia including China, Korea, and Japan. Edwardsiella tarda is a major species of pathogenic bacteria that causes ascites disease and, consequently, a huge economy loss for Japanese flounder farming. After generation selection, traditional breeding methods can hardly improve the E. tarda resistance effectively. Genomic selection is an effective way to predict the breeding potential of parents and has rarely been used in aquatic breeding. In this study, we chose 931 individuals from 90 families, challenged by E. tarda from 2013 to 2015 as a reference population and 71 parents of these families as selection candidates. 1,934,475 markers were detected via genome sequencing and applied in this study. Two different methods, BayesCπ and GBLUP, were used for genomic prediction. In the reference population, two methods led to the same accuracy (0.946) and Pearson's correlation results between phenotype and genomic estimated breeding value (GEBV) of BayesCπ and GBLUP were 0.912 and 0.761, respectively. In selection candidates, GEBVs from two methods were highly similar (0.980). A comparison of GEBV with the survival rate of families that were structured by selection candidates showed correlations of 0.662 and 0.665, respectively. This study established a genomic selection method for the Japanese flounder and for the first time applied this to E. tarda resistance breeding.

Molecular characterization, expression and immune functions of cystatin B in Japanese flounder (Paralichthys olivaceus)

Cystatin B is an intracellular inhibitor that regulates the activities of cysteine proteases. In this study, cystatin B in Japanese flounder (Paralichthys olivaceus) was characterized and its immune function was analyzed. This gene had a high similarity with the sequence of cystatin B in other fish species, and the derived peptide shared typical features of cystatin proteins including the QXVXG motif. The results of quantitative real-time PCR showed that cystatin B mRNA was constitutively expressed in all examined tissues, with the highest level in gill. The stimulations of lipopolysaccharide, peptidoglycan and polyinosinic-polycytidylic acid effectively increased the expression level of cystatin B mRNA. Functional analysis implied that the recombinant P. olivaceus cystatin B purified from Escherichia coli had cysteine protease inhibitory activity and could inhibit bacterial growth by binding to bacteria. Furthermore, we found that P. olivaceus cystatin B had no effects on the expression of inflammatory factors cytokines tumor necrosis factor α, interleukin 10, interleukin 1β and interferon γ. These results indicate that cystatin B of P. olivaceus is potentially involved in immune responses against invading microbial pathogens, and provide a better understanding of the immune mechanisms of cystatins in teleosts.

Genome-wide association and genomic prediction of resistance to viral nervous necrosis in European sea bass (Dicentrarchus labrax) using RAD sequencing

BACKGROUND

European sea bass (Dicentrarchus labrax) is one of the most important species for European aquaculture. Viral nervous necrosis (VNN), commonly caused by the redspotted grouper nervous necrosis virus (RGNNV), can result in high levels of morbidity and mortality, mainly during the larval and juvenile stages of cultured sea bass. In the absence of efficient therapeutic treatments, selective breeding for host resistance offers a promising strategy to control this disease. Our study aimed at investigating genetic resistance to VNN and genomic-based approaches to improve disease resistance by selective breeding. A population of 1538 sea bass juveniles from a factorial cross between 48 sires and 17 dams was challenged with RGNNV with mortalities and survivors being recorded and sampled for genotyping by the RAD sequencing approach.

RESULTS

We used genome-wide genotype data from 9195 single nucleotide polymorphisms (SNPs) for downstream analysis. Estimates of heritability of survival on the underlying scale for the pedigree and genomic relationship matrices were 0.27 (HPD interval 95%: 0.14-0.40) and 0.43 (0.29-0.57), respectively. Classical genome-wide association analysis detected genome-wide significant quantitative trait loci (QTL) for resistance to VNN on chromosomes (unassigned scaffolds in the case of 'chromosome' 25) 3, 20 and 25 (P < 1e06). Weighted genomic best linear unbiased predictor provided additional support for the QTL on chromosome 3 and suggested that it explained 4% of the additive genetic variation. Genomic prediction approaches were tested to investigate the potential of using genome-wide SNP data to estimate breeding values for resistance to VNN and showed that genomic prediction resulted in a 13% increase in successful classification of resistant and susceptible animals compared to pedigree-based methods, with Bayes A and Bayes B giving the highest predictive ability.

CONCLUSIONS

Genome-wide significant QTL were identified but each with relatively small effects on the trait. Tests of genomic prediction suggested that incorporating genome-wide SNP data is likely to result in higher accuracy of estimated breeding values for resistance to VNN. RAD sequencing is an effective method for generating such genome-wide SNPs, and our findings highlight the potential of genomic selection to breed farmed European sea bass with improved resistance to VNN.

SNPs in the Toll1 receptor of Litopenaeus vannamei are associated with immune response

Tolls and Toll-like receptors (TLRs) are important regulators in the innate immune system and their genetic variations usually affect the host's susceptibility/resistance to pathogen infections. In this study, we report on the single nucleotide polymorphisms (SNPs) of Toll1 in Litopenaeus vannamei (LvToll1) and how this is associated with immune response. PCR-DGGE analysis revealed genetic polymorphisms in LvToll1 at both the genomic DNA (gDNA) and cDNA levels. Using high-throughput sequencing, 223 SNPs were identified at the gDNA level, of which 145 were non-synonymous SNP (nsSNP), with 3 nsSNPs having frequency over 1%. On the other hand, 60 SNPs were identified at the cDNA level including 38 nsSNPs and 4 nsSNPs with frequency over 1%. Upon challenging shrimps with Streptococcus iniae, Vibrio parahaemolyticus and white spot syndrome virus (WSSV), LvToll1 was shown to generate 6, 4 and 4 novel bands, respectively when analyzed with PCR-DGGE. Sequencing analysis of these bands showed that they contained 6, 4 and 2 nsSNPs, respectively. Moreover, the nsSNP C1526T was detected in S. iniae-resistant but not in susceptible shrimps. Most significantly, the C1526T mutation could shorten the α-helix of the LRR domain and was predicted to affect the function of LvToll1, indicating that SNP C1526T might be associated with shrimp's resistance to pathogen infections. In sum, our findings here reveal that the genetic polymorphisms of Toll receptor are linked with the immune response to pathogen infections in L. vannamei.

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.

QTL Mapping for Resistance to Iridovirus in Asian Seabass Using Genotyping-by-Sequencing

Identifying quantitative trait loci (QTL) for viral disease resistance is of particular importance in selective breeding programs of fish species. Genetic markers linked to QTL can be useful in marker-assisted selection (MAS) for elites resistant to specific pathogens. Here, we conducted a genome scan for QTL associated with Singapore grouper iridovirus (SGIV) resistance in an Asian seabass (Lates calcarifer) family, using a high-density linkage map generated with genotyping-by-sequencing. One genome-wide significant and three suggestive QTL were detected at LG21, LG6, LG13, and LG15, respectively. The phenotypic variation explained (PVE) by the four QTL ranged from 7.5 to 15.6%. The position of the most significant QTL at LG21 was located between 31.88 and 36.81 cM. The SNP marker (SNP130416) nearest to the peak of this QTL was significantly associated with SGIV resistance in an unrelated multifamily population. One candidate gene, MECOM, close to the peak of this QTL region, was predicted. Evidence of alternative splicing was observed for MECOM and one specific category of splicing variants was differentially expressed at 5 days post-SGIV infection. The QTL detected in this study are valuable resources and can be used in the selective breeding programs of Asian seabass with regard to resistance to SGIV.