Molecular characterization of miiuy croaker CC chemokine gene and its expression following Vibrio anguillarum injection

METTL16, an evolutionarily conserved m6A methyltransferase member, inhibits the antiviral immune response of miiuy croaker (Miichthys miiuy)

Methyltransferase like-16 (METTL16) is an m6A RNA methylation transferase that is known to methylate U6 snRNA and pre-mRNA of S-adenosylmethionine synthase but has been poorly studied in fish. In this study, METTL16 was identified in miiuy croaker (Miichthys miiuy). We first performed bioinformatics analysis of the miiuy croaker METTL16 (mmiMETTL16). MmiMETTL16 and other vertebrates METTL16 have a relatively conserved MTD structural domain and gene structure, suggesting that their methylase activity may also be conservative. In healthy miiuy croaker, mmiMETTL16 was commonly expressed in the tested tissues. Expression of mmiMETTL16 in kidney, liver, and spleen tissues was significantly increased after poly(I:C) stimulation. Consistently, mmiMETTL16 was sensitive to poly(I:C) stimulation in miiuy croaker kidney cell (MKC), suggesting that METTL16 might participate in antiviral immunity. For further functional experiments, immunofluorescence of mmiMETTL16 presents in the nucleus in kidney cells. In addition, the overexpression of mmiMETTL16 could significantly increase the overall m6A level of MKC cells, which shows that the function of METTL16 as methyltransferase is conservative in miiuy croaker. Last, mmiMETTL16 can inhibit the expression of TNF-α, IFN-1, Mx1, and ISG15, suggesting that mmiMETTL16 can suppress the immune response caused by viral stimulation. In summary, studies on mmiMETTL16 will contribute to future studies on the role of METTL16 and potential mechanisms of the m6A regulation network in the teleost immune system.

Rho-GDP-dissociation inhibitor-γ negatively regulates NF-κB signaling by promoting the degradation of TAK1 in miiuy croaker (Miichthys miiuy)

Transforming growth factor-beta activated kinase 1 (TAK1) is an adaptor molecular in TLR-mediated NF-κB signaling pathway and plays indispensable roles in innate immunity. As the most typical innate immune pathway, the strict regulation of NF-κB signaling pathway is particularly important. Rho-GDP-dissociation inhibitor-γ (Rho-GDIγ) is a member of the Rho protein family that regulates many important physiological processes. In this study, we demonstrated the mechanism of suppressing TAK1 expression in the teleost and found that Rho-GDIγ negatively regulated the NF-κB signaling pathway mediated by TAK1. We determined that TAK1 could directly interact with Rho-GDIγ. It is interesting that Rho-GDIγ promotes TAK1 degradation through the ubiquitin proteasome pathway. This study brings a new experimental basis to the teleost fish innate immune signaling pathway. Moreover, this discovery may provide new insights into innate immune regulation mechanism in mammals.

Long non-coding RNA LTCONS7822 positively regulates innate immunity by targeting MITA in miiuy croaker (Miichthys miiuy)

Accumulated studies have shown that long non-coding RNA (lncRNA) is considered a critical regulatory factor in mammals, with a length greater than 200 nucleotides, and it can participate in gene imprinting, dose compensation, transcription enhancement, and antisense regulation. Most of the above studies are carried out in mammals, and there are very few studies on lncRNA of lower vertebrates. Here, we report a novel lncRNA, LTCONS7822, which can play a positive regulatory effect on antiviral immunity in miiuy croaker, Miichthys miiuy. Our results show that the levels of lncRNA LTCONS7822 were significantly increased after poly (I:C) stimulation. Further study, we found that lncRNA LTCONS7822 could positively regulate MITA at the post-transcriptional level. In addition, the dual-luciferase reporter assay analysis showed that the positive regulatory effect of lncRNA LTCONS7822 on NF-κB and IRF3 signaling pathways presented the dose and time-dependent manner. Western blotting experiments proved that lncRNA LTCONS7822 has a positive regulatory effect on MITA. Collectively, our study provided new information to enrich the immune regulation network of lncRNA in teleost fish.

Long Noncoding RNA MIR122HG Inhibits MAVS-Mediated Antiviral Immune Response by Deriving miR-122 in Miiuy Croaker (Miichthys miiuy)

Long noncoding RNAs (lncRNAs) function as micro regulators to impact gene expression after multiple pathogen infections, which have been largely studied in the last few years. Although lncRNA studies on lower vertebrates have received less attention than those on mammals, current studies suggest that lncRNA plays an important role in the immune response of fish to pathogen infections. Here, we studied the effect of MIR122HG as the host gene of miR-122 and indirectly negatively regulate MAVS-mediated antiviral immune responses in miiuy croaker (Miichthysmiiuy). We found that poly(I:C) significantly increases the host MIR122HG expression. The increased MIR122HG expression inhibited the production of the antiviral immune-related genes IFN-1, ISG15 and Viperin upon SCRV treatment. In addition, MIR122HG can act as a pivotally negative regulator involved in the MAVS-mediated NF-κB and IRF3 signaling pathways, which can effectively avoid an excessive immune response. Additionally, we found that MIR122HG can promote the replication of SCRV. Our study provides evidence about the involvement of lncRNAs in the antiviral immune response of fish and broadens the understanding of the function of lncRNAs as a precursor miRNA in teleost fish.

Characterization of MDA5 and microRNA-203 negatively regulates the RLR signaling pathway via targeting MDA5 in miiuy croaker

