Complement system is one of the most important defense systems of innate immunity, which plays a crucial role in disease defense responses in channel catfish. However, inappropriate and excessive complement activation could lead to potential damage to the host cells. Therefore the complement system is controlled by a set of complement regulatory proteins to allow normal defensive functions, but prevent hazardous complement activation to host tissues. In this study, we identified nine complement regulatory protein genes from the channel catfish genome. Phylogenetic and syntenic analyses were conducted to determine their orthology relationships, supporting their correct annotation and potential functional inferences. The expression profiles of the complement regulatory protein genes were determined in channel catfish healthy tissues and after infection with the two main bacterial pathogens, Edwardsiella ictaluri and Flavobacterium columnare. The vast majority of complement regulatory protein genes were significantly regulated after bacterial infections, but interestingly were generally up-regulated after E. ictaluri infection while mostly down-regulated after F. columnare infection, suggesting a pathogen-specific pattern of regulation. Collectively, these findings suggested that complement regulatory protein genes may play complex roles in the host immune responses to bacterial pathogens in channel catfish.

The Rho family GTPases are a group of small monomeric G proteins, which are molecular switches in signaling pathways. They have been known to regulate a diverse range of cellular processes including actin cytoskeleton rearrangement and microtubule dynamics. In particular, their participations in immune responses are also significant. However, little information of the Rho GTPases is available in teleost including channel catfish, an economically important species and one of the best teleost models forimmunological research. In this study, Rho GTPase genes were identified from channel catfish and well annotated by phylogenetic and syntenic analyses. Their expression profiles were determined in channel catfish healthy tissues and infected tissues. Altogether seven Rho GTPase genes were significantly regulated after bacterial infection, with six genes in the gill after Flavobacterium columnare challenge and two genes in the intestine in response to Edwardsiella ictaluri. All the differentially expressed genes were up-regulated soon after bacterial infection. Different expression patterns between the two experiments were observed, which may be attributed to tissue-specific regulation or pathogen-specific regulation. These results suggested that Rho GTPases play important roles in immune responses to bacterial pathogens, setting a foundation for future investigation on Rho GTPases.

Adaptor proteins non-catalytic region of tyrosine kinase (NCK) and Abelson interactor (ABI) are crucial for disease response. NCK1 was identified to be a candidate gene for enteric septicemia of catfish (ESC) disease resistance, and was speculated to play similar roles during ESC and enteropathogenic Escherichia coli (EPEC) pathogenicity. ABI1 was reported as a positional candidate gene for bacterial cold water disease (BCWD) resistance in rainbow trout. In this study, three NCK genes and six ABI genes were identified in the channel catfish (Ictalurus punctatus) genome and blue catfish (I. furcatus) transcriptome, and annotated by domain structures, phylogenetic and syntenic analyses. Their expression patterns were examined in the intestine and liver of catfish after challenge with Edwardsiella ictaluri. In the intestine, NCK1, ABI2a, ABI2b, ABI3a were differentially expressed after E. ictaluri infection. In the liver, NCK2a, NCK2b, ABI1b, ABI2a, ABI2b were significantly upregulated in ESC susceptible fish. In general, the NCK and ABI genes, with exception of ABI3a gene and NCK1 gene, were expressed at higher levels in susceptible fish after infection than in control fish, but were expressed at lower levels in resistant fish than in the control fish. Taken together, these results support the notion that NCK and ABI genes are involved in disease processes facilitating pathogenesis of the E. ictaluri bacteria.

