A genome-wide association study in catfish reveals the presence of functional hubs of related genes within QTLs for columnaris disease resistance

Multi-organ transcriptomic profiles and gene-regulation network underlying vibriosis resistance in tongue sole

Vibrio spp. are major pathogens responsible for mortality and disease in various marine aquaculture organisms. Effective disease control and genetic breeding strategies rely heavily on understanding host vibriosis resistance mechanisms. The Chinese tongue sole (Cynoglossus semilaevis) is economically vital but suffers from substantial mortalities due to vibriosis. Through continuous selective breeding, we have successfully obtained vibriosis-resistant families of this species. In this study, we conducted RNA-seq analysis on three organs, including liver, spleen and intestine from selected resistant and susceptible tongue soles. Additionally, we integrated these data with our previously published RNA-seq datasets of skin and gill, enabling the construction of organ-specific transcriptional profiles and a comprehensive gene co-expression network elucidating the differences in vibriosis resistance. Furthermore, we identified 12 modules with organ-specific functional implications. Overall, our findings provide a valuable resource for investigating the molecular basis of vibriosis resistance in fish, offering insights into target genes and pathways essential for molecular selection and genetic manipulation to enhance vibriosis resistance in fish breeding programs.

Elucidation of the genetic determination of clutch traits in Chinese local chickens of the Laiwu Black breed

BACKGROUND

Egg laying rate (LR) is associated with a clutch, which is defined as consecutive days of oviposition. The clutch trait can be used as a selection indicator to improve egg production in poultry breeding. However, little is known about the genetic basis of clutch traits. In this study, our aim was to estimate genetic parameters and identify quantitative trait single nucleotide polymorphisms for clutch traits in 399 purebred Laiwu Black chickens (a native Chinese breed) using a genome-wide association study (GWAS).

METHODS

In this work, after estimating the genetic parameters of age at first egg, body weight at first egg, LR, longest clutch until 52 week of age, first week when the longest clutch starts, last week when the longest clutch ends, number of clutches, and longest number of days without egg-laying until 52 week of age, we identified single nucleotide polymorphisms (SNPs) and potential candidate genes associated with clutch traits in Laiwu Black chickens. The restricted maximum likelihood method was used to estimate genetic parameters of clutch pattern in 399 Laiwu Black hens, using the GCTA software.

RESULTS

The results showed that SNP-based heritability estimates of clutch traits ranged from 0.06 to 0.59. Genotyping data were obtained from whole genome re-sequencing data. After quality control, a total of 10,810,544 SNPs remained to be analyzed. The GWAS revealed that 421 significant SNPs responsible for clutch traits were scattered on chicken chromosomes 1-14, 17-19, 21-25, 28 and Z. Among the annotated genes, NELL2, SMYD9, SPTLC2, SMYD3 and PLCL1 were the most promising candidates for clutch traits in Laiwu Black chickens.

CONCLUSION

The findings of this research provide critical insight into the genetic basis of clutch traits. These results contribute to the identification of candidate genes and variants. Genes and SNPs potentially provide new avenues for further research and would help to establish a framework for new methods of genomic prediction, and increase the accuracy of estimated genetic merit for egg production and clutch traits.

Dissection of the Genetic Basis of Resistance to Stem Rot in Cultivated Peanuts (Arachis hypogaea L.) through Genome-Wide Association Study

Peanut (Arachis hypogaea) is an important oilseed and cash crop worldwide, contributing an important source of edible oil and protein for human nutrition. However, the incidence of stem rot disease caused by Athelia rolfsii poses a major challenge to peanut cultivation, resulting in significant yield losses. In this study, a panel of 202 peanut accessions was evaluated for their resistance to stem rot by inoculating plants in the field with A. rolfsii-infested oat grains in three environments. The mean disease index value of each environment for accessions in subsp. fasitigiate and subsp. hypogaea showed no significant difference. Accessions from southern China displayed the lowest disease index value compared to those from other ecological regions. We used whole-genome resequencing to analyze the genotypes of the accessions and to identify significant SNPs associated with stem rot resistance through genome-wide association study (GWAS). A total of 121 significant SNPs associated with stem rot resistance in peanut were identified, with phenotypic variation explained (PVE) ranging from 12.23% to 15.51%. A total of 27 candidate genes within 100 kb upstream and downstream of 23 significant SNPs were annotated, which have functions related to recognition, signal transduction, and defense response. These significant SNPs and candidate genes provide valuable information for further validation and molecular breeding to improve stem rot resistance in peanut.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Applying genetic technologies to combat infectious diseases in aquaculture

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

