Identification and functional analysis of circpdlim5a generated from pdlim5a gene splicing in the skeletal muscle of Japanese flounder (Paralichthys olivaceus)

Circular RNAs (circRNAs) are non-coding RNAs with endogenous regulatory functions, including regulating skeletal muscle development. However, its role in the development of skeletal muscle in Japanese flounder (Paralichthys olivaceus) is not clear. Therefore we screened a candidate circpdlim5a, which is derived from the gene pdlim5a, from the skeletal muscle transcriptome of Japanese flounder. We characterized circpdlim5a, which was more stable compared to the linear RNA pdlim5a. Distributional characterization of circpdlim5a showed that circpdlim5a was predominantly distributed in the nucleus and was highly expressed in the skeletal muscle of adult Japanese flounder (24 months). When we further studied the circpdlim5a function, we found that it inhibited the expression of proliferation and differentiation genes according to the over-expression experiment of circpdlim5a in myoblasts. We concluded that circpdlim5a may inhibit the proliferation and differentiation of myoblasts and thereby inhibit skeletal muscle development in Japanese flounder. This experiment provides information for the study of circRNAs by identifying circpdlim5a and exploring its function, and offers clues for molecular breeding from an epigenetic perspective.

Functional analysis of the ube3a response in Japanese flounder (Paralichthys olivaceus) to CSBV infection

Ube3a is a member of the E3 ubiquitin ligase HECTc family, and its role has been established in neurodevelopmental disorders. However, studies on its role in Japanese flounder are scarce. Thus, in this study, the ube3a of Japanese flounder was cloned, and its role in conferring resistance against Chinook salmon bafnivirus (CSBV) was analyzed. Japanese flounder ube3a encoded a protein containing 834 amino acids. Interestingly, its homology with the Atlantic halibut was determined to be 94%. In addition, there were differential expressions of ube3a in different tissues of Japanese flounder, with the highest expression level observed in the fin, followed by the gills and skin (P ≤ 0.05). Subcellular localization analysis revealed that Ube3a is a cytoplasmic protein. We established an in vitro CSBV infection model using Japanese flounder gill cell line (FG). After ube3a overexpression, the viral load was significantly lower than that of the control group (P ≤ 0.05). Contrastingly, after incubation of FG cells with an E3 ubiquitin ligase inhibitor, the viral load was significantly higher than in the control group (P ≤ 0.01). Then, the expression levels of nf-κb, traf3, and tnf-α after incubation with an E3 ubiquitin ligase inhibitor were examined. The results demonstrated that ube3a may exerted a significant antiviral effect in Japanese flounder via the ubiquitination pathway.

Transcriptional regulation of the Japanese flounder Cu,Zn-SOD (Jfsod1) gene in RAW264.7 cells during oxidative stress caused by causative bacteria of edwardsiellosis

Edwardsiellosis is one of the most important bacterial diseases in fish, sometimes causing extensive economic losses in the aquaculture industry. Our previous studies demonstrated that the Cu,Zn-SOD (sod1) activity has significantly increased in Japanese flounder, Paralichthys olivaceus, hepatopancreas infected by causative bacteria of edwardsiellosis Edwardsiella tarda NUF251. In this study, NUF251 stimulated intracellular superoxide radical production in mouse macrophage RAW264.7 cells, which was reduced by N-acetylcysteine. This result suggests that NUF251 infection causes oxidative stress. To evaluate the regulatory mechanism of Jfsod1 at transcriptional levels under oxidative stress induced by NUF251 infection, we cloned and determined the nucleotide sequence (1124 bp) of the 5'-flanking region of the Jfsod1 gene. The sequence analysis demonstrated that the binding sites for the transcription factors C/EBPα and NF-IL6 involved in the transcriptional regulation of the mammalian sod1 gene existed. We constructed a luciferase reporter system with the 5'-flanking region (-1124/-1) of the Jfsod1 gene, and a highly increased transcriptional activity of the region was observed in NUF251-infected RAW264.7 cells. Further studies using several mutants indicated that deletion of the recognition region of NF-IL6 (-272/-132) resulted in a significant decrease in the transcriptional activity of the Jfsod1 gene in NUF251-infected RAW264.7 cells. In particular, the binding site (-202/-194) for NF-IL6 might play a major role in upregulating the transcriptional activity of the 5'-flanking region of the Jfsod1 gene in response to oxidative stress induced by NUF251 infection. These results could be provided a new insight to understand the pathogenic mechanism of causative bacteria of edwardsiellosis.

Molecular characterization and ontogenetic expression profiles of LPXRFa and its receptor in Japanese flounder (Paralichthys olivaceus)

Investigations concerning the LPXRFa system are rarely conducted in flatfish species. Here, we first identified and characterized lpxrfa and its cognate receptor lpxrfa-r genes in the Japanese flounder (Paralichthys olivaceus). The coding DNA sequence of lpxrfa was 579 bp in length, wich encoded a 192-aa preprohormone that can produce three mature LPXRFa peptides. The open reading frame (ORF) of lpxrfa-r was 1446 bp in size, and encoded a 481-aa LPXRFa-R protein that encompassed seven hydrophobic transmembrane domains. Subsequently, tissue distribution expression profiles of lpxrfa and lpxrfa-r transcripts were assayed by quantitative real-time PCR. The results indicated that expressions of lpxrfa transcripts were detected at the highest levels in the brain of both females and males, however, lpxrfa-r transcripts were remarkablely expressed in the brain tissue of female fish and in the testis tissue of male fish. Furthermore, transcript levels of lpxrfa and lpxrfa-r genes were investigated during early ontogenetic development, with the maximum expression levels at 30 days post-hatching. Overall, these data contribute to providing preliminary proof for the existence and structure of the LPXRFa system in Japanese flounder, and the study is just the foundation for researching physiological function of LPXRFa system in this species.

Characterization and function of Japanese flounder (Paralichthys olivaceus) slc2a6 in response to lymphocystis disease virus infection

Slc2a6 is a member of the slc2 family (solute carrier 2 family) and previous reports have indicated its involvement in the inflammatory response. Slc2a6 is regulated by the NF-ĸB signaling pathway. This study investigated the differential expression of slc2a6 in the early embryonic development of Japanese flounder, revealing that the early gastrula stage had the highest level of slc2a6 expression. Moreover, slc2a6 expression was increased in vitro after stimulation by lymphocystis disease virus (LCDV), and in vivo experiments also showed significantly elevated levels in the spleen and muscle tissues following LCDV stimulation. Subcellular localization revealed that Slc2a6 was expressed in both the nucleus and cytoplasm of cells. The pcDNA3.1-slc2a6 overexpression plasmid was successfully constructed; the si-slc2a6 interfering strand was screened and samples were collected. The expression of NF-ĸB signaling pathway-related genes il-1β, il-6, nf-ĸb, and tnf-α was evaluated in overexpressed, silenced, and LCDV-stimulated samples. The results showed that slc2a6 is involved in viral regulation in Japanese flounder by regulating innate immune responses.

Stress responses to warming in Japanese flounder (Paralichthys olivaceus) from different environmental scenarios

Japanese flounder (Paralichthys olivaceus) is one of cold-water species widely farmed in Asia. In recent years, the increased frequency of extreme weather events caused by global warming has led to serious impact on Japanese flounder. Therefore, it is crucial to understand the effects of representative coastal economic fish under increasing water temperature. In this study, we investigated the histological and apoptosis responses, oxidative stress and transcriptomic profile in the liver of Japanese flounder exposed to gradual temperature rise (GTR) and abrupt temperature rise (ATR). The histological results showed liver cells in ATR group were the most serious in all three groups including vacuolar degeneration and inflammatory infiltration, and had more apoptosis cells than GTR group detected by TUNEL staining. These further indicated ATR stress caused more severe damage than GTR stress. Compared with control group, the biochemical analysis showed significantly changes in two kinds of heat stress, including GPT, GOT and D-Glc in serum, ATPase, Glycogen, TG, TC, ROS, SOD and CAT in liver. In addition, the RNA-Seq was used to analyze the response mechanism in Japanese flounder liver after heat stress. A total of 313 and 644 differentially expressed genes (DEGs) were identified in GTR and ATR groups, respectively. Further pathway enrichment of these DEGs revealed that heat stress affected cell cycle, protein processing and transportation, DNA replication and other biological processes. Notably, protein processing pathway in the endoplasmic reticulum (ER) was enriched significantly in KEGG and GSEA enrichment analysis, and the expression of ATF4 and JNK was significantly up-regulated in both GTR and ATR groups, while CHOP and TRAF2 were high expressed in GTR and ATR groups, respectively. In conclusion, heat stress could cause tissue damage, inflammation, oxidative stress and ER stress in the liver of Japanese flounder. The present study would provide insight into the reference for the adaptive mechanisms of economic fish in face of increasing water temperature caused by global warming.

The regulatory function of GATA3 on immune response in Japanese flounder (Paralichthys olivaceus)

