Search

The catfish industry is the largest sector of U.S. aquaculture production. Given its role in food production, the catfish immune response to industry-relevant pathogens has been extensively studied and has provided crucial information on innate and adaptive immune function during disease progression. To further examine the channel catfish immune system, we performed single-cell RNA sequencing on nuclei isolated from whole spleens, a major lymphoid organ in teleost fish. Libraries were prepared using the 10X Genomics Chromium X with the Next GEM Single Cell 3' reagents and sequenced on an Illumina sequencer. Each demultiplexed sample was aligned to the Coco_2.0 channel catfish reference assembly, filtered, and counted to generate feature-barcode matrices. From whole spleen samples, outputs were analyzed both individually and as an integrated dataset. The three splenic transcriptome libraries generated an average of 278,717,872 reads from a mean 8,157 cells. The integrated data included 19,613 cells, counts for 20,121 genes, with a median 665 genes/cell. Cluster analysis of all cells identified 17 clusters which were classified as erythroid, hematopoietic stem cells, B cells, T cells, myeloid cells, and endothelial cells. Subcluster analysis was carried out on the immune cell populations. Here, distinct subclusters such as immature B cells, mature B cells, plasma cells, γδ T cells, dendritic cells, and macrophages were further identified. Differential gene expression analyses allowed for the identification of the most highly expressed genes for each cluster and subcluster. This dataset is a rich cellular gene expression resource for investigation of the channel catfish and teleost splenic immunome., Competing Interests: The authors have declared that no competing interests exist., (Copyright: This is an open access article, free of all copyright, and may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose. The work is made available under the Creative Commons CC0 public domain dedication.)

Several farmed fish species, including carps, tilapia, salmon, and catfish, have experienced significant economic losses in aquaculture due to motile Aeromonas septicemia caused by Aeromonas hydrophila. In the present study, a novel lytic bacteriophage infecting hypervirulent Aeromonas hydrophila (vAh) was isolated and characterized. This is the first report of a phage against vAh. Phage AhFM11 demonstrated lytic activity against both vAh strains and the A. hydrophila reference strain ATCC 35654. The AhFM11 genome was sequenced and assembled, comprising 168,243 bp with an average G + C content of 41.5%. The genome did not harbor any antibiotic resistance genes. Genomic information along with transmission electron microscopy revealed that phage AhFM11 belongs to the Straboviridae family. Therapeutic application of monophage AhFM11 in fish showed 100% survival in injection, 95% in immersion and 93% in oral feeding of phage top-coated feed. Fish and chicken meat spiked with A. hydrophila and phage showed significant reduction of A. hydrophila. These findings support that phage AhFM11 can be used as a biocontrol agent against vAh as an alternative to antibiotics., (© 2024. This is a U.S. Government work and not under copyright protection in the US; foreign copyright protection may apply.)

White bass ( Morone chrysops ) are a popular sportfish throughout the southern United States, and one parent of the commercially-successful hybrid striped bass ( M. chrysops ♂ x M. saxatilis ♀). Currently, white bass are cultured using diets formulated for other carnivorous fish, such as largemouth bass ( Micropterus salmoides ) or hybrid striped bass and contain a significant percentage of marine fish meal. Since there are no studies regarding the utilization of alternative proteins in this species, we evaluated the global gene expression of white bass fed diets in which fish meal was partially or totally replaced by various combinations of soybean meal, poultry by-product meal, canola meal, soy protein concentrate, wheat gluten, or a commercial protein blend (Pro-Cision™). Six isonitrogenous (40% protein), isolipidic (11%), and isocaloric (17.1 kJ/g) diets were formulated to meet the known nutrient and energy requirements of largemouth bass and hybrid striped bass using nutrient availability data for most of the dietary ingredients. One of the test diets consisted exclusively of plant protein sources. Juvenile white bass (40.2 g initial weight) were stocked into a flow-through aquaculture system (three tanks/diet; 10 fish/tank) and fed the test diets twice daily to satiation for 60 days. RNA sequencing and bioinformatic analyses revealed significant differentially expressed genes between all test diets when compared to fish meal control. A total of 1,260 differentially expressed genes were identified, with major ontology relating to cell cycle and metabolic processes as well as immune gene functions. This data will be useful as a resource for future refinements to moronid diet formulation, as marine fish meal becomes limiting and plant ingredients are increasingly added as a reliable protein source., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest. The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision., (Copyright © 2024 Fuller, Abernathy, Sankappa, Beck, Rawles, Green, Rosentrater, McEntire, Huskey and Webster.)

The larval waste, exoskeleton shedding, and leftover feed components of the black soldier fly and its larvae make up the by-product known as frass. In this study, we subjected channel catfish ( Ictalurus punctatus ) to a 10-week feeding trial to assess how different dietary amounts of frass inclusion would affect both systemic and mucosal tissue gene expression, especially in regard to growth and immune-related genes. Fish were divided in quadruplicate aquaria, and five experimental diets comprising 0, 50, 100, 200, and 300 g of frass per kilogram of feed were fed twice daily. At the end of the trial, liver, head kidney, gill, and intestine samples were collected for gene expression analyses. First, liver and intestine samples from fish fed with a no frass inclusion diet (control), low-frass (50 g/kg) inclusion diet, or a high-frass (300 g/kg) inclusion diet were subjected to Illumina RNA sequencing to determine global differential gene expression among diet groups. Differentially expressed genes (DEGs) included the upregulation of growth-related genes such as glucose-6-phosphatase and myostatin, as well as innate immune receptors and effector molecules such as toll-like receptor 5, apolipoprotein A1, C-type lectin, and lysozyme. Based on the initial screenings of low/high frass using RNA sequencing, a more thorough evaluation of immune gene expression of all tissues sampled, and all levels of frass inclusion, was further conducted. Using targeted quantitative PCR panels for both innate and adaptive immune genes from channel catfish, differential expression of genes was identified, which included innate receptors (TLR1, TLR5, TLR9, and TLR20A), proinflammatory cytokines (IL-1β type a, IL-1β type b, IL-17, IFN-γ, and TNFα), chemokines (CFC3 and CFD), and hepcidin in both systemic (liver and head kidney) and mucosal (gill and intestine) tissues. Overall, frass from black soldier fly larvae inclusion in formulated diets was found to alter global gene expression and activate innate and adaptive immunity in channel catfish, which has the potential to support disease resistance in this species in addition to demonstrated growth benefits., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2024 Sankappa, Lange, Yildirim-Aksoy, Eljack, Kucuktas, Beck and Abernathy.)

