Adenylate kinase 1 overexpression increases locomotor activity in medaka fish

Maintenance of the energy balance is indispensable for cell survival and function. Adenylate kinase (Ak) is a ubiquitous enzyme highly conserved among many organisms. Ak plays an essential role in energy regulation by maintaining adenine nucleotide homeostasis in cells. However, its role at the whole organism level, especially in animal behavior, remains unclear. Here, we established a model using medaka fish (Oryzias latipes) to examine the function of Ak in environmental adaptation. Medaka overexpressing the major Ak isoform Ak1 exhibited increased locomotor activity compared to that of the wild type. Interestingly, this increase was temperature dependent. Our findings suggest that cellular energy balance can modulate locomotor activity.

Alternative Fish Species for Nutritional Management of Children with Fish-FPIES—A Clinical Approach

In the Mediterranean region, fish is a common cause of food protein-induced enterocolitis syndrome (FPIES) in children. No laboratory tests specific to FPIES are available, and oral food challenge (OFC) is the gold standard for its diagnosis and testing for achievement of tolerance. Children with FPIES to fish are usually advised to avoid all fish, regardless of the species. Fish are typically classified into bony and cartilaginous, which are phylogenetically distant species and therefore contain less cross-reacting allergens. The protein β-parvalbumin, considered a pan-allergenic, is found in bony fish, while the non-allergenic α-parvalbumin is commonly found in cartilaginous fish. Based on this difference, as a first step in the therapeutic process of children with FPIES caused by a certain fish in the bony fish category (i.e., hake, cod, perch, sardine, gilthead sea bream, red mullet, sole, megrim, sea bass, anchovy, tuna, swordfish, trout, etc.), an OFC to an alternative from the category of cartilaginous fish is suggested (i.e., blue shark, tope shark, dogfish, monkfish, skate, and ray) and vice versa. Regarding the increased mercury content in some sharks and other large species, the maximum limit imposed by the European Food Safety Authority (EFSA) for weekly mercury intake must be considered. An algorithm for the management of fish-FPIES, including alternative fish species, is proposed.

The pink salmon genome: Uncovering the genomic consequences of a two-year life cycle

Pink salmon (Oncorhynchus gorbuscha) adults are the smallest of the five Pacific salmon native to the western Pacific Ocean. Pink salmon are also the most abundant of these species and account for a large proportion of the commercial value of the salmon fishery worldwide. A two-year life history of pink salmon generates temporally isolated populations that spawn either in even-years or odd-years. To uncover the influence of this genetic isolation, reference genome assemblies were generated for each year-class and whole genome re-sequencing data was collected from salmon of both year-classes. The salmon were sampled from six Canadian rivers and one Japanese river. At multiple centromeres we identified peaks of Fst between year-classes that were millions of base-pairs long. The largest Fst peak was also associated with a million base-pair chromosomal polymorphism found in the odd-year genome near a centromere. These Fst peaks may be the result of a centromere drive or a combination of reduced recombination and genetic drift, and they could influence speciation. Other regions of the genome influenced by odd-year and even-year temporal isolation and tentatively under selection were mostly associated with genes related to immune function, organ development/maintenance, and behaviour.

Molecular characterization, expression and function analysis of Epinephelus coioides MKK4 response to SGIV and Vibrio alginolyticus infection

Mitogen-activated protein kinase 4 (MKK4), a member of the MAP kinase family, play important roles in response to many environmental and cellular stresses in mammals. In this study, three MKK4 subtypes, EcMKK4-1, EcMKK4-2 and EcMKK4-3, were obtained from grouper Epinephelus coioides. The open reading frame (ORF) of EcMKK4s are obtained and the EcMKK4s proteins contain highly conserved domains: a S_TKc domain, a canonical diphosphorylation group and two conserved MKKK ATP binding motifs, Asp-Phe-Gly (DFG) and Ala-Pro-Glu (APE). EcMKK4s could be found both in the cytoplasmic and nuclear. The EcMKK4s mRNA were detected in all E. coioides tissues examined with the different expression levels, and the expression were up-regulated during SGIV (Singapore grouper iridescent virus) or Vibrio alginolyticus infection. EcMKK4 could significantly reduce the activation of AP-1 reporter gene. The results suggested that EcMKK4s might play important roles in pathogen-caused inflammation.

Molecular cloning, inducible expression with SGIV and Vibrio alginolyticus challenge, and function analysis of Epinephelus coioides PDCD4

Programmed cell death 4 (PDCD4) in mammals, a gene closely associated with apoptosis, is involved in many biological processes, such as cell aging, differentiation, regulation of cell cycle, and inflammatory response. In this study, grouper Epinephelus coioides PDCD4, EcPDCD4-1 and EcPDCD4-2, were obtained. The open reading frame (ORF) of EcPDCD4-1 is 1413 bp encoding 470 amino acids with a molecular mass of 52.39 kDa and a theoretical pI of 5.33. The ORF of EcPDCD4-2 is 1410 bp encoding 469 amino acids with a molecular mass of 52.29 kDa and a theoretical pI of 5.29. Both EcPDCD4-1 and EcPDCD4-2 proteins contain two conserved MA3 domains, and their mRNA were detected in all eight tissues of E. coioides by quantitative real-time PCR (qRT-PCR) with the highest expression in liver. The expressions of two EcPDCD4s were significantly up-regulated after Singapore grouper iridovirus (SGIV) or Vibrio alginolyticus infection. In addition, over-expression of EcPDCD4-1 or EcPDCD4-2 can inhibit the activity of the nuclear factor-κB (NF-κB) and activator protein-1 (AP-1), and regulate SGIV-induced apoptosis. The results demonstrated that EcPDCD4s might play important roles in E. coioides tissues during pathogen-caused inflammation.

The evolution and functional characterization of CXC chemokines and receptors in lamprey

Chemokines are a large family of soluble peptides guiding cell migration in development and immune defense. They interact with chemokine receptors and are essential for the coordination of cell migration in diverse physiological processes. The CXC subfamily is one of the largest groups in the chemokine family and consists of multiple members. In this study, we identified homologues of three chemokine ligands (CXCL8, CXCL_F5 and CXCL12) and two CXC receptor like molecules (CXCR_L1 and CXCR_L2) in lamprey. Sequence analysis revealed that they share the same genomic organization with their counterparts in jawed vertebrates but synteny was not conserved. Lamprey CXCL8 and CXCL12 have four conserved cysteine residues whilst the CXCL_F5 has two additional cysteine residues. In addition, CXCL_F5 is evolutionarily related to the fish specific CXC chemokine groups previously identified and contains multiple cationic aa residues in the extended C- terminal region. The two CXCRs possess seven transmembrane domains and conserved structural elements for receptor activation and signaling, including the DRYXXI(V)Y motif in TM2, the disulphide bond connecting ECL2 and TM3, the WXP motif in TM6 and NPXXY motif in TM7. The identified CXC chemokines and receptors were constitutively expressed in tissues including the liver, kidney, intestine, heart, gills, supraneural body and primary leukocytes, but exhibited distinct expression patterns. Relatively high expression was detected in the gills for CXCL8, CXCL_F5 and CXCR_L1 and in the supraneural body for CXCL12 and CXCR_L2. All the genes except CXCL12 were upregulated by stimulation with LPS, pokeweed and bacterial infection, and the CXCL8 and CXCL_F5 was induced by poly (I:C). Functional analysis showed that the CXCL8 and CXCL_F5 specifically interacted with CXCR_L1 and CXCR_L2, respectively. Our results demonstrate that the CXC chemokine system had diversified in jawless fish.

Grass carp Mre11A activates IFN 1 response by targeting STING to defend against GCRV infection

Mre11A is considered as a cytosolic DNA receptor in mammals. However, it is rarely known about Mre11A in other vertebrates. Recently, a mammalian ortholog of Mre11A has been identified in grass carp (Ctenopharyngodon idellus) in our lab. Phylogenetic-tree analysis provided evidence for a close genetic relationship between C.idellus Mre11A and Carassius auratus Mre11A. The tissue expression profile of CiMre11A was detected, with a relatively higher level of expression in kidney, intestines, liver and spleen than that in other tissues after grass carp reovirus (GCRV) infection. Similarly, CiMre11A was also up-regulated in CIK cells after treatment with GCRV. Q-PCR and dual-luciferase assays indicated that the transcription levels of IFN1 and ISG15 were inhibited by CiMre11A knockdown, but were gradually augmented after CIK cells were transfected with increasing amounts of CiMre11A. Subcellular localization assays showed that a part of CiMre11A was translocated from the nucleus to the cytoplasm. Co-immunoprecipitation and co-localization assays demonstrated that CiMre11A interacts with CiSTING in response to GCRV infection. In CIK cells, the expressions of both IFN1 and ISG15 were acutely up-regulated by CiMre11A overexpression, as well as by co-overexpression of CiMre11A and CiSTING. CiMre11A and CiSTING induced the phosphorylation and cytoplasmic-to-nuclear translocation of IRF7 in CIK cells. The multiplication of GCRV in CIK cells was inhibited by the overexpression of CiMre11A and CiSTING.

Glutathione-S-transferase alpha-4 in Hippocampus abdominalis (big-belly seahorse): Molecular characterization, antioxidant properties, and its potent immune response

