Transcriptome analysis revealed changes of multiple genes involved in immunity in Cynoglossus semilaevis during Vibrio anguillarum infection

Identification and functional characterization of annexin A2 in half-smooth tongue sole (Cynoglossus semilaevis)

Annexin A2 (AnxA2), belonging to the annexin family, plays a crucial role in immune responses. In this study, the cDNA of the AnxA2 gene was identified in half-smooth tongue sole, Cynoglossus semilaevis. The transcript of AnxA2 gene in C. semilaevis (CsAnxA2) showed broad tissue distribution, with the highest expression level observed in the gut. CsAnxA2 expression was significantly up-regulated in the intestine, spleen, and kidney tissues following exposure to Shewanella algae. Immunohistochemical staining revealed that CsAnxA2 was predominantly expressed in epithelial cells and significantly elevated after S. algae challenge. Subcellular localization showed that CsAnxA2 was primarily localized in the cytoplasmic compartment. Moreover, proinflammatory cytokines (IL-6, IL-8 and IL-1β) exhibited significant upregulation after CsAnxA2 was overexpressed in vivo. One hundred and fifty-eight CsAnxA2-interacting proteins were captured in the intestinal tissue, showing the top two normalized abundance observed for actin beta (ACTB) and protein S100-A10 (p11). Fifty-four high abundance CsAnxA2-interacting proteins (HIPs) were primary enriched in ten pathways, with the top three significantly enriched pathways being Salmonella infection, glycolysis/gluconeogenesis, and peroxisome proliferator-activated receptor (PPAR) signaling pathway. These results provide valuable information for further investigation into the functional mechanism of AnxA2 in C. semilaevis.

Differentially expressed proteins and microbial communities of the skin regulate disease resistance to Chinese tongue sole (Cynoglossus semilaevis)

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

Mechanistic study of C5F10O-induced lung toxicity in rats: An eco-friendly insulating gas alternative to SF6

The global warming and other environmental problems caused by SF6 emissions can be reduced due to the widespread use of eco-friendly insulating gas, perfluoropentanone (C5F10O). However, there is an exposure risk to populations in areas near C5F10O equipment, so it is important to clarify its biosafety and pathogenesis before large-scale application. In this paper, histopathology, transcriptomics, 4D-DIA proteomics, and LC-MS metabolomics of rats exposed to 2000 ppm and 6000 ppm C5F10O are analyzed to reveal the mechanisms of toxicity and health risks. Histopathological shows that inflammatory cell infiltration, epithelial cell hyperplasia, and alveolar atrophy accompanied by alveolar wall thickening are present in both low-dose and high-dose groups. Analysis of transcriptomic and 4D-DIA proteomic show that Cell cycle and DNA replication can be activated by both 2000 ppm and 6000 ppm C5F10O to induce cell proliferation. In addition, it also leads to the activation of pathways such as Antigen processing and presentation, Cell adhesion molecules and Complement and coagulation cascades, T cell receptor signal path, Th1 and T cell receptor signal path, Th1 and Th2 cell differentiation, complement and coagulation cascades. Finally, LC-MS metabolomics analysis confirms that the metabolic pathways associated with glycerophospholipids, arachidonic acid, and linoleic acid are disrupted and become more severe with increasing doses. The mechanism of lung toxicity caused by C5F10O is systematically expounded based on the multi-omics analysis and provided biosafety references for further promotion and application of C5F10O.

Comparative transcriptome analysis reveals potential regulatory mechanisms of genes and immune pathways following Vibrio harveyi infection in red drum (Sciaenops ocellatus)

Red drum (Sciaenops ocellatus), as an important economical marine fish, has been affected by various bacterial diseases in recent years. Vibrio harveyi cause fatal vibriosis in S. ocellatus, leading to massive mortality and causing significant setbacks in aquaculture. However, the regulatory mechanisms of S. ocellatus response to V. harveyi infection are poorly understood. In this regard, we performed transcriptomic analysis with head kidney tissues of S. ocellatus after V. harveyi infection from 12 h to 48 h to reveal genes, gene expression profiles, and pathways involved in immune and inflammation responses. Specifically, a total of 9,599, 5,728, and 7144 differentially expressed genes (DEGs) were identified after V. harveyi infection at 12 h, 24 h, and 48 h, respectively, and 1,848 shared DEGs have been identified from the above three comparison groups. Subsequent pathway analysis revealed that the shared DEGs following V. harveyi were involved in complement and coagulation cascades (C1R, C1QC, C3, C4, C5, C7, C8A, C8B, C8G, C9, CFB, CFH, and CFI), MAPK signaling pathway, chemokine signaling pathway (CCL19, CXCL8, CXCL12, CXCL14, CCR4, CCR7, and CXCR2), PPAR signaling pathway (PPAR-α, PPAR-γ and PPAR-β), and TNF signaling pathway. Finally, the expression patterns of DEGs in head kidney tissues and S. ocellatus macrophages were validated by qRT-PCR, suggesting the reliability of RNA sequencing for gene expression analysis. This dynamic transcriptome analyses provided insights into gene expression regulation and immune related pathways involved in S. ocellatus after V. harveyi infection, and provides useful information for further study on the immune defense mechanisms in S. ocellatus as well as other teleost species.

Regulatory mechanism of miR-722 on C5aR1 and its functions against bacterial inflammation in half-smooth tongue sole (Cynoglossus semilaevis)

MicroRNAs (miRNAs) are small non-coding RNAs involved in various biological processes, including immunity. Previously, we investigated the miRNAs of half-smooth tongue sole (Cynoglossus semilaevis) and found that miR-722 (designated Cse-miR-722) was significantly differentially expressed after infection with Vibrio anguillarum, reflecting its importance in immune response. Our preliminary bioinformatic analysis suggested that Cse-miR-722 could target C5aR1 (designated CsC5aR1), which was known to play crucial roles in complement activation and inflammatory response, as a receptor of C5a. However, the underlying mechanisms of their interactions and specific functions in inflammatory and immune response are still enigmas. In this study, we successfully cloned the precursor sequence of Cse-miR-722 (94 bp) and the full length of CsC5aR1 (1541 bp, protein molecular weight 39 kDa). The target gene of Cse-miR-722 was verified as CsC5aR1 by a dual luciferase reporter assay, and Cse-miR-722 was confirmed to regulate CsC5aR1 at the protein level using quantitative real-time polymerase chain reaction (qRT-PCR) and Western blotting. The expression of CsC5aR1 and Cse-miR-722 in liver cells and four immune tissues of half-smooth tongue sole changed significantly after LPS stimulation and infection with V. anguillarum. To explore the functional role of Cse-miR-722 in half-smooth tongue sole, we performed both in vitro and in vivo experiments. Cse-miR-722 was observed to affect phagocytosis and respiratory burst activity of macrophages by regulating CsC5aR1 in half-smooth tongue sole. Furthermore, we found that Cse-miR-722 regulated the expression of CsC5aR1, CsC5a, and the inflammatory factors CsIL1-β, CsIL6, CsIL8, and CsTNF-α both in vitro and in vivo. In addition, Cse-miR-722 reduced mortality and pathological damage. This study clarified the regulatory mechanism of Cse-miR-722 on CsC5aR1 and provided insight into the regulatory roles of Cse-miR-722 in immune responses, laying a theoretical foundation for the feasibility of using miR-722 to prevent and control bacterial diseases in teleost.

Inflammatory reaction and immune response of half-smooth tongue sole (Cynoglossus semilaevis) after infection with Vibrio anguillarum