The lectin pathway of the complement system is characterized by two groups of soluble pattern recognition molecules, mannose-binding lectins (MBLs) and ficolins. These molecules recognize and bind carbohydrates in pathogens and activate complement leading to opsonization, leukocyte activation, and direct pathogen killing. While MBLs have been reported in many fish species, ficolins do not appear to be present in the teleost lineage, despite their importance in invertebrate and higher vertebrate innate immunity. A protein with a similar fibrinogen-like domain, microfibrillar-associated protein 4, MFAP4, is present in fish, albeit with no described immune function. We examined whether MFAP4 genes in fish may potentially act as pathogen receptors in the absence of ficolin. We isolated and characterized five MFAP4 genes from channel catfish. Linkage mapping and phylogenetic analysis indicated that at least three of the catfish MFAP4 genes are tightly clustered on a single chromosome, suggesting that they may have arisen through tandem duplication. Divergent, duplicated families of MFAP4 genes are also present in other teleost species. Expression analysis of the catfish MFAP4 transcripts revealed unique patterns of homeostatic expression among the genes in gill, spleen, skin, liver, and muscle. Expression of the five MFAP4 transcripts showed significant changes in expression as soon as 4h after infection with either Edwardsiella ictaluri or Flavobacterium columnare with modulation of expression continuing up to 7 d following pathogen exposure. Several different tissues and gene-specific patterns were captured and transcript expression changes of >30-fold were observed over the course of the bacterial challenges. Our results suggest a novel role for MFAP4 in teleost immune responses.

Toll-like receptors (TLRs) were the earliest characterized and the most extensively studied pathogen recognition receptors (PRRs). The majority of tetrapod TLR orthologs have been found in teleost fish. In addition, a group of "fish-specific" TLRs have been identified. In catfish, a number of TLR-related sequences have been reported, but systematic phylogenetic analyses have not been conducted. In this study, we conducted phylogenetic and comparative analysis of 20 catfish TLR genes against their counterparts from various species. TLR25 and TLR26 are TLRs identified only in channel catfish. Phylogenetic analyses suggested that four catfish TLR genes have duplicated copies in the genome, i.e., TLR4, TLR5, TLR8, and TLR20. Six fish-specific TLRs were identified, and the vast majority of these belong to the TLR11 subfamily. In healthy catfish tissues, most of the tested TLR genes were ubiquitously expressed although expression levels varied among the 11 tested tissues. We tested nine TLRs for their expression in response to Edwardsiella ictaluri infection. They were significantly up-regulated in the spleen and liver, but down-regulated in the head kidney, suggesting their involvement in the immune responses against the intracellular bacterial pathogen in a tissue-specific manner in catfish, perhaps through rapid migration of phagocytes to infection sites.

Nitric oxide is well known for its roles in immune responses. As such, its synthesizing enzymes have been extensively studied from various species including some teleost fish species. However, the NOS genes have not been characterized in channel catfish (Ictalurus punctatus). In this study, we identified and characterized three NOS genes including one NOS1 and two NOS2 genes in channel catfish. Comparing with the NOS genes from other fish species, the catfish NOS genes are highly conserved in their structural features. Phylogenetic and syntenic analyses allowed determination of NOS1 and NOS2 genes of channel catfish and their orthology relationships. Syntenic analysis, as well as the phylogenetic analysis, indicated that the two NOS2 genes of catfish were lineage-specific duplication. The NOS genes were broadly expressed in most tested tissues, with NOS1 being expressed at the highest levels in the brain, NOS2b1 highly expressed in the skin and gill, and NOS2b2 lowly expressed in most of the tested tissues. The most striking findings of this study was that the expression of the NOS genes are highly regulated after bacterial infection, with time-dependent expression patterns that parallel the migration of macrophages. After Edwardsiella ictaluri challenge, dramatically different responses among the three NOS genes were observed. NOS1 was only significantly in the skin early after infection, while NOS2b1 was rapidly upregulated in gill, but more up-regulated in trunk kidney with the progression of the disease, suggesting such differences in gene expression may be reflective of the migration of macrophages among various tissues of the infected fish. In contrast to NOS1 and NOS2b1, NOS2b2 was normally expressed at very low levels, but it is induced in the brain and liver while significantly down-regulated in most other tissues.

