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

Unravelling the genomic secrets of bacterial fish pathogens: a roadmap to aquaculture sustainability

In the field of aquaculture, bacterial pathogens pose significant challenges to fish health and production. Advancements in genomic technologies have revolutionized our understanding of bacterial fish pathogens and their interactions with their host species. This review explores the application of genomic approaches in the identification, classification, and characterization of bacterial fish pathogens. Through an extensive analysis of the literature, we have compiled valuable data on 79 bacterial fish pathogens spanning 13 different phyla, encompassing their whole genome sequences. By leveraging high-throughput sequencing techniques, researchers have gained valuable insights into the genomic makeup of these pathogens, enabling a deeper understanding of their virulence factors and mechanisms of host interaction. Furthermore, genomic approaches have facilitated the discovery of potential vaccine and drug targets, opening up new avenues for the development of effective interventions against fish pathogens. Additionally, the utilization of genomics in fish disease resistance and control in aquaculture has shown promising results, enabling the identification of genetic markers associated with disease resistance traits. This review highlights the significant contributions of genomics to the field of fish pathogen research and underscores its potential for improving disease management strategies and enhancing the sustainability of aquaculture practices.

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.

Antimicrobial resistance profile of Aeromonas spp. isolated from asymptomatic Colossoma macropomum cultured in the Amazonas State, Brazil

Bacterial diseases are important factors that limit productivity in aquaculture. To reduce negative economic impacts, fish farmers use antimicrobials, often indiscriminately, and this action has led to bacterial resistance to drugs. The objectives of this study were to isolate and identify the main putative pathogenic bacterial species in tambaqui (Colossoma macropomum), establish the profile of resistance to antimicrobials by the methods of disc diffusion, and determine the minimum inhibitory concentration (MIC) values. Two hundred and ninety asymptomatic fish were collected between March and November 2015 from ten fish farms in the Amazonas state (Brazil). Of the total strains recovered from tambaqui, seven were identified as Aeromonas spp. by sequencing the 16S rRNA gene. These seven isolates showed resistance to ampicillin, 28% to erythromycin, and 28% to sulfonamide. Additionally, the seven isolates showed a MIC higher than the range evaluated for amoxicillin, penicillin, novobiocin, tylosin tartrate, and clindamycin, and 85% showed resistance to erythromycin. The results of this study indicate the need to increase the awareness of fish farmers and, most importantly, the government, about the lack of drug regulations for use in aquaculture, and good management practices, so the indiscriminate prophylactic and systemic use of antimicrobials be inhibited.

The early function of cortisol in liver during Aeromonas hydrophila infection: Dynamics of the transcriptome and accessible chromatin landscapes

In China, channel catfish (Ictalurus punctatus) is an important aquaculture species; however, haemorrhagic disease (Aeromonas hydrophila induced disease) in these fish has caused tremendous economic loss due to high morbidity and mass mortality in the breeding industry. The role of cortisol in bacterial diseases, particularly in the acute phase, remains unclear. In this study, liver transcriptome (RNA-seq) and chromatin accessibility (ATAC-seq) analyses were employed to investigate the early functional role of cortisol in Aeromonas hydrophila-stimulated responses. Our experiments confirmed that A. hydrophila infection can initially significantly increase serum cortisol levels at 1 h after infection. At this time point, the increased serum cortisol levels can significantly regulate A. hydrophila-regulated genes by affecting both transcriptome and chromatin accessibility. Cross-analysis of RNA-seq and ATAC-seq revealed that a certain gene group (92 target_DEGs) was regulated at an early time point by cortisol. KEGG enrichment analysis revealed that the top three pathways according to target_DEGs were cancer, glutathione metabolism, and the Notch signalling pathway. The protein-protein interaction analysis of target_DEGs revealed that they may be primarily involved in cell proliferation, CD8⁺ T cell function, glutathione synthesis, and activation of the NF-κB signalling pathway. This suggests that after the emergence of immune stress, the early regulation of cortisol is positive against the immune response. It is possible that in this situation, the animal is attempting to avoid dangerous situations and risks and then cope with the imbalance produced by the stressor to ultimately restore homeostasis. Our results will contribute to future research on fish and provide valuable insight regarding the mechanism of immune regulation by cortisol and the study of bacterial haemorrhagic disease in channel catfish.

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.