Frequently occurred bacterial diseases have seriously affected the aquaculture industry of half-smooth tongue sole (Cynoglossus semilaevis). Notably, vibriosis, with Vibrio anguillarum as one of the causative pathogens, is the most severe bacterial disease with severe inflammatory response of the host, leading to high mortality rates. In the present study, we explored the relationship between bacterial concentrations and host mortality, inflammatory reaction, and immune response in half-smooth tongue sole after infection with V. anguillarum at different concentrations (Treatment 1, 6.4 × 105 CFU/mL; Treatment 2, 6.4 × 106 CFU/mL). The mortality of Treatment 2 (77.5%) was significantly higher than that of Treatment 1 (10%), corresponding with bacterial concentrations. Although the number of deaths varies, intensive deaths were observed within 24 h post infection (hpi) in both bacterial concentration groups. Histopathological analyses revealed that fish tissues were most severely damaged at 24 or 48 hpi, and Treatment 2 was more severe than Treatment 1. A qRT-PCR-based detection method with virulence factor gene empA was established to quantify the bacterial loads in various tissues, and the bacterial loads were the highest at 24 hpi in Treatment 2, and at 48 hpi in Treatment 1. Additionally, the expression levels of complement genes (C5a, C3, C5, and C6), inflammatory factors (IL-1β, TNF-α, and IL-10), and other immune-related genes (jak2, NF-κB1, stat3, and tlr3) were increased in various tissues after infection in both treatment groups, with most genes being most expressed at 24 or 48 hpi, and expression levels of inflammatory factors in Treatment 2 were higher than those in Treatment 1. Moreover, the expression of C5a was positively correlated with that of proinflammatory cytokines in both bacterial concentration groups. According to the results of this study, 24-48 hpi was a key node for early vibriosis detection and intervention. Compared with the low mortality of Treatment 1, the mass death of fish in Treatment 2 was suggested to be caused by uncontrolled excessive inflammatory reaction induced by the overactivation of complement system, especially C5a. We believe these results could provide theoretical basis for prevention, evaluation, and treatment of vibrio disease in tongue sole aquaculture, and lay a solid foundation for future functional analyses.

Genome-wide characterization of integrin (ITG) gene family and their expression profiling in half-smooth tongue sole (Cynoglossus semilaevis) upon Vibrio anguillarum infection

Integrins (ITGs) are transmembrane heterodimer receptors with ITGα subunit and ITGβ subunit, participating in various physiological processes, including immunity. At present, systematic research on ITGs in teleost is scarce, especially in half-smooth tongue sole (Cynoglossus semilaevis). In this study, a set of 28 ITG genes in half-smooth tongue sole have been identified and characterized. The phylogenetic analysis showed that ITGα and ITGβ subunits were respectively classified into five and two clusters, consistent with previous studies. The selection pressure analysis indicated that most of ITG genes were under purifying selection, except for ITGα11b and ITGαL with positive selection. The expression profiles of eight selected ITG genes, including ITGα1, ITGα5, ITGα8, ITGα11, ITGβ1, ITGβ2, ITGβ3, and ITGβ8, were analyzed in healthy tissues and after infection with Vibrio anguillarum, revealed their implications in immune response. The study provided a comprehensive characterization and expression analysis of ITG genes in half-smooth tongue sole, setting a solid foundation for further functional studies and promising potential in disease control.

Transcriptomic Analysis of Liver Tissue of Black Sea Bass (Centropristis striata) Exposed to High Nitrogen Environment

The black sea bass, Centropristis striata, is a potential candidate for commercial aquaculture. Due to inadequate removal of nitrogen in its breeding environment, C. striata exhibits increased nitrate concentration, which can cause acute toxicity, including energy metabolism damage and tissue damage. Therefore, RNA-seq technology was applied to characterize genes associated with toxicity tolerance under nitrate stress. The nitrate treatment caused significant changes in a total of 8920 genes, of which 2949 genes were up-regulated and 5971 genes were down-regulated. It was found that significantly enriched GO terms and KEGG were associated with blood microparticles, inhibitors of enzyme activity, and complement and coagulation cascade pathways. Furthermore, through bioinformatics analysis, it was found that these different pathways obtained in GO and KEGG enrichment analysis were mostly related to the immune and inflammatory response of fish. This study expands our understanding of the mechanism of nitrate stress affecting the liver function of C. striata.

Comparative proteomics study of exosomes in Vibrio harveyi and Vibrio anguillarum

Exosomes are a class of extracellular vesicles released by bacteria and contain diverse biomolecules. In this study, we isolated exosomes from Vibrio harveyi and Vibrio anguillarum, which are both serious pathogens in mariculture, using a supercentrifugation method, and the proteins in the exosomes of these two vibrios were analyzed by LC-MS/MS proteomics. Exosome proteins released by V. harveyi and V. anguillarum were different; they not only contained virulence factors (such as lipase and phospholipase in V. harveyi, metalloprotease and hemolysin in V. anguillarum), but also participated in the important life activities of bacteria (such as fatty acid biosynthesis, biosynthesis of antibiotics, carbon metabolism). Subsequently, to verify whether the exosomes participated in bacterial toxicity, after Ruditapes philippinarum was challenged with V. harveyi and V. anguillarum, the corresponding genes of virulence factors from exosomes screened by proteomics were tested by quantitative real-time PCR. All the genes detected were upregulated which suggested that exosomes were involved in vibrio toxicity. The results could provide an effective proteome database for decoding the pathogenic mechanism of vibrios from the exosome perspective.

Integrated microRNA-mRNA analysis reveals a possible molecular mechanism of enteritis susceptibility in Litopenaeus vannamei

Enteritis is one of the main diseases affecting Pacific whiteleg shrimp (Litopenaeus vannamei) in recent years, and it has resulted in huge losses to the aquaculture industry. Prior to this study, the molecular mechanism underlying enteritis in L. vannamei was unclear, and comprehensive multi-omics analysis had not been conducted. In this study, 1209 differentially expressed genes (DEGs) were identified from the hepatopancreas of L. vannamei with and without enteritis. Kyoto Encyclopedia of Genes and Genomes analysis showed that genes were significantly enriched in immune, metabolic, and endocrine regulatory pathways. Forty-eight significantly different microRNAs (miRNAs) were identified in the miRNA-Seq analysis. Further functional annotation analysis showed that the regulatory pathway of target gene enrichment of differentially expressed miRNAs was consistent with DEGs. Through miRNA-mRNA integration analysis, 47 meaningful miRNA-mRNA pairs were obtained, of which melanogenesis and pancreatic secretion were considered key pathways. Subsequent miRNA-mRNA interaction network analysis revealed that mja-miR-6493-3p, Mja-miR-6494, novel-198, novel-272, novel-261, novel-200, novel-183, novel-184, novel-237, and novel-192 may be key miRNAs involved in the regulation of these two signaling pathways. Finally, the RAS signaling pathway was found to inhibit the translation level of proteins in the hepatopancreas. These results suggest that target gene integration analysis of mRNA-miRNA can reveal the molecular mechanism underlying enteritis in L. vannamei and also provide valuable new insights for resisting enteritis.

Identification and evaluation of potential probiotics against skin-ulceration disease in the Chinese tongue sole (Cynoglossus semilaevis)

In this study, three highly pathogenic bacterial strains (Vibrio harveyi TB6, Vibrio alginolyticus TN1, and Vibrio parahaemolyticus TN3) were isolated from skin ulcers and intestines of diseased Chinese tongue sole (Cynoglossus semilaevis). The bacteria were investigated using hemolytic activity tests, in vitro co-culture with intestinal epithelial cells, and artificial infection of C. semilaevis. A further 126 strains were isolated from the intestines of healthy C. semilaevis. The three pathogens were used as indicator bacteria, and the antagonistic strains were identified from the 126 strains. The activities of exocrine digestive enzymes in the strains were also tested. Four strains with antibacterial and digestive enzyme activities were obtained and the best strains, Bacillus subtilis Y2 and Bacillus amyloliquefaciens Y9, were selected according to their ability to protect epithelial cells from infection. In addition, the effects of strains Y2 and Y9 at the individual level were investigated, finding that the activities of the immune-related enzymes superoxide dismutase, catalase, acid phosphatase, and peroxidase were significantly increased in the sera of the treatment group compared with the control group (p < 0.05). The specific growth rate (SGR, %) was also increased, especially in the Y2 group, and was significantly higher compared with the controls (p < 0.05). The result of the artificial infection test showed that the cumulative mortality within 72 h in the Y2 group was the lowest (50.5%), and in the Y9 group (68.5%) it was significantly lower than that in the control group (100%) (p < 0.05). Analysis of the intestinal microbial communities indicated that Y2 and Y9 could alter the composition of the intestinal flora, increasing both species richness and evenness, and inhibiting the growth of Vibrio in the intestine. These results suggested food supplemented with Y2 and Y9 could improve both immune function and disease resistance, as well as have a positive effect on the growth performance and the intestinal morphology of C. semilaevis.

Applying genetic technologies to combat infectious diseases in aquaculture

