Search

Complement factor H-related protein (CFHR) plays an important role in regulating complement activation and defensive responses. The function of CFHR2 (complement factor H related 2), a member of the CFHR family, in fish remains unclear. Here, we report the genetic relationship, expression characteristics and regulatory mechanism of cfhl5 (complement factor H like 5) gene, which encodes CFHR2 in Chinese tongue sole. We observed that the cfhl5 gene was widely expressed in several tissues, such as brain, heart and immune organs, and was most abundantly expressed in liver. After injection with Vibrio harveyi, the expression of cfhl5 was up-regulated significantly in liver, spleen and kidney at 12 or 24 hours post infection (hpi), suggesting an involvement of this gene in the acute immune response. Knockdown of cfhl5 in liver cells significantly up-regulated the expression of the pro-inflammatory cytokines tnf-α (tumor necrosis factor-alpha) and il1β (interleukin-1beta), the immunomodulatory factor il10 (interleukin-10) and the lectin complement pathway gene masp1 (MBL-associated serine protease 1), and down-regulated the expression of complement components c3 (complement 3) and cfi (complement factor I). In our previous work, we found that cfhl5 gene was significantly higher methylated and lower expressed in the resistant family compared with the susceptible family. Therefore, we used dual-luciferase reporter system to determine the effect of DNA methylation on this gene and found that DNA methylation could inhibit the promoter activity to reduce its expression. These results demonstrated that the expression of cfhl5 is regulated by DNA methylation, and this gene might play an important role in the immune response by regulating the expression of cytokines and complement components genes in Chinese tongue sole., Competing Interests: Declaration of competing interest All the authors declare no competing interests., (Copyright © 2024 Elsevier Ltd. All rights reserved.)

Vibrio spp. are major pathogens responsible for mortality and disease in various marine aquaculture organisms. Effective disease control and genetic breeding strategies rely heavily on understanding host vibriosis resistance mechanisms. The Chinese tongue sole (Cynoglossus semilaevis) is economically vital but suffers from substantial mortalities due to vibriosis. Through continuous selective breeding, we have successfully obtained vibriosis-resistant families of this species. In this study, we conducted RNA-seq analysis on three organs, including liver, spleen and intestine from selected resistant and susceptible tongue soles. Additionally, we integrated these data with our previously published RNA-seq datasets of skin and gill, enabling the construction of organ-specific transcriptional profiles and a comprehensive gene co-expression network elucidating the differences in vibriosis resistance. Furthermore, we identified 12 modules with organ-specific functional implications. Overall, our findings provide a valuable resource for investigating the molecular basis of vibriosis resistance in fish, offering insights into target genes and pathways essential for molecular selection and genetic manipulation to enhance vibriosis resistance in fish breeding programs., (© 2024. The Author(s).)

Vibriosis, caused by Vibrio , seriously affects the health of fish, shellfish, and shrimps, causing large economic losses. Teleosts are represent the first bony vertebrates with both innate and adaptive immune responses against pathogens. Aquatic animals encounter hydraulic pressure and more pathogens, compared to terrestrial animals. The skin is the first line of defense in fish, constituting the skin-associated lymphoid tissue (SALT), which belongs to the main mucosa-associated lymphoid tissues (MALT). However, little is known about the function of immunity related proteins in fish. Therefore, this study used iTRAQ (isobaric tags for relative and absolute quantitation) to compare the skin proteome between the resistant and susceptible families of Cynoglossus semilaevis . The protein integrin beta-2, the alpha-enolase isoform X1, subunit B of V-type proton ATPase, eukaryotic translation initiation factor 6, and ubiquitin-like protein ISG15, were highly expressed in the resistant family. The 16S sequencing of the skin tissues of the resistant and susceptible families showed significant differences in the microbial communities of the two families. The protein-microbial interaction identified ten proteins associated with skin microbes, including immunoglobulin heavy chain gene (IGH), B-cell lymphoma/leukemia 10 (BCL10) and pre-B-cell leukemia transcription factor 1 isoform X2 (PBX2). This study highlights the interaction between skin proteins and the microbial compositions of C. semilaevis and provides new insights into understanding aquaculture breeding research., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2024 Wang, Tian and Chen.)

The assembly of W and Y chromosomes poses significant challenges in vertebrate genome sequencing and assembly. Here, we successfully assembled the W chromosome of Verasper variegatus with a length of 20.48 Mb by combining population and PacBio HiFi sequencing data. It was identified as a young sex chromosome and showed signs of expansion in repetitive sequences. The major component of the expansion was Ty3/Gypsy. The ancestral Osteichthyes karyotype consists of 24 protochromosomes. The sex chromosomes in four Pleuronectiformes species derived from a pair of homologous protochromosomes resulting from a whole-genome duplication event in teleost fish, yet with different sex-determination systems. V. variegatus and Cynoglossus semilaevis adhere to the ZZ/ZW system, while Hippoglossus stenolepis and H. hippoglossus follow the XX/XY system. Interestingly, V. variegatus and H. hippoglossus derived from one protochromosome, while C. semilaevis and H. stenolepis derived from another protochromosome. Our study provides valuable insights into the evolution of sex chromosomes in flatfish and sheds light on the important role of whole-genome duplication in shaping the evolution of sex chromosomes., (© 2023. The Author(s).)

The Chinese tongue sole (Cynoglossus semilaevis) holds significant economic importance within the fishing industry along the eastern coasts of China. In recent years, the frequent outbreaks of bacterial diseases have become a common concern as the aquaculture scale expands. The majority of the diseased fish exhibit symptoms such as skin congestion, damage and skin ulceration. As the skin serves as the first line of defense against bacterial infections, establishing a skin cell line for immunological research on Chinese tongue sole's response to bacterial infection is of utmost importance. In this study, a cell line named CSS (derived from the skin of the Chinese tongue sole) was successfully established. The cells have demonstrated stability during passages and exhibit a multipolar fibroblast-like morphology. They were cultured in L-15 medium with 20% serum and have been successfully passed through 60 passages over a period of 20 months. The identification of the mitochondrial CO1 gene confirmed that the cell originated from Chinese tongue sole. The karyotype detection revealed that the cell had a chromosome number of 2n = 42. After being stored in liquid nitrogen for 15 months, the cells can maintain more than 75% viability upon recovery. After transfecting with cy3-labeled scramble siRNA and pEGFP-N3 plasmid, clear fluorescence was observed in the transfected cells. We observed that lipopolysaccharide (LPS) from Escherichia coli significantly upregulate the gene expression of various immune-related pathways at 2 h in the CSS cell line. Additionally, the differentially expressed genes showed a higher enrichment in immune-related pathways at 2 and 6 h after stimulation compared to the 24 h point. Moreover, we identified 347 genes that exhibited a gradual increase in expression during the 0-24 h stimulation period. These genes were primarily enriched in pathways related to Autophagy, GABAergic synapse, Apelin signaling and Ferroptosis. In general, the CSS cell line established in this study exhibits stable growth and can serve as a valuable tool for in vitro studies of immunology and other basic biologies of Chinese tongue sole., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023 Elsevier Ltd. All rights reserved.)

Chinese tongue sole (Cynoglossus semilaevis) is an economically important marine fish in China. However, vibriosis has caused huge mortality and economic losses in its culturing industry. To reveal the effect of DNA methylation on the resistance to vibriosis in tongue sole, we conducted RNA sequencing and whole genome bisulfite sequencing (WGBS), and compared the gene expressions and DNA methylation patterns between the resistant and susceptible families. We identified a total of 741 significantly differentially expressed genes (DEGs) in kidney and 17460 differentially methylated genes (DMGs), which were both enriched in immune-related pathways, such as "cAMP signaling pathway" and "NOD-like receptor signaling pathway". Through the correlation analysis of DEGs and DMGs, we identified two important immune pathways, including "complement and coagulation cascades", and "cytokine-cytokine receptor interaction", which played important roles in regulating the inflammation level and immune homeostasis. For example, the expression of proinflammatory cytokine il17c was down-regulated under the regulation of DNA methylation; in addition, the expression of protease-activated receptor 3 (par3) was up-regulated, which could induce the up-expressionof il8. These results demonstrated that the regulation of DNA methylation on the genes involved in immune responses might contribute to the resistance to vibriosis in tongue sole, and provided a basis for the control of diseases in fish aquaculture., Competing Interests: Declaration of competing interest All the authors declare no competing interests., (Copyright © 2023 Elsevier Ltd. All rights reserved.)

