Identification of Quantitative Trait Loci for Resistance to RSIVD in Red Sea Bream (Pagrus major)

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.

Scale Drop Disease Virus Associated Yellowfin Seabream (Acanthopagrus latus) Ascites Diseases, Zhuhai, Guangdong, Southern China: The First Description

Scale drop disease virus (SDDV), an emerging piscine iridovirus prevalent in farmed Asian seabass Lates calcarifer in Southeast Asia, was firstly scientifically descripted in Singapore in 2015. Here, an SDDV isolate ZH-06/20 was isolated by inoculating filtered ascites from diseased juvenile yellowfin seabream into MFF-1 cell. Advanced cytopathic effects were observed 6 days post-inoculation. A transmission electron microscopy examination confirmed that numerous virion particles, about 140 nm in diameter, were observed in infected MFF-1 cell. ZH-06/20 was further purified and both whole genome and virion proteome were determined. The results showed that ZH-06/20 was composed of 131,122 bp with 135 putative viral proteins and 113 of them were further detected by virion proteome. Western blot analysis showed that no (or weak) cross-reaction was observed among several major viral proteins between ZH-06/20 and ISKNV-like megalocytivirus. An artificial challenge showed that ZH-06/20 could cause 100% death to juvenile yellowfin seabream. A typical sign was characterized by severe ascites, but not scale drop, which was considerably different from SDD syndrome in Asian seabass. Collectively, SDDV was confirmed, for the first time, as the causative agent of ascites diseases in farmed yellowfin seabream. Our study offers useful information to better understanding SDDV-associated diseases in farmed fish.

Identification of the causative gene of a transparent phenotype of juvenile red sea bream Pagrus major

Deformities in cultured fish species may be genetic, and identifying causative genes is essential to expand production and maintain farmed animal welfare. We previously reported a genetic deformity in juvenile red sea bream, designated a transparent phenotype. To identify its causative gene, we conducted genome-wide linkage analysis and identified two single nucleotide polymorphisms (SNP) located on LG23 directly linked to the transparent phenotype. The scaffold on which the two SNPs were located contained two candidate genes, duox and duoxa, which are related to thyroid hormone synthesis. Four missense mutations were found in duox and one in duoxa, with that in duoxa showing perfect association with the transparent phenotype. The mutation of duoxa was suggested to affect the transmembrane structure and thyroid-related traits, including an enlarged thyroid gland and immature erythrocytes, and lower thyroxine (T₄) concentrations were observed in the transparent phenotype. The transparent phenotype was rescued by T₄ immersion. Loss-of-function of duoxa by CRISPR-Cas9 induced the transparent phenotype in zebrafish. Evidence suggests that the transparent phenotype of juvenile red sea bream is caused by the missense mutation of duoxa and that this mutation disrupts thyroid hormone synthesis. The newly identified missense mutation will contribute to effective selective breeding of red sea bream to purge the causative gene of the undesirable phenotype and improve seed production of red sea bream as well as provide basic information of the mechanisms of thyroid hormones and its related diseases in fish and humans.

Lack of ∆5 Desaturase Activity Impairs EPA and DHA Synthesis in Fish Cells from Red Sea Bream and Japanese Flounder

Long-chain (≥ C₂₀) polyunsaturated fatty acids (LC-PUFA), such as eicosapentaenoic acid (20:5n-3, EPA) and docosahexaenoic acid (22:6n-3, DHA), are necessary for human health and are obtained from marine fish-derived oils. Marine fish are LC-PUFA-rich animals; however, many of them require LC-PUFA for growth. Therefore, it is suggested that they do not have sufficient ability to biosynthesize LC-PUFA. To evaluate in vivo LC-PUFA synthetic activity in fish cells, fish-derived cell lines from red sea bream (Pagrus major, PMS and PMF), Japanese flounder (Paralichthys olivaceus, HINAE), and zebrafish (Danio rerio, BRF41) were incubated with n-3 fatty acids labeled by radioisotopes or stable isotopes, and then, n-3 PUFA were analyzed by thin-layer chromatography or liquid chromatography-mass spectrometry. Labeled EPA and DHA were biosynthesized from labeled α-linolenic acid (18:3n-3) in BRF41, whereas they were not detected in PMS, PMF, or HINAE cells. We next cloned the fatty acid desaturase 2 (Fads2) cDNAs from PMF cells and zebrafish, expressed in budding yeasts, and then analyzed the substrate specificities of enzymes. As a result, we found that Fads2 from PMF cells was a ∆6/∆8 desaturase. Collectively, our study indicates that cell lines from red sea bream and Japanese flounder were not able to synthesize EPA or DHA by themselves, possibly due to the lack of ∆5 desaturase activity. Furthermore, this study provides a sensitive and reproducible non-radioactive method for evaluating LC-PUFA synthesis in fish cells using a stable isotope and liquid chromatography-mass spectrometry.

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.

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.

PCR Detection and Phylogenetic Analysis of Megalocytivirus Isolates in Farmed Giant Sea Perch Lates calcarifer in Southern Taiwan