Disease and parasitism cause major welfare, environmental and economic concerns for global aquaculture. In this review, we examine the status and potential of technologies that exploit genetic variation in host resistance to tackle this problem. We argue that there is an urgent need to improve understanding of the genetic mechanisms involved, leading to the development of tools that can be applied to boost host resistance and reduce the disease burden. We draw on two pressing global disease problems as case studies-sea lice infestations in salmonids and white spot syndrome in shrimp. We review how the latest genetic technologies can be capitalised upon to determine the mechanisms underlying inter- and intra-species variation in pathogen/parasite resistance, and how the derived knowledge could be applied to boost disease resistance using selective breeding, gene editing and/or with targeted feed treatments and vaccines. Gene editing brings novel opportunities, but also implementation and dissemination challenges, and necessitates new protocols to integrate the technology into aquaculture breeding programmes. There is also an ongoing need to minimise risks of disease agents evolving to overcome genetic improvements to host resistance, and insights from epidemiological and evolutionary models of pathogen infestation in wild and cultured host populations are explored. Ethical issues around the different approaches for achieving genetic resistance are discussed. Application of genetic technologies and approaches has potential to improve fundamental knowledge of mechanisms affecting genetic resistance and provide effective pathways for implementation that could lead to more resistant aquaculture stocks, transforming global aquaculture.

The feasibility of skin mucus replacing exosome as a pool for bacteria-infected markers development via comparative proteomic screening in teleost

In fish, skin mucus forms a protective barrier between the body surface and the external water environment, thus providing the most direct and intuitive clues to monitor the subject's health condition. To explore the impact of the Vibrio harveyi pathogen on teleost, the proteome of epidermal mucus from control and sick Cynoglossus semilaevis were screened through iTRAQ followed with LC-MS/MS. 1531 credible proteins were obtained relating to structural, metabolic and immunological functions. 335 different expressed proteins (DEPs) were identified, with 166 up-regulated and 169 down-regulated in MS. 62 proteins were characterized, including 22 up-regulated proteins and 40 down-regulated proteins. Integrated analysis of DE-miRNAs and DEPs from miRomics and proteomics were conducted to show the indirect regulatory relationship. Comparative analysis of DEPs between mucus and exosomes demonstrated that exosomes contributed the most DEPs of all mucus DEPs. 125 proteins are DEPs only in exosomes, which presented minor difference in total mucus. Expression of Aminopeptidase (anpep), Calcium-transporting ATPase, Histone H2B and H2A confirmed implied fine discriminative power with infected C. semilaevis, among which Calcium-transporting ATPase and H2B also appeared in list of exosomal markers. This study might shed the light on effective biomarker digging at other extended screening scenarios.

Transcriptome analysis of the potential molecular mechanism of enteritis in intestine and hepatopancreas of Litopenaeus vannamei

The prawn, Litopenaeus vannamei (L. vannamei), is the most widely farmed species in the world but the incidence of enteritis in L. vannamei has increased in recent years. However, the pathogenesis of enteritis remains unclear. In this study, high-throughput sequencing was used to analyze the hepatopancreatic and intestinal transcriptome of healthy and enteritis-affected individuals from the same pond. In total, 1209 and 1608 differently-expressed genes (DEGs) were detected in the hepatopancreatic and intestinal transcriptomes, respectively. Significantly changed genes were enriched in the intestinal immune network for IgA Production, Lysosomes, Sphingolipid Metabolism and the Peroxisome Signaling Pathway. Expression of the integrin α4β7 gene was significantly increased in the intestine of L. vannamei with enteritis, while expression of 38 DEGs associated with the lysosome was significantly down-regulated. Furthermore, the expression of sphingolipid metabolism-related enzymes and superoxide dismutase (SOD) genes was also significantly decreased, indicating that abnormal autoimmune function, weak intestinal resistance to external pathogenic microbial invasion, and self-healing ability were important factors associated with enteritis in L. vannamei. In addition, the expression of trypsin and pancreatic lipase was decreased in the hepatopancreas of L. vannamei with enteritis. This study provided new insights into the possible molecular pathogenesis of enteritis in L. vannamei.

Transcriptome analysis in the spleen of Northern Snakehead (Channa argus) challenged with Nocardia seriolae

Northern snakehead (Channa argus) is an indigenous fish species and is one of popularly cultured snakeheads in China and other Asian countries. Unfortunately, Nocardia seriolae infections have caused considerable losses in the snakehead aquaculture industry. However, the infectivity and the immune response induced by N. seriolae in snakehead are unclear. In order to better understand the immune response of Northern snakehead in a series of time points after N. seriolae challenge, we conducted the transcriptomic comparison in snakehead spleen at 48, 96, and 144 h after the challenge of N. seriola against their control counterparts. Gene annotation and pathway analysis of differentially expressed genes (DEGs) were carried out to understand the functions of the DEGs. Additionally, protein-protein interaction networks were conducted to obtain the interaction relationships of immune-related DEGs. These results revealed the expression changes of multiple DEGs and signaling pathways involved in immunity during N. seriolae infection, which will facilitate our comprehensive understanding of the mechanisms involved in the immune response to bacterial infection in the northern snakehead.

Molecular Identification of Nocardia seriolae and Comparative Analysis of Spleen Transcriptomes of Hybrid Snakehead (Channa maculata Female × Channa argus Male) With Nocardiosis Disease

Hybrid snakehead (Channa maculata female × Channa argus male) is a new freshwater aquaculture fish species in southern China. During intensive aquaculture, hybrid snakeheads are often infected by Nocardia seriolae. In this study, hybrid snakehead infected suspiciously by N. seriolae in an artificial breeding pond were examined. Diseased hybrid snakeheads swam slowly without food intake, and the clinical symptoms included skin wound, anal swelling and ascites, and white granulomatous in liver, spleen, and kidney of fish. Through bacterial isolation, 16S rDNA sequencing, fluorescence in situ hybridization (FISH) and artificial infection experiment, the pathogen was identified as N. seriolae. Furthermore, the spleen samples from diseased and healthy male hybrid snakeheads in the same pond were used for RNA-Seq analysis. A total of 3,512 unique transcripts (unigenes) were identified as differentially expressed genes (DEGs), and 1,886 of them were up-regulated in diseased fish. The expression patterns of 20 DEGs were verified by quantitative polymerase chain reaction (qPCR). Several immune-related pathways and many immune-related genes were identified. qPCR results showed that the expression patterns of immune-related genes in the liver and kidney of diseased fish were comparable to that in the spleen. This study provides deep-sequencing data of hybrid snakehead spleen and will help understand the immune response of hybrid snakehead to N. seriolae. It is also helpful for the biomarker screening of fish-borne Nocardia spp. and the breeding of nocardiosis-resistant fish species.

RNA-Sequencing Analysis of the Spleen and Gill of Takifugu rubripes in Response to Vibrio harveyi Infection

Takifugu rubripes is commonly subjected to the disease-causing bacterium, Vibrio harveyi. However, the mechanism involved in the immune response of T. rubripes to V. harveyi infection is unclear. We conducted a transcriptomic analysis of the spleen and gill from T. rubripes infected with V. harveyi. We obtained 60,981,357 and 60,760,550 clean reads from the control and infected spleens, and 57,407,586 and 57,536,651 clean reads from the control and infected gills, respectively. We also identified 1,560 and 1,213 differentially expressed genes in the spleen and gill, respectively. Gene ontology analysis revealed that the most enriched biological process in both the spleen and gill was "immune response". The most enriched Kyoto Encyclopedia of Genes and Genomes immune response-related pathways were the NOD-like receptor signaling pathway in the spleen and cytokine-cytokine receptor interaction in the gill. We found 10 candidate immune-related genes in the spleen and gill. These putative immune pathways and candidate genes will provide insight into the immune response mechanisms of T. rubripes against V. harveyi.

Mucus piRNAs profiles of Vibrio harveyi-infected Cynoglossus semilaevis: A hint for fish disease monitoring

The half-smooth tongue sole, Cynoglossus semilaevis, is an important cultured flatfish species. Vibrio harveyi is a common pathogen to this fish, which may result in great economic loss to C. semilaevis culture industry. piRNAs, a non-coding RNAs with 26-32 nt, have been regarded as promising biomarkers for cancer diagnosis and fish diseases. Here, we extracted the RNA from mucus of C. semilaevis and constructed the differential expression profiles of piRNAs between the sick fish (MS) and healthy fish (MC). We identified 45,696 differentially expressed piRNAs including 22,735 up-regulated piRNAs and 22,961 down-regulated piRNAs in MS group compared with MC group. The GO enrichment and KEGG pathway enrichment analyses of the differential piRNAs were carried out. The result showed immunity-related target genes mainly involved in immune system process, response to stimulus, cell killing, immune system, infectious diseases and cell growth and death. The 10 most differentially expressed piRNAs were chosen to perform the qRT-PCR, while only seven piRNAs were consistent with the sequence result. Compared with MC group, the expression levels of piR-mmu-72173>piR-rno-62831>piR-xtr-704880, piR-dme-15546979, piR-mmu-49941660, piR-mmu-29283297 and piR-mmu-1758399 were significantly lower, and piR-gga-10574 and piR-gga-134812 were significantly higher in MS group. These piRNAs may be potential biomarkers during the V. harveyi infection of C. semilaevis. This study could provide a new method to identify the infection status of C. semilaevis and understand better about the innate and adaptive immune system in C. semilaevis during bacterial infection.