Long non-coding RNAs (lncRNAs) play crucial roles in a variety of biological processes, including stress response. However, the number, characteristics and stress-related expression of lncRNAs in turbot are still largely unknown. In this study, a total of 12,999 lncRNAs were identified at the genome-wide level of turbot for the first time using 24 RNA-seq datasets. Sequence characteristic analyses of transcripts showed that lncRNA transcripts were shorter in average length, lower in average GC content and in average expression level as compared to the coding genes. Expression pattern analyses of lncRNAs in 12 distinct tissues showed that lncRNAs, especially lincRNA, exhibited stronger tissue-specific expression than coding genes. Moreover, 612, 1351, 1060, 875, 420 and 1689 differentially expressed (DE) lncRNAs under Vibrio anguillarum , Enteromyxum scophthalmi , and Megalocytivirus infection and heat, oxygen, and salinity stress conditions were identified, respectively. Among them, 151 and 62 lncRNAs showed differential expression under various abiotic and biotic stresses, respectively, and 11 lncRNAs differentially expressed under both abiotic and biotic stresses were selected as comprehensive stress-responsive lncRNA candidates. Furthermore, expression pattern analysis and qPCR validation both verified the comprehensive stress-responsive functions of these 11 lncRNAs. In addition, 497 significantly co-expressed target genes (correlation coefficient ( R ) > 0.7 and q -value < 0.05) for these 11 comprehensive stress-responsive lncRNA candidates were identified. Finally, GO and KEGG enrichment analyses indicated that these target genes were enriched mainly in molecular function, such as cytokine activity and active transmembrane transporter activity, in biological processes, such as response to stimulus and immune response, and in pathways, such as protein families: signaling and cellular processes, transporters and metabolism. These findings not only provide valuable reference resources for further research on the molecular basis and function of lncRNAs in turbot but also help to accelerate the progress of molecularly selective breeding of stress-resistant turbot strains or varieties.

The Chinese tongue sole ( Cynoglossus semilaevis ) is a traditional, precious fish in China. Due to the large growth difference between males and females, the investigation of their sex determination and differentiation mechanisms receives a great deal of attention. Forkhead Box O ( FoxO ) plays versatile roles in the regulation of sex differentiation and reproduction. Our recent transcriptomic analysis has shown that foxo genes may participate in the male differentiation and spermatogenesis of Chinese tongue sole. In this study, six Csfoxo members ( Csfoxo1a , Csfoxo3a , Csfoxo3b , Csfoxo4 , Csfoxo6-like , and Csfoxo1a-like ) were identified. Phylogenetic analysis indicated that these six members were clustered into four groups corresponding to their denomination. The expression patterns of the gonads at different developmental stages were further analyzed. All members showed high levels of expression in the early stages (before 6 months post-hatching), and this expression was male-biased. In addition, promoter analysis found that the addition of C/EBPα and c-Jun transcription factors enhanced the transcriptional activities of Csfoxo1a , Csfoxo3a , Csfoxo3b , and Csfoxo4. The siRNA-mediated knockdown of the Csfoxo1a , Csfoxo3a , and Csfoxo3b genes in the testicular cell line of Chinese tongue sole affected the expression of genes related to sex differentiation and spermatogenesis. These results have broadened the understanding of foxo 's function and provide valuable data for studying the male differentiation of tongue sole.

Heat shock proteins 70 (HSP70s) are known to play essential roles in organisms' response mechanisms to various environmental stresses. However, no systematic identification and functional analysis has been conducted for HSP70s in the turbot ( Scophthalmus maximus ), a commercially important worldwide flatfish. Herein, 16 HSP70 genes unevenly distributed on nine chromosomes were identified in the turbot at the genome-wide level. Analyses of gene structure, motif composition, and phylogenetic relationships provided valuable data on the HSP70s regarding their evolution, classification, and functional diversity. Expression profiles of the HSP70 genes under five different stresses were investigated by examining multiple RNA-seq datasets. Results showed that 10, 6, 8, 10, and 9 HSP70 genes showed significantly up- or downregulated expression after heat-induced, salinity-induced, and Enteromyxum scophthalmi , Vibrio anguillarum, and Megalocytivirus infection-induced stress, respectively. Among them, hsp70 ( hspa1a ), hspa1b , and hspa5 showed significant responses to each kind of induced stress, and qPCR analyses further validated their involvement in comprehensive anti-stress, indicating their involvement in organisms' anti-stress mechanisms. These findings not only provide new insights into the biological function of HSP70s in turbot adapting to various environmental stresses, but also contribute to the development of molecular-based selective breeding programs for the production of stress-resistant turbot strains in the aquaculture industry.

The tumor necrosis factor receptor-associated factor 6 (TRAF6) can act as a fundamental adaptor protein in a chain reaction of signal transduction and cascade events to finish off immune defenses. However, immunomodulatory research on TRAF6 gene is still limited in fish. In this study, a novel miRNA, Cse-miR-33 was identified from the whole genome of Chinese tongue sole (Cynoglossus semilaevis). After separate infections with three different Vibrio strains (V. harveyi, V. anguillarum, V. parahemolyticus) and one virus (nervous necrosis virus, NNV), the expressions of CsTRAF6 and Cse-miR-33 displayed significant time-dependent changes in immune related tissues and the trends were opposite in general. Through target gene prediction and dual luciferase reporter assay, Cse-miR-33 was proven to regulate CsTRAF6 by combining with 3'-UTR sequence of the gene. The results of qRT-PCR and western blotting (WB) analyses showed that Cse-miR-33 blocked the translation of CsTRAF6 protein at post-transcriptional level, rather than degrading the target mRNA. Further experiment indicated that Cse-miR-33 inhibitor largely reduced the death rate of Chinese tongue sole caused by V. harveyi and NNV. The expressions of CsTRAF6-associated immune genes (such as CsIL-1R, CsMYD88, CsIRAK1, CsTNFα, CsIL6 and CsIL8) were also significantly changed in response to Cse-miR-33 agomir and inhibitor. The study suggested that Cse-miR-33 affected the immune response via targeting CsTRAF6 in C. semilaevis, which would provide us deep insights into miRNA-mediated regulatory network and help improve the immunity in fish., Competing Interests: Declaration of competing interest The authors read and approved this version of the final manuscript, confirmed the integrity of this work and declared no conflicts of interest., (Copyright © 2023. Published by Elsevier Ltd.)

The typical sexual size dimorphism (SSD) phenomenon of Chinese tongue sole (Cynoglossus semilaevis) seriously restricts the sustainable development of the fishing industry. Previous transcriptome analysis has found a close relationship between the steroid biosynthesis and C. semilaevis SSD. The 7-dehydrocholesterol reductase (dhcr7) and lathosterol 5-desaturase (sc5d) are two genes in the steroid biosynthesis pathway, playing important roles in lipid synthesis, cellular metabolism, and growth. The present study assessed their roles in the mechanism of C. semilaevis SSD. The quantitative polymerase chain reaction (qPCR) results showed that C. semilaevis dhcr7 was mainly expressed in female livers, and C. semilaevis sc5d was highly expressed in female livers and gonads. Dual-luciferase experiment showed that dhcr7 and sc5d promoters had strong transcriptional activity. The transcription factors E2F transcription factor 1 (E2F1), and CCAAT enhancer binding protein alpha (C/EBPα) significantly regulated the transcriptional activity of dhcr7 and sc5d promoters, respectively. Furthermore, small interfering RNA (siRNA) knockdown results showed that expression levels of several genes [SREBF chaperone (scap), membrane-bound transcription factor peptidase, site 1 (mbtps1), fatty acid synthase (fasn), sonic hedgehog (shh), bone morphogenetic protein 2b (bmp2b) and AKT serine/threonine kinase 1 (akt1)] were suppressed. Protein subcellular localization results indicated that Dhcr7 and Sc5d were both specifically distributed in the cytoplasm, with co-localization been observed. The present study provides evidence that dhcr7 and sc5d might regulate C. semilaevis sexual size dimorphism by involving in energy homeostasis and cell cycle, or by affecting PI3K-Akt and Shh signaling pathways. The detailed roles of these steroid biosynthesis genes regulating C. semilaevis SSD needed more information., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2022 Elsevier B.V. All rights reserved.)