Glutathione-S-transferase (GST) is a key enzyme in the phase-II detoxification process and is a biomarker of oxidative stress. In this study, we analyzed the molecular, biochemical, and antioxidant properties of GST alpha-4 from Hippocampus abdominalis (HaGSTA-4). Also, the spatial and temporal expression of HaGSTA-4 upon immune challenge with abiotic and biotic stimulants were evaluated. The HaGSTA-4 ORF encodes 223 amino acids with a molecular weight of 25.7 kDa, and an estimated isoelectric point (pI) of 8.47. It consists of the GST_C superfamily and thioredoxin-like superfamily domain. The phylogenetic tree revealed that HaGSTA-4 is evolutionarily conserved with its GST alpha class counterparts. From pairwise alignment, the highest values of identity (78.5%) and similarity (85.7%) were with Parambassis ranga GSTA-4. Protein rHaGSTA-4 exhibited the highest conjugation activity towards 1-chloro-2,4-dinitrobenzene (CDNB) at pH 7 and 20 °C. A disk diffusion assay showed that rHaGSTA-4 significantly protects cells from the stress of exposure to ROS inducers such as CuSO₄, CdCl₂, and ZnCl₂. Furthermore, overexpressed HaGSTA-4 defended cells against oxidative stress caused by H₂O₂; evidence of selenium-independent peroxidase activity. From qPCR, the tissue-specific expression profile demonstrates that HaGSTA-4 is most highly expressed in the kidney, followed by the intestine and stomach, among fourteen different tissues extracted from healthy seahorses. The mRNA expression profile of HaGSTA-4 upon immune challenge varied depending on the tissue and the time after challenge. Altogether, this study suggests that HaGSTA-4 may be involved in protection against oxidative stress, in immune defense regulation, and xenobiotic metabolism.

A Comprehensive Integrated Genetic Map of the Complete Karyotype of Solea senegalensis (Kaup 1858)

Solea senegalensis aquaculture production has experienced a great increase in the last decade and, consequently, the genome knowledge of the species is gaining attention. In this sense, obtaining a high-density genome mapping of the species could offer clues to the aquaculture improvement in those aspects not resolved so far. In the present article, a review and new processed data have allowed to obtain a high-density BAC-based cytogenetic map of S. senegalensis beside the analysis of the sequences of such BAC clones to achieve integrative data. A total of 93 BAC clones were used to localize the chromosome complement of the species and 588 genes were annotated, thus almost reaching the 2.5% of the S. senegalensis genome sequences. As a result, important data about its genome organization and evolution were obtained, such as the lesser gene density of the large metacentric pair compared with the other metacentric chromosomes, which supports the theory of a sex proto-chromosome pair. In addition, chromosomes with a high number of linked genes that are conserved, even in distant species, were detected. This kind of result widens the knowledge of this species’ chromosome dynamics and evolution.

Fingerprinting of hatchery haplotypes and acquisition of genetic information by whole-mitogenome sequencing of masu salmon, Oncorhynchus masou masou, in the Kase River system, Japan

Stocking hatchery fish can lead to disturbance and extinction of the local indigenous population. Masu salmon Oncorhynchus masou masou, which is endemic across Japan, is a commonly stocked fish for recreational fishing in Japan. To conserve the indigenous resource, their genetic information is required, however, especially on Kyushu Island, the paucity of genetic information for this species has hindered proper resource management. Here, to identify hatchery mitogenome haplotypes of this species, stocked in the Kase River system, Kyushu Island, Japan, and to provide mitogenomic information for the resource management of this species, we analyzed the whole-mitogenome of masu salmon in this river system and several hatcheries potentially used for stocking. Whole-mitogenome sequencing clearly identified hatchery haplotypes, like fingerprints: among the 21 whole-mitogenome haplotypes obtained, six were determined to be hatchery haplotypes. These hatchery haplotypes were distributed in 13 out of 17 sites, suggesting that informal stocking of O. m. masou has been performed widely across this river system. The population of no hatchery haplotypes mainly belonged to clade I, a clade not found in Hokkaido Island in previous studies. Sites without hatchery haplotypes, and the non-hatchery haplotypes in clade I might be candidates for conservation as putative indigenous resources. The whole-mitogenome haplotype analysis also clarified that the same reared strain was used in multiple hatcheries. Analysis of molecular variance suggested that stocked hatchery haplotypes reduce the genetic variation among populations in this river system. It will be necessary to pay attention to genetic fluctuations so that the resources of this river system will not deteriorate further. The single nucleotide polymorphism data obtained here could be used for resource management in this and other rivers: e.g., for monitoring of informal stocking and stocked hatchery fishes, and/or putative indigenous resources.

The Evolution of Duplicated Genes of the Cpi-17/Phi-1 (ppp1r14) Family of Protein Phosphatase 1 Inhibitors in Teleosts

The Cpi-17 (ppp1r14) gene family is an evolutionarily conserved, vertebrate specific group of protein phosphatase 1 (PP1) inhibitors. When phosphorylated, Cpi-17 is a potent inhibitor of myosin phosphatase (MP), a holoenzyme complex of the regulatory subunit Mypt1 and the catalytic subunit PP1. Myosin phosphatase dephosphorylates the regulatory myosin light chain (Mlc2) and promotes actomyosin relaxation, which in turn, regulates numerous cellular processes including smooth muscle contraction, cytokinesis, cell motility, and tumor cell invasion. We analyzed zebrafish homologs of the Cpi-17 family, to better understand the mechanisms of myosin phosphatase regulation. We found single homologs of both Kepi (ppp1r14c) and Gbpi (ppp1r14d) in silico, but we detected no expression of these genes during early embryonic development. Cpi-17 (ppp1r14a) and Phi-1 (ppp1r14b) each had two duplicate paralogs, (ppp1r14aa and ppp1r14ab) and (ppp1r14ba and ppp1r14bb), which were each expressed during early development. The spatial expression pattern of these genes has diverged, with ppp1r14aa and ppp1r14bb expressed primarily in smooth muscle and skeletal muscle, respectively, while ppp1r14ab and ppp1r14ba are primarily expressed in neural tissue. We observed that, in in vitro and heterologous cellular systems, the Cpi-17 paralogs both acted as potent myosin phosphatase inhibitors, and were indistinguishable from one another. In contrast, the two Phi-1 paralogs displayed weak myosin phosphatase inhibitory activity in vitro, and did not alter myosin phosphorylation in cells. Through deletion and chimeric analysis, we identified that the difference in specificity for myosin phosphatase between Cpi-17 and Phi-1 was encoded by the highly conserved PHIN (phosphatase holoenzyme inhibitory) domain, and not the more divergent N- and C- termini. We also showed that either Cpi-17 paralog can rescue the knockdown phenotype, but neither Phi-1 paralog could do so. Thus, we provide new evidence about the biochemical and developmental distinctions of the zebrafish Cpi-17 protein family.

Black carp RIPK1 negatively regulates MAVS-mediated antiviral signaling during the innate immune activation

Receptor-interacting serine/threonine protein kinase 1 (RIPK1) is an important regulator of necroptosis and involved in innate immune response in human and mammal; however, its function in teleost fish mains largely unknown. In this paper, the RIPK1 homologue of black carp (Mylopharyngodon piceus) has been cloned and characterized to explore its role in immunity. Black carp RIPK1 (bcRIPK1) possesses the similar structure to its mammalian counterpart, which has been identified as a cytosolic protein by immunofluorescence staining. Overexpressed bcRIPK1 in host cells led to the decreased transcription of interferon (IFN) and interferon stimulated genes, and exogenous bcRIPK1 in EPC cells led to the decreased transcription of interferon promoters in reporter assay. Our previous study has identified that black carp MAVS (bcMAVS) functions as an antiviral adaptor protein against both grass carp reovirus (GCRV) and spring viremia of carp virus (SVCV). The reporter assay showed that the IFN-inducing ability of bcMAVS was dampened by bcRIPK1 and the plaque assay demonstrated that the antiviral activity of bcMAVS was inhibited by bcRIPK1. The immunofluorescent staining and co-immunoprecipitation identified the interaction between these two molecules. Thus, the data generated in this paper support the conclusion that bcRIPK1 interacts with bcMAVS and negatively regulates bcMAVS-mediated antiviral signaling.

Effects of Yersinia ruckeri invasion on the proteome of the Chinook salmon cell line CHSE-214

Yersinia ruckeri is an important bacterial pathogen of fish, in particular salmonids, it has been associated with systemic infections worldwide and, like many enteric bacteria, it is a facultative intracellular pathogen. However, the effect of Y. ruckeri's interactions with the host at the cellular level have received little investigation. In the present study, a culture of Chinook Salmon Embryo (CHSE) cell line was exposed to Y. ruckeri. Afterwards, the proteins were investigated and identified by mass spectrometry and compared to the content of unexposed cultures. The results of this comparison showed that 4.7% of the identified proteins were found at significantly altered concentrations following infection. Interestingly, infection with Y. ruckeri was associated with significant changes in the concentration of surface adhesion proteins, including a significantly decreased presence of β-integrins. These surface adhesion molecules are known to be the target for several adhesion molecules of Yersiniaceae. The concentration of several anti-apoptotic regulators (HSP90 and two DNAj molecules) appeared similarly downregulated. Taken together, these findings suggest that Y. ruckeri affects the proteome of infected cells in a notable manner and our results shed some light on the interaction between this important bacterial pathogen and its host.