A heterophil/lymphocyte-selected population reveals the phosphatase PTPRJ is associated with immune defense in chickens

Quantification of leukocyte profiles is among the simplest measures of animal immune function. However, the relationship between H/L ratio and innate immunity and the measure's utility as an index for heterophil function remains to be analyzed. Variants associated with H/L ratio were fine-mapped based on the resequencing of 249 chickens of different generations and an F2 segregating population generated by crossing selection and control lines. H/L ratio in the selection line was associated with a selective sweep of mutations in protein tyrosine phosphatase, receptor type J (PTPRJ), which affects proliferation and differentiation of heterophils through its downstream regulatory genes. The SNP downstream of PTPRJ (rs736799474) have a universal effect on H/L, with CC homozygotes exhibiting improved heterophil function because of downregulated PTPRJ expression. In short, we systematically elucidated the genetic basis of the change in heterophil function resulting from H/L selection by identifying the regulatory gene (PTPRJ) and causative SNP.

Multi-locus genome-wide association studies reveal genomic regions and putative candidate genes associated with leaf spot diseases in African groundnut (Arachis hypogaea L.) germplasm

Early leaf spot (ELS) and late leaf spot (LLS) diseases are the two most destructive groundnut diseases in Ghana resulting in ≤ 70% yield losses which is controlled largely by chemical method. To develop leaf spot resistant varieties, the present study was undertaken to identify single nucleotide polymorphism (SNP) markers and putative candidate genes underlying both ELS and LLS. In this study, six multi-locus models of genome-wide association study were conducted with the best linear unbiased predictor obtained from 294 African groundnut germplasm screened for ELS and LLS as well as image-based indices of leaf spot diseases severity in 2020 and 2021 and 8,772 high-quality SNPs from a 48 K SNP array Axiom platform. Ninety-seven SNPs associated with ELS, LLS and five image-based indices across the chromosomes in the 2 two sub-genomes. From these, twenty-nine unique SNPs were detected by at least two models for one or more traits across 16 chromosomes with explained phenotypic variation ranging from 0.01 - 62.76%, with exception of chromosome (Chr) 08 (Chr08), Chr10, Chr11, and Chr19. Seventeen potential candidate genes were predicted at ± 300 kbp of the stable/prominent SNP positions (12 and 5, down- and upstream, respectively). The results from this study provide a basis for understanding the genetic architecture of ELS and LLS diseases in African groundnut germplasm, and the associated SNPs and predicted candidate genes would be valuable for breeding leaf spot diseases resistant varieties upon further validation.

Genome-Wide Association Study Reveals SNPs and Candidate Genes Related to Growth and Body Shape in Bighead Carp (Hypophthalmichthys nobilis)

Growth is an economically important trait in bighead carp and other aquaculture species that affects production efficiency. Interestingly, the head of the bighead carp has a high market value in China; therefore, it is important to study the genetic bases of both growth and body shape traits. A genome-wide association study was performed based on 2b-RAD sequencing of 776 individuals to identify SNPs associated with growth and body shape traits, including body weight, body length, body height, and deheaded body length. In total, 26 significant and 19 suggestive SNPs were identified, and more than half of these significant SNPs were clustered in LG16. Two LGs (LG16 and LG21) contained QTLs associated with body weight. Fourteen SNPs of LG16 and two LG21 SNPs were found to be associated with body length. For body height, 12 significantly associated SNPs were identified in LG16. Additionally, 12 SNPs of LG16 and 3 SNPs of LG21 were found to be associated with deheaded body length. Forty-three genes were significantly or suggestively associated with body shape/growth traits based on GWAS results, 18 of which were candidate genes for all BW, BL, BH, and DBL traits. One of these genes, fndc5b, was selected for further analyses. Association analysis revealed that one SNP (g.245 C > T) in the introns of fndc5b was significantly associated with growth-related traits in growth-extreme samples. The mRNA levels of fndc5b in the brains of the lightweight group were significantly higher than those of the heavy-weight group. This study helps to reveal the genetic structure of growth and body development in fish and provides candidate genes for future molecular marker-assisted selection for fast growth and better body conformation in bighead carp.

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

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

Pooled resequencing of larvae and adults reveals genomic variations associated with Ostreid herpesvirus 1 resistance in the Pacific oyster Crassostrea gigas