GATA3 belongs to the GATA family, and it could interact with the target gene promoter. It has been reported to play a central role in regulating lymphocyte differentiation. In this study, the GATA3 cDNA sequence was identified by a homologous clone and the RACE technology from Japanese flounder (Paralichthys olivaceus). The full-length of the GATA3 cDNA sequence was 2904 bp, including 1332 bp open reading frame (ORF), 265 bp 5 '-untranslated region (5' UTR), and 1308 bp 3 '-UTR, encoding 443 amino acids. GATA3 protein sequence was conserved in vertebrates and invertebrates, including two zinc finger domains. qRT-PCR showed that the expression of GATA3 was high in the gill, kidney, and spleen. Expression of GATA3 slowly increased at the earlier stages and culminated at the late gastrula and somatic stages. Immunohistochemistry (IHC) results showed that the GATA3 protein was expressed in lymphocyte cells, undifferentiated basal and pillar cells of the gills, as well as lymphocyte cells and melanin macrophages of the kidney. The expression of GATA3 was significantly regulated in tissues and different types of lymphocytes after stimulation with Edwardsiella tarda. Dual-luciferase reporter assay indicated that the GATA3 protein could directly interact with promoters of target genes involved in the immune response. These findings suggested that GATA3 plays a major role in regulating the immune response. This study provided a theoretical basis for the immune response mechanism of teleost and a useful reference for later research on fish immunology.

Sex-Associated SNP Confirmation of Sex-Reversed Male Farmed Japanese Flounder Paralichthys olivaceus

Female Japanese flounder Paralichthys olivaceus grow more rapidly than the male. The goal of all-female commercial production requires an efficient method of genetic sex identification. We conducted genome-wide association analysis of female and male farmed Japanese flounder (n = 24 per phenotypic sex) and found all regions of chromosome 24 to be significantly associated with phenotypic sex, suggesting it as the sex chromosome. Genetic sex was identified based on single nucleotide polymorphisms (SNP) on chromosome 24 (n = 3568) using multidimensional scaling analysis, and individuals were clearly separated according to sex by the first dimension. The 61 SNPs most highly associated with sex were selected, and an amplicon-based SNP panel was developed. This was used to determine genetic sex of 39 females and 40 males. Eleven phenotypic males were assigned as female with XX genotype, suggesting sex reversal. Genetic sex was also assessed based on the indel of the amh gene promoter, which is the major candidate sex gene of Japanese flounder. We found four SNPs perfectly associated with genotypic sex in the sex-associated SNP panel, one of which was located in exon 2 of the amh gene. Along with the indel of the amh gene promoter, the sex-associated SNP panel will be of value in identifying genetic sex of farmed Japanese flounder. Molecular sexing will facilitate all-female production by breeding sex-reversed males.

Identification of the matrix metalloproteinase (MMP) gene family in Japanese flounder (Paralichthys olivaceus): Involved in immune response regulation to temperature stress and Edwardsiella tarda infection

The Matrix metalloproteinase (MMP) gene family is responsible for regulating the degradation of Extra Cellular Matrix (ECM) proteins, which are important for physiological processes such as wound healing, tissue remodeling, and stress response. Although MMPs have been studied in many species, their role in immune response in Japanese flounder (Paralichthys olivaceus) is still not fully understood. This study conducted a comprehensive analysis of MMPs in flounder, including gene structures, evolutionary relationships, conserved domains, molecular evolution, and expression patterns. Analysis revealed that MMP genes could be grouped into 17 subfamilies and were evolutionarily conserved and functionally-constrained. Meanwhile, MMP genes were found to express in different embryonic and larval stages and might play the role of sentinel in healthy tissues. Furthermore, expression profiling showed that MMPs had diverse functions in environmental stress, with 60% (9/15) and 73% (11/15) of MMPs showing differential expression patterns under temperature stress and Edwardsiella tarda (E. tarda) infection, respectively. These findings provide a useful resource for understanding the immune functions of MMP genes in Japanese flounder.

Genome-wide identification of TNF receptor-associated factors in Japanese flounder (Paralichthys olivaceus) and functional analysis of resistance to temperature and Edwardsiella tarda stress

Tumor necrosis factor receptor-associated factors (TRAFs), as the signaling mediators of the tumor necrosis factor (TNFR) superfamily, toll-like receptors (TLR) and interleukin-1 receptor (IL-1R) superfamily, can activate downstream signal transduction pathways and play an important role in the body's immune process. In this study, six TRAF genes, namely PoTRAF2a, PoTRAF2b, PoTRAF3, PoTRAF4, PoTRAF6 and PoTRAF7, were identified and annotated in Japanese flounder by using bioinformatics methods. Phylogenetic analysis confirmed that TRAF genes can be divided into seven groups. Analysis of motif composition and gene structure demonstrated that all PoTRAF members were evolutionarily conserved. The expression patterns of PoTRAF genes were then further investigated in six different developmental stages and eleven tissues of healthy fish, and it was found that there were spatial and tissue specificities among the members. To investigate the immune response of Japanese flounder to abiotic and biotic stresses, we further analyzed the expression profile of PoTRAFs after temperature stress and pathogen challenge. The result showed that PoTRAF3 and PoTRAF4 were observably differentially expressed under temperature stress, indicating that they were involved in the immune response after temperature stress. The expression of PoTRAF2a, PoTRAF2b and PoTRAF4 was significantly different after E. tarda infection, suggesting that they might have antibacterial effects. These results would help to clarify the molecular roles of PoTRAF genes in the regulation of immune and inflammatory responses in Japanese flounder.

Characterization of Japanese flounder (Paralichthys olivaceus) STAT members: An immune-related gene family involved in Edwardsiella tarda and temperature stress

The signal transducer and activator of transcription (STAT) family members are not only the transcriptional activators, but also play important roles in regulating inflammatory response. Some members have been reported to be involved in innate bacterial and antiviral immunity in aquatic organisms. However, no systematic research on STATs has been found in teleost. In this present study, we characterized six STAT genes in Japanese flounder based on bioinformatics methods, namely PoSTAT1, PoSTAT2, PoSTAT3, PoSTAT4, PoSTAT5 and PoSTAT6. The phylogenetic analysis of STATs in fish indicated that STATs were highly conserved and revealed an absence of STAT5 in a few species. Further analysis of gene structures and motifs showed STAT proteins shared a similar structure and probably had similar functionality in Japanese flounder. The expression profiles of different development stages and tissues demonstrated that PoSTATs exhibited specificity in temporality and spatiality as well as PoSTAT4 was highly expressed in gill. The transcriptome data analysis of E. tarda and temperature stress showed that PoSTAT1 and PoSTAT2 were more respective to these two kinds of stress. In addition, the results also demonstrated that these PoSTATs might regulate immune response in different ways, manifested by up-regulation in E. tarda infection and down-regulation in temperature stress. In a word, this systematic analysis of PoSTATs would provide valuable information about the phylogenetic relationship of STATs in fish species and help understand the role of STAT genes in the immune response of Japanese flounder.

Amh/Amhr2 Signaling Causes Masculinization by Inhibiting Estrogen Synthesis during Gonadal Sex Differentiation in Japanese Flounder (Paralichthys olivaceus)

The anti-Müllerian hormone (Amh) is a protein belonging to the TGF-β superfamily, the function of which has been considered important for male sex differentiation in vertebrates. The Japanese flounder (Paralichthys olivaceus) is a teleost fish that has an XX/XY sex determination system and temperature-dependent sex determination. In this species, amh expression is up-regulated in genetic males and in temperature-induced masculinization during the sex differentiation period. However, to the best of our knowledge, no reports on the Amh receptor (Amhr2) in flounder have been published, and the details of Amh signaling remain unclear. In this study, we produced amhr2-deficient mutants using the CRISPR/Cas9 system and analyzed the gonadal phenotypes and sex-related genes. The results revealed that the gonads of genetically male amhr2 mutants featured typical ovaries, and the sex differentiation-related genes showed a female expression pattern. Thus, the loss of Amhr2 function causes male-to-female sex reversal in Japanese flounder. Moreover, the treatment of genetically male amhr2 mutants with an aromatase inhibitor fadrozole, which inhibits estrogen synthesis, resulted in testicular formation. These results strongly suggest that Amh/Amhr2 signaling causes masculinization by inhibiting estrogen synthesis during gonadal sex differentiation in the flounder.

Identification and functional analysis of dual-specificity phosphatases (DUSP) genes in Japanese flounder (Paralichthys olivaceus) against temperature and Edwardsiella tarda stress

Dual-specificity Phosphatases (DUSPs) are not only the key regulators of dephosphorylating and inactivating mitogen-activated protein kinases (MAPKs), but play a crucial role in the immune response. However, the role of DUSP genes in Japanese flounder (PoDUSPs) is still unclear. In this study, 28 DUSP genes in Japanese flounder were identified and classified based on the whole genome database. Phylogenetic analysis and protein structure analysis revealed that DUSPs had highly conserved domains in teleosts. Molecular evolution analysis indicated that the PoDUSP genes were conservative during evolution and were functional-constrained. Meanwhile, PoDUSP genes were found to express in different embryonic and larval stages which might play the role of sentinel in healthy organisms. Furthermore, PoDUSP genes' expression profiles after temperature stress and Edwardsiella tarda (E. tarda) infection were determined in Japanese flounder without precedent, and the results demonstrated that Podusp1, Podusp2 and Podusp16 were more respective to temperature variation whereas Podusp1 and Podusp6 were more respective to E. tarda infection. In summary, our results provide useful resources for understanding the immune responsibilities of DUSP genes in flatfish.

Genome-wide identification and analysis of scavenger receptors and their expression profiling in response to Edwardsiella tarda infection in Japanese flounder (Paralichthys olivaceus)