Transcriptome analysis revealed multiple immune processes and energy metabolism pathways involved in the defense response of the large yellow croaker Larimichthys crocea against Pseudomonas plecoglossicida

The large yellow croaker (Larimichthys crocea) aquaculture industry is suffering substantial financial losses caused by visceral white nodules disease resulting from Pseudomonas plecoglossicida infection. However, how L. crocea responds to P. plecoglossicida infection remains largely unknown. Here, we characterized the changes in the mRNA profile in the spleen of L. crocea upon P. plecoglossicida infection and explored the related defensive strategies. Sample clustering analysis and qRT-PCR indicated that P. plecoglossicida induced profound and reproducible transcriptome remodeling in the L. crocea spleen. Many innate immune-related genes, such as IL-17 signaling molecules, chemokines and chemokine receptors, complement components, TLR5 signaling molecules, and antimicrobial peptide hepcidins (Hamps), were upregulated by P. plecoglossicida and may play important roles in the L. crocea defense against P. plecoglossicida. The antibacterial activity of Hamp2-5 against P. plecoglossicida was further confirmed by using synthetic mature peptide of Hamp2-5. Additionally, significant enrichment of "Glycolysis/Gluconeogenesis", "Citrate cycle" and "Oxidative phosphorylation" pathways and a significant upregulation of all 6 rate-limiting enzyme genes (HK1, PFK, PKM, CS, IDH2, DLST) in the Glycolysis and Citrate cycle pathways in P. plecoglossicida-infected fish suggested that ATP synthesis may be accelerated to ensure energy supply in response to pathogenic infection. Altogether, our results not only identified the key immune-related genes and immune pathways that participated in the defense response of L. crocea against P. plecoglossicida, but also revealed a novel defensive strategy involving ATP synthesis in this species.

Immune Response and Apoptosis-Related Pathways Induced by Aeromonas schubertii Infection of Hybrid Snakehead (Channa maculata♀ × Channa argus♂)

Aeromonas schubertii is the etiological pathogen of internal organ nodules in snakehead fish. Infections with A. schubertii produce a significant economic loss in aquaculture. Therefore, it is important to examine the immune mechanisms by which snakeheads defend against A. schubertii infection. In this study, we established a hybrid snakehead infection model by intraperitoneal injection of A. schubertii that produced internal organ nodules. The splenic immune response of infected fish was examined at the transcriptome level by Illumina-seq analysis. Results showed 14,796 differentially expressed genes (DEGs) following A. schubertii infection, including 4441 up-regulated unigenes and 10,355 down-regulated unigenes. KEGG analysis showed 2084 DEGs to be involved in 192 pathways, 14 of which were immune-related. Twelve DEGs were used to validate quantitative real-time PCR results with RNA-seq data. Time-course expression analysis of six genes demonstrated modulation of the snakehead immune response by A. schubertii. Furthermore, transcriptome analysis identified a substantial number of DEGs that were involved in the apoptosis signaling pathway. TUNEL analysis of infected spleens confirmed the presence of apoptotic cells. This study provided new information for a further understanding of the pathogenesis of A. schubertii in snakeheads, which can be used to prevent and possibly treat A. schubertii infections.

Analysis of lncRNA and mRNA expression profiles in peripheral blood leukocytes of the half-smooth tongue sole (Cynoglossus semilaevis) treated with chitosan oligosaccharide

Long noncoding RNAs (lncRNAs) play a multifaceted role in transcriptional regulation and are important regulators of immune function. Scarce information is available regarding lncRNAs in fish. Peripheral blood mononuclear cells participate in the immune response of fish and aid resistance to infection with pathogenic microorganisms. Chitosan oligosaccharide can improve cellular and humoral immunity to enhance disease resistance in fish. In this study, we obtained peripheral blood leukocytes from half-smooth tongue sole and studied the effect of chitosan oligosaccharide on the lncRNA-mRNA expression profile of these cells using high-throughput sequencing and bioinformatics techniques. A total of 609 differentially expressed mRNAs and 50 differentially expressed lncRNAs were identified. The GO term enrichment analysis of the differentially expressed genes was annotated by 220 GO terms, 137 biological processes, 18 cellular components, and 65 molecular functions. Sixteen KEGG pathways, including immune signaling pathways, metabolism, and genetic information processing, were significantly enriched in differentially expressed genes. Thirty-six differentially expressed lncRNAs and 32 differentially expressed mRNAs produced a coexpression network containing 90 relationship pairs. The prediction of lncRNA target genes revealed 244 lncRNAs that potentially cis-regulated 294 differentially expressed mRNAs. qPCR verified that the expression levels of 17 differentially expressed lncRNAs and 15 differentially expressed mRNAs were consistent with the RNA-Seq results. Among them, 6 lncRNAs and 7 mRNAs were differentially expressed genes obtained from the prediction and analysis of lncRNA target genes, and 8 lncRNAs and 4 mRNAs were differentially expressed genes that participated in the construction of the coexpression network. In peripheral blood leukocytes after chitosan oligosaccharide treatment, as well as in peripheral blood and spleen after Vibrio anguillarum stimulation, lncRNAs and mRNAs showed significant differential expression. The results indicated that they may be related to the immune response, providing novel reference information for further research on the role of lncRNAs in immune regulation in half-smooth tongue sole.

Combining of transcriptomic and proteomic data to mine immune-related genes and proteins in the liver of Cynoglossus semilaevis challenged with Vibrio anguillarum

The liver is a multi-functional organ including metabolism, substance synthesis, detoxification, and various immune functions, and its role in immunity has attracted more and more attention. However, research on the liver immune response of fish infected by pathogenic bacteria is currently lacking. In this study, the transcriptomics and proteomics of the liver of Cynoglossus semilaevis infected with Vibrio anguillarum were analyzed. A total of 1470 genes and 497 proteins were differentially expressed in the pairwise comparison of obvious symptoms of infection (HOSG), no obvious symptoms of infection (NOSG) and PBS treatment (CG). Gene ontology and KEGG enrichment pathways analysis showed that differentially expressed genes (DEGs) and differentially expressed proteins (DEPs) were mainly enriched in toll-like receptors (TLRs), complement and coagulation cascades, nucleotide oligomerization domain (NOD)-like receptors (NLRs), mitogen-activated protein kinase (MAPK) and phagosome signaling pathways, which suggested the combined action of the five pathways were significant to enhance the liver immune defense. The combination of transcriptomic and proteomic analysis showed that ITGβ1, C3, C5 and MRC1 were significantly up-regulated, which might play an important role in the liver immune response to the recognition of V. anguillarum, inflammatory response and phagocytosis. The transcriptome and proteome data we obtained provide information on some key genes and proteins for further study of the mechanism of liver immune response.

Proteomics of mucosal exosomes of Cynoglossus semilaevis altered when infected by Vibrio harveyi

The cargo of exosomes contains proteins with various functions, which might be promising biomarkers for disease diagnosis and prognosis. To explore the impact of the Vibrio harveyi pathogen on Cynoglossus semilaevis from a different perspective and develop promising biomarkers for infection, the exosomes from epidermal mucus of healthy controls（EC）and sick fish（ES）were extracted and identified, coupled with proteomic screening through iTRAQ followed with LC-MS/MS. 1531 credible proteins were obtained relating to structural, metabolic and immunological functions. 359 different expressed proteins (DEPs) (FC > 2 or FC < 0.5) were found, with 161 up-regulated and 198 down-regulated in ES. Based on the database of C. semilaevis on Uniprot, 71 proteins were characterized as concrete names, including 19 up-regulated proteins and 52 down-regulated proteins, and were selected as subjects for further studies. Ferritin, Toll-like receptor 5S protein and Calcium-transporting ATPase were upregulated, while Histone H2B and Eukaryotic translation initiation factor 5A were downregulated, consistent with the expression levels of related mRNAs in skin tissue verified by qRT-PCR. The integrated analysis between miRomics and proteomics also provided possible regulatory relationships mediated by mucous exosomes during infection. The signature proteins in mucosal exosomes could make sense in the explanation of the infection defending mechanism and the development of biomarkers which can differentiate diseased and healthy C. semilaevis individuals.