Identification, expression profile and analysis of the antimicrobial activity of collectin 11 (CL-11, CL-K1), a novel complement-associated pattern recognition molecule, in half-smooth tongue sole (Cynoglossus semilaevis)

Collectin 11 (CL-11, also known as collectin kidney-1, CL-K1), a new member of the vertebrate C-type lectin superfamily, plays an important role in innate immunity as a pattern recognition molecule of the lectin complement pathway. However, little is known about CL-11 in teleosts. In the present study, a CL-11 homolog was identified and characterized from half-smooth tongue sole (Cynoglossus semilaevis) (designated as CsCL-11). The full-length cDNA of CsCL-11 is 1220 bp long and includes a 5'untranslated region (5'-UTR) of 180 bp, a 3'-UTR of 218 bp and an open reading frame (ORF) of 819 bp encoding 273 amino acids. Multiple sequence alignment revealed that the deduced CsCL-11 protein has the typical modular architecture (EPN and WTD) conserved throughout vertebrates, suggesting a conserved function of CsCL-11. Tissue expression profile analysis by quantitative real-time PCR (qRT-PCR) showed CsCL-11 to be ubiquitously distributed in tissues and highly expressed in the ovary and liver. A pattern of significant upregulation of CsCL-11 expression was observed in the blood, spleen, head kidney and gill at 6 h, 12 h and 24 h after infection with Vibrio anguillarum, and western blotting showed that natural CsCL-11 protein levels in the blood were significantly increased after V. anguillarum infection. Moreover, by binding to various bacteria, recombinant CsCL-11 (rCsCL-11) expressed in HEK-293 T cells displayed strong antibacterial activity. Taken together, these results suggest that CsCL-11 is a unique C-type lectin that is likely involved in host defense against bacterial infection. To our knowledge, this is the first study on CL-11 in marine fish.

Molecular characterization of novel mitochondrial peroxiredoxins from the Antarctic emerald rockcod and their gene expression in response to environmental warming

In the present study we describe the molecular characterization of the two paralogous mitochondrial peroxiredoxins from Trematomus bernacchii, a teleost that plays a pivotal role in the Antarctic food chain. The two putative amino acid sequences were compared with orthologs from other fish, highlighting a high percentage of identity and similarity with the respective variant, in particular for the residues that are essential for the characteristic peroxidase activity of these enzymes. The temporal expression of Prdx3 and Prdx5 mRNAs in response to short-term thermal stress showed a general upregulation of prdx3, suggesting that this isoform is the most affected by temperature increase. These data, together with the peculiar differences between the molecular structures of the two mitochondrial Prdxs in T. bernacchii as well as in the tropical species Stegastes partitus, suggest an adaptation that allowed these poikilothermic aquatic vertebrates to colonize very different environments, characterized by different temperature ranges.

Liver-Specific microRNA Identification in Farmed Carp, Labeo bata (Hamilton, 1822), Fed with Starch Diet Using High-Throughput Sequencing

The liver is an important central organ, which controls carbohydrate metabolism through maintaining glucose homeostasis by a tightly regulated system of genes or enzymes. The microRNAs are small non-coding RNAs playing an important role in the regulation of genes associated with developmental biology, physiology, metabolism, etc. Thus, in this study, we have intended to detect liver-specific microRNAs in farmed carp, Labeo bata, upon being fed a diet with different levels of carbohydrates. Here, we have conducted the experiment for 45 days using fingerlings of farmed carp fed with 20% (control), 40%, and 60% gelatinized starch levels. The liver tissues were collected from each treatment and processed for RNA isolation, small RNA library preparation, and high-throughput sequencing using Illumina NexSeq500. Through sequencing, 15,779,417 reads in 20% CHO, 13,959,039 in 40% CHO, and 13,661,950 in 60% CHO reads were generated for control and treated fishes using three small RNA libraries. We have investigated 445 novel and 231 conserved microRNAs in 20%, 40%, and 60% carbohydrate (CHO), respectively, through computational analysis. The differential expression analysis of miRNAs was carried out between different treatments compared with control and this study depicted 117 known and 114 novel miRNA genes involved in carbohydrate metabolic pathways. Further, target prediction and gene ontology analysis revealed that miRNAs were involved in several pathways such as signaling pathway, G protein pathway, complement receptor-mediated pathway, dopamine receptor signaling pathway, epidermal growth factor pathway, and notch signaling pathway. The predicted miRNA sites in targeted genes were associated with cellular activities, developmental biology, DNA binding, Golgi apparatus, extracellular region, catalytic activity, MAPK cascade, etc. Overall, we have generated a vital resource of liver-specific miRNAs involved in metabolic gene regulation. These studies further will help develop miRNA inhibitors to study their role during carbohydrate metabolism in farmed carp.

Genome-Wide Association Study of Growth and Body-Shape-Related Traits in Large Yellow Croaker (Larimichthys crocea) Using ddRAD Sequencing

Large yellow croaker (Larimichthys crocea) is an economically important marine fish species of China. Due to overfishing and marine pollution, the wild stocks of this croaker have collapsed in the past decades. Meanwhile, the cultured croaker is facing the difficulties of reduced genetic diversity and low growth rate. To explore the molecular markers related to the growth traits of croaker and providing the related SNPs for the marker-assisted selection, we used double-digest restriction-site associated DNA (ddRAD) sequencing to dissect the genetic bases of growth traits in a cultured population and identify the SNPs that associated with important growth traits by GWAS. A total of 220 individuals were genotyped by ddRAD sequencing. After quality control, 27,227 SNPs were identified in 220 samples and used for GWAS analysis. We identified 13 genome-wide significant associated SNPs of growth traits on 8 chromosomes, and the beta P of these SNPs ranged from 0.01 to 0.86. Through the definition of candidate regions and gene annotation, candidate genes related to growth were identified, including important regulators such as fgf18, fgf1, nr3c1, cyp8b1, fabp2, cyp2r1, ppara, and ccm2l. We also identified SNPs and candidate genes that significantly associated with body shape, including bmp7, col1a1, col11a2, and col18a1, which are also economically important traits for large yellow croaker aquaculture. The results provided insights into the genetic basis of growth and body shape in large yellow croaker population and would provide reliable genetic markers for molecular marker-assisted selection in the future. Meanwhile, the result established a basis for our subsequent fine mapping and related gene study.

The RAG2 gene of yellow catfish (Tachysurus fulvidraco) and its immune response against Edwardsiella ictaluri infection

Recombination-activating gene 2 (rag 2) allies with recombination-activating gene 1 (rag 1) and regulates the V(D)J recombination of immunoglobulin (Ig) and T-cell receptor (TCR) genes. Being a key player in the adaptive immune response of vertebrates, functional characterization of rag 2 from yellow catfish is beneficial for understanding the biological response towards the pathogens. In this report, we have cloned and characterized the rag 2 gene of yellow catfish, and a particular pattern of expression was analysed in the major tissues of yellow catfish. The results showed that the open reading frame (ORF) of yellow catfish rag 2 was 1596 bp in length, which encodes a peptide of 531 amino acids. The multiple sequence alignment and phylogenetic analysis of rag 2 of yellow catfish with other species showed the conserved regions and the classical taxonomic evolution among the different vertebrate species. The qRT-PCR and Western blot results revealed that rag 2 transcripts and proteins were present in various tissues of yellow catfish with relatively high expression in the tissues of the thymus, head-kidney, and spleen. The systematic distribution analysis of the rag 2 expression by immunohistochemistry (IHC) using the rabbit polyclonal antibody, exposed relatively high expression in head kidney, spleen and thymus tissues after infected with Edwardsiella ictaluri. Moreover, the temporal expression of rag 2 and pro-inflammatory cytokines (IL-1β and TNF-α) were significantly upregulated at different time points in the specific lymphoid tissues of yellow catfish following E. ictaluri infection. Our findings suggest that rag 2 potentially exhibited the immunological response in primary lymphoid tissues of yellow catfish against bacterial infection. This study will provide an essential source about rag 2 gene and its relationship with the inflammatory cytokines during infection.

Full characterization and transcript expression profiling of the interferon regulatory factor (IRF) gene family in Atlantic cod (Gadus morhua)

Atlantic cod (Gadus morhua) represents a unique immune system among teleost fish, making it a species of interest for immunological studies, and especially for investigating the evolutionary history of immune gene families. The interferon regulatory factor (IRF) gene family encodes transcription factors which function in the interferon pathway, but also in areas including leukocyte differentiation, cell growth, autoimmunity, and development. We previously characterized several IRF family members in Atlantic cod (Irf4a, Irf4b, Irf7, Irf8, and two Irf10 splice variants) at the cDNA and putative amino acid levels, and in the current study we took advantage of the new and improved Atlantic cod genome assembly in combination with rapid amplification of cDNA ends (RACE) to characterize the remaining family members (i.e. Irf3, Irf5, Irf6, Irf9, and two Irf2 splice variants). Real-time quantitative PCR (QPCR) was used to investigate constitutive expression of all IRF transcripts during embryonic development, suggesting several putative maternal transcripts, and potential stage-specific roles. QPCR studies also showed 11 of 13 transcripts were responsive to stimulation with poly(I:C), while 6 of 13 transcripts were responsive to lipopolysaccharide (LPS) in Atlantic cod head kidney macrophages, indicating roles for cod IRF family members in both antiviral and antibacterial responses. This study is the first to investigate expression of the complete IRF family in Atlantic cod, and suggests potential novel roles for several of these transcription factors within immunity as well as in early development of this species.