The phenomenon of sexual size dimorphism (SSD), existing in mammals, birds, reptiles, spiders, amphibians, insects, and fishes, is generally related to feeding efficiency, energy allocation, sex steroids, and somatotropic and reproductive endocrine axes. Recently, positive and negative regulations of sex steroids have been reported on SSD in various species. Chinese tongue soles (Cynoglossus semilaevis) at 4 months were fed with 17β-estradiol (E2) and testosterone (T) supplemented feeds for 8 months to assess the effect of sex steroids on growth traits in different sexes. The potential genetic regulation was examined using several growth-related genes. The results showed that two sex steroid hormones had inhibitory effects on the growth performance of different sexes of C. semilaevis. At the age of 8 months, the expression of insulin-like growth factor 2 gene (igf2), 24-dehydrocholesterol reductase (dhcr24), leptin, and estrogen receptor 2 (esr2) in the liver showed an overall downward trend. The expression of insulin-like growth factor 1 (igf1) was reduced, while thyroid hormone receptor-associated protein 3 (thrap3) expression tended to increase in the gonad after T and E2 treatments. In the brain, somatostatin 1, tandem duplicate 2 (sst1.2) expression increased with the treatment of T and E2 (P < 0.05), while growth hormone-releasing hormone (ghrh) expression decreased. E2 and T had different effects on growth differentiation factor 8 (gdf8) and insulin-like growth factor-binding protein 7 (igfbp7) expression in the muscle. Expression of gdf8 increased in the treated fishes in contrast to the reduction expression of igfbp7. This study provided important clues for understanding the role of sex steroids in flatfish SSD., (© 2022. The Author(s), under exclusive licence to Springer Nature B.V.)

Turbot (Scophthalmus maximus), commercially important flatfish species, is widely cultivated in Europe and China. With the continuous expansion of the intensive breeding scale, turbot is exposed to various stresses, which greatly impedes the healthy development of turbot industry. Here, we present an improved high-quality chromosome-scale genome assembly of turbot using a combination of PacBio long-read and Illumina short-read sequencing technologies. The genome assembly spans 538.22 Mb comprising 27 contigs with a contig N50 size of 25.76 Mb. Annotation of the genome assembly identified 104.45 Mb repetitive sequences, 22,442 protein-coding genes and 3,345 ncRNAs. Moreover, a total of 345 stress responsive candidate genes were identified by gene co-expression network analysis based on 14 published stress-related RNA-seq datasets consisting of 165 samples. Significantly improved genome assembly and stress-related candidate gene pool will provide valuable resources for further research on turbot functional genome and stress response mechanism, as well as theoretical support for the development of molecular breeding technology for resistant turbot varieties., (© 2022. The Author(s).)

Sexual size dimorphism (SSD), whereby females and males exhibit different body sizes, are widely documented in animals. To explore crucial regulators implicated in female-biased SSD of Chinese tongue sole (Cynoglossus semilaevis), GWAS was conducted on 350 females and 59 males. Twenty SNPs and 25 genes including zbed1, nsd3, cdc45, klhl29, and smad4 with -log(p) > 7 were screened, mainly mapping to sex chromosome. The chromosome W-linked gene zbed1 attracted particular attention because it is a master key for cell proliferation. Thus, the regulatory network of zbed1 in C. semilaevis was explored by DAP-seq and 1352 peaks were discovered in the female brain. Moreover, zbed1 potentially regulated hippo signaling pathway, cell cycle, translation, and PI3k-Akt signaling pathway in C. semilaevis. These findings identify crucial SNPs and genes associated with female-biased SSD in C. semilaevis, also provide the first genome-wide survey for the zbed1 regulatory network in fish species., Competing Interests: Declaration of Competing Interest We declare that we do not have any commercial or associative interest that represents a conflict of interest in connection with the work submitted., (Copyright © 2022 The Authors. Published by Elsevier Inc. All rights reserved.)

As intracellular pathogen recognition receptors (PRRs), nucleotide-binding domain, leucine-rich repeat containing receptors (NLRs, NOD-like receptors) are involved in innate immune responses in vertebrates. However, there is no systemic study on NLRs in Chinese tongue sole (Cynoglossus semilaevis), a popular maricultured fish in China. In the present study, a genome-wide survey of NLRs was performed in C. semilaevis, with the identification of 29 NLRs, including five genes from the NLR-A subfamily (referred to as CsNOD1-5), two genes from the NLR-B subfamily, 18 genes from the NLR-C subfamily (referred to as CsNLR-C1 to 18) and four other NLR genes. Phylogenetic analysis implied that CsNOD1-5 contained conserved functional domains and had orthologous relationships with human NOD1-5. Moreover, CsNLR-C genes all possessed the FISHNA domain, which is a fish-specific NACHT subdomain. Expression analysis showed that CsNOD1-5 and CsNLR-C1/2 were ubiquitously expressed in various normal tissues. Bacterial infection with Vibro harveyi revealed distinct expression patterns of all the tested CsNLRs in gill, intestine, trunk kidney, liver and spleen. In particular, CsNOD1-4 and CsNLR-C2 were significantly upregulated in gills at 48 h post bacterial infection. In addition, CsNOD3 and CsNOD4 were significantly elevated in infectious intestine, trunk kidney, liver and spleen, revealing that their expressions were more sensitive to bacterial infection than other CsNLRs. Together with the computational protein-protein interaction network of CsNLRs, it was suggested that individual NLR genes had different roles in the innate immune cascades of C. semilaevi against bacterial infection. This study provides valuable information for further studies on CsNLR immune function., (© 2021 Fisheries Society of the British Isles.)

Sex differentiation in aquatic fish is important both for theoretical study and practical production, as growth dimorphism frequently appears in different sexes, especially in marine fish. The deciphered genome, identification of the male-determining gene dmrt1 and established genotypic sex screening method make Chinese tongue sole (Cynoglossus semilaevis) an ideal model to study sex differentiation in fish. In this study, comparative gonadal transcriptomic analyses were conducted for genetic females and males at 48, 68, and 108 days post hatching (dph), representing pre-, during- and post-gonadal differentiation stages, although the gonad is not completely differentiated and isolable in 48 and 68 dph individuals, while it is in 108 dph individuals. Altogether, 28 libraries were constructed, and a mean of 46.64 M clean reads was obtained. Differentially expressed gene (DEG) analysis revealed that 179 genes had similar expression patterns in males and females in all three stages. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis suggested that the enriched pathways included ubiquitin-mediated proteolysis, lysosomes, and RNA degradation. Moreover, weighted gene coexpression network analyses (WGCNA) identified 14 modules, one of which was closely correlated with female differentiation, exhibiting female-biased expression in all three stages (48, 68, 108 dph). An illustrated core gene interaction network of this module identified 50 genes, most of which are on W chromosomes. Six genes, including two ubiquitin conjugating enzymes, were selected for further investigation, and their female-biased expression was confirmed in even earlier stages, at 10 and 30 dph. These data facilitate our understanding of sex differentiation in fish and provide a genomic rationale for screening candidate genes (preferentially W-linked genes) that could be involved in the female differentiation process., (Copyright © 2021. Published by Elsevier Inc.)

The sexual size dimorphism of the Chinese tongue sole (Cynoglossus semilaevis) has greatly obstructed its sustainable development; however, the underlying mechanism remains unclear. Based on C. semilaevis transcriptomic information, 24-dehydrocholesterol reductase (dhcr24) was identified in steroid biosynthesis, showing female-liver-biased expression. Dhcr24 has been reported to participate in various processes, such as cholesterol synthesis, oxidative stress response, neuroprotection, and cell survival. The present study assessed its role in the sexual size dimorphism in fish. First, detailed expression pattern analysis showed that dhcr24 mRNAs were extensively expressed in tissues and the highest levels were found in the liver and gonads of females. Analysis of the dhcr24 promoter region demonstrated different DNA methylation statuses in female, male, and pseudomale gonads with higher epigenetic modification in males. The confirmation of transcription activity of the dhcr24 promoter and putative transcription factors (e.g., ER, AR, SREBP, and POU1F1a) provides the foundation for studying its regulatory mechanism. Finally, dhcr24-siRNA mediated knock-down assay using C. semilaevis liver cells showed that steroid biosynthesis related genes (e.g., ebp, dhcr7, and sc5d), core component of PI3K/Akt pathway (e.g., pi3k), and igf1r exhibited different expression patterns. Further investigation on the interplay between steroid hormones, dhcr24, PI3K/Akt, and IGF-1 systems will be valuable to better understand the mechanism underlying the sexual size dimorphism in C. semilaevis., (Copyright © 2021. Published by Elsevier Inc.)