Partially replacing dietary fish meal by Saccharomyces cerevisiae culture improve growth performance, immunity, disease resistance, composition and function of intestinal microbiota in channel catfish (Ictalurus punctatus)

The aim of the present study was to investigate the partial replacement of fish meal by Saccharomyces cerevisiae culture on growth performance, immunity, composition and function of intestinal microbiota and disease resistance in channel catfish (Ictalurus punctatus). Two equal nitrogen and energy diets were prepared including a basal diet (containing 10% fish meal, Control) and an experimental diet (replacing 20% of the fish meal of the basal diet with yeast culture, RFM). Channel catfish were fed with the diets for 12 weeks. The results showed that weight gain and condition factor were significantly increased, and FCR was significantly decreased in RFM group (P < 0.05). The gene expression of intestinal HIF1α was significantly increased in RFM group (P < 0.05), while the expressions of NF-κB in the intestine and liver were significantly decreased (P < 0.05). The relative abundance of Firmicutes tended to increase, and the Turicibacter had an upward trend (0.05 < P < 0.2). In addition, the survival rate of channel catfish was significantly increased in RFM group after challenged with Aeromonas veronii Hm091 and Aeromonas hydrophila NJ-1 (P < 0.05). Compared with intestinal microbiota of channel catfish of control group, intestinal microbiota of channel catfish of RFM group significantly increased the expression of HIF1α, and decreased the expression of IL-1β and TNF-α (P < 0.05) in germ-free zebrafish. Intestinal microbiota induced by RFM diet also significantly increased disease resistance to Aeromonas veronii Hm091 and Aeromonas hydrophila NJ-1. In conclusion, replacement of fish meal by the yeast culture improved the growth, immunity and disease resistance of channel catfish, and intestinal microbiota of channel catfish induced by the yeast culture played a critical role in these effects.

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

Background

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

Results

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

Conclusion

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

Supplementary Information

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

Aeromonas sobria regulates proinflammatory immune response in mouse macrophages via activating the MAPK, AKT, and NF-κB pathways

Aeromonas sobria, a Gram-negative bacterium that can colonize both humans and animals, is found in a variety of environments, including water, seafood, meat, and vegetables (Cahill, 1990; Galindo et al., 2004; Song et al., 2019). Aeromonas spp. are conditionally pathogenic bacteria in aquaculture, which can rapidly proliferate, causing disease and even death in fish, especially when the environment is degraded (Neamat-Allah et al., 2020, 2021a, 2021b). In developing countries, Aeromonas spp. have been associated with a wide spectrum of infections in humans, including gastroenteritis, wound infections, septicemia, and lung infections (San Joaquin and Pickett, 1988; Wang et al., 2009; Su et al., 2013). Infections caused by Aeromonas spp. are usually more severe in immunocompromised individuals (Miyamoto et al., 2017). The presence of a plasmid encoding a β‍-lactamase in A. sobria that confers resistance to β-lactam antibiotics poses a huge challenge to the treatment of diseases caused by this microorganism (Lim and Hong, 2020). Consequently, an in-depth understanding of the interaction between A. sobria and its hosts is urgently required to enable the development of effective strategies for the treatment of A. sobria infections.

Exploring genetic resistance to infectious salmon anaemia virus in Atlantic salmon by genome-wide association and RNA sequencing

Background

Infectious Salmonid Anaemia Virus (ISAV) causes a notifiable disease that poses a large threat for Atlantic salmon (Salmo salar) aquaculture worldwide. There is no fully effective treatment or vaccine, and therefore selective breeding to increase resistance to ISAV is a promising avenue for disease prevention. Genomic selection and potentially genome editing can be applied to enhance host resistance, and these approaches benefit from improved knowledge of the genetic and functional basis of the target trait. The aim of this study was to characterise the genetic architecture of resistance to ISAV in a commercial Atlantic salmon population and study its underlying functional genomic basis using RNA Sequencing.

Results