The scavenger receptors (SRs) gene family, as one of pattern recognition receptors, participates in the innate immune response in diverse lineages. However, the systematic identification, characteristics and functions of SRs family are lacking in teleost. Here, we identified all 19 SRs family members in Japanese flounder (Paralichthys olivaceus) based on the genome and transcriptome data. Phylogenetic and Ka/Ks analysis demonstrated that these SRs genes were divided into five classes and all exhibited pronounced purified selection pressures. Whole genome duplication event was found in colec12, scarb2, and lamp1. Gene structure, functional domain and motif distribution analyses indicated that SRs within the different subfamilies are severely conservative. SRs genes showed diverse expression patterns in the embryogenesis and unchanged tissues. The regulations of 14 SRs genes in blood, gill and kidney after E. tarda infection suggested their roles in innate immune response. Meanwhile, ten SRs genes were differentially expressed after E. tarda stimulation in macrophages in vitro. Then we proved that PoSCARA3 could suppress the activity of NF-κB and AP-1 in HEK 293T cells by dual-luciferase assays. In summary, this study provided valuable basis for further functional characterization and immune functions of SRs genes in P. olivaceus.

Identification and functional characterization of Pomstna in Japanese flounder (Paralichthys olivaceus)

Myostatin (MSTN) as a negative regulator of muscle growth has been identified in Japanese flounder. Yet, most fish experienced the teleost specific genome duplication and possess at least two mstn genes. In current study, the second mstn gene named Pomstna is identified in Japanese flounder. Pomstna is clustered with other mstn2 of teleosts and owned highly conserved TGF-beta domain. In addition to muscle, Pomstna also highly expressed in brain and spleen. Using the primarily cultured muscle cells of Japanese flounder, we found that Pomstna could inhibit the proliferation and differentiation of muscle cells in vitro. As a ligand of TGF-beta signaling pathway, Pomstnb could regulate the expression of p21 and myod by activating the TGF-beta signaling pathway. Different from the function of Pomstnb, Pomstna could not activate the TGF-beta signaling pathway in vitro. During the differentiation of PoM cells, the expression of Pomstnb decreased significantly but the expression of Pomstna showed no change. Our study suggests that Pomstna could negatively regulate the growth and differentiation of muscle like Pomstnb yet through a different regulatory mechanism than Pomstnb. The present study suggests that muscle proliferation and differentiation were regulated by mstn not only through the TGF-beta signaling pathway but also other unknown mechanisms.

Acute damage from the degradation of Ulva prolifera on the environmental microbiota, intestinal microbiota and transcriptome of Japanese flounder Paralichthys olivaceus

Green tide outbreaks caused by overgrowth of Ulva prolifera in the Yellow Sea of China can cause serious ecological stress with concomitant economic hardships, especially to marine fisheries. In this study, short-term effects (14 days) were evaluated using fresh algae U. prolifera (FU), and a 7-day assessment of the effects of decomposing U. prolifera (DU) algal effluent was conducted to determine the effects on the environmental and intestinal microbiota, intestinal transcriptome and mortality of the commercial marine benthic fish, Japanese flounder (Paralichthys olivaceus). The results revealed that algal degradation altered the microbial community structure of fish farm water and fish intestines and increased the relative abundance of the pathogens Flavobacteriaceae in water and Vibrio in fish intestines. Fish intestinal tissue structure appeared to be damaged, as indicated in pathological sections, and transcriptome analysis showed intestinal inflammation after exposure, which may have caused an increase in fish mortality. The degradation of U. prolifera led to a bloom of potential pathogenic bacteria and the inflammation of fish intestines, which resulted in disease in the flounder population that reduced fish harvests and might pose a potential health threat.

DNA methylation in promoter region of immune related genes STAT3 and VEGFA and biochemical parameters change in muscle of Japanese flounder under acute hypoxia

Acute hypoxic stress can lead to immune response in fish, but the molecular mechanism of muscle immunity in fish under acute hypoxia are still unclear. In this study, we carried out the effect of signal transducer and activator of transcription3(STAT3) and vascular endothelial growth factor A(VEGFA) on muscle immune responses of Japanese flounder (Paralichthys olivaceus) during acute hypoxic stimulation (1.65 ± 0.28mg/L O₂; 3h, 6h, 12h, 24h) and reoxygenation (7.30 ± 0.40mg/L O₂; R12h, R24h, R48h). In situ hybridization (ISH) showed that STAT3 and VEGFA RNA were co-located in the skeletal muscle of Japanese flounder. Japanese flounder was seriously affected by hypoxia for 3h and 6h. The expression of STAT3 and VEGFA increased significantly. The methylation levels of STAT3 5'UTR region and VEGFA promoter region were significantly lower than those in normoxia group, which was negatively correlated with the expression levels of STAT3 and VEGFA. The enzyme activities (LDH, ALT, AST, ALP) changed significantly. In addition, enzyme-linked immunosorbent assay (ELISA) detected a positive correlation between serum VEGFA concentration and muscle VEGFA mRNA. The current study have shown that Japanese flounder responded to acute hypoxic stress at multiple metabolic levels by changing DNA methylation status and activating transcription factors such as HIF-1α, Nrf2 and STAT3. It is significant for the scientific development of aquaculture through analyzing the effects of hypoxia on biological immunity.

Transcriptome analysis and protein-protein interaction in resistant and susceptible families of Japanese flounder (Paralichthys olivaceus) to understand the mechanism against Edwardsiella tarda

Edwardsiella tarda is one of the most harmful bacterial pathogens for aquaculture flatfish. After artificial infection of 47 Japanese flounder (Paralichthys olivaceus) families, resistant and susceptible families were identified in this study. High-throughput sequencing was performed on the liver transcriptome of uninfected groups (PoRU and PoSU) and infected groups (PoRC and PoSC). Through assembly and annotation, a total of 3012 and 1386 differentially expressed genes (DEGs) were identified in PoRU vs. PoSU and PoRC vs. PoSC. The significant enrichment pathways between PoRU and PoSU were mainly in metabolic and biosynthesis pathways. A total of thirty dominant enrichment pathways between PoRC and PoSC mainly focused on some immune-related pathways, including the hematopoietic cell lineage, cytokine-cytokine receptor interaction, complement and coagulation cascades, antigen processing and presentation, the intestinal immune network for immunoglobulin A (IgA) production and T/B cell receptor signaling pathway. Under the protein-protein interaction (PPI) analysis, hub genes, including CD molecules, complement component factors and chemokines, were identified in the network, and 16 core genes were differentially expressed in resistant and sustainable families in quantitative polymerase chain reaction (qPCR) validation. This study represents the first transcriptome analysis based on resistant and susceptible families and provides resistant genes to understand the potential molecular mechanisms of antibacterial function in marine fish. The results obtained in this study provide crucial information on gene markers for resistant breeding of Japanese flounder.

Molecular Mechanism Based on Histopathology, Antioxidant System and Transcriptomic Profiles in Heat Stress Response in the Gills of Japanese Flounder

As an economically important flatfish in Asia, Japanese flounder is threatened by continuously rising temperatures due to global warming. To understand the molecular responses of this species to temperature stress, adult Japanese flounder individuals were treated with two kinds of heat stress—a gradual temperature rise (GTR) and an abrupt temperature rise (ATR)—in aquaria under experimental conditions. Changes in histopathology, programmed cell death levels and the oxidative stress status of gills were investigated. Histopathology showed that the damage caused by ATR stress was more serious. TUNEL signals confirmed this result, showing more programmed cell death in the ATR group. In addition, reactive oxygen species (ROS) levels and the 8-O-hDG contents of both the GTR and ATR groups increased significantly, and the total superoxide dismutase (T-SOD) activities and total antioxidant capacity (T-AOC) levels decreased in the two stressed groups, which showed damage to antioxidant systems. Meanwhile, RNA-seq was utilized to illustrate the molecular mechanisms underyling gill damage. Compared to the control group of 18 °C, 507 differentially expressed genes (DEGs) were screened in the GTR group; 341 were up-regulated and 166 were down-regulated, and pathway enrichment analysis indicated that they were involved in regulation and adaptation, including chaperone and folding catalyst pathways, the mitogen-activated protein kinase signaling (MAPK) pathway and DNA replication protein pathways. After ATR stress, 1070 DEGs were identified, 627 were up-regulated and 423 were down-regulated, and most DEGs were involved in chaperone and folding catalyst and DNA-related pathways, such as DNA replication proteins and nucleotide excision repair. The annotation of DEGs showed the great importance of heat shock proteins (HSPs) in protecting Japanese flounder from heat stress injury; 12 hsp genes were found after GTR, while 5 hsp genes were found after ATR. In summary, our study records gill dysfunction after heat stress, with different response patterns observed in the two experimental designs; chaperones were activated to defend heat stress after GTR, while replication was almost abandoned due to the severe damage consequent on ATR stress.

[Identification of released population of Japanese flounder based on stable isotopes analysis]