The Ostreid herpesvirus 1 (OsHV-1) is a lethal pathogen of the Pacific oyster (Crassostrea gigas), an important aquaculture species. To understand the genetic architecture of the defense against the pathogen, we studied genomic variations associated with herpesvirus-caused mortalities by pooled whole-genome resequencing of before and after-mortality larval samples as well as dead and surviving adults from a viral challenge. Analysis of the resequencing data identified 5,271 SNPs and 1,883 genomic regions covering 3,111 genes in larvae, and 18,692 SNPs and 28,314 regions covering 4,863 genes in adults that were significantly associated with herpesvirus-caused mortalities. Only 1,653 of the implicated genes were shared by larvae and adults, suggesting that the antiviral response or resistance in larvae and adults involves different sets of genes or differentiated members of expanded gene families. Combined analyses with previous transcriptomic data from challenge experiments revealed that transcription of many mortality-associated genes was also significantly upregulated by herpesvirus infection confirming their importance in antiviral response. Key immune response genes especially those encoding antiviral receptors such as TLRs and RLRs displayed strong association between variation in regulatory region and herpesvirus-caused mortality, suggesting they may confer resistance through transcriptional modulation. These results point to previously undescribed genetic mechanisms for disease resistance at different developmental stages and provide candidate polymorphisms and genes that are valuable for understanding antiviral immune responses and breeding for herpesvirus resistance.

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

Background

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

Results

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

Conclusions

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

Identification of an Epigenetically Marked Locus within the Sex Determination Region of Channel Catfish

Channel catfish has an XY sex determination system. However, the X and Y chromosomes harbor an identical gene content of 950 genes each. In this study, we conducted comparative analyses of methylome and transcriptome of genetic males and genetic females before gonadal differentiation to provide insights into the mechanisms of sex determination. Differentially methylated CpG sites (DMCs) were predominantly identified on the sex chromosome, most notably within the sex determination region (SDR), although the overall methylation profiles across the entire genome were similar between genetic males and females. The drastic differences in methylation were located within the SDR at nucleotide position 14.0–20.3 Mb of the sex chromosome, making this region an epigenetically marked locus within the sex determination region. Most of the differentially methylated CpG sites were hypermethylated in females and hypomethylated in males, suggesting potential involvement of methylation modification in sex determination in channel catfish. Along with the differential methylation in the SDR, a number of differentially expressed genes within the SDR were also identified between genetic males and females, making them potential candidate genes for sex determination and differentiation in channel catfish.

Increased accuracy of genomic predictions for growth under chronic thermal stress in rainbow trout by prioritizing variants from GWAS using imputed sequence data

Through imputation of genotypes, genome-wide association study (GWAS) and genomic prediction (GP) using whole-genome sequencing (WGS) data are cost-efficient and feasible in aquaculture breeding schemes. The objective was to dissect the genetic architecture of growth traits under chronic heat stress in rainbow trout (Oncorhynchus mykiss) and to assess the accuracy of GP based on imputed WGS and different preselected single nucleotide polymorphism (SNP) arrays. A total of 192 and 764 fish challenged to a heat stress experiment for 62 days were genotyped using a customized 1 K and 26 K SNP panels, respectively, and then, genotype imputation was performed from a low-density chip to WGS using 102 parents (36 males and 66 females) as the reference population. Imputed WGS data were used to perform GWAS and test GP accuracy under different preselected SNP scenarios. Heritability was estimated for body weight (BW), body length (BL) and average daily gain (ADG). Estimates using imputed WGS data ranged from 0.33 ± 0.05 to 0.55 ± 0.05 for growth traits under chronic heat stress. GWAS revealed that the top five cumulatively SNPs explained a maximum of 0.94%, 0.86% and 0.51% of genetic variance for BW, BL and ADG, respectively. Some important functional candidate genes associated with growth-related traits were found among the most important SNPs, including signal transducer and activator of transcription 5B and 3 (STAT5B and STAT3, respectively) and cytokine-inducible SH2-containing protein (CISH). WGS data resulted in a slight increase in prediction accuracy compared with pedigree-based method, whereas preselected SNPs based on the top GWAS hits improved prediction accuracies, with values ranging from 1.2 to 13.3%. Our results support the evidence of the polygenic nature of growth traits when measured under heat stress. The accuracies of GP can be improved using preselected variants from GWAS, and the use of WGS marginally increases prediction accuracy.

Omics Strategies in Current Advancements of Infectious Fish Disease Management