MDA5 is a member of retinoic acid-inducible gene I (RIG-I)-like receptors (RLR receptors), which may play a crucial role in the immune regulation process. Recently, microRNAs (miRNAs) have been shown to act as an important regulator in the RLRs signaling pathway. Additionally, the MDA5 gene, as a significant cytosolic pathogen recognition receptor (PRR), its characteristics and functions have been extensively investigated, while less research has been done on the mechanisms of MDA5-miRNA mediated gene regulation. In this study, the evolution and functional characterization of MDA5 from miiuy croaker (mmiMDA5) were characterized. Comparative genomic analysis demonstrated that the ascidiacea and superclass do not have the MDA5 gene in the process of evolution. MDA5 contains four structural domains: CARD, ResIII, Helicase C, and RIG-I C-RD. The MDA5 was ubiquitously expressed in all tested miiuy croaker tissues. Moreover, the expressions were significantly up-regulated after stimulation with poly (I: C), which indicated that MDA5 might be involved in the antiviral immune response. The bioinformatics predicted programs have indicated that miR-203 has a direct negative regulatory effect on MDA5 in miiuy croaker. Furthermore, the dual-luciferase reporter assay have showed that miR-203 was the direct negative regulator of MDA5 in miiuy croaker. This study is the first to demonstrate that miRNA can suppress cytokines by regulating the RLR signaling pathway in teleost fish, providing some new ideas for studying miRNA-mediated regulation of immune responses in mammals.

Genome-wide identification and analysis of chemokine receptor superfamily in miiuy croaker, Miichthys miiuy

The chemokine receptor (ChemR) superfamily, which is divided into 4 subfamilies (CXCR, CCR, XCR, and CX3CR), is the main receptors of chemokines in innate immune responses. In the current study, we have identified 27 ChemRs in miiuy croaker: 13 CCR genes, 11 CXCR genes, and 3 XCR genes. Multiple characteristics of these genes, including phylogeny, gene structures, conserved motifs, chromosome locations, evolutionary mechanism, and expression levels upon the bacterial challenge were analyzed. Gene structure and location analysis showed that all ChemR genes contain fewer introns (≤4) and they are unevenly distributed on the 12 chromosomes. And the XCR subfamily of miiuy croaker don't have the DRY motif of ChemR. Phylogenetic and synteny analysis showed that these genes experienced tandem and segmental duplication event in several species, and tandem duplication might be the main expansion way in miiuy croaker. The major ChemRs of each orthologous group in vertebrates were selected for molecular evolution analysis, the results of which indicated that compared with vertebrates, ChemRs of teleost fishes may have a relatively high evolutionary dynamic. In addition, a total of 21 positively selected codons were detected in vertebrate ChemRs under Model 8. RNA-Seq analysis and qRT-PCR verification demonstrated that CXCR3.2, CXCR5, and XCR1 genes were up-regulated significantly upon the Vibrio harveyi infection. These results provide valuable information for investigating the evolutionary relationships of chemokine receptor superfamily in miiuy croaker and laid the basis for further functional analysis.

Genome-guided transcriptome analysis of miiuy croaker provides insights into pattern recognition receptors and cytokines in response to Vibrio anguillarum

Miiuy croaker (Miichthys miiuy), as an economically important marine fish is affected by numerous bacterial diseases. Infection by bacterial pathogen Vibrio anguillarum causes high mortality and great economic loss in aquaculture. To understand the immune response of the miiuy croaker to V. anguillarum infection, we analyzed the transcriptomic profile of V. anguillarum-challenged M. miiuy, and two cDNA libraries, namely, the normal library (Ctrl) and V. anguillarum challenged library (Van) were constructed. Combined with the whole genome-guided assembly, total 47,971 unigenes were acquired. Moreover, 2482 significantly differentially expressed unigenes were identified based on the expression patterns. Immune-related differentially expressed genes (DEGs) that were significantly up- or down-regulated after V. anguillarum injection were identified via enrichment analysis using the Gene Ontology and Kyoto Encyclopedia of Genes and Genomes databases. To further understand the expression profiles that underlie the response to V. anguillarum infection, we investigated the DEGs that are related to pattern-recognition receptors and cytokines. The expression patterns of relevant DEGs were validated by qRT-PCR. Additionally, the expression profiles of DEGs that are related to pattern-recognition receptors and cytokines were further measured in miiuy croaker macrophages. This study is the first to identify and characterize the transcriptome of miiuy croaker in response to V. anguillarum infection. The results expand the general understanding of the immune defense mechanisms of miiuy croaker and provide helpful information for counteracting V. anguillarum infection.

miR-192 targeting IL-1RI regulates the immune response in miiuy croaker after pathogen infection in vitro and in vivo

Activation of innate and acquired immune responses is regulated by detailed mechanisms to control their onset and termination. MicroRNAs have been implicated as negative regulators controlling the diverse of biophysical and biochemical processes at the post-transcriptional level. However, the physiological roles of miRNAs in aquatic organisms are largely unclear. In this study, we explored the potential roles of mmi-miR-192 in regulating interleukin 1 receptor type I (IL-1RI) involved in immune and inflammatory response in miiuy croakers. This was further evaluated by negative expression profiles in both LPS exposure macrophages and Vibrio anguillarum challenged miiuy croaker. By means of promoter analysis, mmi-miR-192 was found to be an AP-1 dependent gene. Importantly, the dual luciferase reporter assay presented the regulation between mmi-miR-192 and IL-1RI. The result of miiuy croaker miR-192 reduced the wild-type IL-1RI but not the mutant one luciferase levels suggested that mmi-miR-192 modulated IL-1RI expression by directly targeting the 3'UTR of IL-1RI mRNA. Overall, our study revealed the mechanism that the miR-192-IL1RI pathway regulated bacteria infection in miiuy croakers.

Molecular characterization, functional analysis, and defense mechanisms of two CC chemokines in Nile tilapia (Oreochromis niloticus) in response to severely pathogenic bacteria