Japanese flounder (Paralichthys olivaceus) is an important releasing fish in the Yellow and Bohai Seas of China. The identification of wild and releasing population is the premise to evaluate the enhancement effects. In order to study the application of stable isotope in the identification of released P. olivaceus population, captured juveniles in the offshore releasing area of Qinhuangdao were distinguished into wild and released population using previous method (combination of morphology and molecular). Then, the contents of carbon and nitrogen stable isotope in muscle and otoliths (including the whole and the core region) were measured. The cultured population was set as control. The results showed that δ¹³C values (wild population: -17.19‰±0.73‰; released population: -17.10‰±0.61‰; cultured population: -20.75‰±0.07‰) and δ¹⁵N values (wild population: 11.81‰±0.38‰; released population: 11.62‰±0.48‰; cultured population: 8.64‰±0.60‰) of muscle and δ¹³C value (wild popu-lation: -4.47‰±0.35‰; released population: -4.63‰±0.29‰; cultured population: -6.59‰±0.58‰) of the whole otolith could be used to identify the cultured population, but could not be used to distinguish the wild from the released population. The δ¹³C value (wild population: -4.66‰±0.30‰; released population: -5.41‰±0.21‰; cultured population: -5.37‰±0.19‰) of the core region of otolith could be used to identify the wild popu-lation. The δ¹⁸O values of the whole and the core region of otolith from these three groups were overlapped and could not be used to distinguish different populations. Our results indicated that the δ¹³C value of the core area of otolith could be used to identify wild and released population, with application prospect in the identification of broodstocks participating in spawning migration. This study provided basic data and technical methods for evaluating early resources replenishment and the effects of Japanese flounder enhancement.

Porous microplastics enhance polychlorinated biphenyls-induced thyroid disruption in juvenile Japanese flounder (Paralichthys olivaceus)

Microplastics and polychlorinated biphenyls are ubiquitous in the marine environments. To illuminate their combined biological impacts, juvenile Japanese flounder (Paralichthys olivaceus) were exposed to 500 ng/L PCBs alone or 500 ng/L PCBs plus 2, 20, and 200 μg/L 10-μm porous MPs for 21 days. Compared to PCBs alone, co-exposure to PCBs and 20, 200 μg/L MPs reduced fish body length and body weight, and the concurrence of MPs aggravated PCBs-induced thyroid-disrupting effects, including significantly decreased L-thyroxine and L-triiodothyronine levels, more severe damage to the thyroid tissue and gill morphology, and disturbance on the expression of hypothalamus-pituitary-thyroid axis genes. The PCBs concentrations in the seawater were decreased dramatically with the increase of MPs concentrations, confirming that MPs absorbed PCBs from the seawater. Our results demonstrated that MPs enhanced the thyroid disruption of PCBs, suggesting that the risk of MPs and thyroid-disrupting chemicals on marine organisms should be paid more attention.

Molecular identification and function characterization of four alternative splice variants of trim25 in Japanese flounder (Paralichthys olivaceus)

Trim25 is a member of Tripartite Motif (TRIM) family. Previous studies report that trim25 modulates antiviral activity by activating RIG-I. In this study we explored the four alternative splicing (AS) variants X1-X4 of Japanese flounder trim25. The sequences of the AS variants were highly conserved. Expression levels of trim25 X1-X4 were increased after 12 h of poly I:C treatment in vitro. In vivo expression of X2-X4 in liver, kidney (except X2) and blood was significantly up-regulated in early stages of poly I:C treatment. Subcellular localization analysis showed that Trim25 X1-X4 were distributed in different cellular organelles. The recombinant vector pcDNA3.1-Trim25 X1-X4 were successfully overexpressed in Flounder cells and the samples were collected. Expression patterns of RIG-I pathway genes dhx58, traf6, traf2, nfkbia and il-8 were explored in vitro and in vivo after poly I:C treatment, as well as overexpressed samples. The findings of this study imply that AS variants of trim25 confer antiviral activity in Japanese flounder by modulating innate immune response.

Identification, evolution and expression analyses of mapk gene family in Japanese flounder (Paralichthys olivaceus) provide insight into its divergent functions on biotic and abiotic stresses response

Mitogen-activated protein kinase (MAPK) are a series of serine/threonine protein kinases showing evolutionary conservation, which can be activated by many stimulus signals and then transfer them from cell membrane to nucleus. MAPKs regulate a variety of biological processes, such as apoptosis, hormone signaling and immune response. In this study, 14 putative mapk genes in Japanese flounder were identified, and their basic physical and chemical properties were characterized. Phylogenetic analysis showed that mapk genes were divided into three main subfamilies, including ERK, JNK and the p38 MAPK. Selection pressure analysis revealed they were evolutionarily-constrained and undergone strong purifying selection. Gene structure and conserved protein motif comparison suggested high levels of conservation in members of mapk gene family. The expression patterns were further investigated in each embryonic and larval development stages and different tissues. In addition, RNA-seq analyses after bacteria and temperature stresses suggested mapk genes had different expression patterns. Three mapk genes showed significant differences in response to E. tarda challenge and five were induced significantly after temperature stress, indicating their potential functions. This systematic analysis provided valuable information for further understanding of the regulation mechanism of mapk gene family under different stresses in Japanese flounder.

Temporal Transcriptome Analysis Reveals Dynamic Expression Profiles of Gametes and Embryonic Development in Japanese Flounder (Paralichthys olivaceus)

The maternal-to-zygotic transition (MZT) is a crucial event in embryo development. While the features of the MZT across species are shared, the stage of this transition is different among species. We characterized MZT in a flatfish species, Japanese flounder (Paralichthys olivaceus). In this study, we analyzed the 551.57 GB transcriptome data of two types of gametes (sperms and eggs) and 10 embryo developmental stages in Japanese flounder. We identified 2512 maternal factor-related genes and found that most of those maternal factor-related genes expression decreased at the low blastula (LB) stage and remained silent in the subsequent embryonic development period. Meanwhile, we verified that the zygotic genome transcription might occur at the 128-cell stage and large-scale transcription began at the LB stage, which indicates the LB stage is the major wave zygotic genome activation (ZGA) occurs. In addition, we indicated that the Wnt signaling pathway, playing a diverse role in embryonic development, was involved in the ZGA and the axis formation. The results reported the list of the maternal genes in Japanese flounder and defined the stage of MZT, contributing to the understanding of the details of MZT during Japanese flounder embryonic development.

Enhancement of immune proteins expression in skin mucus of Japanese flounder Paralicthys olivaceus upon feeding a diet supplemented with high concentration of ascorbic acid

To search immune defense proteins in skin mucus of Japanese flounder fed with a diet containing high concentration of ascorbic acid, we carried out 2D-PAGE and compared the resolved pattern of proteins between control group that fed commercial diet and ascorbic acid supplemented group (AsA group) fed a diet supplemented with high concentration of ascorbic acid (2,000 mg/kg) for 7 days. The results revealed that there were many proteins exhibited distinct increase in AsA group. Among them, 6 regions that showed a dramatic elevation were chosen for protein identification using LC-MS/MS analysis and Mascot database search. Six proteins were identified, i.e. serotransferrin (Sero), transferrin (Trans), warm temperature acclimation-related 65 kDa protein (Wap65), complement component c3 (C3), hemoglobin beta-A chain (Hbß) and apolipoprotein A-1 (Apo). Quantitative RT-PCR analysis showed that the mRNA level of Hbß in epidermis of AsA group gave much higher increase (11.6 folds) than control group; the levels of Sero/Trans, Wap65, C3 and Apo showed no apparent difference between the two groups. The mRNA levels of wap65 and c3 in the liver and Apo in the kidney of AsA group exhibited significant increase in comparison to control group. In the case of secreted immunoglobulin M (IgM) and lysozyme (lyz), no difference of the mRNA levels of IgM in epidermis, gill, kidney, spleen and intestine, and lyz in epidermis, gill, spleen and intestine, was observed. The results of in situ hybridization confirmed the elevation of Hbß mRNA level in the epidermis tissue of AsA group. Our present study provided additional evidence showing the effectiveness of AsA in activating innate immune defense system in skin mucosal tissue of fish.

Molecular characterization and expression analysis of Japanese flounder (Paralichthys olivaceus) chemokine receptor CXCR2 in comparison with CXCR1

Chemokines are categorized into five families; one of the families is the CXC chemokines, which are critical in the pro-inflammatory process. CXC chemokines transmit signals and mediate a cell's biological activities by binding to cell surface receptors known as chemokine receptors (CXCRs). In this study, the CXCR2 from Japanese flounder (Paralichthys olivaceus) (JfCXCR2) was identified and characterized at the molecular level. The JfCXCR2 gene has a 1077 bp open reading frame that encodes a protein of 359 amino acid residues with seven transmembrane domains. Phylogenetic analysis of JfCXCR2 revealed that it belonged to the fish CXCR2 subfamily. Furthermore, JfCXCR2 was compared with the previously identified Japanese flounder CXCR1 (JfCXCR1). The expression analysis of uninfected Japanese flounder showed that JfCXCR1 and JfCXCR2 were expressed in all the tissues and organs tested but mainly in immune-related organs, including the kidney and spleen. Infection by Streptococcus iniae significantly increased the level of JfCXCR1 and JfCXCR2 mRNA in the kidney at days 1 and 3 post-infection. On the other hand, VHSV (viral hemorrhagic septicemia virus) and Edwardsiella tarda infection significantly increased JfCXCR2 mRNA levels in the kidney at days 3 and 6 post-infection, respectively. Conversely, JfCXCR1 expression was not significantly changed by either E. tarda or VHSV infection. Additionally, the peripheral blood leukocytes (PBLs) stimulated by recombinant proteins rCXCL8_L1a and rCXCL8_L1b were found to have significantly increased levels of JfCXCR1 and JfCXCR2 mRNA. Interestingly, even higher levels of JfCXCR1 and JfCXCR2 expression were observed in PBLs stimulated with rCXCL8_L1a and rCXCL8_L1b than in PBLs stimulated with either recombinant protein. These data suggest that bacterial infections may activate JfCXCR1. By contrast, JfCXCR2 may be activated by both bacterial and viral infection to mediate the immune response. These data can contribute to further understanding the functions of CXCR1 and CXCR2 in the fish immune system.

