Dissecting the dynamic cellular transcriptional atlas of adult teleost testis development throughout the annual reproductive cycle.

Teleost testis development during the annual cycle involves dramatic changes in cellular compositions and molecular events. In this study, the testicular cells derived from adult black rockfish at distinct stages - regressed, regenerating and differentiating - were meticulously dissected via single-cell transcriptome sequencing. A continuous developmental trajectory of spermatogenic cells, from spermatogonia to spermatids, was delineated, elucidating the molecular events involved in spermatogenesis. Subsequently, the dynamic regulation of gene expression associated with spermatogonia proliferation and differentiation was observed across spermatogonia subgroups and developmental stages. A bioenergetic transition from glycolysis to mitochondrial respiration of spermatogonia during the annual developmental cycle was demonstrated, and a deeper level of heterogeneity and molecular characteristics was revealed by re-clustering analysis. Additionally, the developmental trajectory of Sertoli cells was delineated, alongside the divergence of Leydig cells and macrophages. Moreover, the interaction network between testicular micro-environment somatic cells and spermatogenic cells was established. Overall, our study provides detailed information on both germ and somatic cells within teleost testes during the annual reproductive cycle, which lays the foundation for spermatogenesis regulation and germplasm preservation of endangered species.

Deciphering dynamic interactions between spermatozoa and the ovarian microenvironment through integrated multi-omics approaches in viviparous Sebastes schlegelii.

The ovarian microenvironment plays a crucial role in ensuring the reproductive success of viviparous teleosts. However, the molecular mechanism underlying the interaction between spermatozoa and the ovarian microenvironment has remained elusive. This study aimed to contribute to a better understanding of this process in black rockfish (Sebastes schlegelii) using integrated multi-omics approaches. The results demonstrated significant upregulation of ovarian complement-related proteins and pattern recognition receptors, along with remodeling of glycans on the surface of spermatozoa at the early spermatozoa-storage stage (1 month after mating). As spermatozoa were stored over time, ovarian complement proteins were progressively repressed by tryptophan and hippurate, indicating a remarkable adaptation of spermatozoa to the ovarian microenvironment. Before fertilization, a notable upregulation of cellular junction proteins was observed. The study revealed that spermatozoa bind to ZPB2a protein through GSTM3 and that ZPB2a promotes spermatozoa survival and movement in a GSTM3-dependent manner. These findings shed light on a key mechanism that influences the dynamics of spermatozoa in the female reproductive tract, providing valuable insights into the molecular networks regulating spermatozoa adaptation and survival in species with internal fertilization.

Genome-wide identification, characterization and expression profiling of TLR family genes in Chromileptesaltivelis.

Toll-like receptors (TLRs) represent a prominent category of pattern recognition receptors that have been extensively investigated for their pivotal role in combating pathogen incursions. Despite this, there has been a notable absence of comprehensive identification and exploration of the immune response associated with the TLR family genes in C. altivelis. This study successfully identified and named fourteen genes as Catlr1-1, Catlr1-2, Catlr2-1, Catlr2-2, Catlr3, Catlr5, Catlr7, Catlr8, Catlr9, Catlr13-1, Catlr13-2, Catlr18, Catlr21, and Catlr22. A series of bioinformatic analysis were performed, encompassing analysis of protein properties, examination of gene structures, evolutionary assessments, and prediction of protein tertiary structures. The expression patterns of Catlr genes were analyzed in five immune tissues: liver, spleen, kidney, gill, and intestine, in both healthy and bacterial stimulated-fish. The results showed that different tissue and different genes showed differed expression patterns after V. harveyi infection, indicating the involvement of all Catlr members in mounting immune responses following infection in various tissues. Additionally, histological evaluations of immune tissues unveiled varying levels of damage. In conclusion, this investigation into the TLR gene family offers novel information that contribute to a more profound comprehension of the immune response mechanisms in C. altivelis.

A Well-Established Gut Microbiota Enhances the Efficiency of Nutrient Metabolism and Improves the Growth Performance of Trachinotus ovatus.

