Pattern of New Gene Origination in a Special Fish Lineage, the Flatfishes.

Origination of new genes are of inherent interest of evolutionary geneticists for decades, but few studies have addressed the general pattern in a fish lineage. Using our recent released whole genome data of flatfishes, which evolved one of the most specialized body plans in vertebrates, we identified 1541 (6.9% of the starry flounder genes) flatfish-lineage-specific genes. The origination pattern of these flatfish new genes is largely similar to those observed in other vertebrates, as shown by the proportion of DNA-mediated duplication (1317; 85.5%), RNA-mediated duplication (retrogenes; 96; 6.2%), and de novo-origination (128; 8.3%). The emergence rate of species-specific genes is 32.1 per Mya and the whole average level rate for the flatfish-lineage-specific genes is 20.9 per Mya. A large proportion (31.4%) of these new genes have been subjected to selection, in contrast to the 4.0% in primates, while the old genes remain quite similar (66.4% vs. 65.0%). In addition, most of these new genes (70.8%) are found to be expressed, indicating their functionality. This study not only presents one example of systematic new gene identification in a teleost taxon based on comprehensive phylogenomic data, but also shows that new genes may play roles in body planning.

Genome assembly and annotation at the chromosomal level of first Pleuronectidae: Verasper variegatus provides a basis for phylogenetic study of Pleuronectiformes.

High quality genome is of great significance for the mining of biological information resources of species. Up to now, the genomic information of several important economic flatfishes has been well explained. All these fishes are eyes on left side-type, and no high-quality genome of eyes on right side-type species has been reported. In this study, we applied a combined strategy involving stLFR and Hi-C technologies to generate sequencing data for constructing the chromosomal genome of Verasper variegates, which belongs to Pleuronectidae with characteristic of eyes on right side. The size of genome of V. variegatus is 556 Mb. More than 97.2% of BUSCO genes were detected, and N50 lengths of the contigs and scaffolds reached 79.8 Kb and 23.8 Mb, respectively, demonstrating the outstanding completeness and sequence continuity of the genome. A total of 22,199 protein-coding genes were predicted in the assembled genome, and more than 95% of those genes could be functionally annotated. Meanwhile, the genomic collinearity, gene family and phylogenetic analyses of similar species in Pleuronectiformes were also investigated and portrayed for metamorphosis and benthic adaptation. Sex related genes mapping has also been achieved at the chromosome level. This study is the first chromosomal level genome of a Pleuronectidae fish (V. variegatus). The chromosomal genome assembly constructed in this work will not only be valuable for conservation and aquaculture studies of the V. variegatus but will also be of general interest in the phylogenetic and taxonomic studies of Pleuronectiformes.

Population genetic structure of turbot (Scophthalmus maximus L., 1758) in the Black Sea.

Turbot, Scophthalmus maximus, is a commercially important demersal flatfish species distributed throughout the Black Sea. Several studies performed locally with a limited number of specimens using both mitochondrial DNA (mtDNA) and microsatellite markers evidenced notable genetic variation among populations. However, comprehensive population genetic studies are required to help management of the species in the Black Sea. In the present study eight microsatellite loci were used to resolve the population structure of 414 turbot samples collected from 12 sites across the Black Sea. Moreover, two mtDNA genes, COI and Cyt-b, were used for taxonomic identification. Microsatellite markers of Smax-04 and B12-I GT14 were excluded from analysis due to scoring issues. Data analysis was performed with the remaining six loci. Loci were highly polymorphic (average of 17.8 alleles per locus), indicating high genetic variability. Locus 3/20CA17, with high null allele frequency (>30%), significantly deviated from HW equilibrium. Pairwise comparison of the FST index showed significant differences between most of the surveyed sampling sites (P < 0.01). Cluster analysis evidenced the presence of three genetic groups among sampling sites. Significant genetic differentiation between Northern (Sea of Azov and Crimea) and Southern (Turkish Black Sea Coast) Black Sea sampling sites were detected. The Mantel test supported an isolation by distance model of population structure. These findings are vital for long-term sustainable management of the species and development of conservation programs. Moreover, generated mtDNA sequences would be useful for the establishment of a database for S. maximus.

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.

Double control regions of some flatfish mitogenomes evolve in a concerted manner.