FOXO genes in channel catfish and their response after bacterial infection

FOXO proteins are a subgroup of the forkhead family of transcription factors that play crucial roles in lifespan regulation. In addition, FOXO proteins are also involved in immune responses. After a systematic study of FOXO genes in channel catfish, Ictalurus punctatus, seven FOXO genes were identified and characterized, including FOXO1a, FOXO1b, FOXO3a, FOXO3b, FOXO4, FOXO6a and FOXO6b. Through phylogenetic and syntenic analyses, it was found that FOXO1, FOXO3 and FOXO6 were duplicated in the catfish genome, as in the zebrafish genome. Analysis of the relative rates of nonsynonymous (dN) and synonymous (dS) substitutions revealed that the FOXO genes were globally strongly constrained by negative selection. Differential expression patterns were observed in the majority of FOXO genes after Edwardsiella ictaluri and Flavobacterium columnare infections. After E. ictaluri infection, four FOXO genes with orthologs in mammal species were significantly upregulated, where FOXO6b was the most dramatically upregulated. However, after F. columnare infection, the expression levels of almost all FOXO genes were not significantly affected. These results suggested that either a pathogenesis-specific pattern or tissue-specific pattern existed in catfish after these two bacterial infections. Taken together, these findings indicated that FOXO genes may play important roles in immune responses to bacterial infections in catfish.

Differential Transcriptomic and Metabolomic Responses in the Liver of Nile Tilapia (Oreochromis niloticus) Exposed to Acute Ammonia

Ammonia is toxic to aquatic animal. Currently, only limited works were reported on the responses of aquatic animals after ammonia exposure using "omics" technologies. Tilapia suffers from the stress of ammonia-nitrogen during intensive recirculating aquaculture. Optimizing ammonia stress tolerance has become an important issue in tilapia breeding. The molecular and biochemical mechanisms of ammonia-nitrogen toxicity have not been understood comprehensively in tilapia yet. In this study, using RNA-seq and gas chromatograph system coupled with a Pegasus HT time-of-flight mass spectrometer (GC-TOF-MS) techniques, we investigated differential expressed genes (DEGs) and metabolomes in the liver at 6 h post-challenges (6 hpc) and 24 h post-challenges (24 hpc) under high concentration of ammonia-nitrogen treatment. We detected 2258 DEGs at 6 hpc and 315 DEGs at 24 hpc. Functional enrichment analysis indicated that DEGs were significantly associated with cholesterol biosynthesis, steroid and lipid metabolism, energy conservation, and mitochondrial tissue organization. Metabolomic analysis detected 31 and 36 metabolites showing significant responses to ammonia-nitrogen stress at 6 and 24 hpc, respectively. D-(Glycerol 1-phosphate), fumaric acid, and L-malic acid were found significantly down-regulated at both 6 and 24 hpc. The integrative analysis of transcriptomics and metabolomics suggested considerable alterations and precise control of gene expression at both physiological and molecular levels in response to the stress of ammonia-nitrogen in tilapia.

Comparative Transcriptome Analysis Reveals Differentially Expressed Genes and Signaling Pathways Between Male and Female Red-Tail Catfish (Mystus wyckioides)

Sexual dimorphism is widespread in fish species. The red-tail catfish (Mystus wyckioides) is a commercially important catfish in the lower reaches of the Lancang River and the Mekong basin, and it shows a growth advantage in males. Here, RNA-seq was for the first time used to explore the gene expression difference between the sexes in the hypothalamus and pituitary of red-tail catfish, respectively. In the hypothalamus, 5732 and 271 unigenes have significantly higher and lower expressions, respectively, in males compared with females. KEGG analysis showed that 212 DEGs were annotated to 216 signaling pathways, and enrichment analysis suggested different levels of cAMP and glutamatergic synapse signaling between male and female hypothalami and some of the DEGs appear involved in gonad development and growth. In the pituitary, we found only 19 differentially expressed unigenes, which were annotated to 32 signaling pathways, most of which play important roles in gonad development.

Receptor complex and signalling pathway of the two type II IFNs, IFN-γ and IFN-γrel in mandarin fish or the so-called Chinese perch Siniperca chuatsi

IFN-γ, as the sole member of mammalian type II IFN, is a multifunctional cytokine which exerts its effects through two distinct IFN-γ receptors, IFNGR1 and IFNGR2. However, in teleost fish, another IFN-γ homologous gene, namely IFN-γ related gene (IFN-γrel), has been identified. Although IFN-γ and IFN-γrel genes have been described in some fish species, many important aspects remain poorly understood in relation with their signalling and function. In the present study, IFN-γ and IFN-γrel, as well as their receptors, cytokine receptor family B (CRFB) 17, CRFB13, two of which are homologous to IFNGR1 in mammals, and CRFB6, homolomous to IFNGR2, have been characterized in mandarin fish, Siniperca chuatsi. It was revealed that the two type IFN members exhibit antiviral activity, and IFN-γ transduces downstream signalling through CRFB13 and CRFB6, while IFN-γrel interacts with CRFB17 to activate downstream signalling. Moreover, IFN-γ and IFN-γrel have been shown to exert antiviral biological activity in a STAT1-dependent manner. Intracellular domain analysis of CRFB17 and CRFB13 demonstrated that the Y386 tyrosine residue of CRFB13 is required for the activation of the IFN-γ-mediated biologic response, and the Y324 and Y370 residues in CRFB17 are required to activate IFN-γrel signalling.

Transcriptomic Analysis Revealed the Regulatory Mechanisms of Oocyte Maturation and Hydration in Orange-Spotted Grouper (Epinephelus coioides)

Oocyte maturation and hydration are regulated by a complex interplay of various hormones and local factors. We have investigated the morphological changes of follicles and serum steroid levels during the HCG (human choionic gonadotophin)-induced oocyte maturation in the orange-spotted grouper. For the first time, a large-scale transcriptomic analysis of follicles during the maturation has been conducted in a fish species which produce pelagic oocytes. Eight cDNA libraries of follicle samples, from full-grown immature follicles to mature follicles, were constructed. A total of 402,530,284 high-quality clean reads were obtained after filtering, 79.66% of which perfectly mapped to the orange-spotted grouper genome. Real-time PCR results of 12 representative genes related to oocyte maturation and hydration verified the reliability of the RNA-seq data. A large number of genes related to oocyte maturation and hydration were identified in the transcriptome dataset. And the transcriptomic analysis revealed the dynamic changes of the steroid synthesis pathway and the pathway of hydration during oocyte maturation. The present study will facilitate future study on the oocyte maturation and hydration in the orange-spotted grouper and other marine pelagic egg spawner.

Dynamics of MiRNA Transcriptome in Turbot (Scophthalmus maximus L.) Intestine Following Vibrio anguillarum Infection

MicroRNAs (miRNAs) are a group of small non-coding RNAs, which could bind to the 3'-untranslated regions of their target mRNAs to regulate gene expression in various biological processes, including immune-regulated signaling pathways. Turbot (Scophthalmus maximus L.), an important commercial fish species in China, has been suffering with Vibrio anguillarum infection resulted in dramatic economic loss. Therefore, we investigated the expression profiles of miRNAs, as well as the immune-related miRNA-mRNA pairs in turbot intestine at 1 h, 4 h, and 12 h following V. anguillarum infection. As a result, 266 predicted novel miRNAs and 283 conserved miRNAs belonging to 92 miRNA families were detected. A total of 44 miRNAs were differentially expressed in the intestine following V. anguillarum infection. Following prediction, the potential target genes of differentially expressed miRNAs were grouped into a wide range of functional categories, including immune defense/evasion, inflammatory responses, RIG-I signaling pathway, and Toll-like receptor signaling pathway. Moreover, we selected 15 differentially expressed immune genes and their related differentially expressed miRNAs to construct an interaction network for V. anguillarum infection in turbot. These results suggested that in teleost, as in higher vertebrates, miRNAs prominently contribute to immune responses, protecting the host against infection. In addition, this is the first report of comprehensive identification of turbot miRNAs being differentially regulated in the intestine related to V. anguillarum infection. Our results provided an opportunity for further understanding of the molecular mechanisms of miRNA regulation in turbot host-pathogen interactions.

Functional characterization of interleukin (IL)-22 and its inhibitor, IL-22 binding protein (IL-22BP) in Mandarin fish, Siniperca chuatsi