Two full-length cDNAs encoding CC chemokine genes in Nile tilapia (Oreochromis niloticus) (On-CC1 and On-CC2) were cloned and characterized. On-CC1 and On-CC2 showed signature cysteine motifs consisting of four cysteines. The expression levels of On-CC1 and On-CC2 were analyzed by RT-PCR, which showed that low expression of these two genes was only observed in the peripheral blood leukocytes (PBLs) and spleen of normal fish. Expression levels of these two molecules were quantified in 13 tissues of fish infected with virulent strains of Streptococcus agalactiae and Flavobacterium columnare. Most tissues, especially PBLs, the spleen and the liver, expressed significantly higher mRNA levels than the controls, particularly at 12 and 24 h after infection (P < 0.05). The current study strongly indicates that CC chemokine genes in Nile tilapia are crucially involved in the early immune responses to pathogens. Functional analyses clearly demonstrated that 10 and 100 μg/ml of recombinant rOn-CC1 and rOn-CC2 proteins efficiently enhanced the phagocytic activity (in vitro) of Nile tilapia phagocytes. Finally, Southern blot analysis and searching in Ensembl databases demonstrated that two different functional CC chemokine genes and other pseudogene fragments were discovered in the Nile tilapia genome.

Comparative analysis of the small RNA transcriptomes of miiuy croaker revealed microRNA-mediated regulation of TLR signaling pathway response to Vibrio anguillarum infection

MicroRNAs (miRNAs) are an abundant class of endogenous noncoding small RNAs (sRNAs) that are partially complementary to their target messenger RNA (mRNA), which post-transcriptionally regulate various biological processes by repressing mRNA translation or inducing mRNA degradation. MiRNAs have been demonstrated to play crucial roles in the host response to infection by diverse pathogens. As an important bacterial pathogen, Vibrio anguillarum has been caused great economic losses in miiuy croaker aquaculture. To identify immune-related miRNAs of miiuy croaker in response to V. anguillarum infection, we constructed two sRNA libraries with or without bacterial infection. High-throughput deep sequencing and subsequent bioinformatic analysis identified 241 conserved and 137 novel miRNA precursors in miiuy croaker based on its whole genome sequences, encoding 293 and 124 mature miRNAs, respectively. Then we compared the expression patterns of miRNAs in the two libraries. There were significant differences in the expression of 12 miRNAs between the infection group (IG) and control group (CG). Further, the expressions of six miRNAs were validated by real-time quantitative PCR. The target gene prediction and function analysis were conducted for the 12 differential miRNAs. This analysis revealed that these miRNAs participated in the regulation multiple immune-related signaling pathways. Transcription factors in TLR signaling, such as AP-1, IRF5, NF-κB and IRF3, were activated by these miRNAs via post-transcriptionally regulating the expression of TLRs and TLR-associated signaling proteins, inducing effective host immune response to eradicate infectious pathogens. This is the first study of the identification and characterization of miiuy croaker miRNAs in response to V. anguillarum infection. The comprehensive analysis of the expression of miRNAs and the target gene and function prediction of differently expressed miRNAs may help to understand the regulatory mechanisms of miRNA in fish during the interaction between host and bacterial pathogens.

A CCL21 chemokine of tongue sole (Cynoglossus semilaevis) promotes host resistance against bacterial infection

Chemokines are a large family of chemotactic cytokines. Based on the arrangement of the first two cysteine residues, chemokines are divided into four groups, one of which is the CC chemokine group. In this study, we characterized a CC chemokine, CsCCL21, from half-smooth tongue sole (Cynoglossus semilaevis), and analyzed its activity. CsCCL21 contains two conserved N-terminal cysteine residues in a NCCL motif and is phylogenetically related to the CCL19/21/25 subgroup of CC chemokines. CsCCL21 was constitutively expressed in nine tissues and significantly upregulated by bacterial and viral infection. The recombinant CsCCL21 (rCsCCL21) induced migration of peripheral blood leukocytes. When the two conserved cysteine residues in the NCCL motif were mutated, the chemotactic activity of rCsCCL21 was abolished. rCsCCL21 enhanced the resistance of tongue sole against bacterial infection, but the mutant protein with NCCL mutation lacked this antibacterial effect. Taken together, these results suggest that CsCCL21 is a functional CC chemokine with the ability to recruit leukocytes and is involved in antibacterial immunity in a manner that requires the conserved NCCL motif.

CsCCL17, a CC chemokine of Cynoglossus semilaevis, induces leukocyte trafficking and promotes immune defense against viral infection

CC chemokines are the largest subfamily of chemokines, which are important components of the innate immune system. To date, sequences of several CC chemokines have been identified in half-smooth tongue sole (Cynoglossus semilaevis); however, the activities and functions of these putative chemokines remain unknown. Herein, we characterized a CC chemokine, CsCCL17, from tongue sole, and examined its activity. CsCCL17 contains a 303 bp open reading frame, which encodes a polypeptide of 100 amino acids with a molecular mass of 12 kDa CsCCL17 is phylogenetically related to the CCL17/22 group of CC chemokines and possesses the typical arrangement of four cysteines and an SCCR motif found in known CC chemokines. Under normal physiological conditions, CsCCL17 expression was detected in spleen, liver, heart, gill, head kidney, muscle, brain, and intestine. When the fish were infected by bacterial and viral pathogens, CsCCL17 expression was significantly up-regulated in a time-dependent manner. Chemotactic analysis showed that recombinant CsCCL17 (rCsCCL17) induced migration of peripheral blood leukocytes. A mutagenesis study showed that when the two cysteine residues in the SCCR motif were replaced by serine, no apparent chemotactic activity was observed in the mutant protein rCsCCL17M. rCsCCL17 enhanced the resistance of tongue sole against viral infection, but rCsCCL17M lacked this antiviral effect. Taken together, these findings indicate that CsCCL17 is a functional CC chemokine with the ability to recruit leukocytes and enhance host immune defense in a manner that requires the conserved SCCR motif.

Comparative genomic of the BAFF and BAFF-like genes and immune response to bacteria of miiuy croaker (Miichthys miiuy)