Heat shock protein 90 (HSP90) are highly conserved molecular chaperones, playing a pivotal role in cellular progress. In this study, we reported the characterization of the Hsp90α and Hsp90β genes in Megalobrama amblycephala, the expression profiling during early development, in various healthy tissues and in response to bacterial infection, and the assessment of their adaptive evolution. The Hsp90α cDNA contains an open reading frame (ORF) of 2193 bp encoding 731 amino acids and the Hsp90β cDNA has an ORF of 2184 bp encoding 728 amino acids. Using quantitative real-time PCR (qRT-PCR) analysis, the mRNA of both Hsp90α and Hsp90β reached the highest level at 15th day post-hatch. Using qRT-PCR and Western blot, both Hsp90α and Hsp90β were widely expressed in various healthy tissues and significantly higher in blood than in other tissues. Expression of both Hsp90α and Hsp90β were up-regulated upon bacterial infection and reached the peak level at 4 h post infection. Site model analysis indicated that one positive selection site (T717) in Hsp90α was found, while no positive selection site was observed in Hsp90β. Branch-site model test showed that there were adaptively evolutionary evidences in the branches of Salmoniformes and Gasterosteiformes for Hsp90α gene, and in the branch of Salmoniformes for Hsp90β gene.

Bcl-2 proteins are of vital importance in regulation of apoptosis, and are involved in a number of biological processes such as carcinogenesis and immune responses. Bcl-2 genes have been well studied in mammals, while they are not well investigated in teleost fish including channel catfish, the major aquaculture species in the United States. In this study, we identified 34 bcl-2 genes from the channel catfish genome, and verified their identities by conducting phylogenetic and syntenic analyses. The expression profiles of the bcl-2 genes in response to bacterial infections (Edwardsiella ictaluri and Flavobacterium columnare) and hypoxia stress were determined by performing meta-analysis using the existing RNA-Seq datasets. Differential expressions of bcl-2 genes were observed after bacterial infections and hypoxia treatment, including 22 bcl-2 genes after E. ictaluri infection, 22 bcl-2 genes after F. columnare infection, and 19 bcl-2 genes after hypoxia stress. Overall, the expression of the pro-apoptotic bcl-2 genes were repressed after bacterial infection and hypoxia stress, indicating that bcl-2 genes are potentially involved in the stress response by reducing cell apoptosis. Some bcl-2 genes, such as bcl2b, mcl1a, bmf1, and bnip3, showed different expression pattern during the E. ictaluri and F. columnare infection, suggesting the difference in the pathogenicity of diseases. This work presented the first systematic identification and annotation of bcl-2 genes in catfish, providing essential genomic resources for further immune and physiological studies.

The superfamily of serine protease inhibitors (serpins) are broadly distributed in all kingdoms of life. Serpins play critical roles in an array of fundamental biological processes. In this study, we identified a complete set of 25 serpin genes from channel catfish genome by comprehensive data mining of existing genomic resources. Phylogenetic analysis verified their identities and supported the classification of serpins into six families as in mammals. Extensive comparative genomic analyses suggested that most serpins were conserved among vertebrates, while some were lineage-specific. Analysis of serpin gene expression in mucosal tissues after bacterial infections indicated that serpin genes were regulated in a tissue-specific and time-dependent manner. Distinct expression patterns between infections of the two pathogens were observed, indicating that much more rapid host responses of serpin expression were initiated after ESC infection than after columnaris infection. These studies set the foundation for future studies of host-pathogen interactions.

Chemokines are vital regulators of cell mobilization for immune surveillance, inflammation, and development. Chemokines signal through binding to their receptors that are a superfamily of seven-transmembrane domain G-coupled receptors. Recently, a complete repertoire of both CC and CXC chemokines have been identified in channel catfish, but nothing is known about their receptors. In this study, a set of 29 CC chemokine receptor (CCR) genes and 8 CXC chemokine receptor (CXCR) genes were identified and annotated from the channel catfish genome. Extensive phylogenetic and comparative genomic analyses were conducted to annotate these genes, revealing fish-specific CC chemokine receptors, and lineage-specific tandem duplications of chemokine receptors in the teleost genomes. With 29 genes, the channel catfish genome harbors the largest numbers of CC chemokine receptors among all the genomes characterized. Analysis of gene expression after bacterial infections indicated that the chemokine receptors were regulated in a gene-specific manner. Most differentially expressed chemokine receptors were up-regulated after Edwardsiella ictaluri and Flavobacterium columnare infection. Among which, CXCR3 and CXCR4 were observed to participate in immune responses to both bacterial infections, indicating their potential roles in catfish immune activities. In addition, CXCR3.2 was significantly up-regulated in ESC-susceptible fish, and CXCR4b was mildly induced in ESC-resistant fish, further supporting the significant roles of CXCR3 and CXCR4 in catfish immune responses. CXCR4b and CCR9a were both up-regulated not only after bacterial infection, but also after hypoxia stress, providing the linkage between bacterial infection and low oxygen stresses. These results should be valuable for comparative immunological studies and provide insights into their roles in disease and stress responses.