Mitochondrial genomes (mitogenomes) typically contain 13 protein-coding genes, 22 tRNA genes, two rRNA genes, and a single control region (CR). Flatfish mitochondrial genomes (mitogenomes) from three different genera in Bothidae (bothids) contain double CRs that evolved in a concerted manner. How these double CRs maintained identical sequences throughout the evolutionary process is an interesting issue. In the present study, over four hundred arrays of the double CRs of mitogenomes from three bothids (Arnoglossus tenuis, Lophonectes gallus and Psettina iijimae) were performed. Interesting variations between double CRs were observed in P. iijimae mitogenomes, and the networks of CR sequences from P. iijimae indicated a high possibility of genetic information exchange between CRs. No recombination product was detected in our results, indicating that the mechanism of the concerted evolution between the double CRs of P. iijimae was not recombination. We speculate that mismatch repair, a mitochondrial DNA repair mechanism, is a potential explanation for the concerted evolution between these double CRs.

Identification and initial functional characterization of lysosomal integral membrane protein type 2 (LIMP-2) in turbot (Scophthalmus maximus L.).

The immune system protects organism from external pathogens, this progress starts with the pathogen recognition by pattern recognition receptors (PRRs). As a group of PRRs, the class B scavenger receptors showed important roles in phagocytosis. Among three class B scavenger receptors, lysosomal integral membrane protein type 2 (LIMP-2) was reported to present in the limiting membranes of lysosomes and late endosomes, but its immune roles in teleost species are still limited in handful species. Here, we characterized LIMP-2 gene in turbot, and its expression patterns in mucosal barriers following different bacterial infection, as well as ligand binding activities to different microbial ligands and agglutination assay with different bacteria. In our results, one SmLIMP2 gene was identified with a 1,593 bp open reading frame (ORF). The multiple species comparison and phylogenetic analysis showed the closest relationship to Paralichthys olivaceus, the genomic structure analysis and syntenic analysis revealed the conservation of LIMP-2 during evolution. In tissue distribution analysis, SmLIMP-2 was expressed in all the examined turbot tissues, with the highest expression level in brain, and the lowest expression level in liver. In addition, SmLIMP-2 was significantly up-regulated in all the mucosal tissues (skin, gill and intestine) following Gram-negative bacteria Vibrio anguillarum infection, and was only up-regulated in gill following Gram-positive bacteria Streptococcus iniae challenge. Finally, the rSmLIMP-2 showed strong binding ability to all the examined microbial ligands, and strong agglutination with Escherichia coli, Staphylococcus aureus and V. anguillarum. Taken together, our results suggested SmLIMP-2 played important roles in fish immune response to bacterial infection. However, further functional studies should be carried out to better characterize its detailed roles in teleost immunity.

Fish diversity of the largest deltaic formation in the Americas - a description of the fish fauna of the Parnaíba Delta using DNA Barcoding.

Deltas are dynamic and productive systems of enormous ecological significance, encompassing unique and biologically diverse wetland habitats. Here, we present the first data on the molecular diversity of the fish fauna of the Parnaíba Delta, the largest deltaic formation of the Americas. Partial sequences (626 bp) of the mitochondrial COI gene (Cytochrome c oxidase subunit I) were used to barcode 402 individuals, representing 128 species, belonging to 98 genera, 57 families, 17 orders and two classes. The most abundant orders were the Perciformes, Siluriformes, Gobiiformes, and Pleuronectiformes. The Neighbor-Joining (NJ), Bayesian Inference (BI), and BIN analyses produced 103 molecular clusters, while the Automatic Barcode Gap Discovery (ABGD) and Maximum Likelihood (ML) approaches revealed 102 clusters. The mean conspecific, congeneric and confamilial genetic distances were 0.33%, 14.37%, and 18.60%, respectively. Intraspecific divergence ranged from 0.0% to 1.4%, and all species presented barcode gaps, with the exception of two clusters of Cathorops spixii (OTU 96 and OTU 103), which were separated by a low interspecific distance (1.2%), which overlaps the maximum intraspecific genetic distance (1.4%). The barcode data provide new insights into the fish diversity of the Parnaíba Delta, which will be important for the development of further research on this fauna.

Seminal Plasma Exosomes: Promising Biomarkers for Identification of Male and Pseudo-Males in Cynoglossus semilaevis.

