Deciphering transcriptome profile of the yellow catfish (Pelteobagrus fulvidraco) in response to Edwardsiella ictaluri

Toll-like Receptor 3 in the Hybrid Yellow Catfish (Pelteobagrus fulvidraco ♀ × P. vachelli ♂): Protein Structure, Evolution and Immune Response to Exogenous Aeromonas hydrophila and Poly (I:C) Stimuli

As a major mediator of cellular response to viral infection in mammals, Toll-like receptor 3 (TLR3) was proved to respond to double-stranded RNA (dsRNA). However, the molecular mechanism by which TLR3 functions in the viral infection response in teleosts remains to be investigated. In this study, the Toll-like receptor 3 gene of the hybrid yellow catfish was identified and characterized by comparative genomics. Furthermore, multiple sequence alignment, genomic synteny and phylogenetic analysis suggested that the homologous TLR3 genes were unique to teleosts. Gene structure analysis showed that five exons and four introns were common components of TLR3s in the 12 examined species, and interestingly the third exon in teleosts was the same length of 194 bp. Genomic synteny analysis indicated that TLR3s were highly conserved in various teleosts, with similar organizations of gene arrangement. De novo predictions showed that TLR3s were horseshoe-shaped in multiple taxa except for avian (with a round-shaped structure). Phylogenetic topology showed that the evolution of TLR3 was consistent with the evolution of the studied species. Selection analysis showed that the evolution rates of TLR3 proteins were usually higher than those of TLR3-TIR domains, indicating that the latter were more conserved. Tissue distribution analysis showed that TLR3s were widely distributed in the 12 tested tissues, with the highest transcriptions in liver and intestine. In addition, the transcription levels of TLR3 were significantly increased in immune-related tissues after infection of exogenous Aeromonas hydrophila and poly (I:C). Molecular docking showed that TLR3 in teleosts forms a complex with poly (I:C). In summary, our present results suggest that TLR3 is a pattern recognition receptor (PRR) gene in the immune response to pathogen infections in hybrid yellow catfish.

Three IRF4 paralogs act as negative regulators of type Ⅰ IFN responses in yellow catfish (Pelteobagrus fulvidraco)

IRF4 is a master member of the interferon regulatory factor (IRF) family playing vital regulatory roles in immune system development and function. Tetrapods have a single-copy IRF4 gene, while teleosts harbor duplicated IRF4 genes. This work describes three IRF4 paralogs from yellow catfish (Pelteobagrus fulvidraco), designated PfIRF4A, PfIRF4B and PfIRF4B-like. These genes all contain a typical IRF structural architecture. Phylogenic and synteny analyses indicate that they should arise from the teleost-specific whole-genome duplication. PfIRF4 genes are abundantly expressed in the immune-related tissues and upregulated by PolyI:C, LPS, and Edwardsiella ictaluri. Ectopic expression of these genes inhibits the activation of fish type Ⅰ IFN promoters and downregulates the transcription levels of IFN-responsive genes, thus allowing the efficient replication of a fish rhabdovirus, spring viremia of carp virus (SVCV). PfIRF4s possess a repressive effect on MyD88-mediated activation of IFN and NF-κB. Some differences are observed between each individual paralog. PfIRF4B is the main form expressed across the tissues and the most up-regulated one after pathogen induction. It exerts a stronger inhibitory effect on IFN antiviral response than the other two paralogs. PfIRF4A and PfIRF4B-like are primarily present in the nucleus, while PfIRF4B displays colocalization and direct associations with MyD88 in the cytoplasm. Overall, the data demonstrate that three PfIRF4 paralogs show shared and individual functional properties in the negative regulation of type Ⅰ IFN response.

Effects of dietary Bacillus amyloliquefaciens on the growth, immune responses, intestinal microbiota composition and disease resistance of yellow catfish, Pelteobagrus fulvidraco

The influence of dietary probiotic Bacillus amyloliquefaciens on the growth performance, digestive enzyme activity, immune parameters and disease resistance of yellow catfish (Pelteobagrus fulvidraco) was evaluated. Commercial diet (C) or diet containing 106 cfu/g B. amyloliquefaciens (T) was fed for 4 weeks, and final weight (FW), specific growth rate (SGR) and feed conversion ratio (FCR) were improved (p<0.05) in the T group. Dietary B. amyloliquefaciens increased protease and amylase activities in the digestive tract after 2 and 4 weeks, respectively. Respiratory burst (RB), plasma lysozyme (LZM) activity, total antioxidant capacity (T-AOC) and superoxide dismutase (SOD) activity were also elevated (p<0.05). Immune-related genes signal transducer and activator of transcription 1 (STATA-1), immunoglobulin M (IgM) and C-type lectin (CTL) were upregulated (p<0.05), but interleukin-1 β (IL-1β) was not (p >0.05). Intestinal microbiota analysis showed that the community structure was significantly different between the two groups; the relative abundance of Cetobacterium was increased but Plesiomonas was decreased in T. Moreover, challenge tests showed that the resistance of fish fed B. amyloliquefaciens against Aeromonas veronii and Edwardsiella ictaluri was significantly enhanced (p<0.05). In conclusion, dietary supplementation of B. amyloliquefaciens can effectively improve the growth performance, digestive enzyme activity, immune responses, intestinal microbiota composition and disease resistance of yellow catfish.