As a member of the tumor necrosis factor (TNF) family, B cell activating factor (BAFF), also known as TNF ligand superfamily member 13B (TNF13B), playing a critical role in enhancing immune responses. BAFF is a central cytokine for B-cell survival, proliferation, maturation and immunoglobulin secretion. In the present study, we describe the identification of the miiuy croaker BAFF (designed MmBAFF) and BAFF-like (designed MmBAFF-like) genes. The cDNA of MmBAFF contains an open reading frame (ORF) of 795 nucleotides that are translated into a predicted 264 amino acids. The ORF of MmBAFF-like consists of 705 bases encoding 234 amino acids. Amino acid sequence comparison indicated that MmBAFF and MmBAFF-like possessed the TNF signatures, a predicted transmembrane domain, three conserved cysteine residues and a putative furin protease cleavage site, which were the typical characteristics of TNF gene in mammals and birds. The predicted three-dimensional (3D) structure of the MmBAFF and MmBAFF-like monomer analyzed by comparative protein modeling revealed that they were very similar to human counterpart. Comparative genomic analysis revealed that the locations of MmBAFF and MmBAFF-like genes are conserved among the bony fish. Phylogenetic analysis shows the MmBAFF is most closely related to other teleost BAFFs with the highest similarity to Epinephelus awoara. And BAFF-like cluster get together first to BAFF cluster than three closely related TNF superfamily (TNFSF) members. Real-time quantitative PCR analysis shows the MmBAFF and MmBAFF-like transcripts are expressed in a wide range of tissues with the highest expression in skin and lymphoid tissue spleen. Upon induction by Vibrio anguillarum, their expressions are significantly upregulated in liver, spleen and kidney as compared to phosphate-buffered saline injected control fish. The association of increased BAFF expression after bacterial infection suggests that it plays a potentially important role in immune system of fish.

Cloning and expression analysis of three novel CC chemokine genes from Japanese flounder (Paralichthys olivaceus)

Chemokines are small cytokines secreted by various cell types. They not only function in cell activation, differentiation and trafficking, but they also have influences on many biological processes. In this study, three novel CC chemokine genes Paol-SCYA105, 106 and 107 in Japanese flounder (Paralichthys olivaceus) were cloned and characterized. Paol-SCYA105 was mainly detected in gill, kidney and spleen, Paol-SCYA106 was detected in all tissues examined and Paol-SCYA107 was mainly detected in the spleen and kidney. Paol-SCYA105 and Paol-SCYA106 gene expressions peaked in kidney at day 3 after viral hemorrhagic septicemia virus infection and decreased at day 6, but Paol-SCYA106 still remained at a high level at day 6. Paol-SCYA107 gene expression was significantly up-regulated in kidney at day 6 after viral hemorrhagic septicemia virus infection. In response to infection by Gram-negative Edwardsiella tarda and Gram-positive Streptococcus iniae in kidney, only Paol-SCYA106 gene expression significantly increased. Together, these results indicate that these three novel CC chemokines are involved in the immune response against pathogen infections.

Three representative subtypes of caspase in miiuy croaker: genomic organization, evolution and immune responses to bacterial challenge

Caspase proteins are intracellular proteases which function as initiators and effectors of apoptosis. According their difference of functions, the caspases can be divided into apoptosis related caspases and inflammatory mediator, and the former included apoptosis activator and apoptosis executioner. In this study, three different subtype caspases (caspase1, caspase3 and caspase9) from Miichthys miiuy miiuy croaker were analyzed. The caspase1 belongs to the inflammatory mediator, caspase3 belongs to apoptosis executioner, and caspase9 belongs to apoptosis activator. Miichthys miiuy caspase1 (Mmcaspase1) and Mmcaspase9 exhibited three conserved domains including a CARD, a large subunit p20 and a small subunit p10. The Mmcaspase3 encoded two conserved domains including a large subunit p20 and a small subunit p10. Mmcaspase3 and Mmcaspase9 contained a histidine active sequence and a cysteine active sequence. However, Mmcaspase1 only contained the cysteine active sequence. The real-time PCR (RT-PCR) analysis showed these three caspases were expressed constitutively in all examined tissues in miiuy croaker, although the expression levels varied from tissue to tissue. Expression analysis showed that Mmcaspase1 was up-regulated obviously in liver, spleen and kidney, and indicated its positive role in response to Vibrio anguillarum infection, but Mmcaspase3 and Mmcaspase9 showed different expression pattern in liver, spleen and kidney, its showed that different subtypes of caspase having different immune response mechanisms. These results revealed that the organs adjusted the expressions of these three genes with the infected of pathogens, suggesting the immunoprotection of these genes.

Identification of 21 novel immune-type receptors in miiuy croaker and expression pattern of three typical inhibitory members

Novel immune-type receptor (NITR) genes belong to the immunoglobulin superfamily and are encoded by clusters of multigene families. NITRs encode type I transmembrane proteins and are only found in teleosts. In the current study, total 21 NITR genes are identified from miiuy croaker (Miichthys miiuy) and named as MmNITR1 to MmNITR21. Miiuy croaker NITR genes that encoded one or two extracellular immunoglobulin (Ig) domains, a transmembrane (TM) region, an immunoreceptor tyrosine-based inhibitor motif (ITIM) in the cytoplasmic (Cyt) region. The majority of MmNITRs possess cytoplasmic ITIM that can be classified as inhibitory receptors. However, a smaller number of NITRs (MmNITR8, MmNITR15 and MmNITR16) can be classified as activating receptors by the lack of cytoplasmic ITIMs and presence of a positively charged residue within their transmembrane domain. As typical inhibitory receptors, MmNITR1, MmNITR2 and MmNITR3 have different characteristics of the structure. In MmNITR1 gene, variable (V) and intermediate (I) domains are encoded by two separate exons. In contrast to MmNITR1, MmNITR3 gene encode V and I domains in a single exon. And MmNITR2 gene is characterized by the presence of only one Ig-like (V-type) extracellular domain and lack of J or J-like motifs. Also MmNITR2 gene displays an additional exon which is 48bp long between the V domain and the TM region. Two and four potential N-link giycosylation sites (N-X-S/T) are present in the extracellular Ig domains. Real-time RT-PCR results showed that upon induction with Vibrio anguillarum, NITR gene expressions were induced by bacteria in kidney, liver and spleen. Meanwhile, NITRs are also primarily detected in different tissues. Phylogenetic analyses of NITR V domains indicate that MmNITR1 and MmNITR2 are more similar than MmNITR3.