In addition to the typical sexual size dimorphism, considerable size differences within the female population of the Chinese tongue sole (Cynoglossus semilaevis) have become a further bottleneck of the improvement of sole aquaculture. To identify the internal mechanism, transcriptomic analysis and weighted gene co-expression network analysis (WGCNA) were employed simultaneously. Transcriptomic analyses of brain, pituitary gland, liver, gonad, and muscle tissues from two female groups with size differences identified 109, 698, 1325, 2299, and 2141 differentially expressed genes (DEGs), respectively. The results of these enrichment analyses suggest that the up-regulation of neuroactive ligand-receptor interaction, cell cycle, DNA replication, and MAPK signaling pathway in the group with larger females may be involved in the regulation of the observed growth differences. WGCNA of DEGs showed that cell cycle and DNA replication might be crucial pathways for accelerating cell growth in the groups with larger females. Finally, a series of hub genes including 6-phosphofructokinase type C (pfkp), ribosome biogenesis protein (wdr12), bleomycin hydrolase (blmh), and semaphorin-3A (sema3a) were recognized by the illustrated network map of modules. The linkage of cell cycle, DNA replication, and hub genes in the growth regulation of C. semilaevis provides further information for a better understanding of growth differences in fish., (Copyright © 2021 Elsevier Inc. All rights reserved.)

The pore-forming protein perforin is one of the effectors of cell-mediated killing via the granule exocytosis pathway. In this study, a genome-wide association study was conducted in Vibrio harveyi disease-resistant and disease-susceptible families of half smooth tongue sole (Cynoglossus semilaevis) to determine the genes accounting for host resistance, and a perforin homologue was identified, designated perforin-1 like (CsPRF1l). The full-length cDNA of CsPRF1l is 1835 bp, and encodes 514 amino acids. The CsPRF1l gene consists of 10 exons and 9 introns, spanning approximately 7 kb. The amino acid sequence of CsPRF1l shows 60.35, 54.03, 41.92, and 34.17% identities to Morone saxatilis PRF1l, Oryzias melastigma PRF1l, Danio rerio PRF1.5 and Homo sapiens PRF, respectively. Sequence analysis revealed the presence of membrane attack complex/perforin (MACPF) and C2 domains in CsPRF1l. Quantitative real-time PCR showed that CsPRF1l presented a higher intestinal expression level in disease-resistant families than in susceptible families. Tissue expression pattern analysis showed that CsPRF1l is present in most of the tested tissues and highly expressed in the intestine, brain, stomach and gills. After challenge with V. harveyi, CsPRF1l mRNA was markedly upregulated in the liver, spleen, kidney, intestine, gills and skin. In addition, the recombinant CsPRF1l protein exhibited obvious antimicrobial activity against V. harveyi in vitro and in a zebrafish model. Collectively, these data indicate that CsPRF1l modulates host immune defense against V. harveyi invasion and provide clues about the efficacy of rCsPRF1l in fish that will give rise to useful therapeutic applications for V. harveyi infection in C. semilaevis., (Copyright © 2021 Elsevier Ltd. All rights reserved.)

In mammals, Class II, major histocompatibility complex (MHC II) transactivator (CIITA) recognizes microbial pathogens and triggers immune responses. In Chinese tongue sole Cynoglossus semilaevis, Cs-CIITA was prevalently expressed in various tissues. Cs-CIITA, Cs-MHC IIA and Cs-MHC IIB were expressed significantly higher in skin in susceptible families infected with Vibrio harveyi, while higher expression of Cs-CIITA and Cs-MHC IIB was examined in liver in resistant families. In addition, the three genes were up-regulated in gill, skin, intestine, liver, spleen and kidney at 48 h or 72 h after V. harveyi infection. Furthermore, the three genes were co-expressed in the epithelial mucous cells of gill, skin, and intestine. Knockdown of Cs-CIITA regulates the expression of other inflammation-related genes, including CD40, IL-1β, IL-8, RelB, NFκB, and Myd88. These results suggest that CIITA functions in the inflammatory responses of C. semilaevis against V. harveyi, via MHC II transcriptional regulation., (Copyright © 2021 Elsevier Ltd. All rights reserved.)

Based on 1572 re-sequenced Chinese tongue sole (Cynoglossus semilaevis), we investigated the accuracy of four genomic methods at predicting genomic estimated breeding values (GEBVs) of Vibrio harveyi resistance in C. semilaevis when SNPs varying from 500 to 500 k. All methods outperformed the pedigree-based best linear unbiased prediction when SNPs reached 50 k or more. Then, we developed an SNP array "Solechip No.1" for C. semilaevis breeding using the Affymetrix Axiom technology. This array contains 38,295 SNPs with an average of 10.5 kb inter-spacing between two adjacent SNPs. We selected 44 candidates as the parents of 23 families and genotyped them by the array. The challenge survival rates of offspring families had a correlation of 0.706 with the mid-parental GEBVs. This SNP array is a convenient and reliable tool in genotyping, which could be used for improving V. harveyi resistance in C. semilaevis coupled with the genomic selection methods., (Copyright © 2021 Elsevier Inc. All rights reserved.)

Vibrio harveyi is a major bacterial pathogen that causes fatal vibriosis in Chinese tongue sole (Cynoglossus semilaevis), resulting in massive mortality in the farming industry. However, the molecular mechanisms of C. semilaevis response to V. harveyi infection are poorly understood. Here, we performed transcriptomic analysis of C. semilaevis, comparing resistant and susceptible families in response to V. harveyi challenge (CsRC and CsSC) and control conditions (CsRU and CsSU). RNA libraries were constructed using 12 RNA samples isolated from three biological replicates of the four groups. We performed transcriptome sequencing on an Illumina HiSeq platform, and generating a total of 1,095 million paired-end reads, with the number of clean reads per library ranging from 75.27 M to 99.97 M. Through pairwise comparisons among the four groups, we identified 713 genes exhibiting significant differences at the transcript level. Furthermore, the expression levels were validated by real-time qPCR. Our results provide a valuable resource and new insights into the immune response to V. harveyi infection.

Lhx9 is a LIM-homeodomain protein related to gonad development and sex reversal. In this study, we cloned and characterized CS-Lhx9 in the gonads of the Chinese tongue sole (Cynoglossus semilaevis). The full-length cDNA of CS-Lhx9 was 3123 bp, including an ORF of 1149 bp encoding 383 amino acids which contains two LIM domains and a homeobox domain. CS-Lhx9 transcripts were primarily observed in the testis of male and neomale at 1 yah, but nearly undetectable in the ovary. During the development of gonad, CS-Lhx9 exhibited an increasing trend and appeared to reach its peak value of expression in testis at 1 yah. In situ hybridization was performed in male and neomale gonads at 210 dah and 1 yah. The results showed strong CS-Lhx9 signals in the spermatids and spermatozoa (germ cells). Methylation level of CS-Lhx9 promoter was higher in female and lower in male and neomale gonads, showing a negative correlation with CS-Lhx9 expression., (Copyright © 2018. Published by Elsevier B.V.)

Zinc finger GATA like protein-1 (ZGLP1) is a nuclear zinc finger protein that regulates the interaction between somatic cells and germ cells during gonad developmental process in mammals. In this study, the zglp1 of Chinese tongue sole, Cynoglossus semilaevis (cysezglp1), was cloned and characterized for the first time in fish. Cysezglp1 had an open reading frame with five exons and was located to chromosome 9. The open reading frame of cysezglp1 consisted of 1692 nucleotides and encoded a 583 amino acid polypeptide. The predicted protein contained two zinc finger structures (Znf1 and Znf2), one of which was highly homologous to the GATA-type zinc finger domain. Multiple sequence alignment showed that Znf1 was conserved across different species while Znf2 was more divergent. Through quantitative Real-time PCR (qRT-PCR), we found that cysezglp1 was predominantly expressed in gonads, and the expression level of the ovary was significantly higher than that of the testis. We compared expression level in different embryonic stages and found that cysezglp1 mRNAs were mainly expressed in the fertilized egg to the cleavage stage, subsequently declining in the blastula stage. Cysezglp1 expression was not detected from the gastrulation stage onward. In the ovary, cysezglp1 expression was detected at 120 days after hatching and expression gradually increased with the maturation of the ovary. In situ hybridization showed that the cysezglp1 was mainly expressed in oocytes. Taken together, our results suggest that cysezglp1 may play an important role in the process of oogenesis in Chinese tongue sole., (Copyright © 2018 Elsevier B.V. All rights reserved.)