A total of 2833 Atlantic salmon parr belonging to 194 families were exposed to ISAV in a cohabitation challenge in which cumulative mortality reached 63% over 55 days. A total of 1353 animals were genotyped using a 55 K SNP array, and the estimate of heritability for the trait of binary survival was 0.13–0.33 (pedigree-genomic). A genome-wide association analysis confirmed that resistance to ISAV was a polygenic trait, albeit a genomic region in chromosome Ssa13 was significantly associated with resistance and explained 3% of the genetic variance. RNA sequencing of the heart of 16 infected (7 and 14 days post infection) and 8 control fish highlighted 4927 and 2437 differentially expressed genes at 7 and 14 days post infection respectively. The complement and coagulation pathway was down-regulated in infected fish, while several metabolic pathways were up-regulated. The interferon pathway showed little evidence of up-regulation at 7 days post infection but was mildly activated at 14 days, suggesting a potential crosstalk between host and virus. Comparison of the transcriptomic response of fish with high and low breeding values for resistance highlighted TRIM25 as being up-regulated in resistant fish.

Conclusions

ISAV resistance shows moderate heritability with a polygenic architecture, but a significant QTL was detected on chromosome 13. A mild up-regulation of the interferon pathway characterises the response to the virus in heart samples from this population of Atlantic salmon, and candidate genes showing differential expression between samples with high and low breeding values for resistance were identified.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07671-6.

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

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

MiRNA-seq analysis of spleen and head kidney tissue from aquacultured largemouth bass (Micropterus salmoides) in response to Aeromonas hydrophila infection

Recently, the same fish diseases, which have been found in pond farming, have been found in the newly tested largemouth bass (Micropterus salmoides) system. Bacterial septicemia caused by Aeromonas hydrophila occurs frequently in largemouth bass culture leading to significant economic losses. To investigate the role miRNA in the largemouth bass disease resistance, twelve (2 tissues (spleen and head kidney) × 2 experimental groups (infected and control) × three biological replicates) small RNA libraries were constructed and sequenced with miRNA-seq. A total of 26 differentially expressed miRNAs, 8 upregulated and 18 downregulated, were identified in the spleen, and 19 differentially expressed miRNAs, 9 upregulated and 10 downregulated, were identified in head kidney (fold change ≥ 2 or ≤ 0.5 and P ≤ 0.05). The differentially expressed miRNAs with the largest fold change were selected for target gene prediction using GO and KEGG analysis. Six miRNAs in the spleen and 5 miRNAs in the head kidney were selected. Analysis showed that, of all the immune and metabolic pathways, the FoxO signaling pathway was enriched in both the spleen and head kidney groups. Common target genes of the pathway included AMP-activated catalytic subunit alpha 1 (prkaa1), phosphatidylinositol 3-kinase (pik3r3b), serine/threonine-protein kinase (plk2), and forkhead box protein G1 (foxg1a). MiRNAs (such as miR-126-5P, miR-126-3P) are involved in immune response and cell cycle functions as they regulate targeted genes in the FoxO pathway. These results will enhance our understanding of the molecular mechanisms underlying immune responses to bacterial septicemia and facilitate molecular-assisted selection of resistant strains of largemouth bass.

Development of a SNP linkage map and genome-wide association study for resistance to Aeromonas hydrophila in pacu (Piaractus mesopotamicus)

BACKGROUND

Pacu (Piaractus mesopotamicus) is one of the most important Neotropical aquaculture species from South America. Disease outbreaks caused by Aeromonas hydrophila infection have been considered significant contributors to the declining levels of pacu production. The current implementation of genomic selection for disease resistance has been adopted as a powerful strategy for improvement in fish species. This study aimed to investigate the genetic architecture of resistance to A. hydrophila in pacu via Genome-Wide Association Study (GWAS), the identification of suggestive Quantitative Trait Loci (QTLs) and putative genes associated with this trait. The genetic data were obtained from 381 juvenile individuals belonging to 14 full-sibling families. An experimental challenge was performed to gain access to the levels of genetic variation for resistance against the bacteria using the following trait definitions: binary test survival (TS) and time of death (TD).

RESULTS

The analyses of genetic parameters estimated moderate heritability (h2) for both resistance traits: 0.20 (± 0.09) for TS and 0.35 (± 0.15) for TD. A linkage map for pacu was developed to enable the GWAS, resulting in 27 linkage groups (LGs) with 17,453 mapped Single Nucleotide Polymorphisms (SNPs). The length of the LGs varied from 79.95 (LG14) to 137.01 (LG1) cM, with a total map length of 2755.60 cM. GWAS identified 22 putative QTLs associated to A. hydrophila resistance. They were distributed into 17 LGs, and were considered suggestive genomic regions explaining > 1% of the additive genetic variance (AGV) for the trait. Several candidate genes related to immune response were located close to the suggestive QTLs, such as tbk1, trim16, Il12rb2 and lyz2.