Lack of ∆5 Desaturase Activity Impairs EPA and DHA Synthesis in Fish Cells from Red Sea Bream and Japanese Flounder

Long-chain (≥ C₂₀) polyunsaturated fatty acids (LC-PUFA), such as eicosapentaenoic acid (20:5n-3, EPA) and docosahexaenoic acid (22:6n-3, DHA), are necessary for human health and are obtained from marine fish-derived oils. Marine fish are LC-PUFA-rich animals; however, many of them require LC-PUFA for growth. Therefore, it is suggested that they do not have sufficient ability to biosynthesize LC-PUFA. To evaluate in vivo LC-PUFA synthetic activity in fish cells, fish-derived cell lines from red sea bream (Pagrus major, PMS and PMF), Japanese flounder (Paralichthys olivaceus, HINAE), and zebrafish (Danio rerio, BRF41) were incubated with n-3 fatty acids labeled by radioisotopes or stable isotopes, and then, n-3 PUFA were analyzed by thin-layer chromatography or liquid chromatography-mass spectrometry. Labeled EPA and DHA were biosynthesized from labeled α-linolenic acid (18:3n-3) in BRF41, whereas they were not detected in PMS, PMF, or HINAE cells. We next cloned the fatty acid desaturase 2 (Fads2) cDNAs from PMF cells and zebrafish, expressed in budding yeasts, and then analyzed the substrate specificities of enzymes. As a result, we found that Fads2 from PMF cells was a ∆6/∆8 desaturase. Collectively, our study indicates that cell lines from red sea bream and Japanese flounder were not able to synthesize EPA or DHA by themselves, possibly due to the lack of ∆5 desaturase activity. Furthermore, this study provides a sensitive and reproducible non-radioactive method for evaluating LC-PUFA synthesis in fish cells using a stable isotope and liquid chromatography-mass spectrometry.

Molecular characterization and functional study of a tandem-repeat Galectin-9 from Japanese flounder (Paralichthys olivaceus)

Galectin-9 is a β-galactoside-binding lectin which could modulate a variety of biological functions including recognition, aggregation and clearance of pathogen. In this study, one Galectin-9 (named PoGalectin-9) was identified from Japanese flounder Paralichthys olivaceus. PoGalectin-9 belongs to the tandem-repeat type, containing one 127-amino acids CRD domain within N terminal and one 122-amino acids CRD domain within C-terminal. The open reading frame of PoGalectin-9 cDNA was 921 bp encoding 306 amino acids. Sequence similarity comparison confirmed that PoGalectin-9 shared high homology with other Galectin-9. The tissue distribution and expression profiles after bacterial infection were also investigated. PoGalectin-9 was widely distributed in all of the examined tissues of Japanese flounder but was predominantly expressed in the spleen, kidney and intestine. After Edwardsiella tarda challenge, the expression of PoGalectin-9 was up-regulated in spleen and down regulated in kidney. ELISA experiment showed that recombinant PoGalectin-9 (rPoGalectin-9) exhibit binding capacity to lipopolysaccharide (LPS) and peptidoglycan (PGN), which is significantly correlated with the concentration of rPoGalectin-9. Meanwhile, the rPoGalectin-9 protein showed strong agglutinating activities against both Gram-negative bacteria and Gram-positive bacteria. Bacterial binding experiments showed that rPoGalectin-9 could bind all examined bacteria. In conclusion, the present study indicate that PoGalectin-9 might play important roles during the immune responses of Japanese flounder against bacterial pathogens.

Growth differentiation factor 9 (gdf9) and bone morphogenetic protein 15 (bmp15) are potential intraovarian regulators of steroidogenesis in Japanese flounder (Paralichthys olivaceus)

Members of transforming growth factor-β (TGF-β) superfamily are vital regulators during the development of fish ovary. However, its intraovarian functions in teleost are still unclear. As members of the TGF-β superfamily, gdf9 and bmp15 are necessary for follicle formation and granulosa cell proliferation. Here in Japanese flounder, quantitative real-time polymerase chain reaction (qRT-PCR) and in situ hybridization (ISH) analysis showed that gdf9 and bmp15 were mainly expressed in oogonia and oocytes, whereas weakly expressed in non-ovarian tissues. Overexpression of single gdf9 and the co-overexpression with bmp15 could up-regulate the expression of most steroidogenic genes, while the overexpression of single bmp15 could down-regulate the expression of most steroidogenic genes. These findings demonstrate that single gdf9 and the combination with bmp15 may act as "activator", while single bmp15 may act as "inhibitor" in the process of steroidogenesis in flounder. This was also verified in negative feedback regulation of gdf9 and bmp15 during hormone treatment. High concentration of human chorionic gonadotropin (hCG) could down-regulate gdf9 and up-regulate bmp15, which were beneficial for the homeostasis of hCG hormone. Besides, knockdown of either gdf9 or bmp15 could significantly down-regulate most steroidogenic genes. This indicated that heterodimer of GDF9:BMP15 might be the most bioactive ligand in gonad development of flounder. Taken together, our study provided a novel recognition that gdf9 and bmp15 could regulate steroidogenesis in teleost through mechanism different from that in mammals.

Identification and Characterization of lncRNAs Related to the Muscle Growth and Development of Japanese Flounder (Paralichthys olivaceus)

Long noncoding RNAs (lncRNAs) play an important role in many life activities, but the expression pattern and function of lncRNAs in Japanese flounder skeletal muscle are unclear. In this study, 751 lncRNAs were selected from skeletal muscle in different development stages of the Japanese flounder [stage A: larval 7 days post hatching (dph); stage B: juvenile about 90 dph; stage C (female) and stage D (male): adult about 24 months] using coding potential analysis methods. In total, 232, 211, 194, 28, 29, and 14 differentially expressed lncRNAs and 9549, 8673, 9181, 1821, 1080, and 557 differentially expressed mRNAs were identified in comparisons of A versus B, A versus C, A versus D, B versus C, B versus D, and C versus D, respectively. We identified the cis- and trans-regulatory target genes of differentially expressed lncRNAs, and lncRNA-gene interaction networks were constructed using the Cytoscape program. In total, there were 200, 200, 200, 93, 47, and 11 cis-regulation relationships, and 29, 19, 24, 38, 8, and 47 trans-regulation relationships in the comparisons between A versus B, A versus C, A versus D, B versus C, B versus D, and C versus D, respectively. These results indicate that lncRNA may participate in the development of Japanese flounder skeletal muscle through cis- or trans-acting mechanisms, thus providing a scientific basis for further study of the biological function of lncRNA in Japanese flounder skeletal muscle. Based on these relationships, functional annotation of the related lncRNAs was performed by gene ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis. Differentially expressed genes associated with muscle development were enriched in multiple pairs of comparisons (e.g., differentially expressed genes LOC109640370, LOC109634180, LOC109643555, rusc1, and LOC109626999 were enriched in the actin-binding term, and differentially expressed genes LOC109640370, was, LOC109644970, LOC109643555, and itga9 were enriched in the regulation of the actin cytoskeleton pathway in the KEGG pathway analysis in the comparison of stages C and D). We predicted lncRNA-mRNA, miRNA-mRNA, and lncRNA-miRNA regulatory relationships and constructed interactive networks using Cytoscape software. Co-expression networks show that most lncRNAs interact with one or two predicted miRNAs involved in muscle growth and development. These results provide a basis for further study of the function of lncRNAs on skeletal muscle in different developmental stages of Japanese flounder.

Characterization of the Japanese flounder NLRP3 inflammasome in restricting Edwardsiella piscicida colonization in vivo

NLRP3 inflammasome is one of the most well-known inflammasomes in mammals, which plays critical roles in innate immunity. However, knowledge about this inflammasome in non-mammalian species, especially in teleost fish, remains rarely known. Herein, we established an Edwardsiella piscicida-head-kidney macrophages (HKMs) infection model in Japanese flounder, and found a robust caspase-1 activation and IL-1β maturation. To characterize the upstream receptor, we established a bioinformatic screening analysis, and found an NLRP3 homolog (JfNLRP3) from Japanese flounder, which shares an overall conservative structure architecture to human NLRP3. Moreover, the JfNLRP3 can assemble JfASC through PYD-PYD domain interaction and trigger JfCaspase-1 activation and JfIL-1β maturation. Meanwhile, the classical inflammasome activation stimulators, including nigericin, ATP or MSU, can trigger the JfCaspase-1 activation and JfIL-1β maturation in Japanese flounder HKMs. During intraperitoneal infection of E. piscicida in Japanese flounder, we found a dynamic up-regulated transcription of JfNLRP3 and JfCaspase-1 in vivo. Furthermore, knockdown of either JfNLRP3 or JfCaspase-1 reduces the serum JfIL-1β level, and promotes the bacterial colonization in systemic immune organs at 2 day-post infection, while overexpression of JfNLRP3 or JfCaspase-1 hampers the bacterial colonization in these organs of Japanese flounder. Taken together, our results identified the NLRP3 inflammasome paradigm in Japanese flounder, which not only providing new insight into the molecular mechanisms of teleost NLRP3 inflammasome and revealing its role in restricting bacterial infection in vivo, but also shedding light on the evolutionary of NLRP3 inflammasome in teleost.

Characterization and functional study of Galectin3 from Japanese flounder (Paralichthys olivaceus)