Aquaculture is an important industry globally as it remains one of the significant alternatives of animal protein source supplies for humankind. Yet, the progression of this industry is being dampened by the increasing rate of fish mortality, mainly the outbreak of infectious diseases. Consequently, the regress in aquaculture ultimately results in the economy of multiple countries being affected due to the decline of product yields and marketability. By 2025, aquaculture is expected to contribute approximately 57% of fish consumption worldwide. Without a strategic approach to curb infectious diseases, the increasing demands of the aquaculture industry may not be sustainable and hence contributing to the over-fishing of wild fish. Recently, a new holistic approach that utilizes multi-omics platforms including transcriptomics, proteomics, and metabolomics is unraveling the intricate molecular mechanisms of host-pathogen interaction. This approach aims to provide a better understanding of how to improve the resistance of host species. However, no comprehensive review has been published on multi-omics strategies in deciphering fish disease etiology and molecular regulation. Most publications have only covered particular omics and no constructive reviews on various omics findings across fish species, particularly on their immune systems, have been described elsewhere. Our previous publication reviewed the integration of omics application for understanding the mechanism of fish immune response due to microbial infection. Hence, this review provides a thorough compilation of current advancements in omics strategies for fish disease management in the aquaculture industry. The discovery of biomarkers in various fish diseases and their potential advancement to complement the recent progress in combatting fish disease is also discussed in this review.

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.

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.

Toward Genome-Based Selection in Asian Seabass: What Can We Learn From Other Food Fishes and Farm Animals?

Due to the steadily increasing need for seafood and the plateauing output of fisheries, more fish need to be produced by aquaculture production. In parallel with the improvement of farming methods, elite food fish lines with superior traits for production must be generated by selection programs that utilize cutting-edge tools of genomics. The purpose of this review is to provide a historical overview and status report of a selection program performed on a catadromous predator, the Asian seabass (Lates calcarifer, Bloch 1790) that can change its sex during its lifetime. We describe the practices of wet lab, farm and lab in detail by focusing onto the foundations and achievements of the program. In addition to the approaches used for selection, our review also provides an inventory of genetic/genomic platforms and technologies developed to (i) provide current and future support for the selection process; and (ii) improve our understanding of the biology of the species. Approaches used for the improvement of terrestrial farm animals are used as examples and references, as those processes are far ahead of the ones used in aquaculture and thus they might help those working on fish to select the best possible options and avoid potential pitfalls.

A high-density genetic map construction and sex-related loci identification in Chinese Giant salamander

BACKGROUND

The Chinese giant salamander Andrias davidianus is an important amphibian species in China because of its increasing economic value, protection status and special evolutionary position from aquatic to terrestrial animal. Its large genome presents challenges to genetic research. Genetic linkage mapping is an important tool for genome assembly and determination of phenotype-related loci.

RESULTS

In this study, we constructed a high-density genetic linkage map using ddRAD sequencing technology to obtain SNP genotyping data of members from an full-sib family which sex had been determined. A total of 10,896 markers were grouped and oriented into 30 linkage groups, representing 30 chromosomes of A. davidianus. The genetic length of LGs ranged from 17.61 cM (LG30) to 280.81 cM (LG1), with a mean inter-locus distance ranging from 0.11(LG3) to 0.48 cM (LG26). The total genetic map length was 2643.10 cM with an average inter-locus distance of 0.24 cM. Three sex-related loci and four sex-related markers were found on LG6 and LG23, respectively.

CONCLUSION

We constructed the first High-density genetic linkage map and identified three sex-related loci in the Chinese giant salamander. Current results are expected to be a useful tool for future genomic studies aiming at the marker-assisted breeding of the species.

Multivariate linear mixed model enhanced the power of identifying genome-wide association to poplar tree heights in a randomized complete block design

With the advances in high-throughput sequencing technologies, it is not difficult to extract tens of thousands of single-nucleotide polymorphisms (SNPs) across many individuals in a fast and cheap way, making it possible to perform genome-wide association studies (GWAS) of quantitative traits in outbred forest trees. It is very valuable to apply traditional breeding experiments in GWAS for identifying genome variants associated with ecologically and economically important traits in Populus. Here, we reported a GWAS of tree height measured at multiple time points from a randomized complete block design (RCBD), which was established with clones from an F₁ hybrid population of Populus deltoides and Populus simonii. A total of 22,670 SNPs across 172 clones in the RCBD were obtained with restriction site-associated DNA sequencing (RADseq) technology. The multivariate mixed linear model was applied by incorporating the pedigree relationship matrix of individuals to test the association of each SNP to the tree heights over 8 time points. Consequently, 41 SNPs were identified significantly associated with the tree height under the P-value threshold determined by Bonferroni correction at the significant level of 0.01. These SNPs were distributed on all but two chromosomes (Chr02 and Chr18) and explained the phenotypic variance ranged from 0.26% to 2.64%, amounting to 63.68% in total. Comparison with previous mapping studies for poplar height as well as the candidate genes of these detected SNPs were also investigated. We therefore showed that the application of multivariate linear mixed model to the longitudinal phenotypic data from the traditional breeding experimental design facilitated to identify far more genome-wide variants for tree height in poplar. The significant SNPs identified in this study would enhance understanding of molecular mechanism for growth traits and would accelerate marker-assisted breeding programs in Populus.