Granulomatous bacterial diseases in fish: An overview of the host's immune response

Bacterial diseases represent the main impediment to the development of fish aquaculture. Granulomatous diseases caused by bacteria lead to fish culture losses by high mortality rates and slow growth. Bacteria belonging to genera Streptococcus spp., Mycobacterium sp., Nocardia sp., Francisella sp., and Staphylococcus sp. have been implicated in the development of granulomatous processes. The granuloma formation and the fish's immune response continue to be the subject of scientific research. In fish, the first defense line is constituted by non-specific humoral factors through growth-inhibiting substances such as transferrin and antiproteases, or lytic effectors as lysozyme and antimicrobial peptides, and linking with non-specific phagocyte responses. If the first line is breached, fish produce antibody constituents for a specific humoral defense inhibiting bacterial adherence, as well as the mobilization of non-phagocytic host cells and counteracting toxins from bacteria. However, bacteria causing granulomatous diseases can be persistent microorganisms, difficult to eliminate that can cause chronic diseases, even using some immune system components to survive. Understanding the infectious process leading to granulomatosis and how the host's immune system responds against granulomatous diseases is crucial to know more about fish immunology and develop strategies to overcome granulomatous diseases.

The Dynamic Immune Response of Yellow Catfish (Pelteobagrus fulvidraco) Infected With Edwardsiella ictaluri Presenting the Inflammation Process

Edwardsiella ictaluri is a highly destructive pathogen in cultured yellow catfish, thus it was very necessary to study the immune response of yellow catfish against bacterial infection. In this study, RNA-Seq technology was used to study the immune response in two distinct tissues of yellow catfish at eight different time points (h) after E. ictaluri infection. The number of differentially expressed genes (DEGs) in the spleen and liver was low at 3 h and 6 h post-infection, respectively. Afterwards, the most number of DEGs in the spleen was detected at 72 h, while the number of DEGs in the liver maintained a high level from 24 h to 120 h. The GO and KEGG enrichment analyses of DEGs at different time points uncovered that cytokines were continuously transcribed at 6 h to 120 h; whereas the liver is the main organ that secretes the components of the complement system, and metabolic regulation was activated from 12 h to 120 h. Moreover, an overview of the inflammation response of yellow catfish was exhibited including pattern-recognition receptors, inflammatory cytokines, chemokines, complements, and inflammation-related signal pathways. The similar expression tendency of nine genes by qRT-PCR validated the accuracy of transcriptome analyses. The different transcriptomic profiles obtained from the spleen and liver will help to better understand the dynamic immune response of fish against bacterial infection, and will provide basic information for establishing effective measures to prevent and control diseases in fish.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Molecular characterization of yellow catfish (Pelteobagrus fulvidraco) IRF7 suggests involvement in innate immune response

Interferon regulatory factor 7 (IRF7) serves as a critical mediator in the regulation of type Ι interferon (IFN) response to invading pathogens. Here, an ortholog of IRF7 was characterized in yellow catfish (Pelteobagrus fulvidraco). The full-length cDNA of PfIRF7 consisted of 1516 bp encoding a polypeptide of 425 amino acids. PfIRF7 protein comprised a typical IRF structural architecture, including a DNA binding domain (DBD), an IRF association domain (IAD) and a serine-rich domain (SRD). PfIRF7 was expressed predominantly in the immune-related tissues and transcriptionally upregulated by PolyI:C, LPS, and Edwardsiella ictaluri. Ectopic expression of PfIRF7 led to activation of fish type I IFN promoters and induction of IFN and Vig1, thereby conferring a strong antiviral effect against spring viremia of carp virus (SVCV). Overall, the present data suggest that PfIRF7 may play an essential role in type I IFN response of yellow catfish.

Evolution of digestive enzyme genes associated with dietary diversity of crabs