The genomes of seven Aeromonas veronii strains isolated from tissues of healthy or diseased channel catfish obtained from Alabama, USA, fish farms were sequenced and annotated. These genome sequences will enable comparative analyses to determine the roles these bacteria play in catfish aquaculture and the development of new preventative or management strategies.

Vibrio parahaemolyticus is a Gram-negative bacterium commonly found in marine and estuarine environments and is endemic among the global shrimp aquaculture industry. V. parahaemolyticus proteins PirA and PirB have been determined to be major virulence factors that contribute significantly to the development of acute hepatopancreatic necrosis disease. Our previous work had demonstrated the lethality of recombinant PirA and PirB proteins to Pacific white shrimp (Liptopenaeus vannamei). To understand the host response to these proteins, recombinant PirA and PirB proteins were administered using a reverse gavage method and individual shrimp were then sampled over time. Shrimp hepatopancreas libraries were generated and RNA sequencing was performed on the control and recombinant PirA/B-treated samples. Differentially expressed genes were identified among the assayed time points. Differentially expressed genes that were co-expressed at the later time points (2-, 4- and 6-h) were also identified and gene associations were established to predict functional physiological networks. Our analysis reveals that the recombinant PirA and PirB proteins have likely initiated an early host response involving several cell survival signaling and innate immune processes., (Published by Elsevier Ltd.)

Utilizing RNA-seq, this study compared the transcriptomic responses of three improved strains (VSel, PSel, and CSel) of rainbow trout fry during acute stages of challenge with infectious hematopoietic necrosis virus (IHNV). The VSel strain has been selected for resistance against the specific strain of IHNV used in our challenge, PSel has undergone selection for utilization of plant-protein based feeds and previously has shown elevated non-specific disease resistance despite no disease related selection pressures, and the final strain, CSel, is a commercial strain that has been domesticated for several years but has not been selected for specific viral disease resistance. Following a 21-day IHNV challenge, Kaplan-Meier survival estimator curves and cumulative percent mortality (CPM) showed significant differences in IHNV resistance across strains: VSel - 19.3 ± 5.0%, PSel - 67. ± 3.03%, CSel - 94.6 ± 4.1% CPM. To evaluate acute responses to IHNV infection, whole blood, as well as samples from the kidney, liver, and intestine, were collected at 0, 4, 12, 24, and 48 h post infection (hpi). Serum lysozyme activity, a marker of non-specific innate immunity, showed strain and temporal effects during the acute infection phase with PSel showing the highest activity at 0 and 48 hpi. Differential gene expression responses were detected, with varying degrees, in all tissues, both between strains, as well as across acute timepoints within strains. The VSel strain showed upregulation for a particular subset of viral recognition genes during early infection timepoints and rather limited upregulation of immune genes later, while maintaining and reactivating metabolic pathways. The CSel strain showed a downregulation of metabolic related genes and a limited upregulation of immune genes, while the PSel strain showed similar downregulation of metabolic genes during acute infection, yet when compared to the CSel strain, showed a more robust innate immune response. Evaluation of upregulated immune response genes, as well as interferon-related genes showed the PSel strain to have the greatest number of uniquely upregulated immune genes in both the kidney and intestine, with CSel and PSel showing a similar number of such genes upregulated in liver. A moderate number of immune response genes were shared between PSel and CSel in all tissues, though both PSel and VSel showed a high number of uniquely overexpressed immune response genes in the kidney, and PSel showed the highest number of uniquely upregulated interferon related genes in the intestine. Overall, the VSel response was unique from the CSel with very little overlap in activated immune responses. Findings from this study highlight the disparity in IHNV resistance among genetic strains of rainbow trout, while identifying molecular mechanisms underlying differences in disease phenotypes. Furthermore, our results on trout strains with distinct selection backgrounds yields comparative insights into the adaptive gains brought about by selection programs for pathogen-specific disease resistance, as well as the non-specific immune enhancement associated with selection for utilization of plant-based diets., (Published by Elsevier Ltd.)

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

The Gram-negative bacterium, Aeromonas hydrophila, has been responsible for extensive losses in the catfish industry for over a decade. Due to this impact, there are ongoing efforts to understand the basic mechanisms that contribute to virulent A. hydrophila (vAh) outbreaks. Recent challenge models demonstrated that vAh cultured in the presence of the iron chelating agent deferoxamine mesylate (DFO) were more virulent to channel catfish (Ictalurus punctatus). Interestingly, differential gene expression of select iron acquisition genes was unremarkable between DFO and non-DFO cultures, posing the question: why the increased virulence? The current work sought to evaluate growth characteristics and protein expression of vAh after the addition of DFO. A comparative proteome analysis revealed differentially expressed proteins among tryptic soy broth (TSB) and TSB + DFO treatments. Upregulated proteins identified among the TSB + DFO treatment were enriched for gene ontology groups including iron ion transport, siderophore transport and siderophore uptake transport, all iron acquisition pathways. Protein-protein interactions were also evaluated among the differentially expressed proteins and predicted that many of the upregulated iron acquisition proteins likely form functional physiological networks. The proteome analysis of the vAh reveals valuable information about the basic biological processes likely leading to increased virulence during iron restriction in this organism., (© The Author(s) 2020. Published by Oxford University Press on behalf of FEMS.)