Differential susceptibility of white bass (Morone chrysops), striped bass (Morone saxatilis) and hybrid striped bass (M. chrysops × M. saxatilis) to Flavobacterium columnare and effects of mucus on bacterial growth and biofilm development

Columnaris disease generates substantial losses of many freshwater fish species; one is the hybrid striped bass. The ubiquitous aquatic bacterium Flavobacterium columnare can be highly effective in biofilm formation on fish skin and gills. Previous research showed a difference between columnaris disease susceptibility of hybrid striped bass (Morone saxatilis × M. chrysops) and white bass (M. chrysops). To understand these differential susceptibilities and possible mucosal relationship, we assessed total bacterial growth and biofilm formation with mucus derived from each moronid parental species: white bass and striped bass (M. saxatilis). Differential susceptibility was confirmed of the other parent species, the striped bass (M. saxatilis). In addition to intraspecies investigations, individual hybrid striped bass mucosal affects were also studied for deferential responses to bacterial growth and biofilm formation. Species- and concentration-dependent differences were detected in the total growth of the bacteria to host mucus. Our data suggest that bass mucus can significantly affect biofilm formation with the F. columnare isolate tested. There appears to be a correlation between the bacteria's response of growth and biofilms and bass species susceptibility. This study provides insight into our understanding of the host-pathogen interaction between F. columnare and moronids.

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 QTLs for resistance to leaf spots in cultivated peanut (Arachis hypogaea L.) through GWAS analysis

Two QTLs on ChrB09 significantly associated with both early and late leaf spots were identified by genome-wide association study in the US peanut mini-core collection. Early leaf spot (ELS) and late leaf spot (LLS) are two serious peanut diseases in the USA, causing tens of millions of dollars of annual economic losses. However, the genetic factors underlying resistance to those diseases in peanuts have not been well-studied. We conducted a genome-wide association study for the two peanut diseases using Affymetrix version 2.0 SNP array with 120 genotypes mainly coming from the US peanut mini-core collection. A total of 46 quantitative trait loci (QTLs) were identified with phenotypic variation explained (PVE) from 10.19 to 24.11%, in which eighteen QTLs are for resistance to ELS and 28 QTLs for LLS. Among the 46 QTLs, there were four and two major QTLs with PVE higher than 16.99% for resistance ELS and LLS, respectively. Of the six major QTLs, five were located on the B sub-genome and only one was on the A sub-genome, which suggested that the B sub-genome has more potential resistance genomic regions than the A sub-genome. In addition, two genomic regions on chromosome B09 were found to provide significant resistance to both ELS and LLS. A total of 74 non-redundant genes were identified as resistance genes, among which, twelve candidate genes were in significant genomic regions including two candidate genes for both ELS and LLS, and other ten candidate genes for ELS. The QTLs and candidate genes obtained from this study will be useful to breed peanuts for resistances to the diseases.

Genomic Prediction of Columnaris Disease Resistance in Catfish

Catfish is an important aquaculture species in the USA. Columnaris disease is distributed worldwide, affecting a wide variety of fish species including catfish . It leads to huge economic losses each year to the US catfish industry. Channel catfish in general is highly resistant to the disease, while blue catfish is highly susceptible. Genomic selection is an effective and accurate way to predict the breeding values and thus was expected to improve the prediction veracity of columnaris disease resistance in catfish effectively. In this study, two different methods, elastic net genomic best linear unbiased prediction (ENGBLUP) and genomic best linear unbiased prediction (GBLUP), were used to predict the columnaris disease resistance evaluated by binary survival status. Cross-validation showed that the prediction accuracy of ENGBLUP and GBLUP was 0.7347 and 0.4868, respectively, showing that ENGBLUP had a high prediction accuracy. It was shown that fitting QTL and polygenic effect with different distribution will improve genomic prediction accuracy for binary traits. In this study, an accurate and effective genomic selection method was proposed to predict the columnaris resistance in catfish, and its application should be beneficial to catfish breeding.

Genome-Wide Association Mapping and Gene Expression Analyses Reveal Genetic Mechanisms of Disease Resistance Variations in Cynoglossus semilaevis