The gut microbiota has become an essential component of the host organism and plays a crucial role in the host immune system, metabolism, and physiology. Nevertheless, our comprehension of how the fish gut microbiota contributes to enhancing nutrient utilization in the diet and improving host growth performance remains unclear. In this study, we employed a comprehensive analysis of the microbiome, metabolome, and transcriptome to analyze intestines of the normal control group and the antibiotic-treated model group of T. ovatus to investigate how the gut microbiota enhances fish growth performance and uncover the underlying mechanisms. First, we found that the growth performance of the control group was significantly higher than that of the antibiotic-treated model under the same feeding conditions. Subsequent multiomics analyses showed that the gut microbiota can improve its own composition by mediating the colonization of some probiotics represented by Lactobacillus in the intestine, improving host metabolic efficiency with proteins and lipids, and also influencing the expression of genes in signaling pathways related to cell proliferation, which together contribute to the improved growth performance of T. ovatus. Our results demonstrated the important contribution of gut microbiota and its underlying molecular mechanisms on the growth performance of T. ovatus.

Genome-wide identification, characterization and expression profiling of TRAF family genes in Sebastes schlegelii.

Tumor necrosis factor receptor-associated factors (TRAFs) are signaling mediators for Toll-like receptor (TLR) and tumor necrosis factor (TNFR) superfamily that play important roles in organism immune response. However, reports on systematic identification of TRAF gene family in teleost fish and the function of TRAFs in innate immunity of black rockfish (Sebastes schlegelii) are lacked. In our study, eight TRAF genes were identified and characterized, namely, SsTRAF2a, SsTRAF2a-like, SsTRAF2b, SsTRAF3, SsTRAF4, SsTRAF5, SsTRAF6 and SsTRAF7 in S. schegelii. Furthermore, we analyzed their sequences, conserved domains, gene structures, motif compositions, phylogeny, tissue expression patterns in healthy and Vibro. anguillarum challenged individuals. All the SsTRAFs contained typical conserved domain, including C-terminal MATH domain and N-terminal RING finger domain. Analyses of gene structures and motifs showed the distribution of exon-intron and conserved motifs in S. schegelii and serval other teleost fish. We also analyzed the expression file of SsTRAFs in five immune-relate organs, liver, spleen, kidney, gill and intestine in healthy and bacterial challenged fish. The results indicated that all SsTRAF member were widely involved in immune response after pathogenic bacteria infection. In summary, the analyses of TRAFs in S. schegelii will be helpful to better understand the diverse roles of TRAF genes in the innate immune response to bacterial challenge.

Intensive masculinization caused by chronic heat stress in juvenile Cynoglossus semilaevis: Growth performance, gonadal histology and gene responses.

The sea temperature has been observed to chronically increase during the past decades, leaving unpredictable influences to the marine biological resources. Thus, it is of vital significance to study the biological responses of ocean inhabited organisms with the artificially stimulated heat stress environment. Cynoglossus semilaevis provides us with an ideal model to study the influence of chronic heat stress on the sexual differentiation in marine teleosts for its genetic sex determination (GSD) + environmental effected (EE) sex determination system. In this study, the comparative experiment was conducted employing heated seawater (HT group) and ambient seawater (CT group) to cultivate juvenile C. semilaevis respectively. Significant differences were exhibited in growth performance and a delayed germ cell development effect was found in pseudomales formed under chronic heat stress. Using transcriptome analysis, the transcription profile of 55 days post fertilization (dpf) and 100 dpf juveniles' gonads were studied. A total of 47 libraries were constructed with an average mapping rate of 94.63% after assembling. GO and KEGG enrichment were proceeded using DEGs screened out between (1) pseudomale gonads at 55 dpf and 100 dpf in HT and CT group (2) pseudomale and female gonads at 55 dpf and 100 dpf in HT and CT group. Terms and pathways involved in steroid stimulation, reproduction ability, germ cell proliferation et al. were shed light on. The expression pattern of 29 DEGs including amh, hsp90b1, pgr et al. were also provided to supplement the results of functional enrichment. Weighted gene co-expression networks analysis (WGCNA) was constructed and hspb8-like, histone H2A.V were exhibited to play vital roles in the heat-induced masculinization. Our findings facilitate the understanding for transcriptional variations in intensive masculinization cause by chronic heat stress of C. semilaevis and provide referable study of the influences on the teleosts in elevated sea temperature.

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.