Columnaris disease is responsible for substantial losses throughout the production of many freshwater fish species. One of the ways in which the bacterium Flavobacterium columnare is so effective in initiating disease is through the formation of biofilms on fish skin and gills. To further explore the interaction between host factors and bacterial cells, we assayed the ability of vertebrate mucus to enhance F. columnare biofilm development. Different concentrations of catfish, tilapia and pig mucus (5-60 µg/ml) increased biofilm growth at varying degrees among F. columnare isolates. Our data suggest that vertebrate mucus acts as a signalling molecule for the development of F. columnare biofilms; however, there are clear disparities in how individual isolates respond to different mucus fractions to stimulate biofilms. The expression of iron acquisition genes among two genomovar II isolates showed that ferroxidase, TonB receptor and the siderophore synthetase gene were all significantly upregulated among F. columnare biofilms. Interestingly, the siderophore acetyltransferase gene was only shown to be significantly upregulated in one of the genomovar II isolates. This work provides insight into our understanding of the interaction between F. columnare and vertebrate mucus, which likely contributes to the growth of planktonic cells and the transition into biofilms., (© 2019 John Wiley & Sons Ltd.)

Flavobacterium columnare causes substantial losses among cultured finfish species. The Gram-negative bacterium is an opportunistic pathogen that manifests as biofilms on the host's mucosal surfaces as the disease progresses. We previously demonstrated that the dominant mucosal IgM antibody response to F. columnare is to the chaperone protein DnaK that is found in the extracellular fraction. To establish the efficacy of using recombinant protein technology to develop a new vaccine against columnaris disease, we are reporting on two consecutive years of vaccine trials using a recombinant F. columnare DnaK protein (rDnaK). In year one, three groups of channel catfish ( n = 300) were immunized by bath immersion with a live attenuated F. columnare isolate, rDnaK or sham immunized. After 6 weeks, an F. columnare laboratory challenge showed a significant increase in survival (>30%) in both the live attenuated and rDnaK vaccines when compared to the non-immunized control. A rDnaK-specific ELISA revealed significant levels of mucosal IgM antibodies in the skin of catfish immunized with rDnaK at 4- and 6-weeks post immunization. In the second year, three groups of channel catfish ( n = 300) were bath immunized with rDnaK alone or with rDnaK after a brief osmotic shock or sham immunized. After 6 weeks a laboratory challenge with F. columnare was conducted and showed a significant increase in survival in the rDnaK (> 25%) and in rDnaK with osmotic shock (>35%) groups when compared to the non-immunized control. The rDnaK-specific ELISA demonstrated significant levels of mucosal IgM antibodies in the skin of catfish groups immunized with rDnaK at 4- and 6-weeks post immunization. To further understand the processes which have conferred immune protection in the rDnaK group, we conducted RNA sequencing of skin samples from the non-immunized ( n = 6) and rDnaK treated channel catfish at 1-week ( n = 6) and 6 weeks ( n = 6) post immunization. Significantly altered gene expression was identified and results will be discussed. Work to further enhance the catfish immune response to F. columnare rDnaK is underway as this protein remains a promising candidate for additional optimization and experimental trials in a production setting.

Reformulation of aquafeeds in salmonid diets to include more plant proteins is critical for sustainable aquaculture. However, increasing plant proteins can lead to stunted growth and enteritis. Toward an understanding of the regulatory mechanisms behind plant protein utilization, directional RNA sequencing of liver tissues from a rainbow trout strain selected for growth on an all plant-protein diet and a control strain, both fed a plant diet for 12 weeks, were utilized to construct long noncoding RNAs. Antisense long noncoding RNAs were selected for differential expression and functional analyses since they have been shown to have regulatory actions within a genome. A total of 142 unique antisense long noncoding RNAs were differentially expressed between strains, 60 of which could be mapped to a gene. Genes underlying these noncoding RNAs are indicated in lipid metabolism and immunity. Six noncoding transcripts were also found to overlap with differentially expressed protein-coding genes, all of which were co-expressed. Associating variation in regulatory elements between rainbow trout strains with differing tolerance to plant-protein diets will assist in future studies toward increased gains throughout carnivorous aquaculture.

Columnaris disease, caused by Flavobacterium columnare, severely impacts the production of freshwater finfish species. Therefore, efforts to better understand the biological processes of F. columnare, including the formation of biofilms and their contribution to disease, are ongoing. In this study, we incubated F. columnare cultures with channel catfish mucus and used high-throughput RNA sequencing to evaluate global changes in gene expression. Our data show that mucus activates in vitro biofilm formation. The analysis of F. columnare transcriptomes after the addition of mucus revealed significant differentially expressed genes (DEGs) between the planktonic and biofilm states. DEGs common among all biofilms were enriched for gene ontology groups including signal transduction, ligand binding and cellular homeostasis and are likely necessary for biofilm formation. Iron acquisition systems included TonB-dependent receptor and ferroxidase genes were expressed among all biofilms, while siderophore synthesis genes were only expressed in mucus-stimulated biofilms. The current analysis of F. columnare transcriptomes adds valuable information about the basic biological processes that occur during the planktonic and biofilm states. This work serves as a basis for future studies on understanding how biofilms are established and how they contribute to disease progression.

Finding suitable alternative protein sources for diets of carnivorous fish species remains a major concern for sustainable aquaculture. Through genetic selection, we created a strain of rainbow trout that outperforms parental lines in utilizing an all-plant protein diet and does not develop enteritis in the distal intestine, as is typical with salmonids on long-term plant protein-based feeds. By incorporating this strain into functional analyses, we set out to determine which genes are critical to plant protein utilization in the absence of gut inflammation. After a 12-week feeding trial with our selected strain and a control trout strain fed either a fishmeal-based diet or an all-plant protein diet, high-throughput RNA sequencing was completed on both liver and muscle tissues. Differential gene expression analyses, weighted correlation network analyses and further functional characterization were performed. A strain-by-diet design revealed differential expression ranging from a few dozen to over one thousand genes among the various comparisons and tissues. Major gene ontology groups identified between comparisons included those encompassing central, intermediary and foreign molecule metabolism, associated biosynthetic pathways as well as immunity. A systems approach indicated that genes involved in purine metabolism were highly perturbed. Systems analysis among the tissues tested further suggests the interplay between selection for growth, dietary utilization and protein tolerance may also have implications for nonspecific immunity. By combining data from differential gene expression and co-expression networks using selected trout, along with ontology and pathway analyses, a set of 63 candidate genes for plant diet tolerance was found. Risk loci in human inflammatory bowel diseases were also found in our datasets, indicating rainbow trout selected for plant-diet tolerance may have added utility as a potential biomedical model.