Extracellular vesicles piwi-interacting RNAs from skin mucus for identification of infected Cynoglossus semilaevis with Vibrio harveyi

Extracellular vesicles play a regulatory role in intracellular and intercellular transmission through a variety of biological information molecules, including mRNA, small RNAs and proteins. piRNAs are one kind of regulatory small RNAs in the vesicles at the post transcriptional level. Hereby, we isolated the extracellular vesicles from skin mucus and screened the piRNA profiles of these vesicles, aiming at developing biomarkers related to bacterial infections in Cynoglossus semilaevis. The different profilings of piRNAs in mucous extracellular vesicles of C. semilaevis were compared through small RNA sequencing, between fish infected with Vibrio harveyi and healthy ones. The number of clean reads on the alignment of exosome sick (ES) group was 105, 345 and that of exosome control (EC) group was 455, 144. GO and KEGG pathway enrichment analysis showed that most of the target genes were involved in cellular process, response to stimulus, biological regulation, immune system process and signal transduction, signal molecular and interaction, transport and catabolism. The 45 final candidate piRNAs related to immunity or infectious diseases included 20 piRNAs with high expression in the ES group and 25 piRNAs with a low expression in the ES group. After verification by qRT-PCR, there was significant difference of five piRNAs expression level between infected fish and healthy fish, in line with the sequencing. The expression level of piR-mmu-16401212, piR-mmu-26829319 and piR-gga-244092 in infected fish were significantly lower than that of control group, while piR-gga-71717 and piR-gga-99034 were higher, which implying that these piRNAs in mucous extracellular vesicles can be used to identify diseased fish from normal ones. This work supplied a novel class of biomarker for infection diagnosis in fish, and it will be benefit for screening disease resistant breeding of C. semilaevis.

Transcriptome analysis in Takifugu rubripes and Dicentrarchus labrax gills during Cryptocaryon irritans infection

Takifugu rubripes and Dicentrarchus labrax are important commercial fish in China that are under serious threat from Cryptocaryon irritans. C. irritans is a ciliated obligate parasite that causes marine white spot disease and leads to heavy economic losses. We analysed the transcriptome in the gills of T. rubripes and D. labrax to compare differentially expressed genes (DEGs) and pathways during infection with C. irritans. In total, we identified 6,901 and 35,736 DEGs from T. rubripes and D. labrax, respectively. All DEGs were annotated into GO terms; 6,901 DEGs from T. rubripes were assigned into 991 sub-categories, and 35,736 DEGs from D. labrax were assigned into 8,517 sub-categories. We mapped DEGs to the KEGG database and obtained 153 and 350 KEGG signalling pathways from T. rubripes and D. labrax, respectively. Immune-related categories included Toll-like receptors, MAPK, lysosome, C-type lectin receptor and NOD-like receptor signalling pathways were significantly enriched pathways. In immune-related signalling pathways, we found that AP-1, P38, IL-1β, HSP90 and PLA were significantly up-regulated DEGs in T. rubripes, but P38 and PLA were significantly down-regulated in D. labrax. In this study, transcriptome was used to analyse the difference between scaly and non-scaly fish infection by C. irritans, which not only provided a theoretical basis for the infection mechanism of C. irritans, but also laid a foundation for effectively inhibiting the occurrence of this disease. Our work provides further insight into the immune response of host resistance to C. irritans.

The Dynamic Immune Response of Yellow Catfish (Pelteobagrus fulvidraco) Infected With Edwardsiella ictaluri Presenting the Inflammation Process

Edwardsiella ictaluri is a highly destructive pathogen in cultured yellow catfish, thus it was very necessary to study the immune response of yellow catfish against bacterial infection. In this study, RNA-Seq technology was used to study the immune response in two distinct tissues of yellow catfish at eight different time points (h) after E. ictaluri infection. The number of differentially expressed genes (DEGs) in the spleen and liver was low at 3 h and 6 h post-infection, respectively. Afterwards, the most number of DEGs in the spleen was detected at 72 h, while the number of DEGs in the liver maintained a high level from 24 h to 120 h. The GO and KEGG enrichment analyses of DEGs at different time points uncovered that cytokines were continuously transcribed at 6 h to 120 h; whereas the liver is the main organ that secretes the components of the complement system, and metabolic regulation was activated from 12 h to 120 h. Moreover, an overview of the inflammation response of yellow catfish was exhibited including pattern-recognition receptors, inflammatory cytokines, chemokines, complements, and inflammation-related signal pathways. The similar expression tendency of nine genes by qRT-PCR validated the accuracy of transcriptome analyses. The different transcriptomic profiles obtained from the spleen and liver will help to better understand the dynamic immune response of fish against bacterial infection, and will provide basic information for establishing effective measures to prevent and control diseases in fish.

Environmental salinity influences the branchial expression of TCR pathway related genes based on transcriptome of a catadromous fish

Environmental salinity not only affects the physiological processes such as osmoregulation and hormonal control, but also changes the immune system in fishes. Studies are limited in fish on the roles of the T cell receptor (TCR)-related genes in relation to changes in environmental salinity. A large group of salinity-challenged transcripts was obtained in gills of marbled eel (Anguilla marmorata). Moreover, bioinformatic ways were used to identify the enriched TCR pathway related genes which were significantly different expressed in fresh water (FW), brackish water (BW) and seawater (SW). Meanwhile, the RT-qPCR results were validated and consistent with the RNA-seq results. TCR a, TCR b, CD45, CD28, PI3K, LCK and LAT were up-regulated when the salinity increases in BW and SW, which connected with the related signaling pathways (Ras-MAPK and PKC pathway). CD4 and Zap70 were down-regulated when the salinity increases in BW and SW, which connected with the PLC pathway. The research offers a novel viewpoint to explore the immune pathways including the TCR pathway in fish based on transcriptome.

Identification of reproduction-related genes and pathways in the Culter alburnus H-P-G axis and characterization of their expression differences in malformed and normal gynogenetic ovaries

This study applied RNA-seq technology to discover reproduction-related genes and pathways in female topmouth culter brain (including pituitary) and ovarian tissues. In functional analysis, 2479 and 2605 unigenes in the brain and ovary tissue were assigned to the "reproductive process" subcategory in addition to the 2660 and 2845 unigenes assigned to the "reproduction" subcategory. Twenty-three complete cDNA sequences were identified through the different gene expression (DGE) approach from five reproduction-related pathways (MAPK signaling pathway, neuroactive ligand-receptor interaction pathway, gonadotropin-releasing hormone signaling pathway, oocyte meiosis pathway, and steroid biosynthesis pathway). The expression levels of 16 candidate genes using qPCR in this study were in accordance with the results of transcriptome analysis. In addition, the expression levels of the FSH, 3β-HSD, PGR, and NPYR genes in malformed gynogenetic ovaries were considerably low, which was consistent with the progress of oocytogenesis in the ovaries of topmouth culter. The high expression of these four genes in the ovaries of normal topmouth culter suggested they might involve in the preparation for the shift of oogenesis to ovulation. Hence, our work identified a set of annotated gene products that are candidate factors affecting reproduction in the topmouth culter H-P-G axis. These results could be essential for further research in functional genomics and genetic editing for topmouth culter reproduction.

Comparative transcriptome analysis of scaled and scaleless skins in Gymnocypris eckloni provides insights into the molecular mechanism of scale degeneration

Background

The scale degeneration is thought to be related to the adaptation to the extreme environment with cold climate and high-altitude in schizothoracine fishes. Gymnocypris eckloni, a schizothoracine fish living in plateau waters with the elevation above 2500 m, is nearly esquamate and only covered with shoulder scales and anal scales, making it a good model species to study the molecular mechanism of scale degeneration.

Results