Rac1 and Rac2, belonging to the small Rho GTPase family, play an important role during the immune responses. In this study, a Rac1 homolog (CsRac1) and a Rac2 homolog (CsRac2) were cloned from the Cynoglossus semilaevis. The full-length of CsRac1 and CsRac2 cDNA was 1219 bp and 1047 bp, respectively. Both CsRac1 and CsRac2 contain a 579 bp open reading frame (ORF) which encoding a 192 amino acids putative protein. The predicted molecular weight of CsRac1 and CsRac2 was 21.41 kDa and 21.35 kDa, and their theoretical pI was 8.50 and 7.91, respectively. Sequence analysis showed that the conserved RHO domain was detected both from amino acid of CsRac1 and CsRac2. Homologous analysis showed that CsRac1 and CsRac2 share high conservation with other counterparts from different species. The CsRac1 and CsRac2 transcript showed wide tissue distribution, in which CsRac1 and CsRac2 exhibit the highest expression level in liver and gill, respectively. The expression level of CsRac1 and CsRac2 fluctuated in the liver and gill tissues at different time points after challenged by Vibrio harveyi. Specifically, CsRac1 and CsRac2 were significantly up-regulated at 48 h and 96 h post injection. Moreover, the knocking down of CsRac1 and CsRac2 in cell line (TSHKC) reduced the expression of CsPAK1, CsIL1-β and CsTNF-α. The present data suggests that CsRac1 and CsRac2 might play important roles in the innate immunity of half-smooth tongue sole., (Copyright © 2018. Published by Elsevier Ltd.)

Receptor interacting protein kinase 3 (RIP3) has emerged as an essential molecular regulator in human inflammation with accumulating evidence in vertebrates. However, the structure and functions of RIP3 in teleosts remains elusive. Here we identified one RIP3 homologue from half-smooth tongue sole Cynoglossus semilaevis, designated CsRIP3. The open reading frame (ORF) of CsRIP3 contained 1491 nucleotides and encoded a protein with 496 amino acids (Mw = 55.44 kDa). CsRIP3 expression was detected in various tissues in half-smooth tongue soles. CsRIP3 expression was up-regulated after pathogens i.e. Vibrio and poly (i:c) treatment in vivo, indicating its possible role in teleost immune response. Further analysis using human cells demonstrated that CsRIP3 exhibited highly conserved pro-apoptotic and pro-necroptotic functions in comparison with human RIP3. In conclusion, for the first time we reported the homologous structure and functions of CsRIP3 in inflammatory cell death, which provides novel perspectives on fish immunity studies in future., (Copyright © 2018 Elsevier Ltd. All rights reserved.)

Spermatid perinuclear RNA-binding protein (Strbp) is a kind of double-stranded (ds) RNA specific binding protein that plays important roles in mammalian spermatogenesis. In this study, we have isolated and characterized the strbp gene of Chinese tongue sole, Cynoglossus semilaevis, termed as CS-strbp. The CS-strbp genomic sequence contained fifteen exons and fourteen introns. Its cDNA was 2655 bp in length, encoding a 666-amino-acid protein with two conserved ds binding motifs. Using quantitative PCR, we found that CS-strbp mRNA exhibited sex-biased and tissue-specific distribution, predominantly expressed in the fertile male testis, though the expression levels varied throughout different developmental stages. Comparison of methylation profile in different sexual genotypes demonstrated the low methylation level of CS-strbp promoter in male and pseudo-male, which is consistent with the high expression levels in those genotypes. In situ hybridization revealed that CS-strbp mRNA mainly localized in male germ cells, especially in spermatids and spermatozoa. Given these findings, we postulate that CS-strbp might function in spermatogenesis of Chinese tongue sole., (Copyright © 2018 Elsevier Inc. All rights reserved.)

Spexin (SPX), a novel neuropeptide discovered by the bioinformatics approach, has been shown to exert pleiotropic functions in mammals. However, little information regarding the physiological role of SPX is available in teleosts. As a first step, we cloned the spexin gene from a flatfish, the half-smooth tongue sole. The open reading frame (ORF) of tongue sole spexin contained 363 nucleotides encoding a 120 amino acid (aa) preprohormone with a calculated molecular mass and isoelectric point of 14.06 kDa and 5.86, respectively. The tongue sole SPX precursor contained a 27 aa signal peptide and a 14 aa mature peptide flanked by two dibasic protein cleavage sites (RR and GRR). Tissue distribution analysis showed that spexin mRNA could be detected in various tissues, notably in the brain. In addition, fasting stimulated the hypothalamic expression of spexin mRNA. Intraperitoneal injection of SPX increased gnih and gnrh3 mRNA levels in the hypothalamus; however, SPX inhibited the pituitary expression of gh, fshβ, and gthα mRNAs. Overall, our results reveal the existence of a functional SPX in the tongue sole, which could represent an important factor in the neuroendocrine control of flatfish reproduction and growth, and the spexin mRNA expression is regulated by feeding status.

Pod1 encodes a Class II bHLH transcription factor involved in the development of a number of tissues such as gonad, spleen, lungs and heart. However, to date, little is known about its function in teleosts. In this study, we cloned and characterized Pod1 gene from Cynoglossus semilaevis. This gene contains three exons and two introns, with the full-length cDNA of 918 nucleotides that encodes a 183 amino acid protein with a conserved bHLH domain. Realtime quantitative PCR revealed that Pod1 was predominantly expressed in the testes of C. semilaevis. In different stages of testes development, Pod1 expression was undetectable up to 120 days after hatching (dah), and then increased at 210 dah and 1 year after hatching (yah). Furthermore, in situ hybridization (ISH) analysis revealed that Pod1 was mainly localized in the germ cells of testes, but was not detected in ovarian cells; which suggested its possible functions in spermatogenesis of C. semilaevis. The methylation profile analysis of Pod1 genomic sequence in the gonads showed that the differences in their putative promoter regions of Pod1 among ovary, male and pseudo-male testes were not obvious. Thus, further research might be needed to evaluate whether Pod1 expression is regulated by epigenetic level., (Copyright © 2017. Published by Elsevier Inc.)

As an important economic marine species cultured in China, Chinese tongue sole (Cynoglossus semilaevis) has interested us due to its sexual dimorphism and ZW/ZZ sex determination system. In a previous study, dmrt1 was identified as a dosage-dependent male-determining gene. In the present study, a female-specific expressed gene, cse0440, initially annotated as lrp1b-like, was identified from chromosome W of C. semilaevis. In view of the differences between cse0440 and lrp1b in terms of expression pattern, a phylogenetic analysis containing 85 LRP proteins was constructed and provided an evidence to re-annotate cse0440 as cseLRP13. In addition, two orthologues of cseLRP13 were separately identified from W and Z chromosomes: cseLRP13-W and cseLRP13-Z. The subsequent multiple sequence alignment and syntenic arrangements of LRP13 in C. semilaevis, Japanese medaka (Oryzias latipes), large yellow croaker (Larimichthys crocea), striped bass (Morone saxatilis), white perch (Morone americana) and Fugu rubripes (Takifugu rubripes) further supported this re-annotation. RT-PCR and in situ hybridization revealed that cselrp13 was exclusively expressed in the oocytes and follicles of ovaries. These results suggested that lrp13 may play important roles in female reproduction. In future, with the advancement of micromanipulation in flatfish, the detailed function of two lrp13 orthologues in C. semilaevis will be elucidated.

The liver is an important organ for bacterial pathogen attack in fish. The differential proteomic response of the Japanese flounder liver to Edwardsiella tarda infection was examined using isobaric tags for relative and absolute quantitation (iTRAQ) labeling followed by liquid chromatography-tandem mass spectrometry (LC-MS/MS). A total of 3290 proteins were identified and classified into categories related to biological process (51.4%), molecular function (63.6%), and cellular component (57.7%). KEGG enrichment analysis indicated the complement and coagulation cascade pathways and the mineral absorption pathway were significantly enriched. Among the differentially expressed proteins, those involved in mediating complement cascade (e.g. complement component C7, C8, C9, complement factor H, complement factor Bf/C2) and mineral absorption (e.g. ferritin, STEAP-4) were most significantly upregulated during infection. Subsequently, five significantly upregulated (C4, C8beta, ferritin middle subunit, PRDX4-like and KRT18) and one significantly downregulated (transferrin) candidate immune proteins were validated by multiple reaction monitoring (MRM) assays. Furthermore, changes in expression of 15 proteins in the complement cascade and mineral absorption pathways were validated at the transcriptional level using quantitative real-time PCR (qPCR). The transcriptional levels of four transcription factors (p21Ras, Rab-31-like, NF-κB, STAT3) were also investigated by qPCR following infection with E. tarda. This study contributes to understanding the defense mechanisms of the liver in fish., (Copyright © 2017 Elsevier Ltd. All rights reserved.)