The exposure of fish to environmental free-living microbes and its effect on early colonization in the gut have been studied in recent years. However, little is known regarding how the host and environment interact to shape gut communities during early life. Here, we tested whether the early microbial exposure of tilapia larvae affects the gut microbiota at later life stages. The experimental period was divided into three stages: axenic, probiotic and active suspension. Axenic tilapia larvae were reared either under conventional conditions (active suspension systems) or exposed to a single strain probiotic (Bacillus subtilis) added to the water. Microbial characterization by Illumina HiSeq sequencing of 16S rRNA gene amplicons showed the presence of B. subtilis in the gut during the seven days of probiotic application. Although B. subtilis was no longer detected in the guts of fish exposed to the probiotic after day 7, gut microbiota of the exposed tilapia larvae remained significantly different from that of the control treatment. Compared with the control, fish gut microbiota under probiotic treatment was less affected by spatial differences resulting from tank replication, suggesting that the early probiotic contact contributed to the subsequent observation of low inter-individual variation.

RNA sequencing (RNA-Seq) is becoming the standard for transcriptome analysis. Removal of contaminating ribosomal RNA (rRNA) is a priority in the preparation of libraries suitable for sequencing. These methods have been well documented in mammals but typically require some optimization for lower vertebrates. Three commercial kits, including Dynabeads mRNA Purification Kit, RiboMinus Eukaryote System v2, and Ribo-Zero Gold rRNA Removal Kit were examined for the ability to remove rRNAs from rainbow trout (Oncorhynchus mykiss) RNA isolations. Total RNA was isolated from liver and muscle tissue samples (n = 24) and rRNAs removed using one of the three kits. Samples were analyzed visually on the Agilent Bioanalyzer and by Illumina RNA-seq, screening for Oncorhynchus rRNAs. There were significant differences between the kits in regards to their ability to remove rRNA, ranging from 2.74% - 10.94% rRNA sequences left behind per kit on average. Using the Bioanalyzer to evaluate ribosomal contamination in rRNA-depleted samples for RNA-Seq was good for detecting samples with higher concentrations of rRNA (>5%), but not very accurate at lower levels. Although all three kits were able to remove a substantial portion of the rRNA from different fish tissues, the Ribo-Zero Gold rRNA Removal Kit eliminated significantly more contaminating ribosomal RNAs than the others.

Growth hormone (GH) transgenic salmon possesses markedly increased metabolic rate, appetite, and feed conversion efficiency, as well as an increased ability to compete for food resources. Thus, the ability of GH-transgenic fish to withstand periods of food deprivation as occurs in nature is potentially different than that of nontransgenic fish. However, the physiological and genetic effects of transgenic GH production over long periods of food deprivation remain largely unknown. Here, GH-transgenic coho salmon (Oncorhynchus kisutch) and nontransgenic, wild-type coho salmon were subjected to a 3-month food deprivation trial, during which time performance characteristics related to growth were measured along with proximate compositions. To examine potential genetic effects of GH-transgenesis on long-term food deprivation, a group of genes related to muscle development and liver metabolism was selected for quantitative PCR analysis. Results showed that GH-transgenic fish lose weight at an increased rate compared to wild-type even though proximate compositions remained relatively similar between the groups. A total of nine genes related to muscle physiology (cathepsin, cee, insulin-like growth factor, myostatin, murf-1, myosin, myogenin, proteasome delta, tumor necrosis factor) and five genes related to liver metabolism (carnitine palmitoyltransferase, fatty acid synthase, glucose-6-phosphatase, glucose-6-phosphate dehydrogenase, glucokinase) were shown to be differentially regulated between GH-transgenic and wild-type coho salmon over time. These genetic and physiological responses assist in identifying differences between GH-transgenic and wild-type salmon in relation to fitness effects arising from elevated growth hormone during periods of long-term food shortage.

Trisomy mapping is a powerful method for assigning genes to chicken microchromosome 16 (GGA 16). The single chicken nucleolar organizer region (NOR), the 2 major histocompatibility complex regions (MHC-Y and MHC-B), and CD1 genes were all previously assigned to GGA 16 using trisomy mapping. Here, we combined array comparative genomic hybridization with trisomy mapping to screen unassigned genomic scaffolds (consigned temporarily to chrUn_random) for sequences originating from GGA 16. A number of scaffolds mapped to GGA 16. Among these were scaffolds that contain genes for olfactory (OR) and cysteine-rich domain scavenger (SRCR) receptors, along with a number of genes that encode putative immunoglobulin-like receptors and other molecules. We used high-resolution cytogenomic analyses to confirm assignment of OR and SRCR genes to GGA 16 and to pinpoint members of these gene families to the q-arm in partially overlapping regions between the centromere and the NOR. Southern blots revealed sequence polymorphism within the OR/SRCR region and linkage with the MHC-Y region, thereby providing evidence for conserved linkage between OR genes and the MHC within birds. This work localizes OR genes to the vicinity of the chicken MHC and assigns additional genes, including immune defense genes, to GGA 16.

Salmonella enterica serovar enteritidis is an enteric bacterium that can contaminate chicken eggs and meat, resulting in production losses and consumer illness. To provide insight into the systemic metabolic effects of S. enteritidis infection, liver samples were harvested 10-days postinfection from broiler hens. Hepatic global gene expression levels were assessed using a chicken 44K Agilent microarray. Forty-four genes were differentially expressed at a significance level of q value < 0.05. One hundred eighty-three genes were differentially expressed at a suggestive significance level of q value < 0.1. A predominance of downregulation existed among significantly differentially expressed genes. Cell cycle and metabolism networks were created from the differentially expressed genes. Mitochondria-mediated apoptosis, electron transport, peptidase activity, vein constriction, cell differentiation, IL-2 signaling, Jak-Stat signaling, B-cell receptor signaling, GDP/GTP exchange, and protein recycling were among the functions of the differentially expressed genes that were down-regulated in response to S. enteritidis. The effects of S. enteritidis infection on the liver transcriptome profiles of broilers reflect a predominance of downregulation of genes with cell cycle and metabolic functions. The most pronounced response was the downregulation of genes that function in metabolic pathways, inflammation, and mitochondria-mediated apoptosis. These results provide insight into important systemic metabolic mechanisms that are active in the chicken liver in response to S. enteritidis infection at 10-days postinfection., (Copyright © 2012 Wiley Periodicals, Inc.)