CONCLUSION

This study describes the development of the first medium density linkage map for pacu, which will be used as a framework to study relevant traits to the production of this species. In addition, the resistance to A. hydrophila was found to be moderately heritable but with a polygenic architecture suggesting that genomic selection, instead of marker assisted selection, might be useful for efficiently improving resistance to one of the most problematic diseases that affects the South American aquaculture.

First High-Density Linkage Map and QTL Fine Mapping for Growth-Related Traits of Spotted Sea bass (Lateolabrax maculatus)

Possessing powerful adaptive capacity and a pleasant taste, spotted sea bass (Lateolabrax maculatus) has a broad natural distribution and is one of the most popular mariculture fish in China. However, the genetic improvement program for this fish is still in its infancy. Growth is the most economically important trait and is controlled by quantitative trait loci (QTL); thus, the identification of QTLs and genetic markers for growth-related traits is an essential step for the establishment of marker-assisted selection (MAS) breeding programs. In this study, we report the first high-density linkage map of spotted sea bass constructed by sequencing 333 F1 generation individuals in a full-sib family using 2b-RAD technology. A total of 6883 SNP markers were anchored onto 24 linkage groups, spanning 2189.96 cM with an average marker interval of 0.33 cM. Twenty-four growth-related QTLs, including 13 QTLs for body weight and 11 QTLs for body length, were successfully detected, with phenotypic variance explained (PVE) ranging from 5.1 to 8.6%. Thirty potential candidate growth-related genes surrounding the associated SNPs were involved in cell adhesion, cell proliferation, cytoskeleton reorganization, calcium channels, and neuromodulation. Notably, the fgfr4 gene was detected in the most significant QTL; this gene plays a pivotal role in myogenesis and bone growth. The results of this study may facilitate marker-assisted selection for breeding populations and establish the foundation for further genomic and genetic studies investigating spotted sea bass.

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.

Genome-Wide DNA Methylation and RNA Analysis Reveal Potential Mechanism of Resistance to Streptococcus agalactiae in GIFT Strain of Nile Tilapia (Oreochromis niloticus )

Streptococcus agalactiae is an important pathogenic bacterium causing great economic loss in Nile tilapia (Oreochromis niloticus) culture. Resistant and susceptible groups sharing the same genome showed significantly different resistance to S. agalactiae in the genetically improved farmed tilapia strain of Nile tilapia. The resistance mechanism is unclear. We determined genome-wide DNA methylation profiles in spleen of resistant and susceptible O. niloticus at 5 h postinfection with S. agalactiae using whole-genome bisulfite sequencing. The methylation status was higher in the spleen samples from resistant fish than in the susceptible group. A total of 10,177 differentially methylated regions were identified in the two groups, including 3725 differentially methylated genes (DMGs) (3129 hyper-DMGs and 596 hypo-DMGs). The RNA sequencing showed 2374 differentially expressed genes (DEGs), including 1483 upregulated and 891 downregulated. Integrated analysis showed 337 overlapping DEGs and DMGs and 82 overlapping DEGs and differentially methylated region promoters. By integrating promoter DNA methylation with gene expression, we revealed four immune-related genes (Arnt2, Nhr38, Pcdh10, and Ccdc158) as key factors in epigenetic mechanisms contributing to pathogen resistance. Our study provided systematic methylome maps to explore the epigenetic mechanism and reveal the methylation loci of pathogen resistance and identified methylation-regulated genes that are potentially involved in defense against pathogens.

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

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

Construction of a High-Density Genetic Linkage Map and QTL Mapping for Growth-Related Traits in Takifugu bimaculatus