The sustainable development of aquaculture has been impeded by infectious diseases worldwide. However, the genomic architecture and the genetic basis underlying the disease resistance remain poorly understood, which severely hampers both the understanding of the evolution of fish disease resistance traits and the prevention of these diseases in the aquaculture community. Cynoglossus semilaevis is a representative and commercially-important flatfish species. Here we combined genome-wide association study and Fst and nucleotide diversity filtration to identify loci important for the disease resistance. Based on 1,016,774 single-nucleotide polymorphisms (SNPs) identified from 650 Gb genome resequencing data of 505 individuals, we detected 33 SNPs significantly associated with disease resistance and 79 candidate regions after filtration steps. Both the allele frequencies and genotype frequencies of the associated loci were significantly different between the resistant and susceptible fish, suggesting a role in the genetic basis of disease resistance. The SNP with strongest association with disease resistance was located in Chr 17, at 145 bp upstream of fblx19 gene, and overlapped with the major quantitative trait locus previously identified. Several genes, such as plekha7, nucb2, and fgfr2, were also identified to potentially play roles in the disease resistance. Furthermore, the expression of some associating genes were likely under epigenetic regulations between the bacterial resistant and susceptible families. These results provide insights into the mechanism that enable variation of disease resistance to bacterial pathogen infection. The identified polymorphisms and genes are valuable targets and molecular resources for disease resistance and other traits, and for advanced breeding practice for superior germplasm in fish aquaculture.

Genome-Wide Association Study of Growth and Body-Shape-Related Traits in Large Yellow Croaker (Larimichthys crocea) Using ddRAD Sequencing

Large yellow croaker (Larimichthys crocea) is an economically important marine fish species of China. Due to overfishing and marine pollution, the wild stocks of this croaker have collapsed in the past decades. Meanwhile, the cultured croaker is facing the difficulties of reduced genetic diversity and low growth rate. To explore the molecular markers related to the growth traits of croaker and providing the related SNPs for the marker-assisted selection, we used double-digest restriction-site associated DNA (ddRAD) sequencing to dissect the genetic bases of growth traits in a cultured population and identify the SNPs that associated with important growth traits by GWAS. A total of 220 individuals were genotyped by ddRAD sequencing. After quality control, 27,227 SNPs were identified in 220 samples and used for GWAS analysis. We identified 13 genome-wide significant associated SNPs of growth traits on 8 chromosomes, and the beta P of these SNPs ranged from 0.01 to 0.86. Through the definition of candidate regions and gene annotation, candidate genes related to growth were identified, including important regulators such as fgf18, fgf1, nr3c1, cyp8b1, fabp2, cyp2r1, ppara, and ccm2l. We also identified SNPs and candidate genes that significantly associated with body shape, including bmp7, col1a1, col11a2, and col18a1, which are also economically important traits for large yellow croaker aquaculture. The results provided insights into the genetic basis of growth and body shape in large yellow croaker population and would provide reliable genetic markers for molecular marker-assisted selection in the future. Meanwhile, the result established a basis for our subsequent fine mapping and related gene study.

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.

Development of a goosegrass (Eleusine indica) draft genome and application to weed science research

BACKGROUND

Genomes are vital to the study of genomics, population genetics, and evolution of species. To date, only one genome (Echinochloa crus-galli) for C4 annual weedy grass species has been sequenced. Research was conducted to develop a draft genome of goosegrass (Eleusine indica; 2n = 2x = 18), one of the most common and troublesome weeds in the world.

RESULTS

A draft assembly of an approximately 492 Mb whole-genome sequence of goosegrass was obtained by de novo assembly of paired-end and mate-paired reads generated by Illumina sequencing of total genomic DNA. The genome was assembled into 24,072 scaffolds with N50 = 233,459 bp. More than 99% of transcriptome sequences were mapped to the goosegrass draft genome, and 95% of the commonly conserved plant genes were present. The assembled genome contains 25,467 unique protein-coding genes. Genes associated with herbicide resistance were obtained and variant calling allowed the detection of 754,409 single nucleotide polymorphisms. In addition, we also report 115,417 simple sequence repeats which can be deployed in population genetics and phylogenetic analysis.

CONCLUSION

This is the first report of genome sequence of goosegrass. Our assembly was able to identify all major herbicide-resistance related genes and develop a useful tool for other genomic and evolutionary analysis. © 2019 Society of Chemical Industry.

Post-transcriptional regulation through alternative splicing after infection with Flavobacterium columnare in channel catfish (Ictalurus punctatus)