Background: Salmonella enterica serovar Typhimurium is a major foodborne pathogen worldwide. S. Typhimurium encodes type III secretion systems via Salmonella pathogenicity islands (SPI), producing the major effector proteins of virulence. Previously, we identified two genes of Salmonella pyruvate metabolism that were up-regulated during chicken cell infection: pyruvate formate lyase I (pflB) and bifunctional acetaldehyde-CoA/alcohol dehydrogenase (adhE). We were therefore interested in examining the role these genes may play in the transmission of Salmonella to humans., Methods: Mutant strains of Salmonella with single gene deletions for pflB and adhE were created. Invasion and growth in human HCT-8 intestinal epithelial cells and THP-1 macrophages was examined. Quantitative PCR was performed on 19 SPI-1 genes., Results: In HCT-8 cells, both mutant strains had significantly higher intracellular counts than the wild-type from 4 to 48 h post-infection. Various SPI-1 genes in the mutants were up-regulated over the wild-type as early as 1 h and lasting until 24 h post-infection. In THP-1 cells, no significant difference in internal Salmonella counts was observed; however, SPI-1 genes were largely down-regulated in the mutants during the time-course of infection. We also found five SPI-1 genes - hilA, hilC hilD, sicP and rtsA - which were up-regulated in at least one of the mutant strains in log-phase broth cultures alone. We have therefore identified a set of SPI-1 virulence genes whose regulation is effected by the central metabolism of Salmonella.

Construction of high-density genetic linkage maps is crucially important for quantitative trait loci (QTL) studies, and they are more useful when integrated with physical maps. Such integrated maps are valuable genome resources for fine mapping of QTL, comparative genomics, and accurate and efficient whole-genome assembly. Previously, we established both linkage maps and a physical map for channel catfish, Ictalurus punctatus, the dominant aquaculture species in the United States. Here we added 2030 BAC end sequence (BES)-derived microsatellites from 1481 physical map contigs, as well as markers from singleton BES, ESTs, anonymous microsatellites, and SNPs, to construct a second-generation linkage map. Average marker density across the 29 linkage groups reached 1.4 cM/marker. The increased marker density highlighted variations in recombination rates within and among catfish chromosomes. This work effectively anchored 44.8% of the catfish BAC physical map contigs, covering ~52.8% of the genome. The genome size was estimated to be 2546 cM on the linkage map, and the calculated physical distance per centimorgan was 393 Kb. This integrated map should enable comparative studies with teleost model species as well as provide a framework for ordering and assembling whole-genome scaffolds.

The ciliate parasite Ichthyophthirius multifiliis (Ich) infects many freshwater fish, causing white spot disease that leads to heavy economic losses to aquaculture and ornamental industries. Despite its economic importance, molecular studies examining fundamental processes such as life stage regulation and infectivity have been scarce. In this study, we developed an oligo microarray platform using all available I. multifiliis expressed sequence tag (EST) information as well as probes designed through comparative genomics to other protozoa. Gene expression profiling for developmental and virulence factors was conducted using this platform. For the developmental study, the microarray was used to examine gene expression profiles between the three major life stages of Ich: infective theront, parasitic trophont, and reproductive tomont. A total of 135 putative I. multifiliis genes were found to be differentially expressed among all three life-stages. Examples of differentially expressed transcripts among life stages include immobilization antigens and epiplasmin, as well as various other transcripts involved in developmental regulation and host-parasite interactions. I. multifiliis has been shown to lose infectivity at later cell divisions potentially due to cellular senescence. Therefore, the microarray was also used to explore expression of senescence-associated genes as related to the passage number of the parasite. In this regard, comparison between tomont early and late passages yielded 493 differently expressed genes; 1478 differentially expressed genes were identified between trophont early and late passages. The EST-derived oligo microarray represents a first generation array of this ciliate and provided reproducible expression data as validated by quantitative RT-PCR., (Copyright © 2011 Elsevier Inc. All rights reserved.)

Infection and inflammation are often accompanied by oxidative stress caused by the accumulation of reactive oxygen species which can be deleterious to the health of the host. Antioxidant defense mechanisms and components are crucial in limiting cellular and tissue-level damage and restoring homeostasis. In mammals, calreticulin is a 46-kDa multifunctional calcium binding protein of the endoplasmic reticulum that has many critical functions in the eukaryotic cell including regulation of intracellular calcium homoeostasis, lectin binding and chaperoning, and oxidative stress responses. In previous studies from our lab, the calreticulin gene was observed to be strongly upregulated in catfish during challenge with infectious Gram-negative bacteria. However, little is known about the function of this gene in teleost fish. The objective of this study, therefore, was to characterize the calreticulin gene from channel catfish, to determine its genomic organization, to profile its patterns of tissue expression, and to establish its potential for physiological antioxidant and immune responses in catfish after bacterial infection with Edwardsiella ictaluri and iron treatment. Our results indicate that there are at least three calreticulin related genes in the catfish genome. The three calreticulin genes are widely expressed in various tissues under homeostatic conditions and their expression showed significant upregulation following infection and/or iron level changes., (Copyright © 2010 Elsevier Ltd. All rights reserved.)