Apolipoproteins are protein component of plasma lipoproteins. They exert crucial roles in lipoprotein metabolism and serve as enzyme cofactors, receptor ligands, and lipid transfer carriers in mammals. In teleosts, apolipoproteins are also involved in diverse processes including embryonic and ontogenic development, liver and digestive system organogenesis, and innate immunity. In this study, we identified a set of 19 apolipoprotein genes in channel catfish (Ictalurus punctatus). Phylogenetic analysis and syntenic analysis were conducted to determine their identities and evolutionary relationships. The expression signatures of apolipoproteins in channel catfish were determined in healthy tissues and after infections with two major bacterial pathogens, Edwardsiella ictaluri and Flavobacterium columnare. In healthy channel catfish, most apolipoprotein genes exhibited tissue-specific expression patterns in channel catfish. After ESC and columnaris infections, 5 and 7 apolipoprotein genes were differentially expressed respectively, which presented a pathogen-specific and time-dependent pattern of regulation. After ESC infection, three exchangeable apolipoproteins (apoA-IB, apoC-I, and apoE-B) were suppressed in catfish intestine, while two nonexchangeable apolipoproteins (apoB-A and apoB-B) were slightly up-regulated. After columnaris infection, apoB-B, apoD-B, and apoE-A were significantly down-regulated in catfish gill, while apoF, apoL-IV, apoO-like, and apo-14 kDa showed significantly up-regulation. Taken together, these results suggested that apolipoprotein genes may play significant roles in innate immune responses to bacterial pathogens in channel catfish.

Interactions of NF-κB family, IκB family and IKK complex are the key components of NF-κB pathway that is essential for many biological processes including innate and adaptive immunity, inflammation and stress responses. In spite of their importance, systematic analysis of these genes in fish has been lacking. Here we report a systematic study of the NF-κB related genes in channel catfish. Five NF-κB family genes, five IκB family genes and three IKK complex genes were identified in the channel catfish genome. Annotation of these 13 NF-κB related genes was further confirmed by phylogenetic and syntenic analysis. Negative selection was found to play a crucial role in the adaptive evolution of these genes. Expression profiles of NF-κB related genes after Flavobacterium columnare (columnaris) infection were determined by analysis of the existing RNA-Seq dataset. The majority of NF-κB related genes were significantly regulated in mucosal tissues of gill, skin and intestine after columnaris infection, indicating their potential involvement in host defense responses. Distinct expression patterns of NF-κB related genes were observed in susceptible and resistant catfish in response to columnaris infection, suggesting that expression of these genes may contribute to the variations in disease resistance/susceptibility of catfish.

Rab genes, encoding a large family of monomeric G-proteins, contain over 60 members in the human genome. They have been recognized as crucial regulators for membrane trafficking including cargo sorting, vesicle formation, budding, motility, docking, fusion, secretory and endocytic pathway of host immune responses. However, little is known of the Rab gene family in teleost fish species. The development of full-length transcripts and whole genome sequences allow the identification and annotation of Rab GTPase gene family in catfish. In this study, a total of 52 Rab genes were identified from catfish cDNA and genome databases. Phylogenetic analysis assigned them into eleven subfamilies. Most Rab GTPases are conserved among vertebrates, though some of which are absent in fish genomes. Analysis of multiple RNA-seq datasets, along with real time PCR analysis revealed up-regulation of 10 Rab genes after bacterial infection. These included Rab3a, Rab4a, Rab4b, Rab5a, Rab5c, Rab7a, Rab9a, Rab11a, Rab11b, and Rab33a. Their up-regulation are temporally and spatially regulated in various tissues, but mostly induced at early stages after infection and in the gill and liver tissues, with the exception of Rab5c that is mostly up-regulated in the head kidney and trunk kidney. The complex pattern of their induced expression suggested both specific and cooperative actions of a these Rab genes in the acute immune responses to bacterial infection.