Characterization and expression of the CXCR1 and CXCR4 in miiuy croaker and evolutionary analysis shows the strong positive selection pressures imposed in mammal CXCR1

The innate immune system can recognize non-self, danger signals, and pathogen associated molecular patterns and provides a first line of antimicrobial host defense. Therefore, it plays an instructive role and is pretty important in vertebrates. In innate immune responses, CXCRs act as the main receptors of CXC chemokines and play a vital role in host defense and inflammation. In present study, we cloned two cDNA molecules of CXCR1 and CXCR4 in Miichthys miiuy (miiuy croaker). In these two genes, we found the most highly conserved DRY motif in the second intracellular loop adjacent to the third transmembrane domain. The expressions of CXCR1 and CXCR4 showed that they were ubiquitously expressed in ten normal tissues. After infection with Vibrio anguillarum and Vibrio harveyi, the expressions of CXCRs in the immune tissues were significantly regulated in most of tissues except that of CXCR1 in the kidney after V. harveyi injection. Evolutionary analysis showed that only the ancestral lineages of CXCR4 in amphibians underwent positive selection, indicating that the ancestors of amphibians boarded the land and had to further evolve to adapt to terrestrial environments. Multiple ML methods were implemented to detect the robust positively selected candidates for sites. In total, we detected 12 and 3 positively selected sites in the subsets of current mammal and fish CXCR1 genes, and only one site under positive selection was found in mammalian CXCR4 subsets. These positively selected sites were mainly located in the extracellular domains of CXCRs. The sliding window analysis and evolution test tended to favor positive selection acting on the N-terminal domain of CXCR1, which was the critical region for ligand/receptor signaling for neutrophils and receptor-ligand interaction, indicating that the N-terminal of CXCR1 in mammals underwent more positive selection than that of fish.

Characterization of the miiuy croaker (Miichthys miiuy) transcriptome and development of immune-relevant genes and molecular markers

Background

The miiuy croaker (Miichthys miiuy) is an important species of marine fish that supports capture fisheries and aquaculture. At present commercial scale aquaculture of this species is limited due to diseases caused by pathogens and parasites which restrict production and limit commercial value. The lack of transcriptomic and genomic information for the miiuy croaker limits the ability of researchers to study the pathogenesis and immune system of this species. In this study we constructed a cDNA library from liver, spleen and kidney which was sequenced using Illumina paired-end sequencing to enable gene discovery and molecular marker development.

Principal Findings

In our study, a total of 69,071 unigenes with an average length of 572 bp were obtained. Of these, 45,676 (66.13%) were successfully annotated in public databases. The unigenes were also annotated with Gene Ontology, Clusters of Orthologous Groups and KEGG pathways. Additionally, 498 immune-relevant genes were identified and classified. Furthermore, 14,885 putative simple sequence repeats (cSSRs) and 8,510 putative single nucleotide polymorphisms (SNPs) were identified from the 69,071 unigenes.

Conclusion

The miiuy croaker (Miichthys miiuy) transcriptome data provides a large resource to identify new genes involved in many processes including those involved in the response to pathogens and diseases. Furthermore, the thousands of potential cSSR and SNP markers found in this study are important resources with respect to future development of molecular marker assisted breeding programs for the miiuy croaker.

Gene structure, immune response and evolution: comparative analysis of three 2-Cys peroxiredoxin members of miiuy croaker, Miichthys miiuy

Peroxiredoxin family was a superfamily of selenium independent peroxidases. It was divided into six subtypes: Prx1-4 (typical 2-Cys), Prx5 (atypical 2-Cys) and Prx6 (1-Cys). This study reports the isolation and characterization three 2-Cys peroxiredoxin members of full cDNA and genomic clones from miiuy croaker (Miichthys miiuy). The genetic structure analysis showed that the C-terminal catalytic Cys positioned within GEVCPAXW. This sequence was different between Prx3 and Prx4, but was conservative in different species of the same gene, the X base was S in Prx3 but G in Prx4. Tissues expression analysis showed that the expressions of Prx3 in liver and brain were much higher than other tissues; the values of Prx4 in spleen, intestine and kidney were significantly higher than others; and the expression of Prx5 in muscle was higher than that of other tissues. Real-time PCR results showed that there were highest values of these three Prxs emerging with the time post challenge of Vibrio anguillarum in liver, spleen and kidney although the highest value time differed from each other and the expression of these three genes also changed with the change of infection time. These results indicated that expression analysis of these three genes play some positive function against pathogenic bacteria infection in miiuy croaker.

Immune gene discovery in the crucian carp Carassius auratus

The crucian carp Carassius auratus (Cyprinidae) is one of the important fish species in aquaculture. Although the crucian carp has several economic benefits, their immune system and gene information have not been investigated in depth as yet. Here, we performed the transcriptome analysis of C. auratus using the pyrosequencing method and selected several immune-related genes. Of unigenes obtained in this species, we identified a number of immune system-related genes (e.g. adhesive protein, antimicrobial protein, apoptosis- and cell cycle-related protein, cellular defense effector, immune regulator, pattern recognition protein, protease, protease inhibitor, reduction/oxidation-related protein, signal transduction-related protein and stress protein) that are potentially useful for studies on fish immunity. To be of public and practical use, we designed primer pairs of each gene from the crucian carp for real-time RT-PCR application and tested the amplicon identity of entire gene sets with the total RNA sample. For comparative analysis, we measured tissue-preferential transcript profiles of selected genes. This study will be helpful to extend our knowledge on the immune system of the crucian carp in comparative aspects and to develop the crucian carp as a potential model organism for aquatic quality monitoring in fish farming.