The transcriptomes of shoulder scaled skins (SSS), anal scaled skins (ASS) and scaleless skins (NSS) were sequenced and analyzed in G. eckloni at the age of 1 year. Histological examination showed that shoulder scale had completed its differentiation and anal scale just initiated the differentiation. A total of 578,046 unigenes were obtained from the transcriptomes, with 407,799 unigenes annotated in public databases. A total of 428 and 142 differentially expressed unigenes (DEUs) were identified between SSS and NSS, and between ASS and NSS, respectively, with 45 DEUs that were overlapped. Annotation analysis indicated that these DEUs were mainly enriched in Gene Ontology (GO) terms and KEGG pathways associated with bone and muscle formation, such as myofibril, contractile fiber, cytoskeletal protein binding, muscle structure development, cardiac muscle contraction, hypertrophic cardiomyopathy (HCM) and calcium signaling pathway.

Conclusions

Our results would provide insights into the molecular mechanisms of scale degeneration in G. eckloni and other congeneric fishes. In addition, the transcriptome data provides candidate genes and markers for future studies.

Detecting Cynoglossus semilaevis infected with Vibrio harveyi using micro RNAs from mucous exosomes

Exosomes are important mediators of vesicle transportation and contain microRNAs (miRNAs) that mediate transcriptional gene knockout and silencing in biological processes. Moreover, exosomic miRNAs are promising biomarkers for disease diagnosis and physiological status indication in many species, including fish. The impact of the Vibrio harveyi pathogen on Cynoglossus semilaevis aquaculture is becoming more and more serious as the industry expands. To overcome this challenge, miRNAs in mucous exosomes were screened by small RNA sequencing and verified by quantitative real-time PCR to develop biomarkers. This is the first capture of exosomes from flatfish mucus coupled with miRNA profiling. The results revealed significant differences in expression levels of some miRNAs between infected and healthy fish. Three unique miRNAs were identified for V. harveyi infection diagnosis; expression levels of dre-miR-205-5p and dre-miR-205-5p in infected fish were significantly lower than controls, while dre-miR-100-5p expression was higher. These miRNAs in mucous exosomes could be used to differentiate diseased and healthy fish in an early screening method with practical value for breeding disease-resistant C. semilaevis.

Cloning, expression profile of the complement component C9 gene and influence of the recombinant C9 protein on peripheral mononuclear leukocytes transcriptome in half-smooth tongue sole (Cynoglossus semilaevis)

The ninth complement component (C9) is a terminal complement component (TCC) that is involved in creating the membrane attack complex (MAC) on the target cell surface. In this study, the CsC9 (C9 of Cynoglossus semilaevis) cDNA sequence was cloned and characterized. The full-length CsC9 cDNA measured 2,150 bp, containing an open reading frame (ORF) of 1,803 bp, a 5'-untranslated region (UTR) of 24 bp and a 3'-UTR of 323 bp. A domain search revealed that the CsC9 protein contains five domains, including two TSP1s, an LDLRA, an EGF, and a MACPF. Quantitative real-time PCR analysis showed that CsC9 at the mRNA level was expressed in all the tested tissues, with the highest expression being observed in the liver. CsC9 expression is significantly upregulated in the tested tissues after challenge with Vibrio anguillarum. To further characterize the role of CsC9, peripheral blood mononuclear cells of C. semilaevis were used for transcriptome analysis after incubation with recombinant CsC9 (rCsC9) protein. A total of 3,775 significant differentially expressed genes (DEGs) were identified between the control and the rCsC9-treated group, including 2,063 upregulated genes and 1,712 downregulated genes. KEGG analyses revealed that the DEGs were enriched in cell adhesion molecules, cytokine-cytokine receptor interactions, T cell receptor signaling pathways, B cell receptor signaling pathways and Toll-like receptor signaling pathways. The results of this study indicate that in addition to participating in MAC formation, CsC9 might play multiple roles in the innate and adaptive immunity of C. semilaevis.

Differentially expressed proteins in the intestine of Cynoglossus semilaevis Günther following a Shewanella algae challenge

Fish intestine is an important constituent of the mucosal immune system. The gut and gut-associated lymphoid tissue construct a local immune environment. A Shewanella algae strain was previously reported to be a pathogen causing ascitic disease accompanied with intestinal inflammation in Cynoglossus semilaevis. This study aimed to investigate the intestine immune response in C. semilaevis to S. algae infection at the protein level. Two-dimensional electrophoresis coupled with mass spectrometry proteomics was utilized to compare protein expression in the intestines from normal and S. algae-infected C. semilaevis. A total of 70 differentially expressed proteins (DEPs), consisting of 16 upregulated and 54 downregulated proteins, were identified in the intestine tissue of C. Semilaevis. These protein expression changes were further validated using western blot analysis and quantitative real-time PCR. Gene ontology enrichment analysis showed that these 70 DEPs could be assigned across three categories: "cellular components", "molecular function", and "biological process". Forty-one DEPs (six up-regulated and 35 down-regulated proteins) related to metabolic processes were identified. In addition, 20 DEPs (eight up-regulated and 12 down-regulated proteins) related to stress and immune responses were identified. A protein-protein interaction network generated by the STRING (Search Tool for the Retrieval of Interacting Genes/protein) revealed that 30 DEPs interacted with one another to form an integrated network. Among them, 29 DEPs were related to stress, immune, and metabolism processes. In the network, some of the immune related proteins (C9, FGB, KNG1, apolipoprotein A-IV-like, and PDIA3) were up-regulated and most DEPs involved in metabolism processes were down-regulated. These results indicate that the immune defense response of the intestine was activated and the intestinal function associated with metabolism processes was disturbed. This study provides valuable information for further research into the functions of these DEPs in fish.

Integrative Transcriptomic and Small RNA Sequencing Reveals Immune-Related miRNA-mRNA Regulation Network for Soybean Meal-Induced Enteritis in Hybrid Grouper, Epinephelus fuscoguttatus♀ × Epinephelus lanceolatus♂

A 10-week feeding experiment was conducted to reveal the immune mechanism for soybean meal-induced enteritis (SBMIE) in hybrid grouper, Epinephelus fuscoguttatus ♀ × Epinephelus lanceolatus ♂. Four isonitrogenous and isolipidic diets were formulated by replacing 0, 10, 30, and 50% fish meal protein with soybean meal (namely FM, SBM10, SBM30, and SBM50, respectively). The weight gain rate of the SBM50 group was significantly lower than those of the other groups. Plica height, muscular layer thickness, and goblet cells of the distal intestine in the SBM50 group were much lower than those in the FM group. The intestinal transcriptomic data, including the transcriptome and miRNAome, showed that a total of 6,390 differentially expressed genes (DEGs) and 92 DEmiRNAs were identified in the SBM50 and FM groups. DEmiRNAs (10 known and 1 novel miRNAs) and their DE target genes were involved in immune-related phagosome, natural killer cell-mediated cytotoxicity, Fc gamma R-mediated phagocytosis, and the intestinal immune network for IgA production pathways. Our study is the first to offer transcriptomic and small RNA profiling for SBMIE in hybrid grouper. Our findings offer important insights for the understanding of the RNA profile and further elucidation of the underlying molecular immune mechanism for SBMIE in carnivorous fish.

Gene network analysis reveals a core set of genes involved in the immune response of Japanese flounder (Paralichthys olivaceus) against Vibrio anguillarum infection

Japanese flounder (Paralichthys olivaceus) is one of the most economically important marine fish cultured in north Asia. Vibrio anguillarum is a severe bacterial pathogen to Japanese flounder and many other aquaculture species. In order to understand the immune response of flounder during bacterial infection, we systematically examined the transcriptome profiles of flounder spleen at three time points after V. anguillarum challenge. More than one billion high quality reads were obtained, approximately 80.70% of which were successfully mapped to the reference genome of flounder. A total of 6060, 4688 and 4235 differentially expressed genes (DEGs) were captured at 6, 12 and 24-h post-infection, respectively. The DEGs exhibited dynamic changes in expression and were assigned into four different profiles based on expression trend. GO and KEGG analysis showed that the DEGs were enriched in various immune-related terms, including response to stimulation, immune system and pathways of cytokine-cytokine receptor interaction, Jak-STAT signaling and Toll-like receptor signaling. Furthermore, a network of highly interactive DEGs involved in 11 immune-related pathways was detected by utilizing the weighted co-expressing network analysis (WGCNA). Accordingly, 26 hub genes were discovered that constituted an elaborate immune regulatory network and functioned mainly in pathogen recognition, antigen processing, and molecular signaling. The results of this study provided the first systematical transcriptome profile of flounder in association with V. anguillarum infection and can serve as a valuable resource of target genes for future studies on the molecular mechanisms underlying the immune defense of flounder against bacterial infection.