In mammals, small RNAs enclosed in exosomes have been identified as appropriate signatures for disease diagnosis. However, there is limited information on exosomes derived from seminal plasma, and few studies have reported analyzed the composition of exosomes and enclosed small RNAs in fish. The half-smooth tongue sole (Cynoglossus semilaevis) is an economically important fish for aquaculture, and it exhibits sexual dimorphism: the female gender show higher growth rates and larger body sizes than males. Standard karyotype analysis and artificial gynogenesis tests have revealed that this species uses heterogametic sex determination (ZW/ZZ), and so-called sex-reversed pseudo-males exist. In this study, we successfully identified exosomes in the seminal plasma of C. semilaevis; to the best of our knowledge, this is the first report of exosomes in fish seminal plasma. Analysis of the nucleotide composition showed that miRNAs were dominant in the exosomes, and the miRNAs were sequenced and compared to identify signature miRNAs as sexual biomarkers. Moreover, target genes of the signature miRNAs were predicted by sequence matching and annotation. Finally, four miRNAs (dre-miR-141-3P, dre-miR-10d-5p, ssa-miR-27b-3p, and ssa-miR-23a-3p) with significant differential expression in the males and pseudo-males were selected from the signature candidate miRNAs as markers for sex identification, and their expression profiles were verified using real-time quantitative PCR. Our findings could provide an effective detection method for sex differentiation in fish.

Identification and expression analysis of fetuin B (FETUB) in turbot (Scophthalmus maximus L.) mucosal barriers following bacterial challenge.

Fetuin B (FETUB), a recently described cysteine proteinase inhibitor, has numerous conserved N-glycosylation sites, species-specific O-glycosylation sites, and two cystatin (CY) domains. FETUB is likely to play regulatory roles in acute inflammation, female infertility, fish organogenesis and tumor suppression. In the present study, transcript of turbot FETUB gene was captured, its protein structure and expression patterns in different tissues with emphasis on mucosal barriers following different bacterial infection were characterized. Turbot FETUB gene showed the closest relationship with Takifugu rubripes in phylogenetic analysis. In addition, FETUB was ubiquitously expressed in all examined tissues with the highest expression level in skin. Finally, FETUB gene showed different expression patterns following both bacterial challenge. The rapidly and significantly differential expression patterns of FETUB in mucosal surfaces against bacterial infections might indicate its key roles to prevent pathogen attachment and entry in turbot mucosal immunity. Functional studies should be carried out to further characterize the FETUB and avail utilization of its function to increase the disease resistance of turbot in maintaining the integrity of the mucosal barriers against infections and to facilitate selection of the fine family/varieties of disease resistance in turbot.

Can an eel be a flatfish? Observations on enigmatic asymmetrical heterenchelyids from the Guinea coast of West Africa.

Morphological asymmetry is described in the heterenchelyid mud eel Pythonichthys cf. macrurus from inshore coastal waters of Guinea, West Africa. The intensity of asymmetry differs between two examined specimens, with the more extreme case exhibiting strong asymmetry in both external and internal features, including unilateral depigmentation, reductive degeneration and embedding of a blind-side eye, skewed jaws with reduced dentition and tooth loss. The extent and nature of asymmetry suggests that this individual probably lived primarily on its left lateral side, not unlike sinistral pleuronectiform flatfishes.

Phylogeny, expression patterns and regulation of DNA Methyltransferases in early development of the flatfish, Solea senegalensis.