Ceruloplasmin is a serum ferroxidase that carries more than 90% of the copper in plasma and has documented roles in iron homeostasis as well as antioxidative functions. In our previous studies, it has been shown that the ceruloplasmin gene is strongly up-regulated in catfish during challenge with Edwardsiella ictaluri. However, little is known about the function of this gene in teleost fish. The objective of this study, therefore, was to characterize the ceruloplasmin gene from channel catfish, determine its genomic organization, profile its patterns of tissue expression, and establish its potential for physiological antioxidant responses in catfish after bacterial infection with E. ictaluri and iron treatment. The genomic organization suggested that the catfish ceruloplasmin gene had 20 exons and 19 introns, encoding 1074 amino acids. Exon sizes of the catfish ceruloplasmin gene were close to or identical with mammalian and zebrafish homologs. Further phylogenetic analyses suggested that the gene was highly conserved through evolution. The catfish ceruloplasmin gene was mapped to both the catfish physical map and linkage map. The catfish ceruloplasmin gene was mainly expressed in liver with limited expression in other tissues, and it was significantly up-regulated in the liver after bacterial infection alone or after co-injection with bacteria and iron-dextran, while expression was not significantly induced with iron-dextran treatment alone., (Copyright © 2011 Elsevier Ltd. All rights reserved.)

Background: Genome annotation projects, gene functional studies, and phylogenetic analyses for a given organism all greatly benefit from access to a validated full-length cDNA resource. While increasingly common in model species, full-length cDNA resources in aquaculture species are scarce., Methodology and Principal Findings: Through in silico analysis of catfish (Ictalurus spp.) ESTs, a total of 10,037 channel catfish and 7,382 blue catfish cDNA clones were identified as potentially encoding full-length cDNAs. Of this set, a total of 1,169 channel catfish and 933 blue catfish full-length cDNA clones were selected for re-sequencing to provide additional coverage and ensure sequence accuracy. A total of 1,745 unique gene transcripts were identified from the full-length cDNA set, including 1,064 gene transcripts from channel catfish and 681 gene transcripts from blue catfish, with 416 transcripts shared between the two closely related species. Full-length sequence characteristics (ortholog conservation, UTR length, Kozak sequence, and conserved motifs) of the channel and blue catfish were examined in detail. Comparison of gene ontology composition between full-length cDNAs and all catfish ESTs revealed that the full-length cDNA set is representative of the gene diversity encoded in the catfish transcriptome., Conclusions: This study describes the first catfish full-length cDNA set constructed from several cDNA libraries. The catfish full-length cDNA sequences, and data gleaned from sequence characteristics analysis, will be a valuable resource for ongoing catfish whole-genome sequencing and future gene-based studies of function and evolution in teleost fishes.

Ferritins are the major iron storage protein in the cytoplasm of cells, responsible for regulating levels of intracellular iron. Ferritin genes are widely distributed in both prokaryotes and eukaryotes. In mammals, ferritin molecules are composed of heavy- (H) and light- (L) chain subunits; amphibian genomes contain three ferritin-type genes (H; middle, M; and L subunits); and teleost genomes to date contain H and M subunits. The objective of this study was to characterize the ferritin H gene in channel catfish (Ictalurus punctatus) to determine its genomic organization and copy numbers, to determine its patterns of tissue expression, and to establish if it is involved in defense responses of catfish after bacterial infection. The catfish ferritin H gene was completely sequenced and characterized, using both mRNA and genomic DNA. Catfish ferritin H gene has a full-length mRNA sequence of 999 bp, an open reading frame of 534 bp, and 4,704 bp genomic DNA sequence. Catfish ferritin H has a 5 exon and 4 intron genetic organization, containing a long 5'-untranslated region, which shares high similarity with mammalian and zebrafish genes. Based on phylogenetic analyses, the catfish ferritin H gene is highly conserved throughout evolution. Southern blot analysis suggested that the ferritin H gene has only one copy in the catfish genome. The catfish ferritin H gene was widely expressed in various healthy tissues. The catfish ferritin H gene was significantly up-regulated in the liver after intraperitoneal injection of iron dextran and coinjection of Edwardsiella ictaluri and iron dextran treatment, suggesting its role in iron metabolism and immunity., ((c) 2010 Wiley-Liss, Inc.)

Alternative splicing (AS) is a mechanism by which the coding diversity of the genome can be greatly increased. Rates of AS are known to vary according to the complexity of eukaryotic species potentially explaining the tremendous phenotypic diversity among species with similar numbers of coding genes. Little is known, however, about the nature or rate of AS in teleost fish. Here, we report the characteristics of AS in teleost fish and classification and frequency of five canonical AS types. We conducted both same-species and cross-species analysis utilizing the Genome Mapping and Alignment Program (GMAP) and an AS pipeline (ASpipe) to study AS in four genome-enabled species (Danio rerio, Oryzias latipes, Gasterosteus aculeatus, and Takifugu rubripes) and one species lacking a complete genome sequence, Ictalurus punctatus. AS frequency was lowest in the highly duplicated genome of zebrafish (17% of mapped genes). The compact genome of the pufferfish showed the highest occurrence of AS (approximately 43% of mapped genes). An inverse correlation between AS frequency and genome size was consistent across all analyzed species. Cross-species comparisons utilizing zebrafish as the reference genome allowed the identification of additional putative AS genes not revealed by zebrafish transcripts. Approximately, 50% of AS genes identified by same-species comparisons were shared among two or more species. A searchable website, the Teleost Alternative Splicing Database, was created to allow easy identification and visualization of AS transcripts in the studied teleost genomes. Our results and associated database should further our understanding of alternative splicing as an important functional and evolutionary mechanism in the genomes of teleost fish.

Matrix metalloproteinase-13 (MMP-13), referred to as collagenase-3, is a proteolytic enzyme that plays a key role in degradation and remodelling of host extracellular matrix proteins. The objective of this study was to characterize the MMP-13 gene in channel catfish, and to determine its pattern of expression in various healthy tissues and during embryogenesis. Since MMP-13 has been shown to have importance in tissue remodelling and some pathological processes, we further studied its involvement in the defense responses of catfish after bacterial infection. The channel catfish MMP-13 cDNA contains an open reading frame of 1416bp encoding 471 amino acids. Using RT-PCR analysis, MMP-13 was widely expressed in various health tissues. Using quantitative real-time PCR analysis, expression of MMP-13 gene was up-regulated by bacterial infection. During normal embryological development, MMP-13 expression was slightly increased in the first day post-fertilization and sharply up-regulated from 1-day post-fertilization through hatching., (Copyright (c) 2010 Elsevier Ltd. All rights reserved.)