Crabs feed on a wide range of items and display diverse feeding strategies. The primary objective of this study was to investigate 10 digestive enzyme genes in representative crabs to provide insights into the genetic basis of feeding habits among crab functional groups. Crabs were classified into three groups based on their feeding habits: herbivores (HV), omnivores (OV), and carnivores (CV). To test whether crabs' feeding adaptations matched adaptive evolution of digestive enzyme genes, we examined the 10 digestive enzyme genes of 12 crab species based on hepatopancreas transcriptome data. Each of the digestive enzyme genes was compared to orthologous sequences using both nucleotide- (i.e., PAML and Datamonkey) and protein-level (i.e., TreeSAAP) approaches. Positive selection genes were detected in HV crabs (AMYA, APN, and MGAM) and CV crabs (APN, CPB, PNLIP, RISC, TRY, and XPD). Additionally, a series of positive selection sites were localized in important functional regions of these digestive enzyme genes. This is the first study to characterize the molecular basis of crabs' digestive enzyme genes based on functional feeding group. Our data suggest that HV crabs have evolved an enhanced digestion capacity for carbohydrates, and CV crabs have acquired digestion capacity for proteins and lipids.

Identification of Multigene Biomarker for Shrimp White Feces Syndrome by Full-Length Transcriptome Sequencing

The pacific white shrimp, Litopenaeus vannamei, with the largest shrimp industry production in the world, is currently threatened by a severe disease, white feces syndrome (WFS), which cause devastating losses globally, while its causal agents remain largely unknown. Herein, compared to the Control shrimp by metagenomic analysis, we firstly investigated that the altered functions of intestinal microbial community in WFS shrimp were the enrichment of bacterial chemotaxis and flagellar assembly pathways, hinting at a potential role of pathogenic bacteria for growth and development, which might be related to WFS occurrence. Single-molecule real-time (SMRT) sequencing was to further identify the gene structure and gene regulation for more clues in WFS aetiology. Totally 50,049 high quality transcripts were obtained, capturing 39,995 previously mapped and 10,054 newly detected transcripts, which were annotated to 30,554 genes. A total of 158 differentially expressed genes (DEGs) were characterized in WFS shrimp. These DEGs were strongly associated with various immune related genes that regulated the expression of multiple antimicrobial peptides (e.g., antilipopolysaccharide factors, penaeidins, and crustin), which were further experimentally validated using quantitative PCR on transcript level. Collectively, multigene biomarkers were identified to be closely associated with WFS, especially those functional alterations in microbial community and the upregulated immune related gene with antibacterial activities. Our finding not only inspired our cogitation on WFS aetiology from both microbial and host immune response perspectives with combined metagenomic and full-length transcriptome sequencing, but also provided valuable information for enhancing shrimp aquaculture.

New insight into the molecular basis of Fe (III) stress responses of Procambarus clarkii by transcriptome analysis

Iron in excess can have toxic effects on living organisms. In China, the freshwater crayfish Procambarus clarkii is a source of aquatic food with high-quality protein and has significant commercial value. P. clarkii shows oxidative stress on exposure to heavy metals, and antioxidant enzymes, such as ubiquitination enzymes and proteasomes, play important roles in oxidative stress. To understand the antioxidant defense system of P. clarkii, we analyzed the hepatopancreas transcriptomes of P. clarkii after stimulation with FeCl₃. In total, 5199 differentially expressed genes (DEGs) were identified (2747 upregulated and 2452 downregulated). GO analysis revealed that these DEGs belonged to 16 cellular component, 16 molecular function, and 19 biological process subcategories. A total of 1069 DEGs were classified into 25 categories by using COG. Some antioxidant defense pathways, such as "Ubiquitin mediated proteolysis" and "Glutathione metabolism," were identified using KEGG. In addition, quantitative real time-PCR (qRT-PCR) substantiated the up-regulation of a random selection of DEGs including antioxidant and immune defense genes. We obtained information for P. clarkii transcriptome databases and new insights into the responses of P. clarkii hepatopancreas to heavy metals.

Comparative transcriptome analyses reveal changes of gene expression in fresh and cryopreserved yellow catfish (Pelteobagrus fulvidraco) sperm and the effects of Cryoprotectant Me2SO

This study, for the first time in fish, compared the transcriptome of fresh and frozen-thawed sperm, and would help to better understand the effect of cryopreservation on fish sperm and then better preserve the aquatic germplasm resources. Here, we employed high-throughput sequencing technology to obtain the transcriptome of yellow catfish from fresh sperm, cryopreserved sperm with and without cryoprotectant. When cryoprotectant (Me₂SO) was excluded, down-regulated genes were significantly enriched into calcium ion binding, cytoskeletal protein binding, microfilament motor activity, calmodulin binding and carnitine O-acyltransferase activity, which affected Ca²⁺ regulation, cellular morphology, motility and metabolism. Moreover, heat shock proteins and genes associated with regulation of cholesterol, HCO₃^- and protein tyrosine phosphorylation (PTP) were down-regulated, and thus would impair ability against stress, membrane rigidity, pH regulation and signal transduction of cryopreserved sperm. After Me₂SO was added, the amounts of DEGs decreased significantly and down-regulation of genes were found mainly in cytoskeleton and heat shock proteins, thereby suggesting that Me₂SO effectively reduced the impact caused by low temperature on gene expression. Whether adding Me₂SO or not, the up-regulated genes were mainly found in ribosomal proteins genes. However, when Me₂SO was added, over-expression of some genes might contribute to maintain normal function of cryopreserved sperm.