Galectins belong to the β-galactoside binding protein family and participate in both innate and acquired immunity. In this study, we described the molecular characteristics of Galectin3 gene from Japanese flounder (Paralichthys olivaceus), designed as PoGalectin3. Its open reading frame was 1128 bp, encoding a protein composed of 375 amino acids. PoGalectin3 belongs to chimeric galactose agglutinin, which contains a C-terminal carbohydrate recognition domain (CRD) (L₂₅₀-P₃₇₂), and its N-terminal is rich in proline (P) and glycine (G). Multiple sequence alignment and phylogenetic tree showed that PoGalectin3 was conservative in different aquatic animals. Tissue distribution confirmed that PoGalectin3 showed significantly highest expression in brain, moderate expression in liver, intestine and muscle. PoGalectin3 was significantly increased post infection with Edwardsiella tarda from intestine tissue of P. olivaceus. In order to investigate the binding ability of PoGalectin3 to pathogen-associated molecular patterns, the recombinant PoGalectin3 protein (rPoGalectin3) was successfully expressed and purified, and an Enzyme linked immunosorbent assay (ELISA) experiment was performed. ELISA refers to the qualitative and quantitative detection method of immune response by combining soluble antigen or antibody with solid-phase carrier. It was confirmed that rPoGalectin3 exhibited high affinity to lipopolysaccharide and peptidoglycan. The rPoGalectin3 also exhibited a concentration dependent binding capacity with Gram-positive bacteria (Bacillus pumilus, Bacillus subtilis, Bacillus cereus) and Gram-negative bacteria (Aeromonas salmonicida, E. tarda, Vibrio vulnificus). In addition, the results of microbial agglutination experiment showed that rPoGalectin3 could agglutinate Gram-positive bacteria (B. pumilus, B. subtilis) and Gram-negative bacteria (A. salmonicida, E. tarda) in the presence of Ca²⁺. In conclusion, this research laid an important foundation for the specific function analysis of PoGalectin3, which provide theoretical basis for the prevention and control of aquatic diseases.

Roles of two cyp11 genes in sex hormone biosynthesis in Japanese flounder (Paralichthys olivaceus)

The P450 side-chain cleavage enzymes P450scc (Cyp11a) and 11β-hydroxylase (Cyp11b) play important roles in sex steroid and cortisol production. Here, two duplicates of cyp11 genes were identified in Japanese flounder (Paralichthys olivaceus): Pocyp11a and Pocyp11b, respectively. Phylogenetic analysis and amino acid sequence alignment revealed that Pocyp11a and Pocyp11b shared significant identity with sequences of other teleost fish species. The quantitative real-time polymerase chain reaction (qRT-PCR) results indicated that among the studied tissues, brain tissue showed the highest expression of Pocyp11a, followed by kidney and testis tissues, whereas Pocyp11b expression was highest in the testis. The expression patterns of these two genes showed sexual dimorphism, with both genes showing higher expression in the testis than in the ovary. In-situ hybridization analysis demonstrated that Pocyp11a and Pocyp11b mRNA were both detected in oocytes, spermatocytes, and Sertoli cells, indicating that they might be involved in hormone synthesis. The expression levels of Pocyp11a and Pocyp11b were significantly downregulated by treatment with 17α-methyltestosterone (17α-MT) in the testis and ovary in both in vivo and studies. In vivo studies showed that Pocyp11a and Pocyp11b transcripts were suppressed by 17β-estradiol (E₂ ) treatment in both the testis and ovary. In addition, in vitro studies showed that the expression level of Pocyp11b was decreased by treatment with E₂ , whereas that of Pocyp11a was largely unaffected. Moreover, the expression levels of Pocyp11a and Pocyp11b in the testis cell line were significantly upregulated after NR0b1 and NR5a2 (p < .05) treatment. These results indicate that Pocyp11a and Pocyp11b might play important roles in sex hormone biosynthesis. Our research can assist future studies of the mechanisms of steroid biosynthesis and functional differences between cyp11a and cyp11b in Japanese flounder.

Molecular characterization and functional analysis of Japanese flounder (Paralichthys olivaceus) thbs2 in response to lymphocystis disease virus

In mammals, a matricellular protein, thrombospondin 2 (Thbs2) has been reported to play important roles in modulating cell-matrix interactions, vascular integrity and thrombosis formation. However, the role of gene, thbs2 has not yet been studied in teleost. In the present study, this novel fish gene from Japanese flounder was cloned and its function in resistant to lymphocystis disease virus was elucidated. The Japanese flounder thbs2 encoded a 1176-amino acid protein with 91% identity to medaka. Amino acid sequence indicated that Japanese flounder Thbs2 contained 10 typical conserved domains. The thbs2 was expressed in all stages of embryo development, and in hatched larva stage, its expression was significantly higher than that in other stages (P < 0.05). The relative expression level of thbs2 was significantly higher in the head kidney, liver, blood, gill, and heart of the lymphocystis disease virus resistant fish than in sensitive fish (P < 0.05); and in muscle, this difference was at highly significant (P < 0.01). Additionally, the distribution of Thbs2 in tissue was evaluated by immunohistochemical staining. Subcellular localization analysis showed that Thbs2 was distributed throughout the cytoplasm of the cells. Taken together, our results provide new basic data for thbs2 function, especially its role in anti-lymphocystis disease virus immune response.

New methods for purification of Paralichthys olivaceus lipovitellin and immunoassay-based detection of vitellogenin

Increasing levels of estrogenic pollution in marine environments has made the development of reliable biological detection techniques urgently needed. In this study, Japanese flounder (Paralichthys olivaceus) lipovitellin (Lv) was purified and used to establish three immunological methods for the detection of vitellogenin (Vtg), a biomarker for environmental estrogens. Firstly, five different methods were employed to purify Lv, among which water-precipitation was the fastest and easiest way to purify Lv. Japanese flounder Lv was characterized as a phospholipoglycoprotein with a molecular weight of ∼369 kDa. Using purified Lv and its specific polyclonal antibody, a sandwich enzyme-linked immunosorbent assay (ELISA) was developed. This assay had a working range from 7.8 to 250 ng/mL and a detection limit of 3.1 ng/mL. Furthermore, we developed an immunohistochemistry (IHC) and an immunofluorescence (IF) assay, both of which allowed visual detection of liver Vtg. Finally, Vtg induction in plasma and liver of juvenile Japanese flounders exposed to 17β-ethinylestradiol (EE₂) was measured using these three methods. Exposure to 10 and 50 ng/L EE₂ significantly increased plasma Vtg levels, and obvious positive fluorescence signals were observed near the liver sinusoidal vessels. These results confirmed that the methods developed effectively detected estrogenic activity of exogenous chemicals. Therefore, this study provides reliable methodologies for biomonitoring of estrogenic pollution in marine environments.

Methylation status and expression patterns of myomaker gene play important roles in postnatal development in the Japanese flounder (Paralichthys olivaceus)

Myomaker is a membrane protein that plays a crucial role in the fusion of myoblasts during muscle growth. DNA methylation, a significant factor, regulates gene expression. The aim of this study was to examine the methylation and mRNA expression patterns of the myomaker gene during 8 different postnatal developmental stages in the Japanese flounder (L: 7 days post hatch (dph); M1: 21 dph; M2: 28 dph; M3: 35 dph; J1: 90 dph; J2: 180 dph; A1: 24 months; A2: 36 months). Muscle tissue samples were taken from Japanese flounder at different postnatal development stages to measure the extent of DNA methylation and gene expression. Methylation level in the promoter and exon 1 of myomaker was measured using bisulfite sequencing, and the relative expression of myomaker during each developmental stage was measured by quantitative PCR. The relative expression levels of myomaker were up-regulated from stages L to M2, M3 to J2, and methylation of myomaker was negatively correlated with mRNA expression. Furthermore, the CpG site located at -26 bp in the promoter was the lowest methylated region in all developmental stages. These results offer a basis for understanding the mechanism by which myomaker regulates muscle formation during postnatal development.

Two cyp17 genes perform different functions in the sex hormone biosynthesis and gonadal differentiation in Japanese flounder (Paralichthys olivaceus)

P450c17, a key enzyme in the steroid generation pathway, plays an important role in the production of sex steroid and cortisol. In this study, two cyp17 gene isoforms, Pocyp17-I and Pocyp17-II were isolated from Paralichthys olivaceus gonads. Domain architecture analysis of Pocyp17-I and Pocyp17-II revealed that they had three regions important to enzymatic function. Structural analysis showed that Pocyp17-I and Pocyp17-II had 8 and 9 exons respectively, and the difference was caused by the insertion of an extra intron (intron1) in the latter. Quantitative real-time polymerase chain reaction results indicated that the expression of these two genes showed sexually dimorphism that Pocyp17-I and Pocyp17-II were highest expressed in testis and ovary, respectively. The in situ hybridization analysis of gonads indicated that Pocyp17-I and Pocyp17-II mRNA were both detected in oocytes, spermatocytes and Sertoli cells. After injection of androgen and estrogen (17α-methyltestosterone, 17β-estradiol) of different concentrations, the expression level of Pocyp17-I decreased significantly (P < 0.01), whereas estrogen had no influence on Pocyp17-II, but androgen upregulated the expression of Pocyp17-II (P < 0.05). Moreover, Pocyp17-I expression level was down-regulated significantly by NR0b1 but up-regulated by NR5a2 (P < 0.05), whereas Pocyp17-II expression level was down-regulated significantly by NR0b1 and NR5a2 (P < 0.05). All these results demonstrated that there were differences in expression patterns, feedback actions of sex hormones and transcriptional regulations between cyp17-I and cyp17-II, which revealed that cyp17-I and cyp17-II might perform different functions in sex hormones biosynthesis and gonadal differentiation in Japanese flounder.