Background: Through the Community Sequencing Program, a catfish EST sequencing project was carried out through a collaboration between the catfish research community and the Department of Energy's Joint Genome Institute. Prior to this project, only a limited EST resource from catfish was available for the purpose of SNP identification., Results: A total of 438,321 quality ESTs were generated from 8 channel catfish (Ictalurus punctatus) and 4 blue catfish (Ictalurus furcatus) libraries, bringing the number of catfish ESTs to nearly 500,000. Assembly of all catfish ESTs resulted in 45,306 contigs and 66,272 singletons. Over 35% of the unique sequences had significant similarities to known genes, allowing the identification of 14,776 unique genes in catfish. Over 300,000 putative SNPs have been identified, of which approximately 48,000 are high-quality SNPs identified from contigs with at least four sequences and the minor allele presence of at least two sequences in the contig. The EST resource should be valuable for identification of microsatellites, genome annotation, large-scale expression analysis, and comparative genome analysis., Conclusions: This project generated a large EST resource for catfish that captured the majority of the catfish transcriptome. The parallel analysis of ESTs from two closely related Ictalurid catfishes should also provide powerful means for the evaluation of ancient and recent gene duplications, and for the development of high-density microarrays in catfish. The inter- and intra-specific SNPs identified from all catfish EST dataset assembly will greatly benefit the catfish introgression breeding program and whole genome association studies.

Transferrin is important in iron metabolism and has been reported to be involved in disease defence responses after bacterial infection. In this study, we identified, sequenced, and characterized the transferrin gene from channel catfish, Ictalurus punctatus. The catfish transferrin gene was similar to those of other vertebrate species with 17 exons and 16 introns. Sequence analysis indicated the presence of the two duplicated lobes, each containing two sub-domains separated by a cleft harboring the iron-binding site, suggesting their structural conservation. The channel catfish transferrin cDNA encodes 679 amino acids with 42-56% similarity to known transferrin genes from various species. Southern blot analysis suggested the presence of two copies of the transferrin gene in the catfish genome, perhaps arranged in a tandem fashion. The catfish transferrin gene was mapped to a catfish BAC-based physical map. The catfish transferrin gene was highly expressed in the liver, but expression was low in most other tested tissues. Transferrin expression was significantly up-regulated after infection with Edwardsiella ictaluri, the causative agent of enteric septicemia of catfish. Such induction was also found with co-injection of iron-dextran and E. ictaluri, while transferrin expression was not significantly induced with the injection of iron-dextran alone., (Copyright 2009 Elsevier Ltd. All rights reserved.)