Transcriptome and Network Analyses Reveal the Gene Set Involved in PST Accumulation and Responses to Toxic Alexandrium minutum Exposure in the Gills of Chlamys farreri.

Bivalve molluscs are filter-feeding organisms that can accumulate paralytic shellfish toxins (PST) through ingesting toxic marine dinoflagellates. While the effects of PST accumulation upon the physiology of bivalves have been documented, the underlying molecular mechanism remains poorly understood. In this study, transcriptomic analysis was performed in the gills of Zhikong scallop (Chlamys farreri) after 1, 3, 5, 10, and 15 day(s) exposure of PST-producing dinoflagellate Alexandrium minutum. Higher numbers of differentially expressed genes (DEGs) were detected at day 1 (1538) and day 15 (989) than that at day 3 (77), day 5 (82), and day 10 (80) after exposure, and most of the DEGs were only regulated at day 1 or day 15, highlighting different response mechanisms of scallop to PST-producing dinoflagellate at different stages of exposure. Functional enrichment results suggested that PST exposure induced the alterations of nervous system development processes and the activation of xenobiotic metabolism and substance transport processes at the acute and chronic stages of exposure, respectively, while the immune functions were inhibited by PST and might ultimately cause the activation of apoptosis. Furthermore, a weighted gene co-expression network was constructed, and ten responsive modules for toxic algae exposure were identified, among which the yellow module was found to be significantly correlated with PST content. Most of the hub genes in the yellow module were annotated as solute carriers (SLCs) with eight being OCTN1s, implying their dominant roles in regulating PST accumulation in scallop gills. Overall, our results reveal the gene set responding to and involved in PST accumulation in scallop gills, which will deepen our understanding of the molecular mechanism of bivalve resistance to PST.

Pax3 and Pax7 Exhibit Distinct and Overlapping Functions in Marking Muscle Satellite Cells and Muscle Repair in a Marine Teleost, Sebastes schlegelii.

Pax3 and Pax7 are members of the Pax gene family which are essential for embryo and organ development. Both genes have been proved to be markers of muscle satellite cells and play key roles in the process of muscle growth and repair. Here, we identified two Pax3 genes (SsPax3a and SsPax3b) and two Pax7 genes (SsPax7a and SsPax7b) in a marine teleost, black rockfish (Sebastes schlegelii). Our results showed SsPax3 and SsPax7 marked distinct populations of muscle satellite cells, which originated from the multi-cell stage and somite stage, respectively. In addition, we constructed a muscle injury model to explore the function of these four genes during muscle repair. Hematoxylin-eosin (H-E) of injured muscle sections showed new-formed myofibers occurred at 16 days post-injury (dpi). ISH (in situ hybridization) analysis demonstrated that the expression level of SsPax3a and two SsPax7 genes increased gradually during 0-16 dpi and peaked at 16 dpi. Interestingly, SsPax3b showed no significant differences during the injury repair process, indicating that the satellite cells labeled by SsPax3b were not involved in muscle repair. These results imply that the muscle stem cell populations in teleosts are more complicated than in mammals. This lays the foundation for future studies on the molecular mechanism of indeterminant growth and muscle repair of large fish species.

Genome-wide identification of nonvisual opsin family reveals amplification of RPE-retinal G protein receptor gene (RGR) and offers novel insights into functions of RGR(s) in Paralichthys olivaceus (Paralichthyidae, Teleostei).

Opsins play important roles in the image-forming visual pathways and numerous biological systems such as the biological clock and circadian rhythm. However, the nonvisual opsins involved in nonimage forming process are not clear to date. The aim of this study was to characterize nonvisual opsins in Paralichthys olivaceus. A total of 24 nonvisual opsin genes were identified. Expressions of these genes in eye, brain, heart, testis, and fin were investigated by quantitative real-time polymerase chain reaction (qRT-PCR). Testis contained a surprisingly large number of nonvisual opsins including Opn4m2, Tmt2a, Tmt3b, Opn3, RRH, Opn7a, and Opn7b. Syntenic and phylogenetic analyses confirmed that the RGRa and RGRb originated from the teleost-specific genome duplication (TSGD). qRT-PCR results demonstrated high RGRa and RGRb expression in the eye, while the expression levels in the brain, heart, testis, and fin were relatively weak. In situ hybridization results presented here revealed the presence of both RGRa and RGRb mRNA-positive signals in the ganglion cell layer but absence in the intracellular compartment of retinal pigment epithelium (RPE) and Müller glial cells. Therefore, we hypothesized that RGRa and RGRb had undergone subfunctionalization in P. olivaceus after TSGD. In conclusion, this study provides novel insights into the evolutionary fates of the RGR genes, still, further studies need to be done to explore the mechanism about the lack of RGR genes' expression in RPE.

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.

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.