Novel breeding approach for Japanese flounder using atmosphere and room temperature plasma mutagenesis tool

Background

Artificial induction of mutagenesis is effective for genetic resource innovation and breeding. However, the traditional mutation methods for fish breeding are not convenient or safe for daily use. Hence, development of a simple, safe and effective mutagenesis method with a high mutation rate and applicability to multiple fish species, is needed.

Results

We reported the first successful mutagenesis in a marine aquaculture fish species, Japanese flounder, Paralichthys olivaceus, using a novel atmosphere and room temperature plasma (ARTP) mutagenesis tool. ARTP treatment time was optimized for the fertilized eggs and sperm, respectively. Eggs fertilized for 60 min were treated by ARTP with a radio-frequency power input of 120 W, and the ARTP treatment time was 25 min. Under an ARTP radio-frequency power input of 200 W, the optimal treatment time for sperm diluted with Ringer’s solution by 1:40 v/v was 10 min. The ARTP-treated group presented differences in morphological traits such as body height, total length among individuals at day 90 after hatching. Whole-genome sequencing was used to reveal the mutation features of ARTP-treated individuals collected at day 120 after hatching. In total, 69.25Gb clean data were obtained from three controls and eight randomly selected ARTP-treated individuals, revealing 240,722 to 322,978 SNPs and 82,149 to 86,798 InDels located in 17,394~18,457 and 12,907~13,333 genes, respectively. The average mutation rate reached 0.064% at the genome level. Gene ontology clustering indicated that genes associated with cell components, binding function, catalytic activity, cellular process, metabolic process and biological regulation processes had higher mutation rates.

Conclusions

ARTP mutagenesis is a useful method for breeding of fish species to accelerate the selection of economically important traits that would benefit the aquaculture industry, given the variety of mutations detected.

Comparative studies on duplicated tdrd7 paralogs in teleosts: Molecular evolution caused neo-functionalization

The third-round whole genome duplication (3R-WGD) event occurred in the stem lineage of teleost during evolution, and is considered to be responsible for the biological diversification of ray-finned fishes. TUDOR domain containing protein 7 (Tdrd7), which belongs to the Tudor family proteins has been widely discussed in mammals. However, information about this gene in teleost is still lacking. In this study, two teleost tdrd7 genes (tdrd7a and tdrd7b) were identified in the transcriptome of Japanese flounder (Paralichthys olivaceus). Through genomic structure, phylogenetic, synteny analysis and online bioinformatic mining of tdrd7 duplications in other selected species, we confirmed that tdrd7a/7b were originated from the teleost-specific 3R-WGD. The tdrd7a is specific to teleost except for spotted gar. The tdrd7a showed a higher molecular evolution rate than tdrd7b with longer branch-length in the phylogenetic tree and multiple positively selected sites. Interestingly, it showed gonad specific expression pattern in adult tissues and germ cell specific distribution in embryos and gonads. Its 3'-untranslated region (3'UTR) labeled eGFP/DsRED could visualize primordial germ cells (PGCs) in zebrafish embryos. The tdrd7b did not show similar tissue and cell type specificity. These characteristic differences between the duplicated tdrd7 paralogues suggest that tdrd7a and tdrd7b have undergone neofunctionalization in Japanese flounder. Our results provide novel insight into the evolution and functional diversification of teleost tdrd7 genes deserving further investigations.

microRNA-mRNA analysis in pituitary and hypothalamus of sterile Japanese flounder

Double haploidy is an advantageous situation for genetic mapping and genome sequencing studies. In the present study, the hypothalamus and pituitary gland from sterile and fertile double-haploid (DH) Japanese flounders (aged 5 years) were used as experimental materials for studying the expression of genes in individuals with reproductive disorders, using high-throughput sequencing technology. The results revealed abnormal levels of some hormones in sterile DHs during the breeding season. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analysis showed that the significantly different microRNAs and messenger RNAs were related to metabolism, signal transduction, and melanogenesis; those related to steroid hormone synthesis and secretion related pathways were not detected. Our results suggest that the key to sterility in DHs was the arrested ovary development. However, the reason for arrested ovary development was mainly related to the lower levels of expression of genes involved in steroid biosynthesis in gonads, and was not related to the pituitary. For maintaining homeostasis, the hypothalamus and pituitary would have large differences in several processes, including signal transduction, metabolism, and immune response. The present study provides primary data for further studies on sterility in fish, and even in other animals.

High temperature increases the gsdf expression in masculinization of genetically female Japanese flounder (Paralichthys olivaceus)

In teleosts, sex is plastic and is influenced by environmental factors. Elevated temperatures have masculinizing effects on the phenotypic sex of certain sensitive species. In this study, we reared genetic XX Japanese flounder at a high temperature (27.5 ± 0.5 °C) and obtained a population of sex-reversal XX males (male ratio, 95.24%). We comparatively analyzed the dynamic characteristics of germ cells and gsdf (gonadal soma-derived factor) expression during sexual differentiation for the experimental (27.5 ± 0.5 °C) and control (18 °C ± 0.5 °C) groups. The results revealed that the germ cell proliferation inhibited and gsdf expression up-regulated in the experimental group, and the gsdf mRNA and proteins expressed in somatic cells that had direct contact with germline stem cells (with Nanos 2 protein expression) including spermatogonia and oogonia by ISH (in situ hybridization) and IHC (immunohistochemistry). In addition, we also overexpressed the gsdf in XX flounders, and the germ cell number of XX flounders bearing gsdf gene significantly decreased and sometimes disappeared completely, which was consistent with the results from high-temperature induction. Therefore, based on all the results, we speculated that the high expression of gsdf might inhibit germ cell proliferation during sex differentiation, and eventually cause sex reversal in the high-temperature induced masculinization of XX Japanese flounder.

Identification and expression analysis of Fc receptor-like proteins in Japanese flounder (Paralichthys olivaceus)

Fc receptors (FcRs) are specific to the Fc portion of immunoglobulin (Ig) molecules. Here, four Fc receptor-like proteins, JF-FcR-like protein 1-4, were identified in Japanese flounder. Their open reading frames encoded 358, 255, 519 and 441 amino acid residues, respectively. JF-FcR-like protein mRNAs were mainly detected in kidney and spleen of healthy fish. Injection of formalin-killed cells (FKCs) of Edwardsiella tarda significantly increased the spleen mRNA levels of JF-FcR-like protein 1 but not the other JF-FcR-like proteins. Injection of FKC of Streptococcus iniae did not significantly affect any of the JF-FcR-like protein mRNAs. These findings suggest that the FcR-like proteins have different involvements in pathogen responses.

The evolution and possible role of two Sox8 genes during sex differentiation in Japanese flounder (Paralichthys olivaceus)

Sox8 genes, as members of the Sox family, have been studied widely in mammals. However, regulation of sox8 genes in teleosts has rarely been studied, and functional analysis of these genes in teleosts has rarely been performed. Here, two duplicates of sox8 genes were identified in Japanese flounder, Posox8a and Posox8b. The analysis of expression showed that Posox8a and Posox8b were expressed in Sertoli cells of the testis, indicating that they play important roles in development and functional maintenance of the testis. Positive selection and phylogenetic analysis found that both Posox8a and Posox8b underwent the purification selection during evolutionary and that sox8 was most likely to be the ancestor sox8a. These results suggested that both Posox8a and Posox8b had important biological functions after generation from three rounds of whole-genome duplication in Japanese flounder. The functional differentiation of Posox8a and Posox8b was verified using cell transfection and dual-luciferase reporter assays; Posox8a overexpression-promoted 3β-hydroxysteroid dehydrogenase expression and Posox8b overexpression-promoted cytochrome P450 aromatase (cyp19a1; P450arom) expression. Finally, combined with Posox8a and Posox8b expression analysis from 30 to 100 days after hatch, we speculated that Posox8a and Posox8b might participate in the process of sex differentiation and gonadogenesis by regulating sex hormone biosynthesis in the Japanese flounder. Our study is the first to demonstrate the possible mechanism of Posox8a and Posox8b in Japanese flounder sex differentiation and gonadogenesis, laying a solid foundation for functional studies of sox8 genes in teleosts.

Genomic Selection Using BayesCπ and GBLUP for Resistance Against Edwardsiella tarda in Japanese Flounder (Paralichthys olivaceus)

The Japanese flounder is one of the most widely farmed economic flatfish species throughout eastern Asia including China, Korea, and Japan. Edwardsiella tarda is a major species of pathogenic bacteria that causes ascites disease and, consequently, a huge economy loss for Japanese flounder farming. After generation selection, traditional breeding methods can hardly improve the E. tarda resistance effectively. Genomic selection is an effective way to predict the breeding potential of parents and has rarely been used in aquatic breeding. In this study, we chose 931 individuals from 90 families, challenged by E. tarda from 2013 to 2015 as a reference population and 71 parents of these families as selection candidates. 1,934,475 markers were detected via genome sequencing and applied in this study. Two different methods, BayesCπ and GBLUP, were used for genomic prediction. In the reference population, two methods led to the same accuracy (0.946) and Pearson's correlation results between phenotype and genomic estimated breeding value (GEBV) of BayesCπ and GBLUP were 0.912 and 0.761, respectively. In selection candidates, GEBVs from two methods were highly similar (0.980). A comparison of GEBV with the survival rate of families that were structured by selection candidates showed correlations of 0.662 and 0.665, respectively. This study established a genomic selection method for the Japanese flounder and for the first time applied this to E. tarda resistance breeding.

