Transcriptome analysis of Pacific white shrimp (Liptopenaeus vannamei) after exposure to recombinant Vibrio parahaemolyticus PirA and PirB proteins.

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.

Multi-tissue RNAseq reveals genetic and temporal differences in acute response to viral (IHNV) infection among three selected lines of rainbow trout with varying resistance.

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.

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

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

Vertebrate mucus stimulates biofilm development and upregulates iron acquisition genes in Flavobacterium columnare.

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.

Expression of Antisense Long Noncoding RNAs as Potential Regulators in Rainbow Trout with Different Tolerance to Plant-Based Diets.

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.

Comparison of Ribosomal RNA Removal Methods for Transcriptome Sequencing Workflows in Teleost Fish.

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.

Mapping genes to chicken microchromosome 16 and discovery of olfactory and scavenger receptor genes near the major histocompatibility complex.

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.

Novel lytic bacteriophage AhFM11 as an effective therapy against hypervirulent Aeromonas hydrophila.

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.

Impacts of Salmonella enteritidis infection on liver transcriptome in broilers.

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.

Whole-Genome Sequencing and Annotation of Aeromonas veronii Isolates from Channel Catfish.

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.

Transcriptome analysis and immune gene expression of channel catfish (Ictalurus punctatus) fed diets with inclusion of frass from black soldier fly larvae.

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.

Hepatic transcriptome analyses of juvenile white bass (Morone chrysops) when fed diets where fish meal is partially or totally replaced by alternative protein sources.

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.

Molecular responses of ceruloplasmin to Edwardsiella ictaluri infection and iron overload in channel catfish (Ictalurus punctatus).

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.

Alternative splicing in teleost fish genomes: same-species and cross-species analysis and comparisons.

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.

Construction of genetic linkage maps and comparative genome analysis of catfish using gene-associated markers.

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.

Structure and expression of transferrin gene of channel catfish, Ictalurus punctatus.

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.

Proteome analysis of virulent Aeromonas hydrophila reveals the upregulation of iron acquisition systems in the presence of a xenosiderophore.

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.

Comparative analysis of catfish BAC end sequences with the zebrafish genome.

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