BACKGROUND: The identification of DNA methyltransferases (Dnmt) expression patterns during development and their regulation is important to understand the epigenetic mechanisms that modulate larval plasticity in marine fish. In this study, dnmt1 and dnmt3 paralogs were identified in the flatfish Solea senegalensis and expression patterns in early developmental stages and juveniles were determined. Additionally, the regulation of Dnmt transcription by a specific inhibitor (5-aza-2'-deoxycytidine) and temperature was evaluated. RESULTS: Five paralog genes of dnmt3, namely dnmt3aa, dnmt3ab, dnmt3ba, dnmt3bb.1 and dnmt3bb.2 and one gene for dnmt1 were identified. Phylogenetic analysis revealed that the dnmt gene family was highly conserved in teleosts and three fish-specific genes, dnmt3aa, dnmt3ba and dnmt3bb.2 have evolved. The spatio-temporal expression patterns of four dnmts (dnmt1, dnmt3aa, dnmt3ab and dnmt3bb.1) were different in early larval stages although all of them reduced expression with the age and were detected in neural organs and dnmt3aa appeared specific to somites. In juveniles, the four dnmt genes were expressed in brain and hematopoietic tissues such as kidney, spleen and gills. Treatment of sole embryos with 5-aza-2'-deoxycytidine down-regulated dntm1 and up-regulated dntm3aa. Moreover, in lecithotrophic larval stages, dnmt3aa and dnmt3ab were temperature sensitive and their expression was higher in larvae incubated at 16 °C relative to 20 °C. CONCLUSION: Five dnmt3 and one dnmt1 paralog were identified in sole and their distinct developmental and tissue-specific expression patterns indicate that they may have different roles during development. The inhibitor 5-aza-2'-deoxycytidine modified the transcript abundance of dntm1 and dntm3aa in embryos, which suggests that a regulatory feedback mechanism exists for these genes. The impact of thermal regime on expression levels of dnmt3aa and dnmt3ab in lecithotrophic larval stages suggests that these paralogs might be involved in thermal programing.

The involvement of cathepsin F gene (CTSF) in turbot (Scophthalmus maximus L.) mucosal immunity.

Cathepsin F (CTSF) is a recently described papain-like cysteine protease and unique among cathepsins due to an elongated N-terminal pro-region, which contains a cystatin domain. CTSF likely plays a regulatory role in processing the invariant chain which is associated with the major histocompatibility complex (MHC) class II. In this regard, we identified the CTSF gene of turbot as well as its protein structure, phylogenetic relationships, and expression patterns in mucosal tissues following Vibrio anguillarum and Streptococcus iniae challenge. We also determined the expression patterns of CTSF in mucosal tissues after vaccinated with the formalin-inactivated V. vulnificus whole-cell vaccine. Briefly, turbot CTSF gene showed the closest relationship with that of Paralichthys olivaceus in phylogenetic analysis. And CTSF was ubiquitously expressed in all tested tissues with the highest expression level in gill. In addition, CTSF gene showed different expression patterns following different bacterial challenge. The significant quick regulation of CTSF in mucosal surfaces against infection indicated its roles in mucosal immunity. Functional studies should further characterize avail utilization of CTSF function to increase the disease resistance of turbot in maintaining the integrity of the mucosal barriers against infection and to facilitate selection of the disease resistant family/strain in turbot.

Characterization and expression profiling of NOD-like receptor C3 (NLRC3) in mucosal tissues of turbot (Scophthalmus maximus L.) following bacterial challenge.

The mucosal surfaces are important for teleost as they are directly and continuously exposed to pathogen-rich aquatic environments. Scrutinization and recognition of the attached pathogens is the first crucial step of mucosal immunity initiation. Nucleotide oligomerization domain (NOD)-like receptors (NLRs) are a large group of intracellular pathogen recognition receptors (PRRs) which play key roles in pathogen recognition and subsequent immune signaling pathways activation. In this study, we identified two NLRC3 genes (NLRC3a and NLRC3b), a subfamily of NLRs from turbot, and profiled their expression patterns in mucosal tissues following bacterial challenge. NLRC3a transcript contains an open reading frame (ORF) of 3405 bp that encodes a putative peptide of 1134 amino acids. While NLRC3b has an ORF of 3114 bp encoding 1037 amino acids. A caspase recruitment domain (CARD) at N-terminus characterized turbot NLRC3a, while NLRC3b seems to be unique to teleost, containing a fish specific NACHT associated (FISNA) domain and an extra B30.2 (PRY/SPRY) domain at C-terminus. In addition, NLRC3a and NLRC3b were detected in all the examined tissues, with the highest expression levels in kidney and blood, respectively. After bacteria challenge, expression levels of turbot NLRC3 genes were strongly induced in intestine rather than in skin and gill, while NLRC3a had relatively higher expression level than that of NLRC3b. Taken together, NLRC3 genes found in this study were the first NLR members identified in turbot. The different expression signatures of NLRC3a and NLRC3b in mucosal tissues following two bacterial infections indicated they probably have important roles in early response to bacterial infection in the first line of host defense system.

Centromeric Satellite DNA in Flatfish (Order Pleuronectiformes) and Its Relation to Speciation Processes.

Two new centromeric satellite DNAs in flatfish (Order Pleuronectiformes) have been characterized. The SacI-family from Hippoglossus hippoglossus, restricted to this species, had a monomeric size of 334 base pair (bp) and was located in most of the centromeres of its karyotype. The PvuII-family, with a monomeric size of 177 bp, was initially isolated from the genome of Solea senegalensis, and fluorescent in situ hybridization (FISH) localized the repeat to centromeres of most of the chromosomes. This family could only be amplified in 2 other species of the genus Solea (Solea solea and Solea lascaris). Molecular features and chromosomal location indicated a possible structural and/or functional role of these sequence repeats. The presence of species-specific satellite-DNA families in the centromeres and their possible role in the speciation processes in this group of fishes is discussed.

First Report of Sex Chromosomes in Achiridae (Teleostei: Pleuronectiformes) with Inferences About the Origin of the Multiple X1X1X2X2/X1X2Y System and Dispersal of Ribosomal Genes in Achirus achirus.

American soles (family Achiridae) have been characterized by remarkable chromosomal variation even though several species lack basic cytogenetic information. This trend indicates that chromosomal traits can be useful to taxonomy once the morphological identification of some taxa in this family (e.g., Achirus species) is controversial. In this work, we expand the cytogenetic data in Achiridae by providing the first karyotypic analysis of Achirus achirus. An unusual multiple sex chromosome system (X1X1X2X2/X1X2Y) was detected, once males presented 2n = 33 with three unpaired chromosomes (X1X2Y) while females presented 2n = 34 with two acrocentric pairs (X1X1 and X2X2) found in hemizygosis in males. The 18S rDNA clusters were observed interspersed with GC-rich sites in a single pair. However, the 5S rRNA genes were dispersed through the genome of A. achirus in a sex-specific manner (10 clusters in males and 12 in females), as a result of the presence of ribosomal cistrons in X1 and X2 chromosomes. This pattern allowed us to infer that Y chromosome has evolved by partial deletion followed by fusion of proto-X1 and proto-X2 homologous chromosomes. The high rate of genomic evolution in Achiridae could have favored their reproductive isolation and speciation even in sympatric conditions.

The complete mitochondrial genome of the marbled flounder Pseudopleuronectes yokohamae (Pleuronectiformes, Pleuronectidae).

We sequenced the complete mitochondrial genome of the marbled flounder Pseudopleuronectes yokohamae collected from the Yellow Sea off China. The mitogenome comprised a 16 864-bp circular DNA molecule containing 37 genes and an AT-rich control region known as the D-loop. Phylogenetic analysis based on the mitochondrial genomes of marbled flounder indicated that P. yokohamae and Verasper variegatus are the most closely related species, which strongly supports their close phylogenetic affinity.

Redescriptions of two closely related East Asian flatfish species of the genus Pleuronichthys.

Because of the taxonomic confusion, including misapplication of their scientific names, resulting from the morphological similarity of two East Asian flatfish species, Pleuronichthys lighti Wu, 1929 and Pleuronichthys cornutus (Temminck & Schlegel, 1846), both species are redescribed, with particular emphasis on some new key characters. New common names are proposed for each species.

Cloning and expression study of a Toll-like receptor 2 (tlr2) gene from turbot, Scophthalmus maximus.

Toll-like receptor 2 (TLR2) in mammals is a member of the ancient Toll-like family of receptors that predominantly recognizes conserved components of Gram-positive bacteria. In the present study, a tlr2 gene and its 5'-flanking sequence were cloned from turbot, Scophthalmus maximus, its responsive expressions to various immunostimulants were subsequently studied in vivo. The turbot (sm)tlr2 gene spans over 9.0 kb with a structure of 12 exon-11 intron and encodes 816 amino acids. The deduced protein shows the highest sequence identity (76.1%) to Japanese flounder Tlr2 and possesses a signal peptide sequence, a leucine-rich repeat (LRR) domain composed of 19 LRR motifs, a transmembrane region and a Toll/interleukin-1 receptor (TIR) domain. Phylogenetic analysis grouped it with other neoteleostei Tlr2as. A number of transcription factor binding sites known to be important for the basal transcriptional activity of TLR3 and response of TLR2 to lipopolysaccharide (LPS) signalling in mammals were predicted in the 5'-flanking sequence of smtlr2. Quantitative real-time PCR (qPCR) analysis demonstrated the constitutive expression of smtlr2 mRNA in all twelve examined tissues with higher levels in the lymphomyeloid-rich tissues and liver. Further, smtlr2 expression was up-regulated following stimulation with LPS, peptidoglycan (PGN) or polyinosinic: polycytidylic acid [poly(I:C)] in the gills, head kidney, spleen and muscle. Finally, for all three immunostimulants, a two-wave induced smtlr2 expression was observed in the head kidney and spleen in a 7-day time course and the strongest inducibility in the head kidney. These findings suggest a possible role of Smtlr2 in the immune responses to the infections of a broad range of pathogens that include Gram-positive and Gram-negative bacteria and RNA virus.

Phylogenomic analysis of carangimorph fishes reveals flatfish asymmetry arose in a blink of the evolutionary eye.

BACKGROUND: Flatfish cranial asymmetry represents one of the most remarkable morphological innovations among vertebrates, and has fueled vigorous debate on the manner and rate at which strikingly divergent phenotypes evolve. A surprising result of many recent molecular phylogenetic studies is the lack of support for flatfish monophyly, where increasingly larger DNA datasets of up to 23 loci have either yielded a weakly supported flatfish clade or indicated the group is polyphyletic. Lack of resolution for flatfish relationships has been attributed to analytical limitations for dealing with processes such as nucleotide non-stationarity and incomplete lineage sorting (ILS). We tackle this phylogenetic problem using a sequence dataset comprising more than 1,000 ultraconserved DNA element (UCE) loci covering 45 carangimorphs, the broader clade containing flatfishes and several other specialized lineages such as remoras, billfishes, and archerfishes. RESULTS: We present a phylogeny based on UCE loci that unequivocally supports flatfish monophyly and a single origin of asymmetry. We document similar levels of discordance among UCE loci as in previous, smaller molecular datasets. However, relationships among flatfishes and carangimorphs recovered from multilocus concatenated and species tree analyses of our data are robust to the analytical framework applied and size of data matrix used. By integrating the UCE data with a rich fossil record, we find that the most distinctive carangimorph bodyplans arose rapidly during the Paleogene (66.0-23.03 Ma). Flatfish asymmetry, for example, likely evolved over an interval of no more than 2.97 million years. CONCLUSIONS: The longstanding uncertainty in phylogenetic hypotheses for flatfishes and their carangimorph relatives highlights the limitations of smaller molecular datasets when applied to successive, rapid divergences. Here, we recovered significant support for flatfish monophyly and relationships among carangimorphs through analysis of over 1,000 UCE loci. The resulting time-calibrated phylogeny points to phenotypic divergence early within carangimorph history that broadly matches with the predictions of adaptive models of lineage diversification.

Identification and expression analysis of suppressors of cytokine signaling (SOCS) of Japanese flounder Paralichthys olivaceus.

Suppressor of cytokine signaling (SOCS) family members are key regulators of the immune system, particularly cytokine action, and have now been discovered in a number of fish species. Here we identified eight SOCS proteins (CISH, SOCS1a, SOCS1b, SOCS3a, SOCS3b, SOCS5, SOCS6 and SOCS9) in the Japanese flounder and analyzed their mRNA expressions after injection of poly (I:C) and formalin-killed cells (FKC) of Edwardsiella tarda. The expressions of all eight SOCS genes were detected in all the tissues examined. Stimulation of Japanese flounder reared at 15 or 25 °C with poly (I:C) affected the gene expressions of CISH, SOCS1a, SOCS1b and SOCS3a. All SOCS genes mRNA levels were significantly changed after FKC injection. Significant up-regulation of SOCS1a, SOCS1b, SOCS3a and SOCS3b genes was detected at 3, 12 and 24 hpi. SOCS5 and SOCS6 genes were significantly down-regulated at 3 hpi. SOCS9 gene was significantly up-regulated at 12 hpi. These results suggest that all eight of the SOCS genes are involved in immune responses, and that the CISH, SOCS1 and SOCS3 genes have functions distinct from those of the other SOCS members.