In the present study, the peripheral blood cells of half smooth tongue sole (Cynoglossus semilaevis) were examined by blood smear under the light microscopy. The proportion of main types of blood cells are as following: erythrocyte occupied the majority (92.3%), followed by thrombocyte (4.15%), granulocyte (1.7%), lymphocyte (1.5%) and monocyte (0.3%), respectively. Meanwhile, the isolation method of monocytes was established, by density gradient centrifugation to isolate mononuclear leukocytes of peripheral blood. In primary culture, the monocytes were adhered to the bottom of the flask without feeder cells and separated easily with suspended leukocytes in the medium in 3 h. After suspended leukocytes were removed, the monocytes multiplied rapidly with the two doubly during the 24 h, then the cells proliferated and kept stable until 48 h. When co-cultured with suspended leukocytes after three days, the monocytes could derive to typical macrophages, of which the size enlarged significantly and showed various forms such as like fried eggs, and giant irregular shape with pseudopod because cells fusion or deformation occurred until macrophages died in about two weeks. Monocytes showed strong respiratory burst activity after treated with Phorbol ester PMA and challenged by bacteria respectively. In addition, macrophage of half smooth tongue sole had typical macrophage features such as phagocytic capability, positive esterase activity, and the considerable expression of M-CSFR, MHC-II, IL-6, IL-10, TNF and arginase genes. That arginase expression in macrophages (3d and 5d after differentiation) was upregulated fluctuant suggest that the cultivation was mixture of alternatively activated type macrophage (M2) in the majority while the classically activated type (M1) win the minority. Furthermore, MHC-Ⅱ, M-CSFR and IL-6 were significantly induced following LPS challenge. Collectively, the present study will be useful for the study on half smooth tongue sole immune systems and immune function., (Copyright © 2017 Elsevier Ltd. All rights reserved.)

β-catenin is a key signalling molecule in the canonical Wnt pathway, which plays a role in cell adhesion, embryogenesis and sex determination. However, little is known about its function in teleosts. We cloned and characterized the full-length β-catenin1 gene from half-smooth tongue sole (Cynoglossus semilaevis), which was designated CS-β-catenin1. The CS-β-catenin1 cDNA consists of 2,346 nucleotides and encodes a protein with 782 amino acids. Although CS-β-catenin1 was transcribed in the gonads of both sexes, the level was significantly higher in ovaries compared to testes. Furthermore, the mRNA level of CS-β-catenin1 was significantly upregulated at 160 days and constantly increased until 2 years of age. In situ hybridization revealed that CS-β-catenin1 mRNA was mainly localized in oocyte cells, especially in stage I, II and III oocytes. When CS-β-catenin1 expression was inhibited by injection of quercetin in the ovaries, levels of CS-Figla and CS-foxl2 mRNA were significantly down-regulated, and CS-dmrt1 was up-regulated, which suggested that CS-β-catenin1 is a potential upstream gene of CS-Figla and is involved in the development of the ovaries, i.e., folliculogenesis.

Aquaporin 1 (AQP1) is a member of the transmembrane water channel family of proteins with special structural features, and two AQP1 paralogous genes (aqp1aa and aqp1ab) are reported in teleosts. In the present study, the aqp1aa gene of half-smooth tongue sole (Cynoglossus semilaevis) was cloned and characterized. The full-length cDNA of aqp1aa is 1411 bp with a 786 bp open reading frame encoding a 261-amino acid putative protein with a characteristic structure consisting of 6 membrane-spanning α-helical domains and two highly conserved asparagine-proline-alanine motifs. Real-time quantitative PCR revealed that aqp1aa mRNA is expressed predominantly in the testis of males and pseudo-males, while its expression is low in the ovary and lowest in doublesex and mab-3-related transcription factor 1(DMRT1) knock out fish and triploid males. In situ hybridization indicated that aqp1aa mRNA is expressed mainly in the germ cells of males and pseudo-males, especially in spermatozoa and spermatids. These results suggest that the aqp1aa may play a role in spermatogenesis of C. semilaevis.

Mitochondrial genomes can provide basic information for phylogenetic analysis and evolutionary studies. We present here the mitochondrial genome of Cynoglossus gracilis, which is 16,565bp in length. Numerous distinct regions were identified, including 13 protein-coding genes (PCGs), 22 tRNA genes, two rRNA genes, a light-strand replication origin, and a control region. Interestingly, we detected rearrangement of genes in C. gracilis, including a control region translocation, tRNA(Gln) gene inversion, and tRNA(Ile) gene shuffling. Additionally, a phylogenetic analysis based on the nucleotide sequences of the 13 PCGs using maximum likelihood and Bayesian inference methods reveals that C. gracilis is closely related to Cynoglossus semilaevis. This study provides important mitogenomic data for analyzing phylogenetic relationships in the Cynoglossinae., (Copyright © 2016 Elsevier B.V. All rights reserved.)

Thermal stress, encountered frequently in aquaculture, affects diverse physiological processes and behavior of fish; however, the mechanism underlying these effects may vary among species. Half-smooth tongue sole (Cynoglossus semilaevis) widely inhabits Asian coastal waters and has been intensively cultured in China. In this study, the transcriptomes of three tissues (gill, liver and muscle) of half-smooth tongue sole stressed at temperatures up to 35°C were profiled and compared with those of the fish living at a normal temperature of 25°C. Of 105,228 transcripts, 23,213 were annotated in GO terms, and 718, 1236 and 561 were found to express differentially in gill, liver and muscle from control, respectively. Of the differentially expressed genes, 119 were shared by all three tissues. Tissue-overlapping transcripts and the pathways and functions they defined may underline the primary response of fish to thermostressing. High temperature may cause directly protein mis-folding and hypoxia. Cells suffering thermostress either survive hypoxia or are eliminated by immune system. High temperature may have interacted indirectly through HSPs with HIF-1 and JAKs/STATs signaling pathways. The former up-regulates the expression of hypoxia inducible genes while the later up-regulates the expression of genes associating with insulin and NOD-like receptors signaling pathways. Insulin pathway functions to sustain homeostasis of sugar and lipids, aiding to survive cells, while NOD-like receptors signaling pathway strengthens immune, apoptotic and inflammatory responses, helping to survive cells. These understandings may help improve our culturing practice, the culture performance of half-smooth tongue sole and its breeding., (Copyright © 2016 Elsevier B.V. All rights reserved.)

Testis-specific protein kinase 1 (tesk1) represents a conserved gene family functioning in many cellular processes. In this study, we cloned and characterized an autosome-localized tesk1 gene (Altesk1) from Cynoglossus semilaevis. The open reading frame consists of 2088 nucleotides and encodes a 665 amino acid polypeptide. Phylogenetic analyses show that vertebrate Tesk1s are divided into two clusters based on protein length and AlTesk1 belongs to "long-type" group. Semi-quantitative PCR reveals that Altesk1 is predominantly expressed in ovary, despite of relatively low detection in some other tissues. Among different development stages, Altesk1 transcripts are only observed in ovary samples of 210-day and 1-year fish. In situ hybridization analyses have further confirmed its major localization in oocyte cells. Comparison of methylation patterns in different sexual genotypes reveals the low methylation level of Altesk1 promoter in female, which is consistent with Altesk1 high expression level in female. Taken together, this is the first time that tesk1 gene has been found to show female-biased expression and in view of this, we postulate that AlTesk1 might be involved in some cellular processes specific in ovary, e.g. oogenesis., (Copyright © 2016 Elsevier B.V. All rights reserved.)

The pituitary adenylate cyclase activating polypeptide (PACAP) is a new type of hypophysiotropic hormone and plays an important role in regulating the synthesis and secretion of growth hormone and gonadotropin. The research on the relationship between PACAP and different growth traits would contribute to explain its function during the process of growth. Moreover, epigenetic modifications, especially DNA methylation at the CpG sites of the SNPs, play important roles in regulating gene expression. The results suggest that a SNP mutation (c.C151G) in the PACAP gene of male half smooth tongue sole (Cynoglossus semilaevis) is significantly associated with growth traits and serum physiological and biochemical parameters such as inorganic phosphorus (P < 0.05). The SNP is located in a CpG-rich region of exon 1. Intriguingly, the transition (C→G) added a new methylation site of PACAP gene. This SNP was also significantly related to the expression and methylation level of PACAP (P < 0.05). Individuals with GG genotype had faster growth rates than those of CG and CC genotypes. Moreover, GG genotype had significantly higher PACAP expression level and lower methylation level than CG and CC genotypes. In the serum indexes, only inorganic phosphorus content within GG genotypes was significantly higher than CC genotypes. This implied that the mutation and methylation status of PACAP gene could influence growth traits and this locus could be considered as a candidate genetic or epigenetic marker for Cynoglossus semilaevis molecular breeding.