As an important immune regulatory molecule, interleukin (IL)-22 has been reported in several species of fish, but its soluble receptor, IL-22 binding protein (IL-22BP), discovered as a natural antagonist of IL-22 in mammals, has not been functionally characterized in fish to date. In the present study, IL-22 and IL-22BP genes were cloned in mandarin fish Siniperca chuatsi. They all exhibited a high basal expression level in mucosa-enriched tissues, implying their possible roles in mucosal immunity. The IL-22 was found to show a potent response to LPS stimulation, acting as an inducer of antimicrobial peptide (AMP) genes, such as hepcidin and Liver-expressed antimicrobial peptide-2 (LEAP-2) in intestinal cells. IL-22BP, via co-incubation with IL-22, inhibited completely the induction of downstream genes by IL-22. Through a yeast two-hybrid assay, the interaction between IL-22BP and IL-22 was confirmed, which may account for the inhibitory effect of IL-22BP. Moreover, two hot spot residues for IL-22 binding, as reported in mammalian IL-22BP, were found to be conserved both in sequence location and function in mandarin fish IL-22BP, indicating that the interaction mode between IL-22 and IL-22BP may be also conserved in fish and mammals. In conclusion, the mandarin fish IL-22 and IL-22BP are conserved in their interaction and function with their mammalian orthologues, and these findings provide basis for future research on IL-22-IL-22BP axis in fish immunity.

Grass carp (Ctenopharyngodon idella) RSK2 protects cells anti-apoptosis by up-regulating BCL-2

In mammals, toll-like receptor 3 (TLR3) is capable of recognizing double-stranded RNA and then initiates transcription of IFN-β. TLR3 can activate the innate immune system by phosphorylating extracellular signal-regulated kinase 1 (ERK1) in the mitogen-activated protein kinase (MAPK) pathway. As a downstream signaling protein of ERK1, ribosomal protein S6 kinase alpha 3 (RSK2) is activated through the "classical" MAPK pathway. So RSK2 plays a critical role in response to innate immune system induced by TRL3. However, the innate immune mechanism of RSK2 remains indistinct in fish. In this study, we cloned and characterized a full length cDNA sequence of RSK2 from Ctenopharyngodon idella (named CiRSK2, MH844551). The full length cDNA of CiRSK2 is 3930 bp with a coding sequence of 2202 bp encoding a polypeptide of 734 amino acids. The expression of CiRSK2 was ubiquitous and significantly up-regulated under the stimulation of poly (I:C) in eight different tissues of C. idella and C. idella kidney cells (CIK). In addition, poly (I:C) stimulation also up-regulated the expression of CiERK1 mRNA in CIK cells and the phosphorylation of CiERK1. We also demonstrated that the activated CiERK1 interacted with CiRSK2 by CO-IP assay and immunofluorescence assay. To further investigate the relationship between CiRSK2 and CiERK1, we performed subcellular localization of CiRSK2 at different periods of CiERK1 stimulation. The result showed that CiERK1 can make CiRSK2 enter the nucleus. Subsequently, we found that CiRSK2 increased the transcriptional level of CiBCL-2 and protein level of CiBCL-2 significantly. Then cell apoptosis was inhibited to a certain extent. Overall, our results suggested that CiRSK2 plays important roles in fish innate immunity and is able to inhibit cell apoptosis by up-regulating CiBCL-2.

Pushing the limits of photoreception in twilight conditions: The rod-like cone retina of the deep-sea pearlsides

Most vertebrates have a duplex retina comprising two photoreceptor types, rods for dim-light (scotopic) vision and cones for bright-light (photopic) and color vision. However, deep-sea fishes are only active in dim-light conditions; hence, most species have lost their cones in favor of a simplex retina composed exclusively of rods. Although the pearlsides, Maurolicus spp., have such a pure rod retina, their behavior is at odds with this simplex visual system. Contrary to other deep-sea fishes, pearlsides are mostly active during dusk and dawn close to the surface, where light levels are intermediate (twilight or mesopic) and require the use of both rod and cone photoreceptors. This study elucidates this paradox by demonstrating that the pearlside retina does not have rod photoreceptors only; instead, it is composed almost exclusively of transmuted cone photoreceptors. These transmuted cells combine the morphological characteristics of a rod photoreceptor with a cone opsin and a cone phototransduction cascade to form a unique photoreceptor type, a rod-like cone, specifically tuned to the light conditions of the pearlsides' habitat (blue-shifted light at mesopic intensities). Combining properties of both rods and cones into a single cell type, instead of using two photoreceptor types that do not function at their full potential under mesopic conditions, is likely to be the most efficient and economical solution to optimize visual performance. These results challenge the standing paradigm of the function and evolution of the vertebrate duplex retina and emphasize the need for a more comprehensive evaluation of visual systems in general.

A Novel igf3 Gene in Common Carp (Cyprinus carpio): Evidence for Its Role in Regulating Gonadal Development

Since the insulin-like growth factor 3 (igf3) gene was recently discovered in fish ovary, its function in the gonads has received much attention. In this study, we isolated two igf3 subtypes from common carp (Cyprinus carpio), which comprised full-length cDNA of 707 and 1153 nucleotides encoding 205 and 198 amino acids (aa), respectively. The Igf3 aa sequence had the highest gene homology of 72% with the corresponding sequence in zebrafish (Danio rerio). Phylogenetic tree construction revealed that the C. carpio igf3 gene was first clustered with D. rerio and then with other teleost species. Igf3 mRNA was widely expressed, with expression being highest in the gonads and blood. In the gonad development stage, igf3a mRNA expression was highest in the maturity and recession stage of the ovary, and decline phase of the testis, while igf3b was highest in the recession and fully mature periods of the ovaries and testes, respectively. Western blotting of testis protein samples showed two bands of approximately 21 kDa and 34 kDa corresponding to the calculated molecular mass of the two Igf3 subtypes; no signal was detected in the ovary. The Igf3 protein was localized in the ovary granulosa cells and testis spermatogonium and spermatids. 17β-Ethinylestradiol treatment increased both ovary and testis igf3 mRNA expression. These findings suggest that Igf3 may play an important role in C. carpio gonadal development.

Molecular characterization and evolution analysis of five interleukin-17 receptor genes in large yellow croaker Larimichthys crocea

Interleukin-17s (IL-17s) play critical roles in inflammatory response and host defense against extracellular pathogens. IL-17s induce the immune response signaling through the specific IL-17 receptors (IL-17Rs) that consist of five members (IL-17RA to E). In the present work, we have identified the five IL-17R orthologs (LycIL-17Rs) from large yellow croaker Larimichthys crocea. The deduced protein of each LycIL-17R exhibits a typical IL-17R domain architecture, including a signal peptide, the extracellular FNIII domain (IL-17RA/RB/RD) or IL-17_R_N domain (IL-17RC/RE), a transmembrane domain, and a SEFIR domain in cytoplasmic region. In particular, the extracellular regions of teleost IL-17RB are much shorter than those in mammals and lack an FNIII domain (FN2). Phylogenetic tree shows that IL-17Rs are classified into two main groups: IL-17RA/RB/RD group and IL-17RC/RE group, which is distinct from previous proposal that grouped IL-17RB into IL-17RC/RE. The surrounding genes of IL-17Rs are conservatively aligned in genomes between teleosts and mammals. The five LycIL-17Rs were constitutively expressed in all tissues examined, but with different expression patterns. Aeromonas hydrophila infection significantly upregulated LycIL-17RA, RC, RD and RE in both mucosal tissue (gills) and systemic immune tissues (head kidney and spleen), while the increase of LycIL-17RB expression could be detected in gills, indicating that LycIL-17Rs may be involved in host defense against bacterial infection. Thus, these results suggest that teleost IL-17Rs may function in mediating immune response as their mammalian orthologs. To our knowledge, this is the first report of molecular characterization of the five IL-17Rs (IL-17RA/RB/RD and IL-17RC/RE) in teleost fish.

Molecular characterization and expression analysis of Hsp90 in Schizothorax prenanti

Aquatic animals suffer from various environmental stresses because the aquatic environment is a very complex system. To monitor the health status of fish, Hsp90 a potential early warning marker was determined in Schizothorax prenanti after infection with a bacterium. In this study, we cloned Hsp90 from S. prenanti for the first time. The full-length cDNA sequence of SpHsp90 was 2663 bp, contains an open reading frame of 2181 bp, and has a gene encoding 726 amino acids, an estimated molecular mass of 83.38 kDa, and a theoretical isoelectric point of 4.91. The SpHsp90 amino acid sequence has five conserved HSP90 family signatures and shares 87.0-95.5 % identity with other vertebrates. Phylogenetic analysis and structure comparison indicated that SpHsp90 should be a β isoform of the HSP90 family. SpHsp90 was ubiquitously expressed in all examined tissues, and the highest level of expression was in the kidney. After Streptococcus agalactiae infection, the level of SpHsp90 expression had significant changes (P < 0.05) in the hepatopancreas, spleen, kidney, and blood. The expression increased to the highest level at 6 h in the blood and at 24 h in the hepatopancreas, spleen, and kidney. The results suggested that the SpHsp90 gene could be induced by S. agalactiae in S. prenanti and that SpHsp90 may be involved in resistance to bacterial infection and provide an early warning information. The kidney is the most suitable for detecting SpHsp90 after bacterial infection.

Molecular Characterization and Growth Association of Two Apolipoprotein A-Ib Genes in Common Carp (Cyprinus carpio)