Kunitz-type serine protease inhibitor is a novel participator in anti-bacterial and anti-inflammatory responses in Japanese flounder (Paralichthys olivaceus)

Kunitz-type serine protease inhibitor (KSPI) interacts with serine protease (SP) to regulate cascade reactions in vivo and plays essential roles in innate immunity. Theoretical considerations support various functions of kspi, but further studies are required for full characterization of these functions. In this study, a KSPI molecule was identified from Japanese flounder (Paralichthys olivaceus), and was named Pokspi. The full-length cDNA sequence of Pokspi was 2810 nt, containing an open reading frame of 1527 nt, which encoded a polypeptide of 509 amino acid residues. PoKspi protein contained five conversed domains, namely, MANEC, PKD, LDLa and two Kunitz domains. Homology analysis revealed that Pokspi shared the highest similarity (83%) with its homolog in Cynoglossus semilaevis. Phylogenetic analysis indicated that Pokspi clustered with the homologs in other fishes. The mRNA transcripts of Pokspi were detected in all tested tissues, with the highest expression level in gill, followed by kidney and intestine. Its elevated expression in response to the application of Edwardsiella tarda (in vivo) and pathogen-associated molecular pattern (in vitro) suggested the involvement of Pokspi in the essential immune defense against various pathogens. Recombinant PoKspi (rPoKspi) purified from Escherichia coli exhibited not only serine protease inhibitor activities but also a broad spectrum of anti-microbial effect in a manner that was independent of any host factors. In addition, the recombinant PoKspi protein could cause the down-regulation of pro-inflammatory factors TNF-α and IL-1β. In conclusion, Pokspi is a biologically active serine protease inhibitor endowed with anti-bacterial and anti-inflammatory property. This study provides strong evidences for understanding the innate immune defense in Japanese flounder.

The anti-androgenic effect of chronic exposure to semicarbazide on male Japanese flounder (Paralichthys olivaceus) and its potential mechanisms

Semicarbazide (SMC), a new marine pollutant, has anti-estrogenic effects on female Japanese flounder (Paralichthys olivaceus). However, whether SMC also affects the reproductive endocrine system of male marine organisms is currently unclear. In this study, Japanese flounder embryos were exposed to 1, 10, and 100 μg/L SMC for 130 days. Plasma testosterone (T) and 17β-estradiol (E₂) concentrations were significantly decreased in male flounders after SMC exposure. The expression of genes involved in T and E₂ synthesis, including steroidogenic acute regulatory protein, cytochrome P450 11A1, 17α-hydroxylase, 17β-hydroxysteroid dehydrogenase and cytochrome P450 19A, was down-regulated in the gonads, which may explain the decrease in plasma sex hormones levels. Moreover, SMC-mediated changes in the transcription of these steroidogenic genes were associated with reduced levels of follicle-stimulating hormone beta subunit (fshβ), luteinizing hormone beta subunit (lhβ), follicle-stimulating hormone receptor (fshr) and luteinizing hormone receptor (lhr) mRNA. In addition, down-regulated transcription of fshβ and lhβ in the SMC exposure groups was affected by reduced mRNA levels of seabream gonadotropin-releasing hormone (sbgnrh), g-protein-coupled receptor 54 (gpr54) in the kisspeptin/gpr54 system, as well as the gamma-aminobutyric acid (GABA) synthesis enzyme glutamic acid decarboxylase (gad). Overall, our results showed that environmentally relevant concentrations of SMC exerted anti-androgenic effects in male flounders via impacting HPG axis, kiss/gpr54 system and GABA synthesis, providing theoretical support for investigating reproductive toxicity of environmental pollutants that interfere with the neuroendocrine system.

The anti-estrogenicity of chronic exposure to semicarbazide in female Japanese flounders (Paralichthys olivaceus), and its potential mechanisms

This study investigated the anti-estrogenic effects of chronic exposure to a new marine pollutant, semicarbazide (SMC; 1, 10, and 100μg/L), in female Paralichthys olivaceus, as well as the associated mechanism. After 130days of exposure, plasma 17β-estradiol and testosterone concentrations, and hepatic estrogen receptors, vitellogenin, and choriogenin mRNA levels decreased significantly in SMC-exposed groups. Moreover, down-regulation of genes in the hypothalamic-pituitary-gonadal axis, including gonadotropin-releasing hormone, gonadotropic hormones and their receptors, the steroidogenic acute regulatory protein, 17α-hydroxylase, 17β-hydroxysteroid dehydrogenase, and cytochrome P450 19A, was observed after SMC exposure. Furthermore, the kisspeptin/g protein-coupled receptor 54 (kiss/gpr54) system and gamma-aminobutyric acid-ergic (GABAergic) system were also affected by SMC: SMC significantly down-regulated mRNA expression of kiss2, gpr54, and the GABA synthesis enzyme gad67. Our results demonstrated for the first time that environmentally relevant concentrations of SMC exerted anti-estrogenicity in female flounders, providing theoretical support for ecological risk assessments of SMC in marine environments.

The DNA methylation status of MyoD and IGF-I genes are correlated with muscle growth during different developmental stages of Japanese flounder (Paralichthys olivaceus)

Many genes related to muscle growth modulate myoblast proliferation and differentiation and promote muscle hypertrophy. MyoD is a myogenic determinant that contributes to myoblast determination, and insulin-like growth factor 1 (IGF-I) interacts with MyoD to regulate muscle hypertrophy and muscle mass. In this study, we aimed to assess DNA methylation and mRNA expression patterns of MyoD and IGF-I during different developmental stages of Japanese flounder, and to examine the relationship between MyoD and IGF-I gene. DNA and RNA were extracted from muscles, and DNA methylation of MyoD and IGF-I promoter and exons was detected by bisulfite sequencing. The relative expression of MyoD and IGF-I was measured by quantitative polymerase chain reaction. IGF-I was measured by radioimmunoassay. Interestingly, the lowest expression of MyoD and IGF-I emerged at larva stage, and the mRNA expression was negatively associated with methylation. We hypothesized that many skeletal muscle were required to complete metamorphosis; thus, the expression levels of MyoD and IGF-I genes increased from larva stage and then decreased. The relative expression levels of MyoD and IGF-I exhibited similar patterns, suggesting that MyoD and IGF-I regulated muscle growth through combined effects. Changes in the concentrations of IGF-I hormone were similar to those of IGF-I gene expression. Our results the mechanism through which MyoD and IGF-I regulate muscle development and demonstrated that MyoD interacted with IGF-I to regulate muscle growth during different developmental stages.

Roles of Two Sox9 Genes during Gonadal Development in Japanese Flounder: Sex Differentiation, Spermatogenesis and Gonadal Function Maintenance

The transcription factor sox9 has been implicated in cartilage formation and testis determination in mammals. Here, two duplicates of sox9 were found in Japanese flounder (Paralichthys olivaceus) named Posox9a and Posox9b, respectively. Phylogenetic and gene structure analyses revealed that Posox9a and Posox9b were homologous to that of teleosts and tetrapods. Quantitative real-time polymerase chain reaction (qRT-PCR) showed that both Posox9a and Posox9b expressed higher in testis than in ovary of adult tissues. The in situ hybridization (ISH) analysis of gonads showed that Posox9a and Posox9b mRNA were both detected in oocytes, Sertoli cells and spermatocytes. During sex differentiation, the expression of Posox9a exhibited obvious sexual dimorphic expression from 60 days after hatch (dah) with higher expression in male preferred individuals than female preferred individuals and increased gradually from 30 to 100 dah. A similar pattern was detected in Posox9b expression. After injection of androgen (17α-methyltestosterone) of different concentrations, the expression level of Posox9b increased significantly, whereas Posox9a did not change obviously. These results indicated that the two sox9 genes of Japanese flounder had converse functions in sex differentiation, whereas their differences in 17α-methyltestosterone administration were obvious and worthwhile for exploring evolutionary and adaptive significance. This study provided a foundation for further exploration of the roles of Posox9 genes during the sex determination and differentiation, spermatogenesis and gonadal function maintenance of Japanese flounder.

Divalent metal modulation of Japanese flounder (Paralichthys olivaceus) purinergic P2X7 receptor

Paralichthys olivaceus P2X7 receptor (poP2X7R) is a recently identified as a P2X7 purinergic receptor involved in innate immunity of the Japanese flounder Paralichthys olivaceus. Divalent metals are allosteric modulators of mammalian P2XRs, but there is no information for fish P2XRs. Here, we characterized the effects of divalent metals on poP2X7R channel activity by electrophysiology and molecular biology techniques. Copper, zinc and mercury inhibited poP2X7R‐mediated currents with different maximal inhibition potency, while cadmium had no effect on poP2X7R activity. Mercury‐induced inhibition was irreversible, but the inhibitory effects of copper and zinc were reversed after washout. Cooper and zinc also reduced poP2X7R‐mediated interleukin‐1 mRNA production. These findings suggest that divalent metals have potential effects on the Japanese flounder innate immune response through modulation of poP2X7R activity.