Half-smooth tongue sole (Cynoglossus semilaevis) is an important marine commercial fish species in China, which suffers from widespread disease outbreaks. Recently, in this regard, our group identified immune-related microRNAs (miRNAs) of C. semilaevis following Vibrio anguillarum infection. Furthermore, miRNA microarray was utilized to characterize the immune roles of important miRNA candidates in response to bacterial infection. Therefore, in the present study, we characterized miR-221 and miR-222 and profiled their expression after challenge. Here, miR-221 and miR-222 precursors were predicted to have a typical hairpin structure. Both miRNAs were expressed in a broad range of tissues in C. semilaevis, while miR-221 and miR-222 were significantly differentially expressed in the immune tissues of C. semilaevis among three small RNA libraries [control group (CG), bacteria-challenged fish without obvious symptoms of infection (NOSG), and bacteria-challenged fish with obvious symptoms of infection (HOSG)]. In order to further characterize and understand the immune response of miR-221 and miR-222, therefore, we profiled miR-221 and miR-222 expression in selected immune tissues after challenge with V. anguillarum. Both miR-221 and miR-222 were upregulated in the liver and spleen, while different expression patterns were observed in the head kidney. In addition, in half-smooth tongue sole head kidney cell line after challenge with lipopolysaccharide (LPS), polyinosinic:polycytidylic acid (poly I:C), peptidoglycan (PGN), and red-spotted grouper nervous necrosis virus (RGNNV), both miR-221 and miR-222 showed significant difference in expression response to pathogen. Meanwhile, the target gene of miR-221 and miR-222 was predicted, which indicated that tumor necrosis factor receptor-associated factor 6 (TRAF6) and interleukin-1 beta (IL-1β) were the target genes of miR-221 and miR-222, respectively. Collectively, these findings indicated that miR-221 and miR-222 have putative roles in innate immune response during C. semilaevis exposure to pathogens. Our findings could expand the knowledge of immune function of C. semilaevis miRNA and guide future studies on C. semilaevis immunity.

Complement factor I (Cfi) is a soluble serine protease which plays a crucial role in the modulation of complement cascades. In the presence of substrate modulating cofactors (such as complement factor H, C4bp, CR1, etc), Cfi cleaves and inactivates C3b and C4b, thereby controlling the complement-mediated processes. In this study, we sequenced and characterized Cfi gene from Cynoglossus Semilaevis (designated as CsCfi) for the first time. The full-length cDNA of CsCfi was 2230 bp in length, including a 98 bp 5'-untranslated region (UTR), a 164 bp 3'-UTR and a 1968 bp open reading frame (ORF). It encoded a polypeptide of 656 amino acids, with a molecular mass of 72.28 kDa and an isoelectric point of 7.71. A signal peptide was defined at N-terminus, resulting in a 626-residue mature protein. Multiple sequence alignment revealed that Cfi proteins were well conserved with the typical modular architecture and identical active sites throughout the vertebrates, which suggested the conserved function of Cfi. Phylogenetic analysis indicated that CsCfi and the homologous Cfi sequences from teleosts clustered into a clade, separating from another clade from the cartilaginous fish and other vertebrates. Tissue expression profile analysis by quantitative real-time PCR (qRT-PCR) showed that CsCfi mRNA constitutively expressed in all tested tissues, with the predominant expression in liver and the lowest in stomach. Temporal expression levels of CsCfi after challenging with Vibrio anguillarum showed different expression patterns in intestine, spleen, skin, blood, head kidney and liver. The recombinant CsCfi (rCsCfi) protein showed broad-spectrum antimicrobial activities against the Gram-positive bacteria Staphylococcus aureus and the Gram-negative bacteria Escherichia coli, Pseudomonas aeruginosa and Shewanella putrefaciens. The research revealed that CsCfi plays an important role in C. Semilaevis immunity., (Copyright © 2015 Elsevier Ltd. All rights reserved.)

To investigate the roles of microRNAs (miRNA) of Cynoglossus semilaevis in response to Vibrio anguillarum that were previously identified using high-throughput sequencing, microarray analyses was performed on three small RNA libraries (CG, NOSG, and HOSG) prepared from C. semilaevis immune tissues. In total, of 1279 designed probes, 739 (57.78 %) were detectable. The expression levels of these miRNAs were analyzed using pairwise comparisons among the three libraries, and a total of 99 miRNAs were observed to be significantly differentially expressed. The expression patterns of 10 differentially expressed miRNAs were validated by real-time quantitative PCR (RT-qPCR). In addition, expression of miR-142-5p, miR-223, and miR-181a in response to V. anguillarum at numerous time-points in four tissues, as well as the responses to lipopolysaccharide (LPS), polyinosinic:polycytidylic acid (poly I:C), peptidoglycan (PGN), and red-spotted grouper nervous necrosis virus (RGNNV) in head kidney cells, were studied by qRT-PCR. Taken together, all of the expression profiles showed significant differences compared to the control group; both similarities and differences in the expression responses to the same pathogen were observed. Collectively, these findings highlighted the putative roles for miRNAs in the context of the innate immune response of C. semilaevis exposing to pathogens and that further studies are needed to understand the molecular mechanisms of miRNA regulation in C. semilaevis host-pathogen interactions.

Craniofacial remodeling during flatfish metamorphosis, including eye migration, is perhaps the most striking example of asymmetric postembryonic development in the vertebrate world. The asymmetry of the cranium mainly results from distortion of the frontal bones, which depends on eye migration during metamorphosis. However, it is unclear how the up-migrating eye causes distortion of the frontal bones. In this study, we first show that distortion of the frontal bones during metamorphosis in Paralichthys olivaceus is the result of cell apoptosis, rather than cell autophagy or cell proliferation. Secondly, we report that cell apoptosis in the frontal bones is induced by the mechanical force transferred from the up-migrating eye. The mechanical force from the up-migrating eye signals through FAK to downstream molecules that are integrated into the BMP-2 signal pathway. Finally, it is shown that cell apoptosis in the frontal bones is activated by the intrinsic mitochondrial pathway; the extrinsic death receptor is not involved in this process. Moreover, cell apoptosis in frontal bones is not induced directly by thyroid hormones, which are thought to mediate metamorphosis in flatfishes and directly mediate cell apoptosis during amphibian metamorphosis. These findings help identify the major signaling route used during regulation of frontal bone distortion during metamorphosis in flatfish, and indicate that the asymmetry of the cranium, or at least the distortion of frontal bones, is the result of rather than the reason underlying eye migration., (Copyright © 2015. Published by Elsevier Ireland Ltd.)

Half-smooth tongue sole (Cynoglossus semilaevis) is one of the most valuable marine aquatic species in Northern China. Given to the rapid development of aquaculture industry, the C. semilaevis was subjected to disease-causing bacteria Vibrio anguillarum. It therefore is indispensable and urgent to understand the mechanism of C. semilaevis host defense against V. anguillarum infection. In the present study, the extensively analysis at the transcriptome level for V. Anguillarum disease in tongue sole was carried out. In total, 94,716 high quality contigs were generated from 75,884,572 clean reads in three libraries (HOSG, NOSG, and CG). 22,746 unigenes were identified when compared with SwissProt, an NR protein database and NT nucleotide database. 954 genes exhibiting the differentially expression at least one pair of comparison in all three libraries were identified. GO enrichment for these genes revealed gene response to biotic stimulus, immune system regulation, and immune response and cytokine production. Further, the pathways such as complement and coagulation cascades and Vibrio cholerae infection pathways were enriched in defensing of pathogen. Besides, 13,428 SSRs and 118,239 SNPs were detected in tongue sole, providing further support for genetic variation and marker-assisted selection in future. In summary, this study identifies several putative immune pathways and candidate genes deserving further investigation in the context of development of therapeutic regimens and lays the foundation for selecting resistant lines of C. semilaevis against V. anguillarum., (Copyright © 2014 Elsevier Ltd. All rights reserved.)