Molecular characterization, expression and functional analysis of cystatin C in Japanese flounder (Paralichthys olivaceus).

Cystatins are natural tight-binding reversible inhibitors of cysteine proteases found in a wide arrange of organisms. Studies have shown that cystatins play important roles under both physiological and pathological conditions in mammals. However, much less is known about fish cystatins. In this study, we described the identification and analysis of the gene encoding cystatin C in Japanese flounder (Paralichthys olivaceus). This gene had a high homology with the sequence of cystatin C in many fish species and had a signal peptide and three conserved functional sites. The results of qRT-PCR showed that the gene was highly expressed in the liver. Lipopolysaccharide, peptidoglycan and polyinosinic-polycytidylic acid all increased its expression after stimulation. Functional analysis showed that the recombinant P. olivaceus cystatin C purified from Escherichia coli had cysteine protease inhibitory activity and could inhibit bacterial growth by binding to bacteria. Meanwhile, rPocystatin C could up-regulate the expression of cytokines tumor necrosis factor α and interleukin 10. These results indicated that cystatin C of P. olivaceus might be considered to have the similar immunomodulatory function to mammalian cystatin.

Edwardsiella tarda-induced miR-7a functions as a suppressor in PI3K/AKT/GSK3ß signaling pathway by targeting insulin receptor substrate-2 (IRS2a and IRS2b) in Paralichthys olivaceus.

To study the effect of Edwardsiella tarda infection on miRNAs expression profile in Japanese flounder, fish were injected intraperitoneally with E. tarda. The miRNAs involved in regulating immune responses were analyzed by high-throughput sequencing. A total of 164 mature miRNAs were identified, of which 17 miRNAs were differentially expressed (DE miRNAs) after E. tarda infection, indicating that they were immune-related miRNAs. To further examine the relationship between the miRNAs and their predicted target mRNAs, a total of 22 predicted target mRNAs, mainly related to endocytic signaling pathway, NF-κB signaling pathway, and p53 signaling pathway, were detected with miRNA mimics in HEK-293T cells by dual-luciferase reporter experiments. Finally, we confirmed that insulin receptor substrate-2 (IRS2a and IRS2b) were regulated by miR-7a. And the target sites of the 3' untranslated region (UTR) of IRS2a and IRS2b were verified by dual-luciferase reporter experiments. Furthermore, we found that the E. tarda and LPS significantly increased host miR-7a expression. In vivo and in vitro studies revealed that IRS2-mediated PI3K/AKT/GSK3ß signaling pathway was suppressed. Taken together, these results implied that miR-7a might be a key regulator of PI3K/AKT/GSK3ß signaling pathway via suppressing the IRS2a and IRS2b genes.

A functional polymorphism in the promoter of RhoB is associated with susceptibility to Vibrio anguillarum in turbot (Scophthalmus maximus).

As an isoform of Rho family GTPases, RhoB plays a pivotal role in cytoskeletal organization, cell proliferation, apoptosis and immune response. However, the regulatory mechanisms of RhoB expression in aquatic animals are still unknown. In the present study, we first construct Vibrio anguillarum infection model in S. maximus, including susceptible and resistant individuals. Then the temporal expression of RhoB was detected after V. anguillarum challenge using qRT-PCR and found that RhoB transcripts were significantly induced in the liver, gill and blood despite of differential expression levels and responsive time points. In addition, the mRNA levels of RhoB in resistant individuals were significantly higher than in susceptible ones. The length of 2083 bp sequences of RhoB promoter was cloned and characterized. Moreover, DNA methylation of the RhoB promoter was measured by bisulfite sequencing (BSP) and hypo-methylated was detected in the CpG islands. Three SNPs (-1590, -1575 and -1449) and two haplotypes in the promoter region of RhoB were identified to be associated with V. anguillarum resistance in turbot by association analysis in group 17-R and 17-S. Deletion analysis indicated that these SNPs could negatively mediate the activity of RhoB promoter. Site-directed mutagenesis and qRT-PCR of individuals with different genotypes demonstrated that -1575 T/A polymorphism affected promoter activity. Further study showed that this mutation altered the binding site of the transcription factor CREB. Co-transfection of SmCREB and RhoB promoter was performed in HEK293T cells which confirmed the -1575 allelic differences on transcriptional activity, with the susceptibility allele showing reduced activity. Taken together, our findings implicate that losing of binding of CREB to SmRhoB promoter due to -1575T/A polymorphisms enhances SmRhoB expression in resistant turbot, which provide insights into the effect of SmRhoB expression in response to V. anguillarum infection.