Columnaris disease has long been recognized as a serious problem worldwide which affects both wild and cultured freshwater fish including the commercially important channel catfish (Ictalurus punctatus). The fundamental molecular mechanisms of the host immune response to the causative agent Flavobacterium columnare remain unclear, though gene expression analysis after the bacterial infection has been conducted. Alternative splicing, a post-transcriptional regulation process to modulate gene expression and increase the proteomic diversity, has not yet been studied in channel catfish following infection with F. columnare. In this study, genomic information and RNA-Seq datasets of channel catfish were used to characterize the changes of alternative splicing after the infection. Alternative splicing was shown to be induced by F. columnare infection, with 8.0% increase in alternative splicing event at early infection stage. Intriguingly, genes involved in RNA binding and RNA splicing themselves were significantly enriched in differentially alternatively spliced (DAS) gene sets after infection. This finding was consistent with our previous study in channel catfish following infection with Edwardsiella ictaluri. It was suggested to be a universal mechanism that genes involved in RNA binding and splicing were regulated to undergo differential alternative splicing after stresses in channel catfish. Moreover, many immune genes were observed to be differentially alternatively spliced after infection. Further studies need to be performed to get a deeper view of molecular regulation on alternative splicing after stresses, setting a foundation for developing catfish broodstocks with enhanced disease resistance.

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.

A Novel Candidate Gene Associated With Body Weight in the Pacific White Shrimp Litopenaeus vannamei

Improvements of growth traits are always the focus in selective breeding programs for the Pacific white shrimp Litopenaeus vannamei (L. vannamei). Identification of growth-related genes or markers can contribute to the application of modern breeding technologies, and thus accelerate the genetic improvement of growth traits. The aim of this study was to identify the genes and molecular markers associated with the growth traits of L. vannamei. A population of 200 individuals was genotyped using 2b-RAD techniques for genome-wide linkage disequilibrium (LD) analysis and genome-wide association study (GWAS). The results showed that the LD decayed fast in the studied population, which suggest that it is feasible to fine map the growth-related genes with GWAS in L. vannamei. One gene designated as LvSRC, encoding the class C scavenger receptor (SRC), was identified as a growth-related candidate gene by GWAS. Further targeted sequencing of the candidate gene in another population of 322 shrimps revealed that several non-synonymous mutations within LvSRC were significantly associated with the body weight (P < 0.01), and the most significant marker (SRC_24) located in the candidate gene could explain 13% of phenotypic variance. The current results provide not only molecular markers for genetic improvement in L. vannamei, but also new insights for understanding the growth regulation mechanism in penaeid shrimp.

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.

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

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

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

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

pik3r3b, a novel immune-related gene in Nile tilapia (Oreochromis niloticus): Identification, expression and analysis of antibacterial activity

A full-length cDNA encoding phosphatidylinositol 3-kinase regulatory subunit gamma b gene in Nile tilapia (Oreochromis niloticus), termed as On-pik3r3b, was identified and characterized in this study. The sequence analysis demonstrated that the full-length cDNA of On-pik3r3b was 2018 bp, containing a 5' untranslated region (UTR) of 171 bp, an open reading frame (ORF) of 1422 bp and a 3' UTR of 425 bp. Its protein sequence displayed a high degree of identity with other fish. Using qPCR, the expression patterns of On-pik3r3b were investigated. In healthy Nile tilapia, the transcripts of On-pik3r3b were detected in all examined tissues, except the skin. Upon the stimulation with Streptococcus agalactiae, the On-pik3r3b expression level in liver, spleen, kidney and gill were significantly increased at 12 h after infection. The recombinant On-pik3r3b showed in vitro antibacterial activity, against S. agalactiae and E. coli. Our observation strongly indicates that On-pik3r3b is involved in the innate immune response in Nile tilapia.

Netting the Stress Responses in Fish