Takifugu bimaculatus is a euryhaline species, distributed ranging from the southern Yellow Sea to the South China Sea. Their tolerance to a wide range of salinity and temperature, coupled with a desirable firm texture, makes T. bimaculatus a strong candidate for Takifugu aquaculture in subtropics areas. Due to the increasing demand in markets and emerging of the Takifugu aquaculture industry, close attention has been paid to improvement on the T. bimaculatus production. In aquaculture, the great effort has been put into marker-assisted selective breeding, and efficient improvement was realized. However, few genetic resources on T. bimaculatus are provided so far. Aiming at understanding the genetic basis underlying important economic growth traits, facilitating genetic improvement and enriching the genetic resource in T. bimaculatus, we constructed the first genetic linkage map for T. bimaculatus via double digestion restriction-site association DNA sequencing and conducted quantitative traits locus (QTL) mapping for growth-related traits. The map comprised 1976 single nucleotide polymorphism markers distributed on 22 linkage groups (LG), with a total genetic distance of 2039.74 cM. Based on the linkage map, a chromosome-level assembly was constructed whereby we carried out comparative genomics analysis, verifying the high accuracy on contigs ordering of the linkage map. On the other hand, 18 QTLs associated with growth traits were detected on LG6, LG7, LG8, LG10, LG20, and LG21 with phenotypical variance ranging from 15.1 to 56.4%. Candidate genes participating in cartilage development, fat accumulation, and other growth-related regulation activities were identified from these QTLs, including col11a1, foxa2, and thrap3. The linkage map provided a solid foundation for chromosomes assembly and refinement. QTLs reported here unraveled the genomic architecture of some growth traits, which will advance the investigation of aquaculture breeding efforts in T. bimaculatus.

An Update on the Genus Aeromonas: Taxonomy, Epidemiology, and Pathogenicity

The genus Aeromonas belongs to the Aeromonadaceae family and comprises a group of Gram-negative bacteria widely distributed in aquatic environments, with some species able to cause disease in humans, fish, and other aquatic animals. However, bacteria of this genus are isolated from many other habitats, environments, and food products. The taxonomy of this genus is complex when phenotypic identification methods are used because such methods might not correctly identify all the species. On the other hand, molecular methods have proven very reliable, such as using the sequences of concatenated housekeeping genes like gyrB and rpoD or comparing the genomes with the type strains using a genomic index, such as the average nucleotide identity (ANI) or in silico DNA–DNA hybridization (isDDH). So far, 36 species have been described in the genus Aeromonas of which at least 19 are considered emerging pathogens to humans, causing a broad spectrum of infections. Having said that, when classifying 1852 strains that have been reported in various recent clinical cases, 95.4% were identified as only four species: Aeromonas caviae (37.26%), Aeromonas dhakensis (23.49%), Aeromonas veronii (21.54%), and Aeromonas hydrophila (13.07%). Since aeromonads were first associated with human disease, gastroenteritis, bacteremia, and wound infections have dominated. The literature shows that the pathogenic potential of Aeromonas is considered multifactorial and the presence of several virulence factors allows these bacteria to adhere, invade, and destroy the host cells, overcoming the immune host response. Based on current information about the ecology, epidemiology, and pathogenicity of the genus Aeromonas, we should assume that the infections these bacteria produce will remain a great health problem in the future. The ubiquitous distribution of these bacteria and the increasing elderly population, to whom these bacteria are an opportunistic pathogen, will facilitate this problem. In addition, using data from outbreak studies, it has been recognized that in cases of diarrhea, the infective dose of Aeromonas is relatively low. These poorly known bacteria should therefore be considered similarly as enteropathogens like Salmonella and Campylobacter.

Evaluation of Genomic Selection for Seven Economic Traits in Yellow Drum (Nibea albiflora)

Yellow drum (Nibea albiflora) is an important maricultural fish in China, and genetic improvement is necessary for this species. This research evaluated the application of genomic selection methods to predict the genetic values of seven economic traits for yellow drum. Using genome-wide single-nucleotide polymorphisms (SNPs), we estimated the genetic parameters for seven traits, including body length (BL), swimming bladder index (SBI), swimming bladder weight (SBW), body thickness (BT), body height (BH), body length/body height ratio (LHR), and gonad weight index (GWI). The heritability estimates ranged from 0.309 to 0.843. We evaluated the prediction performance of various statistical methods, and no one method provided the highest predictive ability for all traits. We then evaluated and compared the use of genome-wide association study (GWAS)-informative SNPs and random SNPs for prediction and found that GWAS-informative SNPs obviously increased. It only needed 5 and 100 informative SNPs for LHR and BT to achieve almost the same predictive abilities as using genome-wide SNPs, and for BL, SBI, SBW, BH, and GWI, about 1000 to 3000 informative SNPs were needed to achieve whole-genome level predictive abilities. It can be concluded from the test results that breeders can use fewer SNPs to save the breeding costs of genomic selection for some traits.

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.