Spotted knifejaw (Oplegnathus punctatus) MyD88: Intracellular localization, signal transduction function and immune responses to bacterial infection.

Myeloid differentiation factor 88 (MyD88) links members of the toll-like receptor (TLR) and interleukin-1 receptor (IL-1R) superfamily to the downstream activation of NF-κB as a "bridge" molecular in response to exogenous pathogen, but the function in spotted knifejaw (Oplegnathus. punctatus), a commercial fish in China, is still unknown. We present a functional analysis of spotted knifejaw MyD88 (OppMyD88) with a typical death domain (DD) at the N-terminus and a conservative Toll/IL-1R (TIR) domain at the C-terminus and suggest that MyD88 is important for the activation of TLR-mediated NF-κB with the synergy between domains. Subcellular localization showed that OppMyD88 was distributed in the cytoplasm in a condensed form. Tissues expression profiling analysis showed that OppMyD88 ubiquitously expressed in all tested tissues with the highest expression in the liver, as determined by real-time PCR. The expression of OppMyD88 significantly upregulated in the liver, spleen, kidney and gills within 120 h post Vibrio anguillarum infection. Moreover, we further confirmed that over-expressed OppMyD88 could also induce apoptosis. These results indicate that OppMyD88 might possess important roles in defense against microbial infection and other biological processes in spotted knifejaw similar to those in mammals, which will deepen our understandings in innate immunity of spotted knifejaw.

Whole-genome resequencing from bulked-segregant analysis reveals gene set based association analyses for the Vibrio anguillarum resistance of turbot (Scophthalmus maximus).

Many achievements have been made to develop quantitative trait loci (QTLs) and gene-associated single nucleotide polymorphisms (SNPs) to facilitate practical marker-assisted selection (MAS) in aquatic animals. However, the systematic studies of SNPs associated with extreme threshold traits were poor in populations lacking of parental genomic information. Coupling next generation sequencing with bulked segregant analysis (BSA) should allow identification of numerous associated SNPs with extreme phenotypes. In the present study, using combination of SNP frequency difference and Euclidean distance, we conducted linkage analysis of SNPs located in genes involved in immune responses, and identified markers associated with Vibrio anguillarum resistance in turbot (Scophthalmus maximus). A total of 221 SNPs was found as candidate SNPs between resistant and susceptible individuals. Among these SNPs, 35 loci located in immune related genes were genotyped in verification population and 7 of them showed significant association with V. anguillarum resistance in both alleles and genotypes (P < 0.05). Among these 7 genes, PIK3CA-like, CYLD, VCAM1, RhoB and RhoGEF are involved in PI3K/Akt/mTOR pathway and NF-κB pathway, which influence the efficiency of bacteria entering the host and inflammation. SNP-SNP interaction analysis was performed by generalized multifactor dimensionality reduction (GMDR). The combination of SNP loci in RhoB, PIK3CA-like and ADCY3 showed a significant effect on V. anguillarum resistance with the verification rate in the sequencing population up to 70.8%. Taken all, our findings demonstrated the feasibility of BSA-seq approach in identifying genes responsible for the extreme phenotypes and will aid in performing MAS in turbot.

Functional differentiation of three phosphatidylinositol 3-kinase catalytic subunit alpha (PIK3CA) in response to Vibrio anguillarum infection in turbot (Scophthalmus maximus).