Transcriptomic Changes in Liver of Juvenile Cynoglossus semilaevis following Perfluorooctane Sulfonate Exposure

Perfluorooctane sulfonate (PFOS) is an increasingly important environmental pollutant that is pervasive in the environment. A number of studies have focused on the toxicological effects of PFOS on model fish species (zebrafish and medaka), but little is known about the impact of PFOS on commercially important marine fish. Thus, the present study examined transcriptome responses to PFOS exposure in the liver of juvenile Cynoglossus semilaevis, an important farmed flatfish in China. Then, in response to PFOS challenges, 1695 and 5244 genes were identified as significantly increased and depressed, respectively. Significant expression changes were observed in immune-related genes (cytokine-cytokine receptor interaction, T-helper [Th]17 cell differentiation, and the chemokine nuclear factor-kappa B and T-cell receptor signaling pathways), indicating that immunotoxicity is a key aspect of the effects of PFOS on C. semilaevis. Exposure to PFOS also altered the gene expression levels of hormones (inhibin, insulin, somatostatin, and glucagon), which could lead to severe metabolic and endocrine dysfunction. As expected from previous studies, several phase I and phase II detoxification enzymes were significantly up-regulated, which could facilitate the biotransformation and detoxification of PFOS in C. semilaevis. The present study provides new insights into the molecular toxicology of PFOS in a commercially important fish species. Environ Toxicol Chem 2020;39:556-564. © 2019 SETAC.

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

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

Genome-Wide Association Mapping and Gene Expression Analyses Reveal Genetic Mechanisms of Disease Resistance Variations in Cynoglossus semilaevis

The sustainable development of aquaculture has been impeded by infectious diseases worldwide. However, the genomic architecture and the genetic basis underlying the disease resistance remain poorly understood, which severely hampers both the understanding of the evolution of fish disease resistance traits and the prevention of these diseases in the aquaculture community. Cynoglossus semilaevis is a representative and commercially-important flatfish species. Here we combined genome-wide association study and Fst and nucleotide diversity filtration to identify loci important for the disease resistance. Based on 1,016,774 single-nucleotide polymorphisms (SNPs) identified from 650 Gb genome resequencing data of 505 individuals, we detected 33 SNPs significantly associated with disease resistance and 79 candidate regions after filtration steps. Both the allele frequencies and genotype frequencies of the associated loci were significantly different between the resistant and susceptible fish, suggesting a role in the genetic basis of disease resistance. The SNP with strongest association with disease resistance was located in Chr 17, at 145 bp upstream of fblx19 gene, and overlapped with the major quantitative trait locus previously identified. Several genes, such as plekha7, nucb2, and fgfr2, were also identified to potentially play roles in the disease resistance. Furthermore, the expression of some associating genes were likely under epigenetic regulations between the bacterial resistant and susceptible families. These results provide insights into the mechanism that enable variation of disease resistance to bacterial pathogen infection. The identified polymorphisms and genes are valuable targets and molecular resources for disease resistance and other traits, and for advanced breeding practice for superior germplasm in fish aquaculture.

Transcriptome analysis and the effects of polyunsaturated fatty acids on the immune responses of the critically endangered angtze sturgeon (Acipenser dabryanus)

The poor understanding of nutrition needed has become a significant obstruction to artificial conservation of Yangtze sturgeon (Acipenser dabryanus) and the relationship between ployunsaturated fatty acid nutrition and the immune response of Yangtze sturgeon is remains unclear. To explore this relationship, the immune response was determined by the activities of serum immune-related enzymes and the transcriptome pattern in the spleen after feeding different fat source diets for 7 weeks. In addition, the gene expression pattern after a lipopolysaccharide (LPS) challenge was investigated in the presence of docosahexaenoic acid (DHA) and eicosapentaenoic acid (EPA). Long-term feeding of the fish oil diets increased the serum immune-related enzyme activities, including lysozyme, acid phosphatase, and alkaline phosphatase of Yangtze sturgeon. More than 653,999 transcripts with an N50 length of 1047 bp were obtained and a final set of 280,408 unigenes was generated. After annotating the unigenes, 3549 genes were assigned to the immune system and 2839 were identified to participate in the response to the different fat sources. A transcriptome assay showed the fish oil diets moderately upregulated immune-related signaling pathways in the spleen of Yangtze sturgeon, including NLR signaling, platelet activation, Fc gamma R-mediated phagocytosis, Th17 cell differentiation, and Th1 and Th2 cell differentiation. The quantitative polymerase chain reaction (qPCR) results of candidate genes for these pathways showed similar results. The LPS challenge study revealed that DHA and EPA moderately upregulated the candidate immune-related genes and modulated excessive activation of the immune pathway by the pathogen. This study confirmed the immunomodulatory function of unsaturated fatty acids in Yangtze sturgeon. This research will provide a reference for the preparation of artificial diets for Yangtze sturgeon.

Comparative transcriptome profiling of immune response against Vibrio harveyi infection in Chinese tongue sole

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

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.

T-box transcription factor eomesodermin/Tbr2 in Atlantic cod (Gadus morhua L.): Molecular characterization, promoter structure and function analysis

Eomesodermin (Eomes) is a member of T-box transcription factor family and plays an important role in the regulation of a wide variety of developmental processes and immune response in animals. Here we report cloning and characterization of the full-length cDNA of Atlantic cod Eomes (GmEomes), which possesses a TBOX_3 domain similar to its counterpart in mammals. The regulated expression was observed in head kidney and spleen in response to live Vibrio anguillarum infection in vivo, and spleen leukocytes in vitro after PMA and poly I:C stimulation. Furthermore, we determined a 694 bp sequence, upstream of the transcriptional start site (TSS), to contain a number of sequence motifs that matched known transcription factor-binding sites. Activities of the presumptive regulatory gene were assessed by transfecting different 5'-deletion constructs in CHSE-214 cells. The results showed that the basal promoters and positive transcriptional regulator activities of GmEomes were dependent by sequences located from -694 to -376 bp upstream of TSS. Furthermore, we found that some Eomes binding sites were present in the 5'-flanking regions of the cod IFNγ gene predicted by bioinformatics. However, Co-transfection of eomesodermin overexpression plasmids with INFγ reporter vector into CHSE-214 cells determined that Atlantic cod eomesodermin played a minor role in activation of the INFγ promoter.

Differences in DNA Methylation Between Disease-Resistant and Disease-Susceptible Chinese Tongue Sole (Cynoglossus semilaevis) Families

DNA methylation, the most widely studied and most well-understood epigenetic modification, has been reported to play crucial roles in diverse processes. Although it has been found that DNA methylation can modulate the expression of immune-related genes in teleosts, a systemic analysis of epigenetic regulation on teleost immunity has rarely been performed. In this research, we employed whole-genome bisulfite sequencing to investigate the genome-wide DNA methylation profiles in select disease-resistant Cynoglossus semilaevis (DR-CS, family 14L006) and disease-susceptible C. semilaevis (DS-CS, family 14L104) against Vibrio harveyi infection. The results showed that following selective breeding, DR-CS had higher DNA methylation levels and different DNA methylation patterns, with 3,311 differentially methylated regions and 6,456 differentially methylated genes. Combining these data with the corresponding transcriptome data, we identified several immune-related genes that exhibited differential expression levels that were modulated by DNA methylation. Specifically, DNA methylation of tumor necrosis factor–like and lipopolysaccharide-binding protein-like was significantly correlated with their expression and significantly contributed to the disease resistance of the selected C. semilaevis family. In conclusion, we suggest that artificial selection for disease resistance in Chinese tongue sole causes changes in DNA methylation levels in important immune-related genes and that these epigenetic changes are potentially involved in multiple immune responses in Chinese tongue sole.

Comparative transcriptomic analysis reveals the gene expression profiles in the liver and spleen of Japanese pufferfish (Takifugu rubripes) in response to Vibrio harveyi infection

Japanese pufferfish (Takifugu rubripes) is one of the main marine aquatic fish species cultured in Asia due to its high nutritional value. In recent years, disease caused by Vibrio harveyi infections have led to serious mortality in Japanese pufferfish industry. To understand the complex molecular mechanisms between V. harveyi and Japanese pufferfish, we performed a transcriptome analysis of liver and spleen samples from Japanese pufferfish at 1 and 2 day post-infection. Between-group comparisons revealed 922 genes that were significantly differentially expressed. The altered genes emphasized the function in several immune related pathways including MAPK signaling pathway, JAK-STAT signaling pathway, toll-like receptor signaling pathway, cytokine-cytokine receptor interaction and lysosomal pathway. The data generated in this study provided insight into the responses of Japanese pufferfish against V. harveyi at the transcriptome level, promoting our comprehensive understanding of immune responses for aquatic animal against V. harveyi.