Apolipoprotein A-I (ApoA-I) is functionally involved in the transportation and metabolism of lipids in vertebrates. In this study, two isoforms of apoA-Ib in common carp (Cyprinus carpio L.) were characterized. Sequence comparison and phylogenetic analysis showed that C. carpio ApoA-Ib is relatively conserved within cyprinid fishes. During embryonic development, C. carpioapoA-Ib was first expressed at the stage of multi-cells, and the highest mRNA level was observed at the stage of optic vesicle. A ubiquitous expression pattern was detected in various tissues with extreme predominance in the liver. Significantly different expression levels were observed between light and heavy body weight groups and also in the compensatory growth test. Seventeen and eight single-nucleotide polymorphisms (SNPs) were identified in matured mRNA of the C. carpioapoA-Ib.1 and apoA-Ib.2, respectively. Two of these SNPs (apoA-Ib.2-g.183A>T and apoA-Ib.2-g.1753C>T) were significantly associated with body weight and body length in two populations of common carp. These results indicate that apoA-Ib may play an important role in the modulation of growth and development in common carp.

Genome-Wide Identification, Characterization and Phylogenetic Analysis of ATP-Binding Cassette (ABC) Transporter Genes in Common Carp (Cyprinus carpio)

The ATP-binding cassette (ABC) gene family is considered to be one of the largest gene families in all forms of prokaryotic and eukaryotic life. Although the ABC transporter genes have been annotated in some species, detailed information about the ABC superfamily and the evolutionary characterization of ABC genes in common carp (Cyprinus carpio) are still unclear. In this research, we identified 61 ABC transporter genes in the common carp genome. Phylogenetic analysis revealed that they could be classified into seven subfamilies, namely 11 ABCAs, six ABCBs, 19 ABCCs, eight ABCDs, two ABCEs, four ABCFs, and 11 ABCGs. Comparative analysis of the ABC genes in seven vertebrate species including common carp, showed that at least 10 common carp genes were retained from the third round of whole genome duplication, while 12 duplicated ABC genes may have come from the fourth round of whole genome duplication. Gene losses were also observed for 14 ABC genes. Expression profiles of the 61 ABC genes in six common carp tissues (brain, heart, spleen, kidney, intestine, and gill) revealed extensive functional divergence among the ABC genes. Different copies of some genes had tissue-specific expression patterns, which may indicate some gene function specialization. This study provides essential genomic resources for future studies in common carp.

Chemokine receptors in Atlantic salmon

Teleost sequence data have revealed that many immune genes have evolved differently when compared to other vertebrates. Thus, each gene family needs functional studies to define the biological role of individual members within major species groups. Chemokine receptors, being excellent markers for various leukocyte subpopulations, are one such example where studies are needed to decipher individual gene function. The unique salmonid whole genome duplication that occurred approximately 95 million years ago has provided salmonids with many additional duplicates further adding to the complexity and diversity. Here we have performed a systematic study of these receptors in Atlantic salmon with particular focus on potential inflammatory receptors. Using the preliminary salmon genome data we identified 48 chemokine or chemokine-like receptors including orthologues to the ten receptors previously published in trout. We found expressed support for 40 of the bona fide salmon receptors. Eighteen of the chemokine receptors are duplicated, and when tested against a diploid sister group the majority were shown to be remnants of the 4R whole genome duplication with subsequent high sequence identity. The salmon chemokine receptor repertoire of 40 expressed bona fide genes is comparably larger than that found in humans with 23 receptors. Diversification has been a major driving force for these duplicate genes with the main variability residing in ligand binding and signalling domains.

Ctenopharyngodon idella IRF2 plays an antagonistic role to IRF1 in transcriptional regulation of IFN and ISG genes

Interferon Regulatory Factors (IRFs) make up a family of transcription factors involved in transcriptional regulation of type I IFN and IFN-stimulated genes (ISG) in cells. In the present study, an IRF2 gene (termed CiIRF2, JX628585) was cloned and characterized from grass carp (Ctenopharyngodon idella). The full-length cDNA of CiIRF2 is 1809 bp in length, with the largest open reading frame (ORF) of 981 bp encoding a putative protein of 326 amino acids. CiIRF2 contains a conserved DNA-binding domain (DBD) in N-terminal and a non-conserved C-terminal region. Protein sequence analysis revealed that CiIRF2 shares significant homology to the known IRF2 counterparts. Phylogenetic reconstruction confirmed its closer evolutionary relationship with other fish counterparts, especially with zebra fish IRF2. CiIRF2 was ubiquitously expressed at low level in all tested grass carp tissues and significantly up-regulated except in brain following poly I:C 6-12 h post stimulation. In order to understand fish innate immune and resistance to virus diseases, recombinant CiIRF2 with His-tag was over-expressed in BL21 Escherichia coli, and the expressed protein was purified by affinity chromatography with Ni-NTA His-Bind Resin. Promoter sequences of grass carp type I IFN gene (CiIFN) and two ISG genes (CiPKR and CiPKZ) were amplified and cloned. In vitro, gel mobility shift assays were employed to analyze the interaction of CiIRF2 protein with promoters of CiIFN, CiPKR and CiPKZ respectively. The results showed that CiIRF2 bound to these promoters with high affinity by means of its DBD. Afterwards, recombinant plasmids of pGL3-CiIFN, pGL3-CiPKR and pGL3-CiPKZ were constructed and transiently co-transfected with pcDNA3.1-CiIRF2 or pcDNA3.1-CiIRF1 respectively into C. idella kidney (CIK) cells. Dual-luciferase reporter assays demonstrated that CiIRF2 down-regulates the transcription activity of CiIFN, CiPKR and CiPKZ genes in CIK cells. To further understand the function of fish IRF2, expression plasmids (pcDNA3.1-IRF2 and pcDNA3.1-IRF1) were transiently co-transfected with pGL3-IFN or pGL3-CiPKZ into CIK cells, respectively. The results revealed that CiIRF2 plays an antagonistic role to CiIRF1 in transcriptional regulation of IFN and ISG genes.

Identification of interferon-γ-inducible-lysosomal thiol reductase (GILT) gene from Mefugu (Takifugu obscures) and its immune response to LPS challenge

Interferon-γ-inducible-lysosomal thiol reductase (GILT) plays a key role in the processing and presentation of MHC class II-restricted antigen (Ag) by catalyzing disulfide bond reduction. In this study, a Mefugu cDNA (ToGILT) encodes a deduced protein of 242 amino acids with a putative molecular weight of 28.6 kDa. It contains typical features of GILT proteins including the signature sequence CQHGX2ECX2NX4C, CXXC motif and other five cysteines. Genomic analysis revealed that ToGILT gene exhibited a similar exon-intron organization to human and mouse GILT. Phylogenetic analysis showed that ToGILT derived from a common ancestor with other vertebrate GILT proteins. The ToGILT mRNA was expressed in a tissue-specific manner and obviously up-regulated in spleen and kidney after LPS induction. These results suggest that ToGILT may be involved in the immune response to bacteria challenge in Takifugu obscurus.

Characterizations of two grass carp Ctenopharyngodon idella HMGB2 genes and potential roles in innate immunity

High-mobility group box 2 (HMGB2) protein is a chromatin-associated nonhistone protein, involved in transcriptional regulation and nucleic-acid-mediated innate immune responses in mammalian. However, the function of piscine HMGB2 in innate immune responses is still unknown. In the present study, two HMGB2 homologue genes (CiHMGB2a, CiHMGB2b) were identified and characterized in grass carp (Ctenopharyngodon idella). Both CiHMGB2a and CiHMGB2b genes encode proteins with 213 amino acids, sharing 71.4% identities and containing two basic HMG boxes and an acidic tail. The deduced protein sequences showed the most identities to HMGB2a (93%) and HMGB2b (86.4%) of zebrafish (Danio rerio), respectively. Quantitative real-time RT-PCR (qRT-PCR) analysis showed that CiHMGB2a and CiHMGB2b were constitutively expressed in all the 15 tested tissues. Post grass carp reovirus (GCRV) infection, mRNA levels of CiHMGB2a and CiHMGB2b were strongly up-regulated in spleen and head kidney and mildly modulated in C. idella kidney (CIK) cells. Meanwhile, mRNA expressions of CiHMGB2a and CiHMGB2b were significantly regulated by viral pathogen associated molecular patterns (PAMPs) polyinosinic-polycytidylic potassium salt (poly(I:C)) and bacterial PAMPs lipopolysaccharide (LPS), peptidoglycan (PGN) challenge in CIK cells. In CiHMGB2a and CiHMGB2b over-expression cells, expressions of CiHMGB2a and CiHMGB2b facilitated each other; transcription levels of CiTRIF, CiMyD88, CiIPS-1 and CiMx1 were remarkably enhanced, whereas CiIFN-I was inhibited, compared with those in cells transfected with pCMV (control plasmid); after GCRV challenge, all those tested genes were up-regulated with divergent expression profiles. Antiviral activities of CiHMGB2a and CiHMGB2b were manifested by the delayed appearance of cytopathic effect (CPE) and inhibition of GCRV yield. All those results demonstrate that CiHMGB2a and CiHMGB2b not only mediate antiviral immune responses but also involve in responding to viral/bacterial PAMPs challenge, which provides novel insights into the essential role of HMGB2 in innate immunity.

Phylogenetic analysis of vertebrate CXC chemokines reveals novel lineage specific groups in teleost fish