Characterization of the CCR3 and CCR9 genes in miiuy croaker and different selection pressures imposed on different domains between mammals and teleosts

The innate immune system can recognize non-self through pattern recognition receptors and provides a first line of antimicrobial host defense. Thus innate immunity plays a very important role in resistance against major bacterial disease in vertebrates. In the innate immune responses, the chemokine receptors act as the main receptors of chemokines which are released at the sites of infection, inflammation and injury. In this study, the Miichthys miiuy CCR3 and CCR9 genes were cloned and characterized, showing that MIMI-CCR3 possesses a highly conserved DRYLA motif similar to that of other fishes. MIMI-CCR3 and CCR9 were ubiquitously expressed in all tested tissues and the expressions were significantly up-regulated after infection with Vibrio anguillarum except that of CCR9 in spleen. Evolutionary analysis showed that all the ancestral lineages of CCR3 and CCR9 in mammals and teleosts underwent positive selection, indicating that the ancestor of terrestrial animals further evolved to adapt to terrestrial environments and the continuous intrusion of microbes stimulated the evolution of CCR genes in the ancestor of teleost. Multiple ML methods were used to detect the robust candidates for sites under positive selection. In total, 11 and 8 positively selected sites were found in the subsets of current mammal and teleost CCR3 genes, and 3 and 2 sites were detected in the subsets of current mammals and teleosts in CCR9. Interestingly, for mammal CCR3 and CCR9 genes, the robust candidates of positively selected sites were mainly located in the extracellular domains which were the ligand binding and pathogen interaction regions. For teleost CCR3 and CCR9 genes, the positively selected sites were not only located in the extracellular domains but also in the C-terminal and intracellular domains, indicating mammals and teleosts experienced different selection pressures upon their N-terminus, C-terminus and intracellular loops of CCRs.

Immunological role of C4 CC chemokine-1 from snakehead murrel Channa striatus

In this study, we have reported a cDNA sequence of C4 CC chemokine identified from snakehead murrel (also known as striped murrel) Channa striatus (named as CsCC-Chem-1) normalized cDNA library constructed by Genome Sequencing FLX™ Technology (GS-FLX™). CsCC-Chem-1 is 641 base pairs (bp) long that contain 438 bp open reading frame (ORF). The ORF encodes a polypeptide of 146 amino acids with a molecular mass of 15 kDa. The polypeptide contains a small cytokine domain at 30-88. The domain carries the CC motif at Cys(33)-Cys(34). In addition, CsCC-Chem-1 consists of another two cysteine residues at C(59) and C(73), which, together with C(33) and C(34), make CsCC-Chem-1 as a C4-CC chemokine. CsCC-Chem-1 also contains a 'TCCT' motif at 32-35 as CC signature motif; this new motif may represent new characteristic features, which may lead to some unknown function that needs to be further focused on. Phylogenitically, CsCC-Chem-1 clustered together with CC-Chem-1 from rock bream Oplegnathus fasciatus and European sea bass Dicentrarchus labrax. Significantly (P<0.05) highest gene expression was noticed in spleen and is up-regulated upon fungus (Aphanomyces invadans), bacteria (Aeromonas hydrophila) and virus (poly I:C) infection at various time points. The gene expression results indicate the influence of CsCC-Chem-1 in the immune system of murrel. Overall, the gene expression study showed that the CsCC-Chem-1 is a capable gene to increase the cellular response against various microbial infections. Further, we cloned the coding sequence of CsCC-Chem-1 in pMAL vector and purified the recombinant protein to study the functional properties. The cell proliferation activity of recombinant CsCC-Chem-1 protein showed a significant metabolic activity in a concentration dependent manner. Moreover, the chemotaxis assay showed the capability of recombinant CsCC-Chem-1 protein which can induce the migration of spleen leukocytes in C. striatus. However, this remains to be verified further at molecular and proteomic level.

A heat shock protein 90 β isoform involved in immune response to bacteria challenge and heat shock from Miichthys miiuy

Heat shock protein 90 (HSP90) is highly conserved molecular chaperone that plays a critical role in cellular stress response. In this study, we reported the identification and functional analysis of a heat shock protein 90 gene from miiuy croaker (designated Mimi-HSP90). Mimi-HSP90 contained five conserved HSP90 protein family signatures and shared 89.6%-99.5% similarity with other known HSP90 β isoform. Homology analysis and structure comparison further indicated that Mimi-HSP90 should be β isoform member of the HSP90 family. The molecular evolutionary analysis showed that HSP90 was under an overall strong purifying select pressure among fish species. Mimi-HSP90 gene was constitutively expressed in ten examined tissues, and the expression level of liver was higher than in other tissues. The expression level of Mimi-HSP90 gene under bacterial infection and heat shock were analyzed by real-time quantitative RT-PCR, resulted in significant changes in liver, spleen, and kidney tissues. The purified recombinant pET-HSP90 protein was used to produce the polyclonal antibody in mice. The specificity of the antibody was determined by Western blot analysis. All results suggested that Mimi-HSP90 was involved in thermal stress and immune response in miiuy croaker.

Molecular characterisation and biological activity of a novel CXC chemokine gene in rock bream (Oplegnathus fasciatus)