Tumor suppressor genes are negative regulators of tumor formation. While their anti-tumor functions have been well studied, they have been found to be also involved in immune responses and innate immunity. In this study, 21 tumor suppressor genes in channel catfish (Ictalurus punctatus) were characterized. Phylogenetic and syntenic analyses allowed annotation of all 21 catfish tumor suppressor genes. The expression profiles of the 21 catfish tumor suppressor genes were determined using the RNA-Seq datasets. After Edwardsiella ictaluri infection, expression of five of the 21 tumor suppressor genes was up-regulated at 3 days in the intestine, and four of the 21 genes were up-regulated in the liver 14 days post-infection. With Flavobacterium columnare infection, seven genes were up-regulated in the gill at 48 h post-infection. These results expanded our knowledge on the tumor suppressor genes in teleosts, setting a foundation for future studies to unravel functions of tumor suppressor genes in response to stresses, particularly after bacterial disease infections.

Peptidoglycan recognition proteins (PGRPs) can recognize bacterial cell wall (peptidoglycan) and activate innate immune system. In addition to its function as pathogen recognition receptors (PRRs), PGRPs are also involved in directly killing bacteria, and regulating multiple signaling pathways. Recently, we have reported catfish PRRs including nucleotide-binding domain, leucine-rich repeat containing receptors (NLRs), retinoic acid inducible gene I (RIG-I) like receptors (RLRs), and Toll-like receptors (TLRs). In this study, we identified and characterized the PGRP gene family in channel catfish which included two members, PGLYRP-5 and PGLYRP-6. Phylogenetic analysis, syntenic analysis and protein structural analysis were conducted to determine their identities and evolutionary relationships. In order to gain insight into the roles of PGRPs in catfish innate immune responses, quantitative real-time PCR was used to investigate the expression profiles in catfish healthy tissues and after bacterial infection. Both PGLYRP-5 and PGLYRP-6 were ubiquitously expressed in all 12 healthy tissues, and most highly expressed in gill and spleen, respectively. Distinct expression patterns were observed for PGRPs after infection with Edwardsiella ictaluri and Flavobacterium columnare, both Gram-negative bacteria. After infection with E. ictaluri, both PGLYRP-5 and PGLYRP-6 were significantly down-regulated at a certain time-point, while both genes were generally up-regulated in the gill after infection with F. columnare. Collectively, these findings suggested that PGRPs may play complex roles in the host immune response to bacterial pathogens in catfish.

Outbreaks of columnaris disease (Flavobacterium columnare) are common in wild and cultured freshwater fish worldwide. Disease occurrences, particularly those caused by virulent genomovar II isolates, in aquaculture species such as channel catfish can be devastating. In contrast to other important aquaculture pathogens, little is known about host immune responses to columnaris. Adhesion of F. columnare to gill tissue has been correlated in some previous studies to virulence and host susceptibility. Here, therefore, we conducted the first transcriptomic profiling of host responses to columnaris following an experimental challenge. We utilized Illumina-based RNA-seq expression profiling to examine transcript profiles at three timepoints (4h, 24h, and 48h) in catfish gill after bath immersion infection. Enrichment and pathway analyses of the differentially expressed genes revealed several central signatures following infection. These included the dramatic upregulation of a rhamnose-binding lectin, with putative roles in bacterial attachment and aggregation, suppression of NF-κB signalling via IκBs, BCL-3, TAX1BP1, and olfactomedin 4, and strong induction of IFN-inducible responses including iNOS2b, IFI44, and VHSV genes. Fifteen differentially expressed genes with varying expression profiles by RNA-seq, were validated by QPCR (correlation coefficients 0.85-0.94, p-value <0.001). Our results highlight several putative immune pathways and individual candidate genes deserving of further investigation in the context of development of therapeutic regimens and laying the foundation for selection of resistant catfish lines against columnaris.