In the last decade, the concept of animal stress has been stressed thin to accommodate the effects of short-term changes in cell and tissue physiology, major behavioral syndromes in individuals and ecological disturbances in populations. Seyle's definition of stress as "the nonspecific (common) result of any demand upon the body" now encompasses homeostasis in a broader sense, including all the hierarchical levels in a networked biological system. The heterogeneity of stress responses thus varies within individuals, and stressors become multimodal in terms of typology, source and effects, as well as the responses that each individual elicits to cope with the disturbance. In fish, the time course of changes after stress strongly depends on several factors, including the stressful experiences in early life, the vertical transmission of stressful-prone phenotypes, the degree of individual phenotypic plasticity, the robustness and variety of the epigenetic network related to environmentally induced changes, and the intrinsic behavioral responses (individuality/personality) of each individual. The hierarchical heterogeneity of stress responses demands a code that may decrypt and simplify the analysis of both proximate and evolutionary causes of a particular stress phenotype. We propose an analytical framework, the stressotope, defined as an adaptive scenario dominated by common environmental selective pressures that elicit common multilevel acute stress-induced responses and produce a measurable allostatic load in the organism. The stressotope may constitute a blueprint of embedded interactions between stress-related variations in cell states, molecular mediators and systemic networks, a map of circuits that reflect the inherited and acquired stress responses in an ever-changing, microorganismal-loaded medium. Several features of the proposed model are discussed as a starting point to pin down the maximum common stress responses across immune-neuroendocrine relevant physiological levels and scenarios, including the characterization of behavioral responses, in fish.

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

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

Genetic Mapping of Head Size Related Traits in Common Carp (Cyprinus carpio)

Head size is important economic trait for many aquaculture fish which is directly linked to their carcass yield. The genetic basis of head size trait remains unclear in many widely cultured fish species. Common carp (Cyprinus carpio) is one of the most widely studied fish due to its importance on both economic and environmental aspects. In this study, we performed genome-wide association study using 433 Yellow River carp individuals from multiple families to identify loci and genes potentially associated with head size related traits including head length (HL), head length/body length ratio (HBR), eye diameter (ED), and eye cross (EC). QTL mapping was utilized to filter the effects of population stratification and improve power for the candidates identification in the largest surveyed family with a published genetic linkage map. Twelve SNPs showed significant for head size traits in GWAS and 18 QTLs were identified in QTL mapping. Our study combining both GWAS and QTL mapping could compensate the deficiency from each other and advance our understanding of head size traits in common carp. To acquire a better understanding of the correlation between head size and body growth, we also performed comparisons between QTLs of head size traits and growth-related traits. Candidate genes underlying head size traits were identified surrounding the significant SNPs, including parvalbumin, srpk2, fsrp5, igf1, igf3, grb10, igf1r, notch2, sfrp2. Many of these genes have been identified with potential functions on bone formation and growth. Igf1 was a putative gene associated with both head size and body growth in Yellow River carp. The teleost-specific igf3 was a candidate head size related gene, related to both HL and HBR. Our study also indicated the importance of Igf signaling pathway for both growth and head size determination in common carp, which could be potentially used in future selective breeding in common carp as well as other species.

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.

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

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

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

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

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

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

Teleosts Genomics: Progress and Prospects in Disease Prevention and Control

Genome wide studies based on conventional molecular tools and upcoming omics technologies are beginning to gain functional applications in the control and prevention of diseases in teleosts fish. Herein, we provide insights into current progress and prospects in the use genomics studies for the control and prevention of fish diseases. Metagenomics has emerged to be an important tool used to identify emerging infectious diseases for the timely design of rational disease control strategies, determining microbial compositions in different aquatic environments used for fish farming and the use of host microbiota to monitor the health status of fish. Expounding the use of antimicrobial peptides (AMPs) as therapeutic agents against different pathogens as well as elucidating their role in tissue regeneration is another vital aspect of genomics studies that had taken precedent in recent years. In vaccine development, prospects made include the identification of highly immunogenic proteins for use in recombinant vaccine designs as well as identifying gene signatures that correlate with protective immunity for use as benchmarks in optimizing vaccine efficacy. Progress in quantitative trait loci (QTL) mapping is beginning to yield considerable success in identifying resistant traits against some of the highly infectious diseases that have previously ravaged the aquaculture industry. Altogether, the synopsis put forth shows that genomics studies are beginning to yield positive contribution in the prevention and control of fish diseases in aquaculture.

Catfish Biology and Farming

This article summarizes the biology and culture of ictalurid catfish, an important commercial, aquaculture, and sport fish family in the United States. The history of the propagation as well as spawning of common catfish species in this family is reviewed, with special emphasis on channel catfish and its hybridization with blue catfish. The importance of the channel catfish female×blue catfish male hybrid, including current and future methods of hybrid catfish production, and the potential role it plays in the recovery of the US catfish industry are discussed. Recent advances in catfish culture elements, including environment, management, nutrition, feeding, disease control, culture systems, genetic improvement programs, transgenics, and the application of genome-based approaches in catfish production and welfare, are reviewed. The current status, needs, and future projections are discussed, as well as genetically modified organism developments that are changing the future.

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.

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

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

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

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