Comparative Study of Immune Reaction Against Bacterial Infection From Transcriptome Analysis

Transcriptome analysis is a powerful tool that enables a deep understanding of complicated physiological pathways, including immune responses. RNA sequencing (RNA-Seq)-based transcriptome analysis and various bioinformatics tools have also been used to study non-model animals, including aquaculture species for which reference genomes are not available. Rapid developments in these techniques have not only accelerated investigations into the process of pathogenic infection and defense strategies in fish, but also used to identify immunity-related genes in fish. These findings will contribute to fish immunotherapy for the prevention and treatment of bacterial infections through the design of more specific and effective immune stimulants, adjuvants, and vaccines. Until now, there has been little information regarding the universality and diversity of immune reactions against pathogenic infection in fish. Therefore, one of the aims of this paper is to introduce the RNA-Seq technique for examination of immune responses in pathogen-infected fish. This review also aims to highlight comparative studies of immune responses against bacteria, based on our previous findings in largemouth bass (Micropterus salmoides) against Nocardia seriolae, gray mullet (Mugil cephalus) against Lactococcus garvieae, orange-spotted grouper (Epinephelus coioides) against Vibrio harveyi, and koi carp (Cyprinus carpio) against Aeromonas sobria, using RNA-seq techniques. We demonstrated that only 39 differentially expressed genes (DEGs) were present in all species. However, the number of specific DEGs in each species was relatively higher than that of common DEGs; 493 DEGs in largemouth bass against N. seriolae, 819 DEGs in mullets against L. garvieae, 909 in groupers against V. harveyi, and 1471 in carps against A. sobria. The DEGs in different fish species were also representative of specific immune-related pathways. The results of this study will enhance our understanding of the immune responses of fish, and will aid in the development of effective vaccines, therapies, and disease-resistant strains.

Transcriptome analysis of Artemia sinica in response to Micrococcus lysodeikticus infection

To enhance genomic resources and understand the molecular immune mechanisms underlying the response topathogens, we first performed a comparative gene transcription analysis from Micrococcus lysodeikticus-immunized Artemia sinica and from a control group through RNA-Seq technology, meanwhile the differentially expressed genes (DEGs) were investigated. In total, 80, 113, 984 clean reads were obtained and then assembled into 71,536 unigenes with an average length of 1115 bp and an N50 of 1783 bp. Unigenes were annotated by comparing against nr, Swiss-Prot\KEGG\ COG\ KOG\ GO and Pfam databases, and 27,689 unigenes (38.7%) were annotated in at least one database. After bacterial challenge, 183 and 298 genes were identified as remarkably up-regulated or down-regulated, respectively, amongst 481 were associated with 168 pathways, including classical immune-related pathways, such as 'Toll-like receptor signaling', 'the complement cascades', 'MAPK signaling pathway' and 'Apoptosis'. Besides, eight genes which were differently expressed immune-related were confirmed by using quantitative real-time PCR. This study characterized a gene expression pattern for normal and M. lysodeikticus -immunized A. sinica for the first time and sheds new light on the molecular mechanisms thus enabling future efforts on disease control programs in this valuable aquaculture species.

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.

Transcriptome Analysis Based on RNA-Seq in Understanding Pathogenic Mechanisms of Diseases and the Immune System of Fish: A Comprehensive Review

In recent years, with the advent of next-generation sequencing along with the development of various bioinformatics tools, RNA sequencing (RNA-Seq)-based transcriptome analysis has become much more affordable in the field of biological research. This technique has even opened up avenues to explore the transcriptome of non-model organisms for which a reference genome is not available. This has made fish health researchers march towards this technology to understand pathogenic processes and immune reactions in fish during the event of infection. Recent studies using this technology have altered and updated the previous understanding of many diseases in fish. RNA-Seq has been employed in the understanding of fish pathogens like bacteria, virus, parasites, and oomycetes. Also, it has been helpful in unraveling the immune mechanisms in fish. Additionally, RNA-Seq technology has made its way for future works, such as genetic linkage mapping, quantitative trait analysis, disease-resistant strain or broodstock selection, and the development of effective vaccines and therapies. Until now, there are no reviews that comprehensively summarize the studies which made use of RNA-Seq to explore the mechanisms of infection of pathogens and the defense strategies of fish hosts. This review aims to summarize the contemporary understanding and findings with regard to infectious pathogens and the immune system of fish that have been achieved through RNA-Seq technology.