In this study, we have identified 421 molecules across the vertebrate spectrum and propose a unified nomenclature for CXC chemokines in fish, amphibians and reptiles based on phylogenetic analysis. Expanding on earlier studies in teleost fish, lineage specific CXC chemokines that have no apparent homologues in mammals were confirmed. Furthermore, in addition to the two subgroups of the CXCL8 homologues known in teleost fish, a third group was identified (termed CXCL8_L3), as was a further subgroup of the fish CXC genes related to CXCL11. Expression of the CXC chemokines found in rainbow trout, Oncorhynchus mykiss, was studied in response to stimulation with inflammatory and antiviral cytokines, and bacterial. Tissue distribution analysis revealed distinct expression profiles for these trout CXC chemokines. Lastly three of the trout chemokines, including two novel fish specific CXC chemokines containing three pairs of cysteines, were produced as recombinant proteins and their effect on trout leucocyte migration studied. These molecules increased the relative expression of CD4 and MCSFR in migrated cells in an in vitro chemotaxis assay.

Genomic characterization and transcriptional evidence for the involvement of complement component 7 in immune response of rock bream (Oplegnathus fasciatus)

The complement component 7 (C7) is the central mediator of pathogenic attack at the membrane surface and its binding to the C5b-7 complex triggers cytolytic signaling. In this study, C7 of rock bream (Oplegnathus fasciatus) was identified (Rb-C7) and characterized at the genomic level. The Rb-C7 gene contains 18 exons and 17 introns and is composed of a 2490 bp complete open reading frame (ORF). The encoded polypeptide (830 amino acids) contains a number of well-conserved C7 signature domains. Important putative transcription factor binding sites, including those for NF-κB, SP-1, C/EBP, AP-1 and OCT-1, are present in the 5'-flanking region of Rb-C7. Phylogenetic analysis revealed a close proximity of Rb-C7 with the orthologues in tilapia and Japanese flounder. Quantitative real-time PCR (qPCR) analysis confirmed constitutive Rb-C7 expression throughout all the examined tissue of healthy rock bream, with highest expression in liver. In immune challenge experiment, Rb-C7 expression was up-regulated in head kidney and liver in response to Edwardsiella tarda, Streptococcus iniae, lipopolysaccharide and rock bream iridovirus (RBIV). Furthermore, significant increases of both intracellular expression level and the number of Rb-C7-expressing cells were detected by in situ hybridization assay in head kidney and liver tissues upon E. tarda infection. These results suggested that Rb-C7 is lytic pathway gene in complement system and its transcriptional regulation may be an important immune response in pathogenic defense mechanism of rock bream.

Bisphenol A affects gene expression of gonadotropin-releasing hormones and type I GnRH receptors in brains of adult rare minnow Gobiocypris rarus

Recent studies support the notion that endocrine disrupting chemicals (EDCs) could affect the reproductive regulations of the neuroendocrine system. The objectives of the present study were to determine whether the weak estrogenic chemical, bisphenol A (BPA), disrupts gonadotropin-releasing hormone (GnRH) system by altering the transcription of GnRHs and GnRH receptor (GnRHR) genes in adult rare minnow Gobiocypris rarus. In the present study, the histological examination of the ovary after 35-day BPA exposure at 15 μg/L demonstrated the perturbing effects of environmentally relevant BPA on the ovarian development in G. rarus. In addition mRNA expression of ovarian P450 aromatase in both ovaries and testes were significantly down-regulated by 15 μg/L BPA. GnRH2, GnRH3, GnRHR1A and GnRHR1B gene were identified in G. rarus. The expression patterns of GnRHs and GnRHR1s were analyzed in various tissues of G. rarus by quantitative real-time PCR. GnRHs and GnRHR1s were all predominantly expressed in the brains. Both GnRH3 and GnRHR1A were significantly upregulated in the brains of female exposed to 15 μg/L BPA for 35 days. It would suggest a potential negative feedback in the GnRH system in response to the disturbance of downstream of the brain-pituitary-gonadal axis. Collectively, the present findings suggest that the transcripts of some key genes in the neuroendocrine system can be used as critical biomarkers in endocrine disruption assays of teleost fish.

The rad1 gene in Rainbow Trout (Oncorhynchus mykiss) is highly conserved and may express proteins from non-canonical spliced isoforms

Cell-cycle checkpoint proteins maintain genomic integrity by sensing damaged DNA and initiating DNA repair or apoptosis. RAD1 is a checkpoint protein involved in the sensing of damaged DNA and is a part of the 9-1-1 complex. In this project rainbow trout rad1 (rtrad1) was cloned, sequenced, expressed as a recombinant protein and anti-rtRAD1 antibodies were developed. RAD1 protein levels were characterized in various rainbow trout tissues. It was determined that an 840 bp open-reading frame encodes 279 aa with a predicted protein size of 31 kDa. The rtRAD1 amino-acid sequence is highly conserved and contains conserved exonuclease and leucine zipper domains. RT-PCR was used to identify three non-canonical splice variants of rtrad1, two of which are capable of forming functional proteins. The rad1 splice variant that encodes an 18 kDa protein appears to be abundant in rainbow trout spleen, heart and gill tissue and in the RTgill-W1 cell-line. Based on the genomic rtrad1 sequence the splice variants contain only partial exons which are consistent with the splicing of rad1 variants in mammals. This is the first time that rad1 has been fully characterized in a fish species.

Dynamic transcriptional regulation of autocrine/paracrine igfbp1, 2, 3, 4, 5, and 6 in the skeletal muscle of the fine flounder during different nutritional statuses

The IGF-binding proteins (IGFBPs) play a dual role in the regulation of the activity and bioavailability of IGFs in different tissues. Diverse evidence has shown that IGFBPs can inhibit and/or potentiate IGF actions. In this study, igfbp1, 2, 3, 4, 5, and 6 were isolated in the fine flounder, a flat fish species that shows slow growth and inherent Gh resistance in muscle. Subsequently, the expression of all igfbps was assessed in the skeletal muscle of flounder that underwent different nutritional statuses. igfbp1 was not expressed in muscle during any of the nutritional conditions, whereas igfbp3 and igfbp5 were the lowest and the highest igfbps expressed respectively. A dynamic expression pattern was found in all the igfbps expressed in skeletal muscle, which depended on the nutritional status and sampling period. During the fasting period, igfbp2, 4, and 5 were downregulated, whereas igfbp3 was upregulated during part of the fasting period. The restoration of food modulated the expression of the igfbps dynamically, showing significant changes during both the long- and short-term refeeding. igfbp3 and igfbp6 were downregulated during short-term refeeding, whereas igfbp5 was upregulated, and igfbp2 and igfbp4 remained stable. During long-term refeeding, the expression of igfbp2, 4, 5, and 6 increased, while igfbp3 remained unchanged. In conclusion, this study shows for the first time the isolation of all igfbps in a single fish species, in addition to describing a dynamic nutritional and time-dependent response in the expression of igfbps in the skeletal muscle of a nonmammalian species.

The complement system in teleost fish: progress of post-homolog-hunting researches

Studies on the complement system of bony fish are now finishing a stage of homologue-hunting identification of the components, unveiling existence of almost all the orthologues of mammalian complement components in teleost. Genomic and transcriptomic data for several teleost species have contributed much for the homologue-hunting research progress. Only an exception is identification of orthologues of mammalian complement regulatory proteins and complement receptors. It is of particular interest that teleost complement components often exist as multiple isoforms with possible functional divergence. This review summarizes research progress of teleost complement system following the molecular identification and sequence analysis of the components. The findings of extensive expression analyses of the complement components with special emphasis of their prominent extrahepatic expression, acute-phase response to immunostimulation and various microbial infections, and ontogenic development including maternal transfer are discussed to infer teleost-specific functions of the complement system. Importance of the protein level characterization of the complement components is also emphasized, especially for understanding of the isotypic diversity of the components, a unique feature of teleost complement system.

Structural and functional diversity of the lectin repertoire in teleost fish: relevance to innate and adaptive immunity

Protein-carbohydrate interactions mediated by lectins have been recognized as key components of innate immunity in vertebrates and invertebrates, not only for recognition of potential pathogens, but also for participating in downstream effector functions, such as their agglutination, immobilization, and complement-mediated opsonization and killing. More recently, lectins have been identified as critical regulators of mammalian adaptive immune responses. Fish are endowed with virtually all components of the mammalian adaptive immunity, and are equipped with a complex lectin repertoire. In this review, we discuss evidence suggesting that: (a) lectin repertoires in teleost fish are highly diversified, and include not only representatives of the lectin families described in mammals, but also members of lectin families described for the first time in fish species; (b) the tissue-specific expression and localization of the diverse lectin repertoires and their molecular partners is consistent with their distinct biological roles in innate and adaptive immunity; (c) although some lectins may bind endogenous ligands, others bind sugars on the surface of potential pathogens; (d) in addition to pathogen recognition and opsonization, some lectins display additional effector roles, such as complement activation and regulation of immune functions; (e) some lectins that recognize exogenous ligands mediate processes unrelated to immunity: they may act as anti-freeze proteins or prevent polyspermia during fertilization.

Identification and characterization of a novel cathelicidin from ayu, Plecoglossus altivelis