piwi is an important regulator gene in germ cell division during spermatogenesis. piwi homologous genes are involved in gametogenesis and germline specification, and knocking down these genes could affect germ cell meiotic progression. To understand the function of piwi-related genes in spermatogenesis, we cloned a Piwi-subfamily member (piwil2 gene) from the gonad of Cynoglossus semilaevis. The full-length of piwil2 cDNA was 3314bp, including a 3162bp open reading frame (ORF), a 60bp 5'-UTR, along with a 92bp 3'-UTR, and encoded a predicted protein of 1053 amino acid residues. Phylogenetic analysis showed that the PIWIL2 putative protein belonged to the Argonaute protein family, and Piwi-subfamily, with typical PAZ and Piwi domains. Ultrathin sections of different gonadal stages, and real-time quantitative PCR showed that the relative expression of the piwil2 gene couldn't be detected until day 95 (95days) larvae, when germ cell divided rapidly in C. semilaevis. The piwil2 transcripts were more abundant in gonads of males and neo-males than in females, and weakly expressed in other tissues and organs. Compared with the relative expression of piwil2 in gonads, the CpG methylation levels were significantly higher in females. Chromosomal fluorescence in situ hybridization (FISH) showed that the piwil2 gene was located on the Z sex chromosome of C. semilaevis. These results suggest that piwil2 plays an important role in spermatogenesis of C. semilaevis., (Copyright © 2014 Elsevier B.V. All rights reserved.)

The csdazl gene is a sex related gene of half-smooth tongue sole (Cynoglossus semilaevis). Our research group have cloned full length cDNA of csdazl, and studied its expression pattern. To get the further information of csdazl, we constructed a prokaryotic expression plasmid, pET-32a-CSDAZL, expressed in Escherichia coli BL21 and purified the fusion protein by His Trap. In order to detect the biological activity of the fusion protein, we injected the protein with liposome into fish, and detected other sex-related genes' mRNA expression. The results showed that the expression levels of half-smooth tongue sole female-related genes Cyp19a and Foxl2 significantly decreased between 6 and 24 h; however, both genes' expressions returned to their normal levels 72 h after injection, indicating that recombinant CSDAZL protein could down-regulated the expression of female-related genes, Foxl2 and Cyp19a genes, implying that the fusion protein has biological activity and csdazl plays a role in sex differentiation by regulating sex related genes' expression., (Copyright © 2014 Elsevier Inc. All rights reserved.)

Vasa is a highly conserved ATP-dependent RNA helicase expressed mainly in germ cells. The vasa gene plays a crucial role in the development of germ cell lineage and has become an excellent molecular marker in identifying germ cells in teleosts. However, little is known about the structure and function of the vasa gene in flatfish. In this study, the vasa gene (Csvasa) was isolated and characterized in half-smooth tongue sole (Cynoglossus semilaevis), an economically important flatfish in China. In the obtained 6425-bp genomic sequence, 23 exons and 22 introns were identified. The Csvasa gene encodes a 663-amino acid protein, including highly conserved domains of the DEAD-box protein family. The amino acid sequence also shared a high homology with other teleosts. Csvasa expression was mainly restricted to the gonads, with little or no expression in other tissues. Real-time quantitative polymerase chain reaction analysis revealed that Csvasa expression levels decreased during embryonic and early developmental stages and increased with the primordial germ cell proliferation. A typical sexually dimorphic expression pattern of Csvasa was observed during early development and sex differentiation, suggesting that the Csvasa gene might play a differential role in the proliferation and differentiation of male and female primordial germ cells (PGCs). Csvasa mRNA expression levels in neomales were significantly lower than those in normal males and females, indicating that the Csvasa gene might be implicated in germ cell development after sex reversal by temperature treatment. In addition, medaka (Oryzias latipes) PGCs could be transiently labeled by microinjection of synthesized mRNA containing the green fluorescence protein gene and 3'-untranslated region of Csvasa, which confirmed that the Csvasa gene has the potential to be used as a visual molecular marker of germ cells and laid a foundation for manipulation of PGCs in tongue sole reproduction., (Copyright © 2014 Elsevier Inc. All rights reserved.)

Environmental sex determination (ESD) occurs in divergent, phylogenetically unrelated taxa, and in some species, co-occurs with genetic sex determination (GSD) mechanisms. Although epigenetic regulation in response to environmental effects has long been proposed to be associated with ESD, a systemic analysis on epigenetic regulation of ESD is still lacking. Using half-smooth tongue sole (Cynoglossus semilaevis) as a model-a marine fish that has both ZW chromosomal GSD and temperature-dependent ESD-we investigated the role of DNA methylation in transition from GSD to ESD. Comparative analysis of the gonadal DNA methylomes of pseudomale, female, and normal male fish revealed that genes in the sex determination pathways are the major targets of substantial methylation modification during sexual reversal. Methylation modification in pseudomales is globally inherited in their ZW offspring, which can naturally develop into pseudomales without temperature incubation. Transcriptome analysis revealed that dosage compensation occurs in a restricted, methylated cytosine enriched Z chromosomal region in pseudomale testes, achieving equal expression level in normal male testes. In contrast, female-specific W chromosomal genes are suppressed in pseudomales by methylation regulation. We conclude that epigenetic regulation plays multiple crucial roles in sexual reversal of tongue sole fish. We also offer the first clues on the mechanisms behind gene dosage balancing in an organism that undergoes sexual reversal. Finally, we suggest a causal link between the bias sex chromosome assortment in the offspring of a pseudomale family and the transgenerational epigenetic inheritance of sexual reversal in tongue sole fish.

Genetic sex determination by W and Z chromosomes has developed independently in different groups of organisms. To better understand the evolution of sex chromosomes and the plasticity of sex-determination mechanisms, we sequenced the whole genomes of a male (ZZ) and a female (ZW) half-smooth tongue sole (Cynoglossus semilaevis). In addition to insights into adaptation to a benthic lifestyle, we find that the sex chromosomes of these fish are derived from the same ancestral vertebrate protochromosome as the avian W and Z chromosomes. Notably, the same gene on the Z chromosome, dmrt1, which is the male-determining gene in birds, showed convergent evolution of features that are compatible with a similar function in tongue sole. Comparison of the relatively young tongue sole sex chromosomes with those of mammals and birds identified events that occurred during the early phase of sex-chromosome evolution. Pertinent to the current debate about heterogametic sex-chromosome decay, we find that massive gene loss occurred in the wake of sex-chromosome 'birth'.

GRIM-19 (gene associated with retinoid-interferon-induced mortality 19), a novel cell death regulatory gene, plays important roles in cell apoptosis, embryogenesis, mitochondrial respiratory chain and immune response. To date, little information is known about fish GRIM-19 characteristics except orange-spotted grouper (Epinephelus coioides). Here a new GRIM-19 gene is identified and characterized from turbot (Scophthalmus maximus), an economic marine fish in China and Europe. Briefly, turbot GRIM-19 is a 595-bp gene encoding a 144 amino acids protein, which shares the closest relationship with Atlantic halibut (Hippoglossus hippoglossus). The expression of turbot grim-19 in liver, spleen and kidney is up-regulated by the infection of Vibrio anguillarum and LCDV (lymphocystis disease virus). Subsequently, a recombinant protein of turbot GRIM-19 is acquired and the anti-bacterial function is proved by liquid culture inhibition experiment. The subcellular location indicates that turbot GRIM-19 is co-localized with STAT3 in the cytoplasm, which is mainly determined by GRIM-19 41-84 amino acids and STAT3 1-321 amino acids. Finally, the involvements of turbot GRIM-19 in cell apoptosis and NF-κB pathway are investigated. All these data help to understand GRIM-19 function in fish, as well as provide the application possibility of GRIM-19 in fish disease resistance breeding., (Copyright © 2013 Elsevier Ltd. All rights reserved.)

The small ubiquitin-like modifier (SUMO) pathway is an essential biological process in eukaryote, and Ubc9 is an important E2 conjugating enzyme (UBE2) for SUMO pathway and plays a critical role in cellular differentiation, development and sex modification in various species. However, the relationship between Ubc9 and sex modification and development in fish remains elusive. To elucidate the impact of Ubc9 on sex modification and development, the full length of the cDNA and genomic sequence was cloned from half-smooth tongue sole, Cynoglossus semilaevis. Real-time quantitative RT-PCR demonstrated that ubc9 was ubiquitously expressed in different tissues, and the expression levels varied in the different stages of embryonic and gonadal development. In addition, the expression level was significantly higher in the temperature-treated females than the normal females and males. Moreover, the PET-32-Ubc9 plasmid was constructed and the recombinant protein was expressed in Escherichia coli. Follistatin gene expression was initially up-regulated and FSE genes (cyp19a1a, ctnnb1, foxl2) were initially down-regulated after the injection of Ubc9 protein, prior to 96 h eventually recovered to normal levels. Taken together, the results show that Ubc9 is involved in embryogenesis, gametogenesis and sex modification, and exerts an effect on gene expression., (Copyright © 2013 Elsevier Inc. All rights reserved.)