The Megalocytivirus genus includes three genotypes, red sea bream iridovirus (RSIV), infectious spleen and kidney necrosis virus (ISKNV), and turbot reddish body iridovirus (TRBIV), and has caused mass mortalities in various marine and freshwater fish species in East and Southeast Asia. Of the three genotypes, TRBIV-like megalocytivirus is not included in the World Organization for Animal Health (OIE)-reportable virus list because of its geographic restriction and narrow host range. In 2017, 39 cases of suspected iridovirus infection were isolated from fingerlings of giant sea perch (Lates calcarifer) cultured in southern Taiwan during megalocytivirus epizootics. Polymerase chain reaction (PCR) with different specific primer sets was undertaken to identify the causative agent. Our results revealed that 35 out of the 39 giant sea perch iridovirus (GSPIV) isolates were TRBIV-like megalocytiviruses. To further evaluate the genetic variation, the nucleotide sequences of major capsid protein (MCP) gene (1348 bp) from 12 of the 35 TRBIV-like megalocytivirus isolates were compared to those of other known. High nucleotide sequence identity showed that these 12 TRBIV-like GSPIV isolates are the same species. Phylogenetic analysis based on the MCP gene demonstrated that these 12 isolates belong to the clade II of TRBIV megalocytiviruses, and are distinct from RSIV and ISKNV. In conclusion, the GSPIV isolates belonging to TRBIV clade II megalocytiviruses have been introduced into Taiwan and caused a severe impact on the giant sea perch aquaculture industry.

Identification of Single Nucleotide Polymorphisms Related to the Resistance Against Acute Hepatopancreatic Necrosis Disease in the Pacific White Shrimp Litopenaeus vannamei by Target Sequencing Approach

Acute hepatopancreatic necrosis disease (AHPND) is a major bacterial disease in Pacific white shrimp Litopenaeus vannamei farming, which is caused by Vibrio parahaemolyticus. AHPND has led to a significant reduction of shrimp output since its outbreak. Selective breeding of disease-resistant broodstock is regarded as a key strategy in solving the disease problem. Understanding the relationship between genetic variance and AHPND resistance is the basis for marker-assisted selection in shrimp. The purpose of this study was to identify single nucleotide polymorphisms (SNPs) associated with the resistance against AHPND in L. vannamei. In this work, two independent populations were used for V. parahaemolyticus challenge and the resistant or susceptible shrimp were evaluated according to the survival time after Vibrio infection. The above two populations were genotyped separately by a SNP panel designed based on the target sequencing platform using a pooling strategy. The SNP panel contained 508 amplicons from DNA fragments distributed evenly along the genome and some immune-related genes of L. vannamei. By analyzing the allele frequency in the resistant and susceptible groups, 30 SNPs were found to be significantly associated with the resistance of the shrimp against V. parahaemolyticus infection (false discovery rate corrected at P < 0.05). Three SNPs were further validated by individual genotyping in all samples of population 1. Our study illustrated that target sequencing and pooling sequencing were effective in identifying the markers associated with economic traits, and the SNPs identified in this study could be used as molecular markers for breeding disease-resistant shrimp.

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.

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.

Major Histocompatibility Complex (MHC) Genes and Disease Resistance in Fish

Fascinating about classical major histocompatibility complex (MHC) molecules is their polymorphism. The present study is a review and discussion of the fish MHC situation. The basic pattern of MHC variation in fish is similar to mammals, with MHC class I versus class II, and polymorphic classical versus nonpolymorphic nonclassical. However, in many or all teleost fishes, important differences with mammalian or human MHC were observed: (1) The allelic/haplotype diversification levels of classical MHC class I tend to be much higher than in mammals and involve structural positions within but also outside the peptide binding groove; (2) Teleost fish classical MHC class I and class II loci are not linked. The present article summarizes previous studies that performed quantitative trait loci (QTL) analysis for mapping differences in teleost fish disease resistance, and discusses them from MHC point of view. Overall, those QTL studies suggest the possible importance of genomic regions including classical MHC class II and nonclassical MHC class I genes, whereas similar observations were not made for the genomic regions with the highly diversified classical MHC class I alleles. It must be concluded that despite decades of knowing MHC polymorphism in jawed vertebrate species including fish, firm conclusions (as opposed to appealing hypotheses) on the reasons for MHC polymorphism cannot be made, and that the types of polymorphism observed in fish may not be explained by disease-resistance models alone.

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.

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.

Identification, Characterization, and Mapping of a Novel SNP Associated with Body Color Transparency in Juvenile Red Sea Bream (Pagrus major)

We previously reported a body color deformity in juvenile red sea bream, which shows transparency in the juvenile stage because of delayed chromatophore development compared with normal individuals, and this finding suggested a genetic cause based on parentage assessments. To conduct marker-assisted selection to eliminate broodstock inheriting the causative gene, developing DNA markers associated with the phenotype was needed. We first conducted SNP mining based on AFLP analysis using bulked-DNA from normal and transparent individuals. One SNP was identified from a transparent-specific AFLP fragment, which significantly associated with transparent individuals. Two alleles (A/G) were observed in this locus, and the genotype G/G was dominantly observed in the transparent groups (97.1%) collected from several production lots produced from different broodstock populations. A few normal individuals inherited the G/G genotype (5.0%), but the A/A and A/G genotypes were dominantly observed in the normal groups. The homologs region of the SNP was searched using a medaka genome database, and intron 12 of the Nell2a gene (located on chromosome 6 of the medaka genome) was highly matched. We also mapped the red sea bream Nell2a gene on the previously developed linkage maps, and this gene was mapped on a male linkage group, LG4-M. The newly found SNP was useful in eliminating broodstock possessing the causative gene of the body color transparency observed in juvenile stage of red sea bream.