- Research article
- Open Access
Transcriptome analysis of the effect of Vibrio alginolyticus infection on the innate immunity-related complement pathway in Epinephelus coioides
BMC Genomics volume 15, Article number: 1102 (2014)
Orange-spotted grouper (Epinephelus coioides) with protogynous hermaphroditic features are one of the most economically important aquaculture species in Taiwan. However, larvae stage grouper are susceptible to infection by the bacterial pathogen Vibrio alginolyticus. To better understand the molecular mechanisms of the immune response to V. alginolyticus in Epinephelus coioides larvae, we used high-throughput deep sequencing technology to study the effect of infection on gene expression.
A total of 114,851,002 reads were assembled, consisting of 9,687,355,560 nucleotides; these were further assembled into 209,082 contigs with a mean length of 372 bp. Gene ontology (GO) analysis of the transcriptome revealed 12 cellular component subcategories, 16 molecular function subcategories, and 42 biological process subcategories (P value <0.05). A total of 32664 Epinephelus coioides genes were mapped to the Kyoto Encyclopedia of Genes and Genomes (KEGG); 1504 differentially expressed genes (DEGs) were subsequently identified, in 12 categories (P value <0.05). Vibrio infection affected the expression of genes involved in complementation, coagulation cascades, pathogen (Staphylococcus aureus) infection, phagosome activity, antigen processing, and the antigen presentation pathway.
We conclude that the complement pathway of innate immunity and the hepicidin antimicrobial peptide may play important roles in the defense of Epinephelus coioides larvae against V. alginolyticus, and the immune response may activate at 4 h after bacterial infection. These results implicate the complement pathway signal pathway in immunity during V. alginolyticus infection at early developmental stages, enhancing our understanding of the mechanisms underlying the immune response to Vibrio infection in Epinephelus coioides.
Groupers are an economically important aquaculture species in Southeast Asian countries, with a high market demand in several locales, including Hong Kong, Taiwan, China, Mexico, Japan, and the USA. However, intensive culture of grouper can result in outbreaks of infectious disease, caused by viral pathogens, such as nodaviruses and iridoviruses, or bacterial pathogens, such as Vibrio carchariae and V. alginolyticus. Grouper zygotes develop by 24 h post-fertilization, and their yolk sacs disappear by 72 h. By 10 days, grouper larvae begin to expand their dorsal and ventral fin rays to form an inverted triangle. At this stage, grouper larvae are easily infected by pathogens. After one month, the long fin reaches its full size, and body shape begins to resemble that of the mature grouper . The orange-spotted grouper (Epinephelus coioides) possesses protogynous hermaphroditic features, and is easier to culture than certain other grouper species, such as Epinephelus lanceolatus. Furthermore, Epinephelus coioides is an excellent source of nutrients; consequently, this species is widely cultured in Taiwan. An earlier study reported that expression of C3 mRNA in Epinephelus coioides is influenced by pH and temperature stress, and may play an important role in antioxidation mechanisms .
Vibrio alginolyticus is a Gram-negative bacterial species with a straight rod shape, and is positive for oxidase and catalase. This bacterial species can cause gastroenteritis, with swelling of the intestine. Furthermore, it is a potential pathogen of marine fish and shrimp . V. alginolyticus strain S3y has been isolated from grouper (Epinephelus malabaricus) larvae with vibriosis in Taiwan; this strain is a pathogen of particular concern, as it causes enormous economic losses in the aquaculture industry . Next-generation high-throughput DNA sequencing techniques, including that provided by the Illumina Genome Analyzer, provide high speed and throughput (gigabase level). Such techniques can be used to identify and quantify rare transcripts without prior knowledge of gene sequence, and provide information regarding alternative splicing and sequence variation in identified genes; as a result, high-throughput sequencing is more effective at detecting genes than microarrays [6, 7].
In this study, we used high throughput sequencing to identify differentially expressed genes (DEGs) between normal grouper larvae and larvae with vibriosis. The DEGs were classified based on their Gene Ontology (GO) categories and the Kyoto Encyclopedia of Genes and Genomes (KEGG). Furthermore, comparative RT-PCR was used to confirm the observed effects on genes involved in likely pathways affected by infection, to elucidate the molecular mechanisms underlying vibriosis in grouper larvae.
De novo sequencing and read assembly of the Epinephelus coioides transcriptome following V. alginolyticusinfection
The assembled transcriptome consisted of a total of 114,851,002 reads of 9,687,355,560 nucleotides; these were further assembled into 209,082 contigs with a mean length of 372 bp. The total contig length (nt) was 77,845,532, and 116,678 unigenes were identified, with an average length of 685 bp; the total unigene length (nt) was 79,966,605. Of these unigenes, 69334 unigenes were present in all databases (Table 1). We identified 3977 DEGs (FDR ≤ 0.001 and |log2Ratio| ≥ 1) between the control and Vibrio challenge group. These include 1104 up-regulated unigenes, and 2873 down-regulated unigenes (>2 fold-change in value). A linear ratio was observed between the RPKM of the V. alginolyticus-infected group and the TSB (tryptic soy broth)-injected group after 24 h (10,2857 unigenes) (Figure 1A). After selection of genes with FDR ≤ 0.001 and |log2Ratio| ≥ 1 (3,977 unigenes), we observed that gene expression was considerably higher in the TSB group than in the V. alginolyticus group (Figure 1B).
Identification of differentially expressed genes (DEGs) via GO and KEGG analysis
Gene ontology (GO) analysis of the 3,977 unigenes was performed using open source clustering software for Annotation , and cluster analysis was performed using cluster software and Java treeview software. GO analysis of the transcriptome revealed 12 cellular component subcategories, 16 molecular function subcategories, and 42 biological process subcategories (P value <0.05) (Additional file 1: Table S2). The largest subcategory in the molecular function group was ‘hydrolase activity’, which was represented by 24.9% of the clustered genes. In the cellular component and biological process categories, ‘cytoskeleton’ and ‘small molecule metabolic process’ were the most abundant GO terms, making up 15.3% and 19.5% of each cluster, respectively (Table 2).
KEGG mapping identified a total of 32,664 genes, including 1,504 DEG in 12 categories (P value <0.05) (Additional file 1: Table S3). The complement and coagulation cascade signaling pathways and Staphylococcus aureus infection signaling pathway were significantly affected by infection; expression of 60 genes related to the complement and coagulation cascades and 44 genes related to the Staphylococcus aureus infection signaling pathway was altered. Other affected immune-related pathways included the phagosome signaling pathway and the antigen processing and presentation signaling pathway, which included 55 and 20 genes, respectively (Table 3).
Bacterial numbers in Epinephelus coioides larvae infected with V. alginolyticus
To study V. alginolyticus infection in grouper larvae, we calculated the colony-forming units (CFUs) in whole fish over time (Figure 2). CFUs were significantly greater in infected larvae than in control fish between 2 and 7 h post-injection. By 8 h, no significant difference was detected, and by 24 h, CFUs had returned to baseline levels.
Analysis of immune-related signal transduction pathways in infected fish
GO (Gene Ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) analyses of immune factors yielded similar findings; for example, both analyses revealed that complement activation was affected by infection. However, KEGG analysis uncovered additional immune responses and clearly disrupted pathways. We thus subsequently focus on the findings of KEGG analysis. Based on KEGG analysis, we selected the following pathways for analysis: complement and coagulation cascades with a p-value of 3.13 × 10−19 (Additional file 2: Figure S1), the Staphylococcus aureus infection pathway with a p-value of 1 × 10−17 (Additional file 3: Figure S2), the phagosome pathway with a p-value of 1.27 × 10−5 (Additional file 4: Figure S3), and the antigen processing and presentation pathway with a p-value of 2.35 × 10−5 (Additional file 5: Figure S4). These pathways are part of the teleost immune response. We combined the Staphylococcus aureus infection pathway with the complement and coagulation cascades (henceforth referred to as the complement-related pathway), and the phagosome pathway with the antigen processing and presentation pathway (henceforth referred to as the phagocytosis-related pathway), designed primers against unigenes involved in these processes based on the transcriptome sequences, and examined RNA expression by real-time qPCR, as described below.
Analysis of gene expression in the complement-related pathway
As shown in Figures 3 and 4 and Additional file 1: Table S4, we examined the effect of infection on the expression of genes involved in the complement-related pathway by qPCR; the gene names given in the figures are the abbreviations used in KEGG. The complement factor B-like (BF) gene was significantly up-regulated between 3 h and 12 h as compared to the control (Figure 3B), while the C2r subcomponent-like (C2r) gene was significantly down-regulated at 2 h and 3 h (Figure 3C). The complement C1q-like protein 2 (C1q) (Figure 3D) and complement C1r subcomponent-like (C1r) (Figure 3E) genes were significantly elevated in infected larvae at 8 h and 4 h, respectively. Both the haptoglobin-like (C1s) (Figure 3F) and complement component C3 (C3) (Figure 3G) genes were up-regulated between 4 h and 8 h in infected fish. The complement C4-like (C4) gene exhibited an erratic pattern, with up-regulation observed in infected fish at 4 h, 5 h, 7 h, 8 h, and 24 h (Figure 3H). The beta-2-glycoprotein 1 precursor (C4BP) gene was also up-regulated at various time points (5 h, 8 h, 12 h, and 24 h) (Figure 3I). Expression levels of the complement C5 (C5) (Figure 3J) and complement component C6-like (C6) (Figure 3K) genes were significantly increased between 8 h and 24 h, and 4 h to 24 h, respectively. Levels of the Spondin-2-like (C6-d) gene were significantly increased at 1 h, 5 h, 6 h, and 7 h (Figure 3L).The complement component C7-like (C7) gene was significantly elevated between 4 h and 16 h (Figure 4A), and the complement component C8 beta (C8) gene was up-regulated at 5 h, 8 h, 12 h, 16 h, and 24 h in infected fish (Figure 4B). Both the complement component C8 alpha chain-like (C8-d) (Figure 4C) and the complement component C9 (C9) (Figure 4D) genes were significantly increased between 4 h and 24 h. On the other hand, the mannan-binding lectin serine protease 1 (MASP1/2) gene was significantly down-regulated between 2 h and 4 h (Figure 4F). The complement regulatory plasma protein (HF) (Figure 4G) and complement factor I-like (IF) (Figure 4H) genes were significantly up-regulated between 5 h and 8 h, and between 2 h and 16 h, respectively. Finally, the Ca2 + −dependent complex C1r/C1s subunit (PLG) gene exhibited both up- and down-regulation, at 1 h and 7 h, respectively (Figure 4I). Expression levels of alpha-2-macroglobulin-like (A2M) (Figure 3A), decay accelerating factor (DAF) (Figure 4E), and minus strand C1 inhibitor (SERPING1) (Figure 4J) were unaffected by infection. Certain genes were affected in a time-dependent manner, i.e., BF, C1s, C3, C5, C6, C7, C8-d, C9, HF, and IF.
Analysis of gene expression in the phagocytosis-related pathway
The effect of infection on the expression of genes in the phagocytosis-related pathway was more erratic than that on genes of the complement-related pathway (Additional file 1: Table S5). The protein HEG-like (αVβ5) (Figure 5C) and Calnexin-like (Calnexin) (Figure 5D) genes were significantly up-regulated by infection at 5 h and 7 h, respectively. Cathepsin L precursor (Cathepsin) gene was up-regulated between 5 h and 7 h (Figure 5E), while the lactose-binding lectin l-2-like (Collectins) gene was up-regulated at 6 h only (Figure 4F). Expression of the type II antifreeze protein I (DCSIGN) gene was elevated at 24 h and 48 h (Figure 5G), that of the early endosome antigen 1 (EEA1) gene at 16 h and 24 h (Figure 5J), and that of the CDH1-D (F-actin) gene at 12 h (Figure 5K).The heat shock 70 kDa protein 14-like (HSP70) gene was significantly up-regulated by infection at 6 h (Figure 6A), and the invariant chain-like protein (Ii-d) gene at 3 h (Figure 6B). The lysosomal membrane glycoprotein 2 precursor (LAMP) gene was up-regulated at both 7 h and 8 h (Figure 6C). The beta-centractin-like (MHCI) gene was also up-regulated at 3 h (Figure 6D), while the eosinophil peroxidase-like (MPO-d) gene was down-regulated at 16 h (Figure 6F). Infection increased expression of the macrophage mannose receptor 1-like (MR) gene between 5 h and 8 h (Figure 6G), the tapasin-like (TAPBP) gene at 5 h and 6 h (Figure 6J), and the transferrin receptor 1a (TfR) gene at 4 h (Figure 6K). Finally, the tubulin, beta 5 (TUBB) gene was significantly down-regulated at 8 h and 16 h (Figure 6L). Infection did not affect the expression levels of Integrin, alpha V (αVβ3) (Figure 5A), Integrin beta-3-like (αVβ3-2) (Figure 5B), Nattectin (DCSIGN2) (Figure 5H), Cytoplasmic dynein 1 heavy chain 1-like isoform 2 (Dynein) (Figure 5I), Beta-centractin-like (F-actin-d) (Figure 5L), COP9 signalosome complex subunit 7a-like (MHCII) (Figure 6E), L-rhamnose-binding lectin CSL2-like (MR-d) (Figure 6H), or Nuclear transcription factor Y subunit alpha-like (NFY) genes (Figure 6I).
Analysis of antimicrobial peptide gene expression
The innate immune response includes antimicrobial peptides, which can damage bacterial membranes. We identified three antimicrobial peptide genes in the transcriptome library, and analyzed their expression in infected larvae over time (Additional file 1: Table S6). Interestingly, expression of the antimicrobial peptide Epinecidin-1, which our laboratory previously isolated from Epinephelus coioides, was unaffected (data not shown). The sequence of epinecidin-1 is similar to that of piscidin (54% identity), and these two peptides were found to have similar effects. However, here we focused on the three antimicrobial peptide genes in our transcriptome database.
The hepcidin-like antimicrobial peptide precursor (Hepcidin) gene was significantly up-regulated by infection from 4 h to 8 h (Figure 7A), while the liver-expressed antimicrobial peptide 2-like (LAP2) was unaffected (Figure 7B). Conversely, the Piscidin-like antimicrobial peptide precursor (Piscidin) gene was significantly down-regulated at 8 h (Figure 7C). These findings suggest that Epinephelus coioides is dependent on innate immunity to defend against V. alginolyticus infection.
Analysis of the complement- and phagocytosis-related pathways by KEGG
As described above, we combined the complement and coagulation cascades with the Staphylococcus aureus infection pathway to form the complement-related pathway (Additional file 6: Figure S5A), and the phagosome pathway with the antigen processing and presentation pathway to form the phagocytosis-related pathway (Additional file 6: Figure S5B). The complement-related pathway transduces signals to classical pathways, the alternative pathway, and lectin pathway. All complement-related pathways involve cleavage of C3 to C3b and then to C5, in order to form a membrane attack complex (MAC) for bacterial lysis. The phagocytosis pathway uses nitric oxide synthase to produce NO via the phagosome, which then digests the bacterium and presents the antigen fragment, thereby activating MHCI and MHCII to stimulate the adaptive immune system.
Here, we describe the use of NGS technology to uncover the response of the transcriptome of Epinephelus coioides larvae to infection by V. alginolyticus. Few studies have focused on infection of the grouper larvae stage; however, because the long dorsal fin of the larvae stage is yet to completely disappear, it makes fish of this developmental stage prone to getting trapped in nets. Furthermore, they are susceptible to dying for many reasons, including changes in water temperature, aeration rate, salinity, and illumination . We constructed a transcriptome library from Epinephelus coioides larvae, as transcriptome profiling is a powerful method for evaluating the relative importance of gene products in a given tissue , and it enabled us to determine the effects of infection on gene expression at the larval stage. We deduced the immune-related signal transduction pathway from KEGG enrichment analysis. This pathway was predicted to consist of the complement and coagulation cascades, Staphylococcus aureus infection pathway, phagosome pathway, and antigen processing and presentation pathway. The complement system is an ancient mechanism, found in both protostome and invertebrate deuterostome species [12–14]. Like other higher vertebrates, teleost fish contain three complement pathways. Earlier studies have demonstrated that the alternative and classical pathways have a significant effect in teleost fish . However, very little is known about the molecules involved in the lectin pathway in fish . The complement components of fish are different to those of mammals, and some consist of multiple isoforms [17–19]. The complement pathway is known to be one of the key mechanisms for bacterial clearance in teleost fish [20, 21]; the alternative pathway can be activated by the lipopolysaccharides (LPS) of Gram-negative bacteria, enabling lysis of the bacterial cell . Here, we observed that infection had time-dependent effects on several genes related to the complement pathway; however, of the genes of the phagocytosis-related pathway, only the MR and cathepsin genes were affected in a time-dependent manner. The complement factor B-like (BF) gene cleaves C3 and acts as a convertase in rainbow trout . Although the C1r/C1s/MASP-like genes of grouper have not been functionally characterized, previous findings suggest that the C1s-like molecule may cleave C4 to C4a and C4b fragments in rainbow trout [15, 24].
In teleost fish, C4 plays an important role in activation of the classical pathways . However, our current results indicate that C4 gene expression is not affected in a time-dependent manner by V. alginolyticus infection in grouper larvae. Mammalian C3 is encoded by a single gene, but almost all teleost fish studied produce multiple forms of C3 encoded by different genes . In trout, carp, and seabream, these C3 isoforms exhibit different binding efficiencies to several active complements; as such, these isoforms may perform separate roles in the destruction of microbes and innate recognition . In Epinephelus coioides, C3 may be inducible and involved in stress responses . C5 is a part of the membrane attack complex (MAC), which cleaves C5 into C5a and C5b fragments . The C5 gene has been partially cloned and purified from trout and seabream [27, 28], and MAC has been shown to consist of C5b, C6, C7, and the beta chains of C8 and C9 in these species [29, 30]. It has been hypothesized that complement cofactor protein (SBP1) regulates both C4b binding protein and factor H  in barred sand bass (Paralabrax nebulifer), and two factor I isotypes have been identified in carp . Based on our results, we hypothesize that V. alginolyticus infection of orange-spotted grouper initially results in the activation of genes such as BF and IF, which mediate C3 production between 2 and 3 h post-infection. C3 is then cleaved by C3 convertase through the alternative or classical pathway at 4 h (via C1s or C4). Formation of MAC occurs between 4 and 5 h (although C5 expression is not significantly up-regulated by infection until 8 h, a non-significant increase can be observed at 4 h). Finally, MAC may clear V. alginolyticus before 48 h.
The phagosome and antigen processing system include phagocytic cells (granulocytes, monocytes, and macrophages), non-specific cytotoxic cells, and dendritic cells [33, 34]. From our gene expression results, it is difficult to determine how infection affects the phagocytosis-related pathway (Figure 8). However, we did observe an increase in macrophage mannose receptor (MR) expression between 5 and 8 h, and thus conclude it may be involved in the complement lectin pathway and initiated by the binding of a lectin, such as C-type lectin; it may then activate the complement pathway upon binding of a collectin to a microbial surface . The MR may also be involved in the phagocytosis of yeast cells by head-kidney leucocytes in seabream (Sparus aurata L.) . Complement component C3 has been identified in fertilized cod eggs , and phagocytic activity has been detected in zebrafish embryo and 2-day-post-fertilization carp embryo [38, 39].
Antimicrobial peptides play important roles in the innate immune response to bacterial infection. Antimicrobial peptides range in size from 6 amino acid residues for anionic peptides, to as many as 59 amino acid residues; larger proteins possess several features of secondary structure, including α -helices, relaxed coils, and antiparallel β-sheet structures. Such features are hydrophobic, which enables water-soluble antimicrobial peptides to pass through the membrane lipid bilayer. Antimicrobial peptides kill bacteria via one of three mechanisms, known as the barrel-stave model, carpet model, and toroidal model. The barrel-stave model involves aggregation of peptides at the membrane bilayer; the hydrophobic peptide regions align to form a lipid hole like a barrel, disrupting the osmotic balance of the bacterial inner membrane. In the carpet model, the peptides are oriented parallel to the surface of the lipid bilayer, like a carpet. In the toroidal model, the peptides aggregate, and form a pore in the lipid monolayers, and the inserted peptides and the lipid head groups induce a water core line .
Hepcidin was first identified as a protein involved in innate immunity and iron regulation in the human liver. Teleost fish hepcidin has been previously demonstrated to be involved in both iron regulation and immunity. Seabream hepcidin is abundant in liver, skin, head-kidney, and peritoneal exudate leucocytes, and flounder hepcidin-like is distributed in liver, esophagus, and cardiac stomach. Gene expression is up-regulated by poly I:C, iron dextran, bacteria, or LPS [41, 42]. In a recent study, a four-cysteine hepcidin isoform gene, EC-hepcidin3, was cloned from Epinephelus coioides, and was reported to be effective against Staphylococcus aureus and Pseudomonas stutzeri. Human hepcidin is induced by IL-6 , and we previously observed that transgenic zebrafish expressing tilapia hepcidin 2–3 had higher transcript levels of IL-10, IL-26, TLR4a, and TNF-α as compared with wild-type zebrafish . The gene encoding the hepcidin-like antimicrobial peptide precursor was strongly induced by infection at 7 and 8 h, and signs of increased expression suggest that grouper larvae may be dependent on hepcidin function from 4 to 8 h.
In conclusion, the present study suggests that the response of Epinephelus coioides larvae to V. alginolyticus infection is dependent on the complement pathway and antimicrobial peptides. Hepcidin, which plays important roles during the larvae stages of grouper, may also be involved in the defense against bacterial infection (Figure 8). These results may be useful to research on fish, as they suggest that the complement pathway and antimicrobial peptides may be beneficial in terms of enhancing grouper anti-bacterial defenses.
Fish and bacteria
Epinephelus coioides larvae were purchased from the Institute of Biotechnology, National Cheng Kung University Core Facility. The Epinephelus coioides larvae were kept in a 2 tonne FRP tank prior to bacterial infection. Vibrio alginolyticus were cultured as previously described . The animal experimental and ethics protocol (12-12-447) was approved by the Academia Sinica Institutional Animal Care and Use Committee (IACUC) of the Institute of Cellular and Organismic Biology, Academia Sinica, Taiwan.
Bacterial infection and colony-forming unit (CFU) counts
At 30-days-old, the whole bodies of Epinephelus coioides larvae with an average body length of about 1.6 ± 0.2 cm were injected with 20 μl (1.3 × 106 colony-forming units (CFU)/ml) of V. alginolyticus in Tryptic soy broth (NaCl 1.5%). Vibrio were detected by spreading the culture onto thiosulfate–citrate–bile salt–sucrose agar (TCBS) plates, and incubating the plates at 28°C for 16 h; V. alginolyticus forms yellow colonies . Control grouper were injected with TSB (NaCl 1.5%) without V. alginolyticus. Fish were sacrificed at 1, 2, 3, 4, 5, 6, 7, 8, 12, 24, and 48 h after infection (five fish were sacrificed at each time point for each group). Whole fish were homogenized using Lyser II (Qiagen, USA) solution, and CFU were determined as previously described .
Epinephelus coioides larval RNA preparation for next generation sequencing
Epinephelus coioides larvae (30-days-old) were injected with V. alginolyticus as described in the preceding section, and fish were sacrificed at 24 h after infection. Ten wild-type Epinephelus coioides larvae were used to build an EST library for transcriptome analysis. RNA was extracted as previously described . RNA concentrations were quantified using a Nano-Drop spectrophotometer (Thermo, USA), and quality was determined with an RNA gel, as described in Molecular Cloning. RNA samples for real-time qPCR were extracted from fish sacrificed at 1, 2, 3, 4, 5, 6, 7, 8, 12, 24, and 48 h after infection (n = 5); three independent trials were performed.
Oligo (dT) magnetic beads were used to enrich mRNA, which were then broken into short fragments (about 200 bp) in fragmentation buffer; these fragments were then reverse transcribed into first strand cDNA using random hexamer primers. The appropriate buffer, dNTPs, RNase H, and DNA polymerase I were added to synthesize second-strand cDNA. Double-stranded cDNA was purified with the QiaQuick PCR extraction kit, and washed with EB buffer. Sequencing adaptors were ligated to the fragments. The required fragments were purified by agarose gel electrophoresis and copied by PCR amplification. The library products were prepared for sequencing analysis using an Illumina HiSeq™ 2000. Raw data were saved as fastq files. The following were removed: reads with adaptors, reads in which over 10% bases were unknown, and low quality reads (i.e., the percentage of bases of quality value ≤5 exceeds 50% in the read). Clean reads were mapped to reference sequences using SOAP aligner/soap2 . The randomness of RNA fragmentation was used to construct the library, and the numbers of reads mapped to the reference sequence were calculated. The RPKM method (Reads Per kb per Million reads) was used to calculate gene expression level , and differentially expressed genes (DEGs) were subsequently screened for as previously described . DEGs were subjected to GO function and KEGG pathway analysis, as described below. All transcriptome databases used in our study can be downloaded from the foot of GSE63148 used by NCBI Gene Expression Omnibus (http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE63148); this includes bacteria-infected group gene rpkm (GSE63148_B4_1.Gene.rpkm.txt.gz), control group gene rpkm (GSE63148_MRS17.Gene.rpkm.txt.gz), transcriptome EST library (GSE63148_MRS-16-Unigene.fa.gz), gene differential expression (GSE63148_MRS17-VS-B4_1.GeneDiffExp.txt.gz), and gene annotation (GSE63148_annotation.txt.gz).
Gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) analysis
GO enrichment analysis was performed by collating all the GO terms that were significantly enriched in the identified DEG, and then filtering the DEGs based on these biological functions. First, all DEGs were mapped to GO terms in the database (http://www.geneontology.org/), and then gene numbers were calculated for every term using the hypergeometric test in order to obtain significantly enriched GO terms for DEGs; these were compared to the genomic background, as described in a previous study . Pathway enrichment analysis was performed using KEGG (http://www.genome.jp/kegg/), which is a public pathway-related database. Such analysis was used to identify significant enrichment of genes involved in metabolic or signal transduction pathways. DEGs were compared with the genomic background, and the formula was calculated as for GO analysis.
Analysis of gene expression
Harikrishnan R, Balasundaram C, Heo M-S: Molecular studies, disease status and prophylactic measures in grouper aquaculture: Economic importance, diseases and immunology. Aquaculture. 2010, 309 (1): 1-14.
Kohno H, Diani S, Supriatna A: Morphological development of larval and juvenile grouper, Epinephelus fuscoguttatus. Jap J Ichthyol. 1993, 40: 307-316.
Qi Z-H, Liu Y-F, Wang W-N, Wu X, Xin Y, Lu Y-F, Wang A-L: Molecular characterization and functional analysis of a complement C3 molecule in the orange-spotted grouper (Epinephelus coioides). Fish Shellfish Immunol. 2011, 31 (6): 1284-1290. 10.1016/j.fsi.2011.09.018.
Cheng A-C, Tu C-W, Chen Y-Y, Nan F-H, Chen J-C: The immunostimulatory effects of sodium alginate and iota-carrageenan on orange-spotted grouper Epinephelus coicoides and its resistance against Vibrio alginolyticus. Fish Shellfish Immunol. 2007, 22 (3): 197-205. 10.1016/j.fsi.2006.04.009.
Lee K-K: Pathogenesis studies on Vibrio alginolyticus in the grouper, Epinephelus malabaricus, Bloch et Schneider. Microb Pathog. 1995, 19 (1): 39-48. 10.1016/S0882-4010(85)90000-2.
Metzker ML: Sequencing technologies—the next generation. Nat Rev Genet. 2010, 11 (1): 31-46. 10.1038/nrg2626.
Ansorge WJ: Next-generation DNA sequencing techniques. New Biotechnol. 2009, 25 (4): 195-203. 10.1016/j.nbt.2008.12.009.
de Hoon MJ, Imoto S, Nolan J, Miyano S: Open source clustering software. Bioinformatics. 2004, 20 (9): 1453-1454. 10.1093/bioinformatics/bth078.
Pan C-Y, Chen J-Y, Cheng Y-SE, Chen C-Y, Ni I-H, Sheen J-F, Pan Y-L, Kuo C-M: Gene expression and localization of the epinecidin-1 antimicrobial peptide in the grouper (Epinephelus coioides), and its role in protecting fish against pathogenic infection. DNA Cell Biol. 2007, 26 (6): 403-413. 10.1089/dna.2006.0564.
Sakakura Y, Shiotani S, Chuda H, Hagiwara A: Improvement of the survival in the seven‒band grouper Epinephelus septemfasciatus larvae by optimizing aeration and water inlet in the mass‒scale rearing tank. Fisheries Sci. 2006, 72 (5): 939-947. 10.1111/j.1444-2906.2006.01241.x.
Mu Y, Ding F, Cui P, Ao J, Hu S, Chen X: Transcriptome and expression profiling analysis revealed changes of multiple signaling pathways involved in immunity in the large yellow croaker during Aeromonas hydrophila infection. BMC Genomics. 2010, 11 (1): 506-10.1186/1471-2164-11-506.
Zhu Y, Thangamani S, Ho B, Ding JL: The ancient origin of the complement system. EMBO J. 2005, 24 (2): 382-394. 10.1038/sj.emboj.7600533.
Al-Sharif WZ, Sunyer JO, Lambris JD, Smith LC: Sea urchin coelomocytes specifically express a homologue of the complement component C3. The J Immunol. 1998, 160 (6): 2983-2997.
Nonaka M: Evolution of the complement system. Subcell Biochem. 2014, 80: 31-43. 10.1007/978-94-017-8881-6_3. doi:10.1007/978-94-017-8881-6_3
Boshra H, Gelman AE, Sunyer JO: Structural and functional characterization of complement C4 and C1s-like molecules in teleost fish: insights into the evolution of classical and alternative pathways. J Immunol. 2004, 173 (1): 349-359. 10.4049/jimmunol.173.1.349.
Boshra H, Li J, Sunyer J: Recent advances on the complement system of teleost fish. Fish Shellfish Immunol. 2006, 20 (2): 239-262. 10.1016/j.fsi.2005.04.004.
Sunyer JO, Zarkadis IK, Sahu A, Lambris JD: Multiple forms of complement C3 in trout that differ in binding to complement activators. Proc Natl Acad Sci U S A. 1996, 93 (16): 8546-8551. 10.1073/pnas.93.16.8546.
Sunyer J, Tort L, Lambris JD: Structural C3 diversity in fish: characterization of five forms of C3 in the diploid fish Sparus aurata. J Immunol. 1997, 158 (6): 2813-2821.
Nakao M, Mutsuro J, Obo R, Fujiki K, Nonaka M, Yano T: Molecular cloning and protein analysis of divergent forms of the complement component C3 from a bony fish, the common carp (Cyprinus carpio): presence of variants lacking the catalytic histidine. Eur J Immunol. 2000, 30 (3): 858-866. 10.1002/1521-4141(200003)30:3<858::AID-IMMU858>3.0.CO;2-M.
Holland MCH, Lambris JD: The complement system in teleosts. Fish Shellfish Immunol. 2002, 12 (5): 399-420. 10.1006/fsim.2001.0408.
Ellis A: Innate host defense mechanisms of fish against viruses and bacteria. Dev Comp Immunol. 2001, 25 (8): 827-839.
Oriol Sunyer J, Tort L: Natural hemolytic and bactericidal activities of sea bream Sparus aurata serum are effected by the alternative complement pathway. Vet Immunol Immunopathol. 1995, 45 (3): 333-345.
Sunyer JO, Zarkadis I, Sarrias MR, Hansen JD, Lambris JD: Cloning, structure, and function of two rainbow trout Bf molecules. J Immunol. 1998, 161 (8): 4106-4114.
Wang T, Secombes CJ: Complete sequencing and expression of three complement components, C1r, C4 and C1 inhibitor, of the classical activation pathway of the complement system in rainbow trout Oncorhynchus mykiss. Immunogenetics. 2003, 55 (9): 615-628. 10.1007/s00251-003-0622-5.
Oriol Sunyer J, Zarkadis IK, Lambris JD: Complement diversity: a mechanism for generating immune diversity?. Immunol Today. 1998, 19 (11): 519-523. 10.1016/S0167-5699(98)01341-3.
Müller-Eberhard HJ: The killer molecule of complement. Br J Dermatol. 1985, 85 (1 Suppl): 47s-52s.
Franchini S, Zarkadis IK, Sfyroera G, Sahu A, Moore WT, Mastellos D, LaPatra SE, Lambris JD: Cloning and purification of the rainbow trout fifth component of complement (C5). Dev Comp Immunol. 2001, 25 (5): 419-430.
Sunyer J, TORT L, Lambris J: Diversity of the third form of complement, C3, in fish: functional characterization of five forms of C3 in the diploid fish Sparus aurata. Biochem J. 1997, 326: 877-881.
Nonaka M, Yamaguchi N, Natsuume-Sakai S, Takahashi M: The complement system of rainbow trout (Salmo gairdneri). I. Identification of the serum lytic system homologous to mammalian complement. J Immunol. 1981, 126 (4): 1489-1494.
Nakao M, Uemura T, Yano T: Terminal components of carp complement constituting a membrane attack complex. Mol Immunol. 1996, 33 (11): 933-937.
Kemper C, Zipfel PF, Gigli I: The complement cofactor protein (SBP1) from the barred sand bass (Paralabrax nebulifer) mediates overlapping regulatory activities of both human C4b binding protein and factor H. J Biol Chem. 1998, 273 (31): 19398-19404. 10.1074/jbc.273.31.19398.
Nakao M, Hisamatsu S, Nakahara M, Kato Y, Smith SL, Yano T: Molecular cloning of the complement regulatory factor I isotypes from the common carp (Cyprinus carpio). Immunogenetics. 2003, 54 (11): 801-806.
Neumann NF, Stafford JL, Barreda D, Ainsworth AJ, Belosevic M: Antimicrobial mechanisms of fish phagocytes and their role in host defense. Dev Comp Immunol. 2001, 25 (8): 807-825.
GANASSIN R, BOLS N: Development of long-term rainbow trout spleen cultures that are haemopoietic and produce dendritic cells. Fish Shellfish Immunol. 1996, 6 (1): 17-34. 10.1006/fsim.1996.0003.
Vasta GR, Nita-Lazar M, Giomarelli B, Ahmed H, Du S, Cammarata M, Parrinello N, Bianchet MA, Amzel LM: Structural and functional diversity of the lectin repertoire in teleost fish: relevance to innate and adaptive immunity. Dev Comp Immunol. 2011, 35 (12): 1388-1399. 10.1016/j.dci.2011.08.011.
Rodrı́guez A, Esteban M, Meseguer J: A mannose-receptor is possibly involved in the phagocytosis of Saccharomyces cerevisiae by seabream (Sparus aurata L.) leucocytes. Fish Shellfish Immunol. 2003, 14 (5): 375-388. 10.1006/fsim.2002.0446.
Magnadóttir B, Lange S, Steinarsson A, Gudmundsdóttir S: The ontogenic development of innate immune parameters of cod (Gadus morhua L.). Comp Biochem Physiol B Biochem Mol Biol. 2004, 139 (2): 217-224. 10.1016/j.cbpc.2004.07.009.
Herbomel P, Thisse B, Thisse C: Ontogeny and behaviour of early macrophages in the zebrafish embryo. Development. 1999, 126 (17): 3735-3745.
Romano N, Picchietti S, Taverne-Thiele J, Taverne N, Abelli L, Mastrolia L, Verburg-van Kemenade B, Rombout J: Distribution of macrophages during fish development: an immunohistochemical study in carp (Cyprinus carpio, L.). Anat Embryol (Berl). 1998, 198 (1): 31-41. 10.1007/s004290050162.
Brogden KA: Antimicrobial peptides: pore formers or metabolic inhibitors in bacteria?. Nat Rev Microbiol. 2005, 3 (3): 238-250. 10.1038/nrmicro1098.
Cuesta A, Meseguer J, Esteban MÁ: The antimicrobial peptide hepcidin exerts an important role in the innate immunity against bacteria in the bony fish gilthead seabream. Mol Immunol. 2008, 45 (8): 2333-2342. 10.1016/j.molimm.2007.11.007.
Douglas SE, Gallant JW, Liebscher RS, Dacanay A, Tsoi S: Identification and expression analysis of hepcidin-like antimicrobial peptides in bony fish. Dev Comp Immunol. 2003, 27 (6): 589-601.
Qu H, Chen B, Peng H, Wang K: Molecular cloning, recombinant expression, and antimicrobial activity of EC-hepcidin3, a New four-cysteine hepcidin isoform from Epinephelus coioides. Biosci Biotechnol Biochem. 2013, 77 (1): 103-110. 10.1271/bbb.120600.
Nemeth E, Valore EV, Territo M, Schiller G, Lichtenstein A, Ganz T: Hepcidin, a putative mediator of anemia of inflammation, is a type II acute-phase protein. Blood. 2003, 101 (7): 2461-2463. 10.1182/blood-2002-10-3235.
Hsieh J-C, Pan C-Y, Chen J-Y: Tilapia hepcidin (TH) 2–3 as a transgene in transgenic fish enhances resistance to Vibrio vulnificus infection and causes variations in immune-related genes after infection by different bacterial species. Fish Shellfish Immunol. 2010, 29 (3): 430-439. 10.1016/j.fsi.2010.05.001.
Wang Y-D, Peng K-C, Wu J-L, Chen J-Y: Transgenic expression of salmon delta-5 and delta-6 desaturase in zebrafish muscle inhibits the growth of Vibrio alginolyticus and affects fish immunomodulatory activity. Fish Shellfish Immunol. 2014, 39 (2): 223-230. 10.1016/j.fsi.2014.04.021.
Wang Y-D, Kung C-W, Chi S-C, Chen J-Y: Inactivation of nervous necrosis virus infecting grouper (Epinephelus coioides) by epinecidin-1 and hepcidin 1–5 antimicrobial peptides, and downregulation of Mx2 and Mx3 gene expressions. Fish Shellfish Immunol. 2010, 28 (1): 113-120. 10.1016/j.fsi.2009.10.001.
Sambrook J, Fritsch EF, Maniatis T: Molecular Cloning, Volume 2. 1989, New York: Cold spring harbor laboratory press
Li R, Yu C, Li Y, Lam T-W, Yiu S-M, Kristiansen K, Wang J: SOAP2: an improved ultrafast tool for short read alignment. Bioinformatics. 2009, 25 (15): 1966-1967. 10.1093/bioinformatics/btp336.
Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B: Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008, 5 (7): 621-628. 10.1038/nmeth.1226.
Audic S, Claverie J-M: The significance of digital gene expression profiles. Genome Res. 1997, 7 (10): 986-995.
Zhang B, Schmoyer D, Kirov S, Snoddy J: GOTree machine (GOTM): a web-based platform for interpreting sets of interesting genes using gene ontology hierarchies. BMC bioinformatics. 2004, 5 (1): 16-10.1186/1471-2105-5-16.
This study was partially supported by funding from the National Science Council (102-2313-B-001-003-) and partially by funding from an ITAR Project Application (Transcriptomic and immunological responses to Vibrio infection during different developmental stages in Epinephelus lanceolatus), Academia Sinica, to Dr. Jyh-Yih Chen. Research funding was also received from the Marine Research Station (Jiaushi, Ilan), Institute of Cellular and Organismic Biology, Academia Sinica, Taiwan, to Dr. Jyh-Yih Chen. We thank Dr. Duncan Wright at the editorial office of the Institute of Cellular and Organismic Biology for manuscript editing and revision.
The authors declare that they have no competing interests.
YDW and JYC carried out the V. alginolyticus infection experiments, counted the CFUs, and extracted RNA. YDW and SJH carried out the transcriptome and gene expression experiments. HNC, WLL, and HYG provided experimental materials for the transcriptome and gene expression experiments. YDW and JYC designed the experiments and wrote the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Table S1.: Real-Time PCR primers used in this study. Gene ontology analysis of Epinephelus coioides larvae. Table S3. KEGG pathway enrichment analysis of Epinephelus coioides larvae. Table S4. Summary of the variations in complement-related gene expression between controls and grouper challenged with Vibrio alginolyticus. Table S5. Summary of the variations in phagocytosis-related gene expression between controls and grouper challenged with Vibrio alginolyticus. Table S6. Summary of the variations in antimicrobial peptide gene expression between controls and grouper challenged with Vibrio alginolyticus. (PDF 360 KB)
Additional file 2: Figure S1.: Complement and coagulation cascades signal pathway. Enrichment analysis of DEGs from the KEGG database; red borders indicate up-regulated genes, green borders indicate down-regulated genes, and red/green borders indicate genes that are both up- and down-regulated at different times. (TIFF 233 KB)
Additional file 3: Figure S2.: Staphylococcus aureus infection signal pathway. Enrichment analysis of DEGs from the KEGG database; red borders indicate up-regulated genes, green borders indicate down-regulated genes, and red/green borders indicate genes that are both up- and down-regulated at different times. (TIFF 280 KB)
Additional file 4: Figure S3.: Phagosome signal pathway. Enrichment analysis of DEGs from the KEGG database; red borders indicate up-regulated genes, green borders indicate down-regulated genes, and red/green borders indicate genes that are both up- and down-regulated at different times. (TIFF 259 KB)
Additional file 5: Figure S4.: Antigen processing and presentation signal pathway. Enrichment analysis of DEGs from the KEGG database; red borders indicate up-regulated genes, green borders indicate down-regulated genes, and red/green borders indicate genes that are both up- and down-regulated at different times. (TIFF 177 KB)
Additional file 6: Figure S5.: Hypothetical model for the response of grouper larvae to Vibrio alginolyticus infection, as predicted by KEGG analysis. (A) Complement pathway. (B) Phagocytosis pathway. Red borders indicate up-regulated genes, green borders indicate down-regulated genes, and purple borders indicate genes that are both up- and down-regulated at different times. Red arrows indicate increased RNA expression and green arrows indicate decreased RNA expression. Numbers adjacent to borders are the log2 ratio of significantly affected genes. (TIFF 287 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Wang, YD., Huang, SJ., Chou, HN. et al. Transcriptome analysis of the effect of Vibrio alginolyticus infection on the innate immunity-related complement pathway in Epinephelus coioides. BMC Genomics 15, 1102 (2014). https://doi.org/10.1186/1471-2164-15-1102
- Gene Ontology
- Antimicrobial Peptide
- Complement Pathway
- Mannose Receptor
- Membrane Attack Complex