Background: Comparative mapping is a powerful tool to transfer genomic information from sequenced genomes to closely related species for which whole genome sequence data are not yet available. However, such an approach is still very limited in catfish, the most important aquaculture species in the United States. This project was initiated to generate additional BAC end sequences and demonstrate their applications in comparative mapping in catfish., Results: We reported the generation of 43,000 BAC end sequences and their applications for comparative genome analysis in catfish. Using these and the additional 20,000 existing BAC end sequences as a resource along with linkage mapping and existing physical map, conserved syntenic regions were identified between the catfish and zebrafish genomes. A total of 10,943 catfish BAC end sequences (17.3%) had significant BLAST hits to the zebrafish genome (cutoff value

The complement system in vertebrates plays a crucial role in immune defense via recognition and removal of pathogens. Complement is tightly regulated by a group of both soluble and cell-associated proteins. Complement factor I is a soluble serine protease that regulates multiple pathways in complement activation. In this work, a complement factor I transcript was isolated and sequenced from channel catfish (Ictalurus punctatus) liver after screening expressed sequence tags. The full-length cDNA is comprised of 2284bp in length, encoding a polypeptide of 668 amino acids. The complement factor I protein was found to be well conserved, with similar domain structures and architecture from fish to mammals. The catfish complement factor I exists as a single-copied gene in the catfish genome. Expression analysis revealed that the catfish complement factor I is constitutively expressed in all tissues and leukocyte cell lines tested, indicating its importance as a regulatory enzyme throughout channel catfish. While expression of complement factor I is often found to be in the liver in mammals, it is constitutively expressed in channel catfish and carp throughout in various tissues and organs., (2009 Elsevier Ltd.)

Polyadenylation of eukaryotic transcripts is usually restricted to mRNA, providing transcripts with stability from degradation by nucleases. Conversely, an RNA degradation pathway can be signaled through poly (A) tailing in prokaryotic, archeal, and organellar biology. Recently polyadenylated transcripts have also been discovered in rRNA in some eukaryotes including humans and yeast. Here we report the discovery of polyadenylated rRNAs in the ciliate teleost parasite Ichthyophthirius multifiliis, an important fish pathogen. Through large-scale analysis of ESTs, a large contig composed of the 28S rRNA with poly (A) tails was identified. Analysis using multiple sequence alignments revealed four potential polyadenylation sites including three internal regions and the 3? end of the rRNA. Further analysis using a polyadenylation test, re-sequencing, and gene-specific PCR using primers flanking the presumed poly (A) sites confirmed the presence of polyadenylated rRNA in this parasite. The functions of polyadenylation of rRNA in this organism are largely unknown at present, but the presence of internal polyadenylation sites, along with the presence of truncated segments of the rRNA, may suggest a role of the polyadenylation in the degradation pathway, a function typical of prokaryotes, archaea, and organelles. These results are in congruence with reports of a similar phenomenon in humans and yeast.

The NLRs (nucleotide-binding domain and leucine-rich repeat containing family receptors) are a recently identified family of pattern recognition receptors in vertebrates. Several subfamilies of NLRs have been characterized in human, mouse, and zebrafish, but studies of NLRs in other species, especially teleost species, have been lacking. Here we report characterization of five NLRs from channel catfish: NOD1, NOD2, NLRC3, NLRC5, and NLRX1. Structural analysis indicated that the genes were organized in a similar fashion as in the mammals and in zebrafish. Phylogenetic analysis suggested that they were orthologous to the NOD-like subfamily of NLRs. All five NOD-like genes exist as a single copy gene in the catfish genome. Hybridization of gene-specific probes allowed mapping of three NLR genes to the catfish physical map, laying a foundation for genome characterization and for establishing orthologies with NLR genes from other species. These genes are widely expressed in various tissues and leukocyte cell lines. While the majority of the NLR genes appeared to be constitutively expressed, NOD1 was induced after infection with a bacterial pathogen, Edwardsiella ictaluri, the causative agent of enteric septicemia of catfish (ESC), suggesting its involvement in immunity against the intracellular pathogen.

A genetic linkage map of the channel catfish genome (N=29) was constructed using EST-based microsatellite and single nucleotide polymorphism (SNP) markers in an interspecific reference family. A total of 413 microsatellites and 125 SNP markers were polymorphic in the reference family. Linkage analysis using JoinMap 4.0 allowed mapping of 331 markers (259 microsatellites and 72 SNPs) to 29 linkage groups. Each linkage group contained 3-18 markers. The largest linkage group contained 18 markers and spanned 131.2 cM, while the smallest linkage group contained 14 markers and spanned only 7.9 cM. The linkage map covered a genetic distance of 1811 cM with an average marker interval of 6.0 cM. Sex-specific maps were also constructed; the recombination rate for females was 1.6 times higher than that for males. Putative conserved syntenies between catfish and zebrafish, medaka, and Tetraodon were established, but the overall levels of genome rearrangements were high among the teleost genomes. This study represents a first-generation linkage map constructed by using EST-derived microsatellites and SNPs, laying a framework for large-scale comparative genome analysis in catfish. The conserved syntenies identified here between the catfish and the three model fish species should facilitate structural genome analysis and evolutionary studies, but more importantly should facilitate functional inference of catfish genes. Given that determination of gene functions is difficult in nonmodel species such as catfish, functional genome analysis will have to rely heavily on the establishment of orthologies from model species.

Background: The ciliate protozoan Ichthyophthirius multifiliis (Ich) is an important parasite of freshwater fish that causes 'white spot disease' leading to significant losses. A genomic resource for large-scale studies of this parasite has been lacking. To study gene expression involved in Ich pathogenesis and virulence, our goal was to generate expressed sequence tags (ESTs) for the development of a powerful microarray platform for the analysis of global gene expression in this species. Here, we initiated a project to sequence and analyze over 10,000 ESTs., Results: We sequenced 10,368 EST clones using a normalized cDNA library made from pooled samples of the trophont, tomont, and theront life-cycle stages, and generated 9,769 sequences (94.2% success rate). Post-sequencing processing led to 8,432 high quality sequences. Clustering analysis of these ESTs allowed identification of 4,706 unique sequences containing 976 contigs and 3,730 singletons. These unique sequences represent over two million base pairs (~10% of Plasmodium falciparum genome, a phylogenetically related protozoan). BLASTX searches produced 2,518 significant (E-value < 10-5) hits and further Gene Ontology (GO) analysis annotated 1,008 of these genes. The ESTs were analyzed comparatively against the genomes of the related protozoa Tetrahymena thermophila and P. falciparum, allowing putative identification of additional genes. All the EST sequences were deposited by dbEST in GenBank (GenBank: EG957858-EG966289). Gene discovery and annotations are presented and discussed., Conclusion: This set of ESTs represents a significant proportion of the Ich transcriptome, and provides a material basis for the development of microarrays useful for gene expression studies concerning Ich development, pathogenesis, and virulence.

Background: The eastern oyster, Crassostrea virginica (Gmelin 1791), is an economically important species cultured in many areas in North America. It is also ecologically important because of the impact of its filter feeding behaviour on water quality. Populations of C. virginica have been threatened by overfishing, habitat degradation, and diseases. Through genome research, strategies are being developed to reverse its population decline. However, large-scale expressed sequence tag (EST) resources have been lacking for this species. Efficient generation of EST resources from this species has been hindered by a high redundancy of transcripts. The objectives of this study were to construct a normalized cDNA library for efficient EST analysis, to generate thousands of ESTs, and to analyze the ESTs for microsatellites and potential single nucleotide polymorphisms (SNPs)., Results: A normalized and subtracted C. virginica cDNA library was constructed from pooled RNA isolated from hemocytes, mantle, gill, gonad and digestive tract, muscle, and a whole juvenile oyster. A total of 6,528 clones were sequenced from this library generating 5,542 high-quality EST sequences. Cluster analysis indicated the presence of 635 contigs and 4,053 singletons, generating a total of 4,688 unique sequences. About 46% (2,174) of the unique ESTs had significant hits (E-value

Toll-like receptors (TLRs) are important components of innate immunity. They were found to recognise specific structures on pathogens termed pathogen-associated molecular patterns (PAMPs) and utilise conserved signaling pathways to activate pro-inflammatory cytokines and type-1 interferons. In spite of much understanding gained from the mammalian systems, many fish TLRs are unknown. Recent studies in Japanese flounder as well as in zebrafish suggested that the ligand binding and activation of inflammatory responses in fish may be different from and more complex than those found in mammals. In channel catfish, the major aquaculture species in the United States, only partial sequences of TLR3 and TLR5 were reported. As a part of efforts to characterise the innate immune components in channel catfish, here we cloned and sequenced both the cDNA and the gene for TLR2, a receptor believed mostly responsible for recognition of lipopeptides on the surface of most Gram-positive bacteria. However, expression analysis after infection with a Gram-negative bacterium, Edwardsiella ictaluri indicated that TLR2 was modestly down-regulated in the head kidney tissue of blue catfish, and with a similar pattern in the head kidney of channel catfish though the down-regulation in channel catfish was not statistically significant. In the spleen, an insignificant down-regulation was initially observed early after infection, with an increase of TLR expression later after infection. These results suggest the involvement of TLR2 in the responses after the bacterial infection. As LPS is believed to be the major PAMP for Gram-negative bacteria, additional research is warranted to determine the functions and mechanisms of TLR2 in infections of Gram-negative bacteria.