Integrated transcriptome analysis of immunological responses in the pearl sac of the triangle sail mussel (Hyriopsis cumingii) after mantle implantation

For pearl culture of bivalve Hyriopsis cumingii, implantation of the sabio may cause nucleus discharge and increased host death rates. We performed a transcriptome analysis of the pearl sac of H. cumingii for 30 days after mantle implantation; 293863 unigenes were obtained, and 27176 unigenes were identified using nr, nt, KO, Swiss-Prot, Pfam, GO, and KOG databases. We detected 4878 differentially expressed genes (DEGs) through pairwise comparisons. We speculated that the physical condition of the recipient mussels returned to normal in about one month; the period was divided into six vital phases (0, 2 h-6 h, 12 h-24 h, 48 h to 7 days, 14 days and 30 days) on the basis of the overall similarities in DEGs. We compared the DEGs between time points and identified key immune-related genes. Our findings provide information on the immunological reactions induced by implantation in pearl mussels.

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.

Transcriptome profiles in the spleen of African catfish (Clarias gariepinus) challenged with Aeromonas veronii

The African catfish, Clarias gariepinus, an important cultured freshwater species in many countries, possess the characteristic of high disease resistance. However, little genomic information for this character of the fish is available up to now. To address the shortfall and to better understand C. gariepinus immune response to pathogen infection at molecular level, C. gariepinus were challenged with potent A. veronii and the high-throughput RNA sequencing (RNA-seq) technology were employed to produce transcriptomes from spleen. In total, an average of 46,073,372 clean reads obtained were de novo assembled into 156,955 unigenes with an average length of 1082 bp. All of unigenes were annotated to seven public databases. Three comparisons were separately conducted between the infected groups at 3 h, 24 h, 48 h post-challenge and control group. A total of 2482 differentially expressed unigenes (DEGs) were identified. Among these, 114 immune-related DEGs were captured, including 88, 42, and 31 genes at 3 h, 24 h and 48 h after infection respectively, for analysis of expression pattern and enrichment. The 114 DEGs displayed four expression patterns by cluster analysis and they were significantly enriched in 38 pathways (q < 0.01) related to the immune or disease, five of which were NF-kappa B, TNF, NLR, TLR and RLR pathways. Finally, the expression levels of twelve selected immune-related DEGs involved in above five pathways were scrutinized. Seven of which were up-regulated at 3 h after infection, afterward, their expression dropped to control level. In summary, this study provides valuable transcriptome resource for understanding the defense mechanisms of C. gariepinus in resistance to pathogens from the gene expression viewpoint, which also open up the possibility to study the immune complexity and to better comprehend the interrelationships between some immune pathways in C. gariepinus.

Comparative Study of Immune Reaction Against Bacterial Infection From Transcriptome Analysis

Transcriptome analysis is a powerful tool that enables a deep understanding of complicated physiological pathways, including immune responses. RNA sequencing (RNA-Seq)-based transcriptome analysis and various bioinformatics tools have also been used to study non-model animals, including aquaculture species for which reference genomes are not available. Rapid developments in these techniques have not only accelerated investigations into the process of pathogenic infection and defense strategies in fish, but also used to identify immunity-related genes in fish. These findings will contribute to fish immunotherapy for the prevention and treatment of bacterial infections through the design of more specific and effective immune stimulants, adjuvants, and vaccines. Until now, there has been little information regarding the universality and diversity of immune reactions against pathogenic infection in fish. Therefore, one of the aims of this paper is to introduce the RNA-Seq technique for examination of immune responses in pathogen-infected fish. This review also aims to highlight comparative studies of immune responses against bacteria, based on our previous findings in largemouth bass (Micropterus salmoides) against Nocardia seriolae, gray mullet (Mugil cephalus) against Lactococcus garvieae, orange-spotted grouper (Epinephelus coioides) against Vibrio harveyi, and koi carp (Cyprinus carpio) against Aeromonas sobria, using RNA-seq techniques. We demonstrated that only 39 differentially expressed genes (DEGs) were present in all species. However, the number of specific DEGs in each species was relatively higher than that of common DEGs; 493 DEGs in largemouth bass against N. seriolae, 819 DEGs in mullets against L. garvieae, 909 in groupers against V. harveyi, and 1471 in carps against A. sobria. The DEGs in different fish species were also representative of specific immune-related pathways. The results of this study will enhance our understanding of the immune responses of fish, and will aid in the development of effective vaccines, therapies, and disease-resistant strains.

Liver transcriptome analysis of the Sparus macrocephlus in response to Vibrio parahaemolyticus infection

The black seabream (Sparus macrocephlus) is an economically pivotal aquaculture species cultured in China and Southeast Asian countries. To understand the molecular immune mechanisms underlying the response to Vibrio parahaemolyticus, a comparative gene transcription analysis were performed with utilized fresh livers of V. parahaemolyticus-immunized Sparus macrocephlus with a control group through RNA-Seq technology. A total of 256663 contigs were obtained after excluded the low-quality sequences and assembly. The average length of contigs collected from this research is 1066.93 bp. Furthermore, blast analysis indicates 30747 contigs were annotated based on homology with matches in the NT, NR, gene, and string databases. A gene ontology analysis was employed to classify 21598 genes according to three major functional categories: molecular function, cellular component, and biological process. A total of 14470 genes were discovered in 303 KEGG pathways. RSEM and EdgeR were introduced to estimate 3841 genes significantly different expressed (False Discovery Rate<0.001) which includes 4072 up-regulated genes and 3771 down-regulated genes. A significant enrichment analysis of these differentially expressed genes and isogenes were conducted to reveal the major immune-related pathways which refer to the toll-like receptor, complement, coagulation cascades, and chemokine signaling pathways. In addition, 92175 potential simple sequence repeats (SSRs) and 121912 candidate single nucleotide polymorphisms (SNPs) were detected and identified sequencely in the Sparus macrocephlus liver transcriptome. This research characterized a gene expression pattern for normal and the V. parahaemolyticus -immunized Sparus macrocephlus for the first time and not only sheds new light on the molecular mechanisms underlying the host-V. parahaemolyticus interaction but contribute to facilitate future studies on Sparus macrocephlus gene expression and functional genomics.

Transcriptome analysis of the critically endangered Dabry's sturgeon (Acipenser dabryanus) head kidney response to Aeromonas hydrophila

Dabry's sturgeon (Acipenser dabryanus), as a living fossil, is considered a critically endangered aquatic animal in China. To date, the immune system of this species remains largely unknown, with limited available sequence information. In addition, increasing incidence of bacterial pathogenic diseases has been reported. Hence, the present study aimed to characterize comprehensively transcriptome profile of the head kidney from Dabry's sturgeon infected with Aeromonas hydrophila using Illumina platform. Over 42 million high-quality reads were obtained and de novo assembled into a final set of 195240 unique transcript fragments (unigenes), with an average length of 564 bp. Approximately 41702 unigenes were annotated in the NR NCBI database. Dabry's sturgeon unigenes had the highest number of hits with 14365 (34.45%) to Lepisosteus oculatus. The 195240 unigenes were assigned to three Gene Ontology (GO) categories: biological process, cellular component, and molecular function. Among them, 27770 unigenes were clustered into 26 Eukaryotic Orthologous Group (KOG) functional categories, and 36031 unigenes were mapped to 335 known Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. After A. hydrophila administration, 1728 differentially expressed unigenes were identified, including 980 upregulated and 748 downregulated unigenes. Further KEGG enrichment analysis of these unigenes identified 16 immune-related pathways, including the Toll-like receptor signaling pathway, chemokine signaling pathway, complement and coagulation pathway, RIG-I-like receptor signaling pathway, and NOD-like receptor signaling pathway. 20 DEGs were selected and their expression patterns are largely consistent with the transcriptome profile analysis, which clearly validated the reliability of the DEGs in transcriptome analysis. This work revealed novel gene expression patterns of Dabry's sturgeon host defense and contributes to a better understanding of the immune system and defense mechanisms of Dabry's sturgeon in response to bacterial infection. The results provide valuable references for studies in sturgeons that lack complete genomic sequences, and could also be helpful for the analyzing evolution among cartilaginous and teleost fish.