PIK3CA has been extensively investigated from its molecular mechanism perspective and association with its mutations in different types of cancers. However, little has been reported regarding the pathological significance of PIK3CA expression in teleost. Here, in our present study, three PIK3CA genes termed SmPIK3CAa, SmPIK3CAb and SmPIK3CA-like were firstly identified in the genome of turbot S. maximus. Although these three genes located in different chromosomes, all of them share the same five domains. Phylogenetic and synteny analysis indicated that SmPIK3CAa, SmPIK3CAb and SmPIK3CA-like were three paralogs that may originate from duplication of the same ancestral PIK3CA gene. Subcellular localization analysis confirmed the cytoplasm distribution of these three paralogs. All three SmPIK3CA were ubiquitously expressed in examined tissues in turbot, with the higher expression levels in immune-related tissues such as blood, spleen, kidney, gills and intestines. Upon Vibrio anguillarum challenge, SmPIK3CAa and SmPIK3CA-like transcripts were significantly induced in spleen, intestine and blood despite of differential expression levels and responsive time points. Additionally, individuals in resistant group showed significantly higher expression level of both two genes than in the susceptible group. Moreover, four SNPs (102, 2530, 3027 and 3060) and one haplotype (Hap2) located in exon region of SmPIK3CA-like were identified and confirmed to be associated with V. anguillarum resistance in turbot by association analysis in different populations. Taken together, these results suggested that functional differentiation occurred in three SmPIK3CA paralogs with Vibrio anguillarum resistance and SmPIK3CAa and SmPIK3CA-like probable play potential roles in innate immune response to pathogenic invasions in turbot.

A novel C-type lectin from spotted knifejaw, Oplegnathus punctatus possesses antibacterial and anti-inflammatory activity.

C-type lectin is a type of carbohydrate-binding protein and plays significant roles in innate immune response against pathogen infection. To date, thousands of C-type lectin had been identified in teleost. In the present study, we isolated a novel isoform of C-type lectin (OppCTL) from spotted knifejaw (Oplegnathus punctatus). The OppCTL encoded a typical Ca2+-dependent carbohydrate-binding protein, and was mainly expressed in liver in a tissue specific fashion. The expression of OppCTL was significantly up-regulated following Vibrio anguillarum infection in vivo, suggesting involvement in immune response. Hemagglutination analysis showed that the recombinant OppCTL (rOppCTL) could agglutinate erythrocyte from Mus musculus, Oplegnathus punctatus, Sebastes schlegelii and Paralichthys olivaceus. The rOppCTL could bind and agglutinate all tested bacteria. The rOppCTL possessed capacities of calcium-dependent agglutination to all tested bacteria. Sugar binding assay revealed that rOppCTL could also bind to the glycoconjugates of the bacterial surface, including lipopolysaccharide and peptidoglycan. Interestingly, Dual-luciferase analysis revealed that OppCTL could inhibit the activity of NF-κB in HEK-293T cells after OppCTL overexpression. Taken together, these results indicate that OppCTL has immune activity capable of defending invading pathogens and possesses potential immunoregulatory activity, enriching our understanding of the function of C-type lectin.

Molecular characterization, expression and immune functions of cystatin B in Japanese flounder (Paralichthys olivaceus).

Cystatin B is an intracellular inhibitor that regulates the activities of cysteine proteases. In this study, cystatin B in Japanese flounder (Paralichthys olivaceus) was characterized and its immune function was analyzed. This gene had a high similarity with the sequence of cystatin B in other fish species, and the derived peptide shared typical features of cystatin proteins including the QXVXG motif. The results of quantitative real-time PCR showed that cystatin B mRNA was constitutively expressed in all examined tissues, with the highest level in gill. The stimulations of lipopolysaccharide, peptidoglycan and polyinosinic-polycytidylic acid effectively increased the expression level of cystatin B mRNA. Functional analysis implied that the recombinant P. olivaceus cystatin B purified from Escherichia coli had cysteine protease inhibitory activity and could inhibit bacterial growth by binding to bacteria. Furthermore, we found that P. olivaceus cystatin B had no effects on the expression of inflammatory factors cytokines tumor necrosis factor α, interleukin 10, interleukin 1ß and interferon γ. These results indicate that cystatin B of P. olivaceus is potentially involved in immune responses against invading microbial pathogens, and provide a better understanding of the immune mechanisms of cystatins in teleosts.