Chemokines are chemoattractant cytokines defined by the presence of four conserved cysteine residues. In mammals, these cytokines can be divided into four subfamilies depending on the arrangement of the first two conserved cysteines in the sequence, and include the CXC(α), CC(β), C(γ), and CX3C(δ) classes. We identified CXC chemokine cDNA, designated RbCXC, isolated using expressed sequence tag analysis of a lipopolysaccharide (LPS)-stimulated rock bream liver cDNA library. The full-length RbCXC cDNA (742 bp) contained an open reading frame of 342 bp encoding 114 amino acids. Results from phylogenetic analysis showed that RbCXC was strictly separated into a distinct clade compared to other known CXC chemokine subgroups. RbCXC was significantly expressed in the trunk kidney, liver, spleen, gill, peripheral blood leukocytes (PBLs), and head kidney. Rock bream PBLs were stimulated with several mitogens, including LPS and polyinosinic-polycytidylic acid (poly I:C), which significantly induced the expression of RbCXC mRNA. RbCXC mRNA expression was examined in several tissues under conditions of bacterial and viral challenge. Experimental challenges revealed that all examined tissues from fish infected with Edwardsiella tarda and red sea bream iridovirus showed significant increases in RbCXC expression compared to the control. In the case of Streptococcus iniae infection, RbCXC mRNA expression was markedly upregulated in the kidney, spleen, and liver. In addition, a maltose binding protein fusion recombinant RbCXC (~53 kDa) was produced in an Escherichia coli expression system and purified. Subsequently, the addition of purified recombinant RbCXC (rRbCXC) to kidney leukocytes was examined to investigate the impact of proliferative and chemotactic activity. The rRbCXC induced significant kidney leukocyte proliferation and attraction at concentrations ranging from 10 to 300 μg/mL, suggesting that it can be utilised as an immune stimulant and/or molecular adjuvant to enhance the immunological effects of vaccines.

Molecular identification and expression analysis of the CC chemokine gene in rock bream (Oplegnathus fasciatus) and the biological activity of the recombinant protein

We identified the CC chemokine cDNA designated as RbCC1 (CC chemokine 1 in rock bream, Oplegnathus fasciatus), which was isolated using expressed sequence tag (EST) analysis of a lipopolysaccharide (LPS)-stimulated rock bream liver cDNA library. The full-length RbCC1 cDNA (850 bp) contained an open reading frame (ORF) of 366 bp encoding 122 amino acids. Results from our phylogenetic analysis demonstrated that the RbCC1 was closest relationship to the orange-spotted grouper and Mi-iyu croaker CC chemokines located within the fish CC chemokine group. RbCC1 was significantly expressed in the intestine, spleen, liver, and PBLs (peripheral blood leukocytes). Rock bream PBLs were stimulated with several mitogens, LPS and Con A/PMA which significantly induced the expression of RbCC1 mRNA in the PBLs. The RbCC1 mRNA expression in several tissues under conditions of bacterial and viral challenge was examined. The experimental challenge revealed that the kidney and spleen of fish infected with Streptococcus iniae showed the most significant increases in RbCC1 expression compared to the control. In the case of RSIV infection, the RbCC1 mRNA expression was markedly up-regulated in the liver. In this study, recombinant RbCC1 (approximately 53 kDa) was produced using an Escherichia coli expression system followed by purification. Subsequently, the addition of purified rRbCC1 was examined to investigate the impact on the proliferative and chemotactic activity on kidney leukocytes from rock bream. The results demonstrated that the rRbCC1 induces significant biological activity on kidney leukocyte proliferation and attraction at concentrations in the range of 10-300 μg/mL and suggests that rRbCC1 could be utilized as an immune-stimulant and/or molecular adjuvant to enhance the immune effects of vaccines.

Expression profile of immune-related genes in Lates calcarifer infected by Cryptocaryon irritans

Cryptocaryon irritans causes Cyptocaryonosis or white spot disease in a wide range of marine fish including Lates calcarifer (Asian seabass). However, the immune response of this fish to the parasite is still poorly understood. In this study, quantitative polymerase chain reaction (qPCR) was performed to assess the expression profile of immune-related genes in L. calcarifer infected by C. irritans. A total of 21 immune-related genes encoding various functions in the fish immune system were utilized for the qPCR analysis. The experiment was initiated with the infection of juvenile fish by exposure to theronts from 200 C. irritans cysts, and non-infected juvenile fish were used as controls. Spleen, liver, gills and kidney tissues were harvested at three days post-infection from control and infected fish. In addition, organs were also harvested on day-10 post-infection from fish that had been allowed to recover from day-4 up to day-10 post-infection. L. calcarifer exhibited pathological changes on day-3 post-infection with the characteristic presence of white spots on the entire fish body, excessive mucus production and formation of a flap over the fish eye. High quality total RNA was extracted from all tissues and qPCR was performed. The qPCR analysis on the cohort of 21 immune-related genes of the various organs harvested on day-3 post-infection demonstrated that most genes were induced significantly (p < 0.05) in all tissues, particularly liver (11/21 genes) and kidney (11/21). The expression profile demonstrated that induction of the MHC Class IIα gene was the highest compared to the other genes followed by serum amyloid A, CC chemokine and hepcidin-2 precursor genes. In fish that were allowed to recover from the C. irritans infection (10 days post-infection), expression of the immune-related genes was down-regulated to levels similar to the control fish. These results provide insights into the interaction between C. irritans and L. calcarifer and suggest that the innate immune system plays an important role in early defence against parasite infection allowing the fish to eventually recover from the infection.

Miiuy croaker (Miichthys miiuy) Peroxiredoxin2: molecular characterization, genomic structure and immune response against bacterial infection