Vertebrates including teleost fish have evolved an array of pathogen recognition receptors (PRRs) for detecting and responding to various pathogen-associated molecular patterns (PAMPs), including Toll-like receptors (TLRs), nucleotide-binding domain, leucine-rich repeat containing receptors (NLRs), and the retinoic acid inducible gene I (RIG-I) like receptors (RLRs). As a part of the series of studies targeted to characterize catfish PRRs, we described 22 NLR receptors in the sister contribution. Here in this study, we focused on cytosolic PRRs recognizing nucleotide pathogen-associated molecular patterns (PAMPs) of invading viruses, the retinoic acid-inducible gene I (RIG-I)-like receptors (RLR receptors). Three RLRs with DExD/H domain containing RNA helicases, retinoic acid inducible gene-I (RIG-I), melanoma differentiation-associated gene 5 (MDA5) and laboratory of genetics and physiology 2 (LGP2), were identified from channel catfish, Ictalurus punctatus. The catfish RIG-I encodes 937 amino acids that contains two CARDs, a DExDc, a HELICc and a RD domains. MDA5 encodes 1005 amino acids with all the domains identified for RIG-I. LGP2 encodes 677 amino acids that contain other domains but not the CARD domain at the N-terminus. Phylogenetic analyses of the three genes of catfish showed close clustering with their counterparts from other teleost fish. All the genes were found to be constitutively expressed in various tissues of catfish with minor variations. Channel catfish ovarian cells when infected with channel catfish virus showed significant increase in the transcript abundance of all the three genes. Further, RLR genes showed significant increases in expression in the liver tissue collected at different time-points after bacterial infection as well. The results indicate that the catfish RLRs may play important roles in antiviral and anti-bacterial immune responses.

Chemokines are a superfamily of structurally related chemotactic cytokines exerting significant roles in regulating cell migration and activation. Based on the arrangement of the first four cysteine residues, they are classified into CC, CXC, C and CX3C subfamilies. In this study, a complete set of 64 CC chemokine ligand (CCL) genes was systematically identified, annotated, and characterized from the channel catfish genome. Extensive phylogenetic and comparative genomic analyses supported their annotations, allowing establishment of their orthologies, revealing fish-specific CC chemokines and the expansion of CC chemokines in the teleost genomes through lineage-specific tandem duplications. With 64 genes, the channel catfish genome harbors the largest numbers of CC chemokines among all the genomes characterized to date, however, they fall into 11 distinct CC chemokine groups. Analysis of gene expression after bacterial infections indicated that the CC chemokines were regulated in a gene-specific and time-dependent manner. While only one member of CCL19 (CCL19a.1) was significantly up-regulated after Edwardsiella ictaluri infection, all CCL19 members (CCL19a.1, CCL19a.2 and CCL19b) were significantly induced after Flavobacterium columnare infection. In addition, CCL19a.1, CCL19a.2 and CCL19b were also drastically up-regulated in ESC-susceptible fish, but not in resistant fish, suggesting potential significant roles of CCL19 in catfish immune responses. High expression levels of certain CC appeared to be correlated with susceptibility to diseases and intolerance to hypoxia.

The Janus kinase/signal transducers and activators of transcription (JAK/STAT) signaling pathway is one of the main pleiotropic cascades used to transmit information from extracellular receptors to the nucleus, which results in DNA transcription and expression of genes involved in immunity, proliferation, differentiation, migration, apoptosis, and cell survival. Members of JAK family and STAT family have been extensively studied in different mammalian species because of their important roles in innate and adaptive immune responses. However, they have not been systematically studied among teleost fish species. In this study, five JAK family members and eight STAT family members were identified and characterized from channel catfish. Phylogenetic analysis was conducted to properly annotate these genes. Syntenic analysis was also conducted to establish orthology, and confirm the results from phylogenetic analysis. Compared to mammals, more members of the JAK and STAT family were identified in channel catfish genome. Expression of JAK and STAT family members was detected in healthy catfish tissues, but was induced in gill, liver, and intestine after bacterial challenge. Notably, the significant upregulation of STAT1b gene in catfish liver, gill and intestine after Edwardsiella ictaluri infection supported the notion that high STAT1 expression are involved in defense against pathogens. Collectively, the increased expression of JAK and STAT members in tested tissues suggested their crucial function in defending the host against pathogen invasion.