Cathelicidins, one family of antimicrobial peptides, play important roles against infections in animals. In this study, a cDNA sequence coding for cathelicidin was cloned from constructed liver cDNA library of ayu, Plecoglossus altivelis. The deduced ayu cathelicidin (aCATH) has a 20 amino acid residue signal peptide, a conserved cathelin domain of 110 amino acid residues and a mature antimicrobial peptide of 61 amino acid residues. Sequence comparison and phylogenetic tree analysis confirmed aCATH as a distinct member of fish cathelicidins. Real-time quantitative PCR revealed that the aCATH transcripts dramatically increased in various tissues after bacterial infection. Subsequently, aCATH was prokaryotic expressed and purified. Western blot and mass spectrometry revealed that aCATH was cleaved at residue Ile130-Arg131 by human neutrophil elastase to release the mature antimicrobial peptide. The mature peptide of aCATH was chemically synthesized and exhibited potent antimicrobial activity. Thus, aCATH may play an important role in the innate immunity of ayu, and this work enriches our knowledge in fish antimicrobial peptides.

Circulating salmon 41-kDa insulin-like growth factor binding protein (IGFBP) is not IGFBP-3 but an IGFBP-2 subtype

In vertebrates, most circulating insulin-like growth factor (IGF) is bound to multiple forms of IGF-binding proteins (IGFBPs) that differ both structurally and functionally. In mammals, the largest reservoir of IGF in the circulation comes from a large (150kDa) ternary complex comprised of IGF bound to IGFBP-3, which is bound to an acid label subunit (ALS), and this variant of IGFBP is regulated by growth hormone (GH) and feed intake. Although multiple variants of IGFBPs ranging from 20 to 50kDa have been found in fishes, no ternary complex is present and it has been assumed that the majority of circulating IGF is bound to fish IGFBP-3. Consistent with this assumption is previous work in salmon showing the presence of a 41-kDa IGFBP that is stimulated by GH, decreases with fasting and increases with feeding. However, the hypothesis that the salmon 41-kDa IGFBP is structurally homologous to mammalian IGFBP-3 has not been directly tested. To address this issue, we cloned cDNAs for several Chinook salmon IGFBPs, and found that the cDNA sequence of the 41-kDa IGFBP is most similar to that of mammalian IGFBP-2 and dissimilar to IGFBP-3. We found an additional IGFBP (termed IGFBP-2a) with high homology to mammalian IGFBP-2. These results demonstrate that salmon 41-kDa IGFBP is not IGFBP-3, but a paralog of IGFBP-2 (termed IGFBP-2b). Salmon IGFBP-2s are also unique in terms of having potential N-glycosylation sites and splice variants. Additional research on non-mammalian IGFBPs is needed to fully understand the molecular/functional evolution of the IGFBP family and the significance of the ternary complex in vertebrates.

Isolation and characterization of koi herpesvirus (KHV) from Indonesia: identification of a new genetic lineage

Koi herpesvirus (KHV) is the aetiological agent of an emerging disease (KHVD) associated with mass mortalities in koi and common carp and reported from at least 30 countries. We report the first isolation of KHV from koi and common carp in Indonesia and initial characterization of the isolates. Clinical signs, histopathology and virion morphology are similar to those of isolates from other countries. Phylogenetic analyses using the thymidine kinase gene amplified from each isolate and from carp tissue samples collected from KHVD outbreaks throughout Indonesia indicated that the Indonesian isolates are more closely related to the Asian than the European KHV lineage. Sequence analysis of two other variable regions between ORF29 and ORF31 (marker I) and near the start of ORF 133 (marker II) indicated that all Indonesian isolates displayed a marker I allele (I(++)) previously identified only in isolates of the Asian lineage. However, in the marker II region, all Indonesian isolates displayed the II(-) allele, which has been reported previously only amongst isolates of the European lineage, and nine of these displayed a mixed genotype (II(+)II(-)). The I(++)II(-) genotype has not been reported previously and appears to represent a new intermediate lineage that may have emerged in Indonesia.

Molecular characterization and gene expression of thyrotropin receptor (TSHR) and a truncated TSHR-like in Senegalese sole

Thyroid hormones (THs) play a key role in larval development, growth and metamorphosis in flatfish. Their synthesis is tightly regulated by the hypothalamic-pituitary-thyroid axis. Thyroid-stimulating hormone receptor (TSHR) is a key protein in the control of thyroid function stimulating TH synthesis after binding its ligand, the thyrotropin. In teleost fish, numerous reports have associated the TSHR with gametogenesis. However, little information about its role during larval development is available. In this study, we report the cloning of two different cDNAs with high similarity to TSHR. Phylogenetic analysis clustered both cDNAs separately. One of them (referred to TSHR) grouped with TSHR orthologs in tetrapods and teleost fish and possessed the three typical conserved domains and regulatory motifs. The second receptor (referred to as TSHRtr-like) represented a novel truncated cDNA bearing the extracellular and part of the transmembrane domain. TSHRtr-like orthologs were only found in teleosts, which suggests that it could have appeared after fish-specific 3R genome duplication. Expression profiles of both genes are analyzed in juvenile tissues and during larval development using a real-time PCR approach. In juvenile fish, TSHR and TSHRtr-like are expressed ubiquitously although transcript levels varied between organs. In both cases, the highest mRNAs levels are detected in brain. During larval development, both genes are expressed to a high level during the first stages (2-3days after hatching) reducing progressively their abundance in the whole larvae during metamorphosis. This reduction in mRNA abundance is more accentuated for the TSHRtr-like gene. To evaluate the possible regulation of both receptors by T4 during sole metamorphosis, larvae are exposed to the goitrogen thiourea (TU). Only TSHRtr-like modifies its expression, increasing its transcripts at 11days after treatment. Moreover, adding exogenous T4 hormone to TU-treated larvae restores the TSHRtr-like steady-state levels similar to the untreated control. Overall, these results demonstrate the existence of two thyrotropin receptors differentially regulated by THs in teleosts.

Differential expressions of melanocortin receptor subtypes in melanophores and xanthophores of barfin flounder

alpha-Melanocyte-stimulating hormone (alpha-MSH) is a member of the melanocortin (MC) family, and the MC receptor (MCR) is a member of the G protein-coupled receptor (GPCR) superfamily. We previously found that in barfin flounder, a flatfish, alpha-MSH with an acetyl group at the N-terminus stimulated pigment dispersion in xanthophores; however, this effect was not observed in melanophores. Therefore, the present study was undertaken to find which MCR subtypes are expressed in these pigment cells in order to elucidate how acetylation regulates activities of alpha-MSH-related peptides. Here, we also report the cloning of Mc1r and Mc5r from barfin flounder. Three types of cells-melanophores, xanthophores, and nonchromatophoric dermal cells-were isolated from the skin samples collected from the dorsal fin. These cells were then tested for the expression of Mc1r and Mc5r as well as Mc2r and Mc4r that we had previously cloned. Mc1r and Mc5r transcripts were detected in melanophores, and a sole Mc5r transcript was detected in xanthophores. We had previously found that the efficiency of alpha-MSH was higher than that of desacetyl-alpha-MSH for pigment dispersion in xanthophores. Acetylated MSH peptide may have increased binding affinity to MC5R, whereas alpha-MSH lacks melanin-dispersing activity. Increasing evidences indicate that many GPCRs form heterodimers, and this may affect the affinity of the ligand for the corresponding GPCR. Taken together, the expression of two different Mcr subtypes in melanophores may suggest that a heterodimer consisting of MC1R and MC5R may have a low binding affinity toward alpha-MSH. The present results clarify the types of MCRs that are expressed in melanophores and xanthophores of barfin flounder; furthermore, the results provide important clues about the functional regulation of alpha-MSH-related peptides.

Salmonid genomes have a remarkably expanded akirin family, coexpressed with genes from conserved pathways governing skeletal muscle growth and catabolism

Metazoan akirin genes regulate innate immunity, myogenesis, and carcinogenesis. Invertebrates typically have one family member, while most tetrapod and teleost vertebrates have one to three. We demonstrate an expanded repertoire of eight family members in genomes of four salmonid fishes, owing to paralog preservation after three tetraploidization events. Retention of paralogs secondarily lost in other teleosts may be related to functional diversification and posttranslational regulation. We hypothesized that salmonid akirins would be transcriptionally regulated in fast-twitch skeletal muscle during activation of conserved pathways governing catabolism and growth. The in vivo nutritional state of Arctic charr (Salvelinus alpinus L.) was experimentally manipulated, and transcript levels for akirin family members and 26 other genes were measured by quantitative real-time PCR (qPCR), allowing the establishment of a similarity network of expression profiles. In fasted muscle, a class of akirins was upregulated, with one family member showing high coexpression with catabolic genes coding the NF-kappaB p65 subunit, E2 ubiquitin-conjugating enzymes, E3 ubiquitin ligases, and IGF-I receptors. Another class of akirin was upregulated with subsequent feeding, coexpressed with 14-3-3 protein genes. There was no similarity between expression profiles of akirins with IGF hormones or binding protein genes. The level of phylogenetic relatedness of akirin family members was not a strong predictor of transcriptional responses to nutritional state, or differences in transcript abundance levels, indicating a complex pattern of regulatory evolution. The salmonid akirins epitomize the complexity linking the genome to physiological phenotypes of vertebrates with a history of tetraploidization.