Peroxiredoxin2 (Prx2) protein is an important member in cellular antioxidant protein superfamily. Prx2 exists widely in prokaryotes and eukaryotes, it not only plays a part in eliminate reactive oxygen, but also takes effect in many other metabolic activities, such as stimulate epithelial cell proliferation, participate in the signal transduction in cells and so on. After molecular cloning we got that the complete cDNA sequence of Prx2 consists 882 bp, including a 5'-UTR of 46 bp, an open reading frame (ORF) of 591 bp, and a 3'-UTR of 245 bp. The complete gene of miiuy croaker Prx2 has 5 exons and 4 introns. The deduced 197 amino acid residues of miiuy croaker Prx2 consists a Val-Cys-Pro (VCP) motifs. In order to better elucidate the immune mechanisms of the Prx2 in the lower vertebrates, we conducted a research about the Prx2 gene of miiuy croaker and its expression pattern after bacterial infection. Real-time PCR (RT-PCR) results showed that expression of Prx2 was up-regulated in kidney, liver and spleen under infection with Vibrio anguillarum, and expressed level differently in ten different uninjected tissues. Our results suggested that Prx2 might be involved in immune defence in miiuy croaker.

Molecular cloning, characterization and expression analysis of a CC chemokine gene from miiuy croaker (Miichthys miiuy)

Chemokines are a family of structurally related chemotactic cytokines that regulate the migration of leukocytes, under both physiological and inflammatory conditions. A partial cDNA of CC chemokine gene designed as Mimi-CC3 was isolated from miiuy croaker (Miichthys miiuy) spleen cDNA library. Unknown 3' part of the cDNA was amplified by 3'-RACE. The complete cDNA of Mimi-CC3 contains an 89-nt 5'-UTR, a 303-nt open reading frame and a 441-nt 3'-UTR. Three exons and two introns were identified in Mimi-CC3. The deduced Mimi-CC3 protein sequences contain a 22 amino acids signal peptide and a 78 amino acids mature polypeptide, which possesses the typical arrangement of four cysteines as found in other known CC chemokines. It shares low amino acid sequence identities with most other fish and mammalian CC chemokines (less than 54.1 %), but shares very high identities with large yellow croaker CC chemokine (94.6 %). Phylogenetic analysis showed that Mimi-CC3 gene may have an orthologous relationship with mammalian/amphibian CCL25 gene. Tissue expression distributed analysis showed that Mimi-CC3 gene was constitutively expressed in all nine tissues examined, although at different levels. Upon stimulated with Vibrio anguillarum, the time-course analysis using a real-time PCR showed that Mimi-CC3 transcript in kidney and liver was obviously up-regulated and reached the peak levels, followed by a recovery. Mimi-CC3 expression in kidney was more strongly increased than in liver. However, down-regulation was observed in spleen. These results indicated that Mimi-CC3 plays important roles in miiuy croaker immune response as well as in homeostatic mechanisms.

Characterization, genomic organization, and expression profiles of MyD88, a key adaptor molecule in the TLR signaling pathways in miiuy croaker (Miichthys miiuy)

Myeloid differentiation factor 88 (MyD88) is an important adaptor protein in the TLR signaling pathways. In the present study, we firstly cloned and characterized Miichthys miiuy MyD88 (Mimi-MyD88) cDNA and gene. The Mimi-MyD88 gene was 3,470 bp consisting of five exons and four introns. The cDNA was composed of 1,627 bp with an 867-bp open reading frame encoding a polypeptide of 288 amino acid residues. The theoretical molecular mass and isoelectric point of this polypeptide were 33.25 and 4.96 kDa. Comparison of the deduced amino acid sequence showed that the conserved death domain and the typical Toll/IL-1 receptor domain are very similar to those presented in other mammals, amphibians, and fishes. To identify potential role of MyD88 in fish innate immunological surveillance, the constitutive Mimi-MyD88 mRNA is detected by real-time quantitative polymerase chain reaction. Results demonstrated that Mimi-MyD88 is broadly expressed in ten normal tissues, with the lowest expression was observed in kidney and the highest was in liver. The transcriptional expression also revealed that Mimi-MyD88 was significantly up-regulated in liver, kidney, and spleen after challenge by Gram-negative bacteria, Vibrio anguillarum. Via contrasted the expression of MyD88 and TLR2 in kidney, we evaluated TLR2 plays an indispensable role in MyD88 transcription, but not absolutely dominant. The combined expression still indicated that MyD88 plays a universal role in keeping immune surveillance for pathogens. Phylogenetic analysis suggested that Mimi-MyD88 gene is classified into the piscine cluster and most closely related to large yellow croaker Larimichthys crocea.

Characterization and SNP variation analysis of a HSP70 gene from miiuy croaker and its expression as related to bacterial challenge and heat shock

Heat shock proteins (HSPs) play crucial roles in the immune response of vertebrates. In order to study immune defense mechanism of heat shock protein gene in miiuy croaker (Miichthys miiuy), a cDNA encoding heat shock protein 70 (designated Mimi-HSP70) gene was cloned from miiuy croaker. The cDNA was 2195 bp in length, consisting of an open reading frame (ORF) of 1917 bp encoding a polypeptide of 638 amino acids with estimated molecular mass of 70.3 kDa and theoretical isoelectric point of 5.55. Genomic DNA structure analysis revealed that the Mimi-HSP70 gene contain no introns in coding region and four SNPs with 373 C/T, 789 G/A, 1005 C/T, and 1185 G/A were detected by direct sequencing of 20 samples from six different populations. BLAST analysis, structure comparison and phylogenetic analysis indicated that Mimi-HSP70 should be an inducible cytosolic member of the HSP70 family. The deduced amino acid sequence of Mimi-HSP70 had 82.4%-92.2% identity with those of vertebrate. A real-time quantitative RT-PCR demonstrated that the HSP70 gene was ubiquitously expressed in ten normal tissues. Under different temperature shock stress, the expression of Mimi-HSP70 gene in miiuy croaker increased at first and then decreased with the rise of temperature, finally, reached a maximum level in liver, spleen and kidney tissues. Infection of miiuy croaker with Vibrio anguillarum resulted in significant changes expression of Mimi-HSP70 gene in the immune-related tissues. These results indicated that expression analysis of Mimi-HSP70 gene provide theoretical basis to further study the mechanism of anti-adverseness in the miiuy croaker.