Assembly of the poorly differentiated Verasper variegatus W chromosome by different sequencing technologies.

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

The evolution of NLRC3 subfamily genes in Sebastidae teleost fishes.

BACKGROUND: With more than 36,000 valid fish species, teleost fishes constitute the most species-rich vertebrate clade and exhibit extensive genetic and phenotypic variation, including diverse immune defense strategies. NLRC3 subfamily genes, which are specific to fishes, play vital roles in the immune system of teleosts. The evolution of teleosts has been impacted by several whole-genome duplication (WGD) events, which might be a key reason for the expansions of the NLRC3 subfamily, but detailed knowledge of NLRC3 subfamily evolution in the family Sebastidae is still limited. RESULTS: Phylogenetic inference of NLRC3 subfamily protein sequences were conducted to evaluate the orthology of NLRC3 subfamily genes in black rockfish (Sebastes schlegilii), 13 other fish species from the families Sebastidae, Serranidae, Gasterosteidae and Cyclopteridae, and three species of high vertebrates (bird, reptile and amphibian). WGD analyses were used to estimate expansions and contractions of the NLRC3 subfamily, and patterns of expression of NLRC3 subfamily genes in black rockfish following bacterial infections were used to investigate the functional roles of these genes in the traditional and mucosal immune system of the Sebastidae. Different patterns of gene expansions and contractions were observed in 17 fish and other species examined, and one and two whole-genome duplication events were observed in two members of family Sebastidae (black rockfish and honeycomb rockfish, Sebastes umbrosus), respectively. Subsequently, 179 copy numbers of NLRC3 genes were found in black rockfish and 166 in honeycomb rockfish. Phylogenetic analyses corroborated the conservation and evolution of NLRC3 orthologues between Sebastidae and other fish species. Finally, differential expression analyses provided evidence of the immune roles of NLRC3 genes in black rockfish during bacterial infections and gene ontology analysis also indicated other functional roles. CONCLUSIONS: We hypothesize that NLRC3 genes have evolved a variety of different functions, in addition to their role in the immune response, as a result of whole genome duplication events during teleost diversification. Importantly, this study had underscored the importance of sampling across taxonomic groups, to better understand the evolutionary patterns of the innate immunity system on which complex immunological novelties arose. Moreover, the results in this study could extend current knowledge of the plasticity of the immune system.

Genome-wide identification, evolution and expression analysis of tight junction gene family and the immune roles of claudin5 gene in turbot (Scophthalmus maximus L.).

Tight junction proteins (TJs) are important component proteins that maintaining the structure and function of TJs, connecting to each other to form a TJ complex between cells, maintaining the biological homeostasis of the internal environment. In this study, a total of 103 TJ genes were identified in turbot according to our whole-transcriptome database. Transmembrane TJs were divided into seven subfamilies, including claudin (CLDN), occludin (OCLD), tricellulin (MARVELD2), MARVEL domain containing 3 (MARVELD3), junctional adhesion molecules (JAM), immunoglobulin superfamily member 5 (IGSF5/JAM4), blood vessel epicardial substance (BVEs). Moreover, the majority of homologous pairs of TJ genes showed highly conserved alongside length, exon/intron number and motifs. As for phylogenetic analysis for 103 TJ genes, eight of them have undergone a positive selection and JAMB-like has undergone the most neutral evolution. The expression patterns of several TJ genes showed the lowest expression levels in blood, while the highest expression levels were detected in intestine, gill and skin, which all belong to mucosal tissues. Meanwhile, most examined TJ genes showed down-regulated expression patterns during bacterial infection, while several TJ genes exhibited up-regulated expression patterns at a later stage (24 h). At the same time, several potential candidate genes (such as CLDN-15, CLDN-3, CLDN-12, CLDN-5 and OCLD) were significantly down-regulated, which may indicate their important functions that involved in the regulation of bacterial infection. Currently, there is little research on CLDN5 in the intestine, but it is highly expressed in the intestine and has significant changes in intestinal expression after bacterial infection. Thus, we knocked down CLDN5 by the method of lentiviral infection. The result showed CLDN5 was related to cell migration (wound healing) and apoptosis, and the method of dualluciferasereporterassay showed that the functions of CLDN5 could be regulated by miR-24. The study of TJs may lead to a better understanding of the function of TJs in teleost.

Identification and functional characterization of caspases in turbot (Scophthalmus maximus) in response to bacterial infection.

Apoptosis is the autonomous and orderly death of cells under genetic control to maintain the stability of the internal environment, and is a programmed cell death process with unique morphological and biochemical properties that is regulated by a variety of factors. Caspase gene family has a significant function in the process of apoptosis. However, the knowledge of caspases in turbot remains largely unknown. In present study, a total of nine turbot caspase genes were identified. The mRNA length of these caspase genes was ranged from 1225 bp (caspase-7) to 3216 bp (caspase-2), and the protein length was ranged from 281 aa (caspase-3a) to 507 aa (caspase-10). Phylogenetic analysis showed these caspase genes were divided into three subfamilies. The qRT-PCR results showed that turbot caspase genes were expressed in all the examined organs, especially the intestine, kidney, blood and gills. Meanwhile, we explored the expression patterns of caspase genes in the intestine, skin and gills after Vibrio anguillarum and Aeromonas salmonids infections. The results showed that caspase genes showed different expression patterns in mucosal tissues after bacterial infection, demonstrating the critical role of caspase genes in mucosal immune responses. In addition, protein-protein interaction analysis showed that caspase proteins interacted with immune molecules such as NLR, IL-1ß, and birc. The results of interference and overexpression experiments showed that caspase-1 might play key roles in the regulation of the IL-1ß production, but the detailed mechanism needs to be further studied. The results of this study provide valuable information for further study the roles of caspase genes in turbot, which could help us to further understand the inflammatory pathways in teleost.

Molecular characterization, expression analysis and function identification of TNFα in black rockfish (Sebastes schlegelii).

TNFα, as a pro-inflammatory cytokine, plays an important role in inflammation and immune homeostasis maintaining. However, the knowledge about the immune functions of teleost TNFα against bacterial infections is still limited. In this study, the TNFα was characterized from black rockfish (Sebastes schlegelii). The bioinformatics analyses showed the evolutionary conservations in sequence and structure. The expression levels of Ss_TNFα mRNA were significantly up-regulated in the spleen and intestine after Aeromonas salmonicides and Edwardsiella tarda infections, and dramatically down-regulated in PBLs after LPS and poly I:C stimulations. Meanwhile, the extremely up-regulated expressions of other inflammatory cytokines (especially for IL-1ß and IL17C) were observed in the intestine and spleen after bacterial infection and down-regulations were obtained in PBLs. The significant regulation with expression patterns of Ss_TNFα and other inflammatory cytokine mRNAs illustrated the variations of immunity in different tissues and cells of black rockfish. The regulated functions of Ss_TNFα in the up/downstream signaling pathways were preliminarily verified on the transcription and translation levels. Subsequently, in vitro knockdown of Ss_TNFα in the intestine cells of black rockfish confirmed the important immune roles of Ss_TNFα. Finally, the apoptotic analyses were conducted in PBLs and intestine cells of black rockfish. The rapid increases of the apoptotic rates were obtained in both PBLs and intestine cells after treatment with rSs_TNFα, but distinct apoptotic rates at the early and late stages of apoptosis were observed between these two types of cells. The results of apoptotic analyses suggested that Ss_TNFα could trigger apoptosis of different cells in different strategies in black rockfish. Overall, the findings in this study indicated the important roles of Ss_TNFα in the immune system of black rockfish during pathogenic infection, as well as the potential function on biomarker for monitoring the health status.

Genome-wide identification and functional characterization of inhibitor of nuclear factor-κB (IκB) kinase (IKK) in turbot (Scophthalmus maximus).

The inhibitor of nuclear factor-κB (IκB) kinase (IKK) is involved in a variety of intracellular cell signaling pathways and is an important component of the NF-κB signaling pathway. IKK genes have been suggested to play important roles in the innate immune response to pathogen infection in both vertebrates and invertebrates. However, little information is available about IKK genes in turbot (Scophthalmus maximus). In this study, six IKK genes were identified including SmIKKα, SmIKKα2, SmIKKß, SmIKKε, SmIKKγ, and SmTBK1. The IKK genes of turbot showed the highest identity and similarity with Cynoglossus semilaevis. Then, phylogenetic analysis showed that the IKK genes of turbot were most closely related to C. semilaevis. In addition, IKK genes were widely expressed in all the examined tissues. Meanwhile, the expression patterns of IKK genes were investigated by QRT-PCR after Vibrio anguillarum and Aeromonas salmonicida infection. The results showed that IKK genes had varying expression patterns in mucosal tissues after bacteria infection, indicating that they may play key roles in maintaining the integrity of the mucosal barrier. Subsequently, protein and protein interaction (PPI) network analysis showed that most proteins interacting with IKK genes were located in the NF-κB signaling pathway. Finally, the double luciferase report and overexpression experiments showed that SmIKKα/SmIKKα2/SmIKKß involved in the activation of NF-κB in turbot. In summary, our results suggested that IKK genes of turbot played important roles in the innate immune response of teleost, and provide valuable information for further study of the function of IKK genes.

Systematic analysis of circRNA-related ceRNA networks of black rockfish (Sebastes schlegelii) in response to Aeromonas salmonicides infection.

Aeromonas salmonicides is a type of Gram-negative bacteria and has become the main fish pathogen in aquaculture because of its characteristics of worldwide distribution, broad host range and potentially devastating impacts. In the past years, studies have been focused to explore the regulatory roles of circRNA-miRNA-mRNA network in fish diseases. However, there are only few systematic studies linked to the anti-bacterial roles of circRNA-related ceRNA networks in the spleen immune system of black rockfish (Sebastes schlegelii). In this study, the whole-transcriptome sequencing (RNA-seq) was conducted in the black rockfish spleen with A. salmonicida challenging. The differentially expressed (DE) circRNAs were identified comprehensively for the following enrichment analysis. Interactions of miRNA-circRNA pairs and miRNA-mRNA pairs were predicted for the construction of circRNA-related ceRNA regulatory networks. Then, protein-protein interaction (PPI) analysis of mRNAs from these ceRNA networks were conducted. Finally, a total number of 39 circRNAs exhibited significantly differential expressions during A. salmonicida infection in the black rockfish spleen in 4338 identified circRNAs from 12 samples in 4 libraries. Functional enrichment analysis suggested that they were significantly enriched in several immune-related pathways, including Endocytosis, FoxO signaling pathway, Jak-STST signaling pathway, Herpes simplex infection, etc. Subsequently, 290 circRNA-miRNA-mRNA pathways (91 at 2 h, 142 at 12 h and 65 at 24 h) were constructed including 31 circRNAs, 50 miRNAs, and 156 mRNAs. In conclusion, the circRNA-related ceRNA networks were established, which will provide some novel insights in molecular mechanistic investigations of anti-bacterial immune response in teleost. Also, these findings will propose significant predictive values for the development of methods of treatment and prevention in black rockfish after bacterial infection in the future.

The immune-related circRNA-miRNA-mRNA ceRNA regulatory network in the liver of turbot (Scophthalmus maximus L.) induced by Vibrio anguillarum.

Recently, Vibrio anguillarum, a Gram-negative pathogenic bacterium, has been becoming a major constraint on the development of the turbot aquaculture industry because of its characteristics of worldwide distribution, broad host range and potentially devastating impacts. Although the functions of protein-coding mRNAs in the immune response against bacterial infection have been reported, as well as several non-coding RNAs (ncRNAs), such as circular RNAs (circRNAs) and microRNAs (miRNAs), the relationships between mRNAs and ncRNAs in the immune system of turbot liver are still limited during bacterial infection. In present study, the comprehensive analyses of whole-transcriptome sequencing were conducted in turbot liver infected by V. anguillarum. The differential expression was analyzed in the data of circRNAs, miRNAs, and mRNAs. The interactions of miRNA-circRNA pairs and miRNA-mRNA pairs were predicted basing on the negative regulatory relationships between miRNAs and their target circRNAs\mRNAs. The circRNA-related ceRNA regulatory networks were constructed for the analyses of regulated mechanism in turbot immune system. Subsequently, the RT-qPCR was carried out to verify the results of sequencing. Finally, we identified 31 circRNAs, 53 miRNAs and 948 mRNAs with differential expression. Gene set enrichment analyses using Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways showed that innate immunity was principally activated at the early stages of infection, while adaptive immunity was activated after 24 h. Finally, 65 circRNA-miRNA-mRNA pathways were constructed, based on the hypothesis of ceRNA regulatory networks. In conclusion, our findings provide new insights on the underlying immune response to bacterial infection and identify novel target genes for the prevention and control of disease in turbot.

Systematic identification and characterization of lncRNAs and lncRNA-miRNA-mRNA networks in the liver of turbot (Scophthalmus maximus L.) induced with Vibrio anguillarum.

Long noncoding RNAs (lncRNAs), can regulate mRNA by targeting miRNA in a competing endogenous RNA network, have become a hot topic in the research of fish immune mechanism recent years. While in turbot (Scophthalmus maximus L.), an economically important marine fish, there are limited researches about the role of lncRNAs in its immune response to bacterial infection. In this study, a total of 184 differentially expressed lncRNAs (DElncRNAs) were systematically identified and characterized using whole-transcriptome sequencing of the liver of turbot challenged with Vibrioanguillarum at 0 h (control) and three different time points post infection (2 h, 12 h and 24 h, respectively). Subsequently, GO and KEGG signaling pathways of differentially expressed lncRNAs were analyzed to predict their function. We found that lncRNAs in our results were significantly enriched in several immune-related signaling pathways, including the NOD-like receptor signaling pathway, Toll-like receptor signaling pathway, Cytokine-cytokine receptor, MAPK signaling pathway, phagosome, PPAR signaling pathway and the regulation of autophagy. In addition, a total of 492 DE lncRNA - DE miRNA -DE mRNA networks were identified at three different time points post infection, which were consisted of 102 networks at 2 h, 122 networks at 12 h and 81 networks at 24 h post infection, respectively. Noticeably, 92 of these regulated networks were immune-related. These observations suggested that lncRNAs can regulate the expression of immune-related genes in the response to bacterial infection in turbot. Moreover, our findings would provide a new insight into the immune response of turbot to pathogen infection and lay a foundation for future study.

Transcriptomic characterization of Atlantic salmon (Salmo salar) head kidney following administration of Aeromonas salmonicida subsp. masoucida vaccine.

Atlantic salmon is one of the most famous and economically important fish species globally. However, bacterial diseases constantly constrain salmon aquaculture. Thereinto, Aeromonas salmonicida subsp. masoucida (ASM), classified as atypical A. salmonicida, caused huge losses to salmonid industry in China. In this regard, we conducted transcriptome analysis in Atlantic salmon head kidney following the administration of ASM vaccination to reveal genes, their expression patterns, and pathways involved in immune responses. A total of 448.71 million clean reads were obtained, and 397.69 million reads were mapped onto the Atlantic salmon reference genome. In addition, 117, 1891, 741, 207, and 377 genes were significantly up-regulated, and 183, 1920, 695, 83, and 539 genes were significantly down-regulated post ASM vaccination at 12 h, 24 h, 1 m, 2 m, and 3 m, respectively. Furthermore, KEGG pathway analysis revealed that many differentially expressed genes (DEGs) following ASM vaccination were involved in cell adhesion molecules (H2-Aa-l and CD28-l),cytokine-cytokine receptor interaction (IL10, CXCL9, CXCL11, CXCR3, and CCL19), herpes simplex infection (IL1B, SOCS3-l, and C3-l), HTLV-I infection (Il1r2 and BCL2L1), influenza A (CXCL8 and Il12b), and PI3K-Akt signaling pathway (PIK3R3-l and Ddit4-l). Finally, the results of qRT-PCR showed a significant correlation with RNA-Seq results, suggesting the reliability of RNA-Seq for gene expression analysis. This study sets the foundation for further study on the vaccine protective mechanism in Atlantic salmon as well as other teleost species.

Comparative Transcriptome Analysis of Spleen Reveals Potential Regulation of Genes and Immune Pathways Following Administration of Aeromonas salmonicida subsp. masoucida Vaccine in Atlantic Salmon (Salmo salar).

Aeromonas salmonicida is a global fish pathogen. Aeromonas salmonicida subsp. masoucida (ASM) is classified as atypical A. salmonicida and caused huge losses to salmonid industry in China. Hence, it is of great significance to develop ASM vaccine and explore its protection mechanism in salmonids. In this regard, we conducted RNA-seq analysis with spleen tissue of Atlantic salmon after ASM vaccination to reveal genes, their expression patterns, and pathways involved in immune protections. In our results, a total of 441.63 million clean reads were obtained, and 389.37 million reads were mapped onto the Atlantic salmon reference genome. In addition, 1125, 2126, 1098, 820, and 1351 genes were significantly up-regulated, and 747, 2626, 818, 254, and 908 genes were significantly down-regulated post-ASM vaccination at 12 h, 24 h, 1 month, 2 months, and 3 months, respectively. Subsequent pathway analysis revealed that many differentially expressed genes (DEGs) following ASM vaccination were involved in cytokine-cytokine receptor interaction (TNFRSF11b, IL-17RA, CCR9, and CXCL11), HTLV-I infection (MR1 and HTLV-1), MAPK signaling pathway (MAPK, IL8, and TNF-α-1), PI3K-Akt signaling pathway (PIK3R3, THBS4, and COL2A1), and TNF signaling pathway (PTGS2, TNFRSF21-l, and CXCL10). Finally, the results of qRT-PCR showed a significant correlation with RNA-seq results, suggesting the reliability of RNA-seq for gene expression analysis. This study provided insights into regulation of gene expression and their involved pathways in Atlantic salmon spleen in responses to vaccine, and set the foundation for further study on the vaccine protective mechanism in Atlantic salmon as well as other teleost species.

Identification of mRNA-miRNA-lncRNA regulatory network associated with the immune response to Aeromonas salmonicides infection in the black rockfish (Sebastes schlegelii).

In aquaculture, Aeromonas salmonicides (A. salmonicida) is a main fish pathogen because of its nearly worldwide distribution, and broad host range. Recently, an increasing number of evidences have uncovered the roles of mRNA-miRNA-lncRNA network in fish diseases. In current study, RNA-seq was conducted in the black rockfish spleen following A. salmonicida infection at 0 h (Sp0 or control) and three different post-infection time-points (2 h: Sp2, 12 h: Sp12 and 24 h: Sp24, respectively) to comprehensively identify differentially expressed (DE) mRNAs, miRNAs and lncRNAs. Enrichment analysis and protein-protein interaction (PPI) analysis of DE mRNAs were performed. Then, expression and correlation analysis for mRNAs and their upstream miRNAs and lncRNAs were conducted. Finally, a total of 1364 mRNAs, 17 miRNAs and 1584 lncRNAs exhibited significantly differential expressions during bacterial infection in the black rockfish spleen. Functional enrichment analysis suggested that they were significantly enriched in several immune-related pathways, including Amino sugar and nucleotide sugar metabolism, Cell adhesion molecules (CAMs), Neuroactive ligand-receptor interaction, Nicotinate and nicotinamide metabolism, Pentose and glucuronate interconversions, Phagosome, Proteasome, etc. Subsequently, 1091 lncRNA-miRNA-mRNA pathways (323 in Sp2, 609 in Sp12 and 207 in Sp24) were constructed including 400 lncRNAs, 69 miRNAs, and 70 mRNAs. Meanwhile, NLRC3/novel-264/LNC_00116154 pathway demonstrated important immune modulating function in the black rockfish against A. salmonicida infection. Finally, the novel mRNA-miRNA-lncRNA sub-networks were established, among which all mRNAs and ncRNAs possessed significant predictive values for further studies for immune responses in the black rockfish.

Characterization and the potential immune role of class A scavenger receptor member 4 (SCARA4) in bacterial infection in turbot (Scophthalmus maximus L.).

The class A scavenger receptors play important roles in innate immunity and are distributed on plasma membrane of macrophages and other cell types. Notably, the class A scavenger receptor 4 (SCARA4) contains a typical C-type (calcium-dependent) lectin domain, which belongs to the collectin family of pattern recognition receptors and is involved in the immune response against infection. Here, one turbot SCARA4 gene was identified with a 2,292 bp open reading frame (ORF) encoding 763 amino acid residues. Multiple sequence analysis and phylogenetic analysis confirmed that SmSCARA4 gene was more close to that of P. olivaceus. Gene structure and syntenic analysis showed conserved exon/intron organization pattern and syntenic pattern across selected vertebrate species. Tissue distribution analysis showed SmSCARA4 was expressed in all the tested healthy tissues with the relative high expression levels in skin, gill and spleen. Following both E. tarda and V. anguillarum challenge in vivo, SmSCARA4 was significantly repressed in gill and intestine. Remarkably, SmSCARA4 showed the strongest binding ability to LPS and strongest upregulation in turbot head kidney macrophages in response to LPS. Knockdown and overexpression of SmSCARA4 revealed its interactions with the two pro-inflammatory cytokines, TNF-α and IL-1ß. Finally, repression of SmSCARA4 via combined treatment of LPS and overexpression of SmSCARA4 construct in turbot head kidney macrophages further indicated an inhibitory role of SmSCARA4 in LPS-stimulated inflammation. Taken together, turbot SmSCARA4 plays an important role in turbot immunity, especially in the mucosa-related systems; SmSCARA4 possesses strong binding specificity to LPS, and exerts protective roles in response to LPS infection by reducing the release of pro-inflammatory cytokines. The mechanisms of inhibitory role of SmSCARA4 in LPS-elicited inflammation await further investigation.

Genome-wide characterization of gap junction (connexins and pannexins) genes in turbot (Scophthalmus maximus L.): evolution and immune response following Vibrio anguillarum infection.

Gap junction (GJ), a special intercellular junction between different cell types, directly connects the cytoplasm of adjacent cells, allows various molecules, ions and electrical impulses to pass through the intercellular regulatory gate, and plays vital roles in response to bacterial infection. Up to date, the information about the GJ in turbot (Scophthalmus maximus L.) is still limited. In current study, 43 gap junction genes were identified in turbot, phylogeny analysis suggested that gap junctions from turbot and other species were clustered into six groups, GJA, GJB, GJC, GJD, GJE and PANX, and turbot GJs together with respective GJs from Japanese flounder, half-smooth tongue sole and large yellow croaker, sharing same ancestors. In addition, these 43 GJ genes distributed in different chromosomes unevenly. According to gene structure and domain analysis, these genes (in GJA-GJE group) were highly conserved in that most of them contain the transmembrane area, connexin domain (CNX) and cysteine-rich domain (connexin CCC), while PANXs contain Pfam Innexin. Although only one tandem duplication was identified in turbot gap junction gene, 235 pairs of segmental duplications were identified in the turbot genome. To further investigate their evolutionary relationships, Ka/Ks was calculated, and results showed that most ratios were lower than 1, indicating they had undergone negative selection. Finally, expression analysis showed that gap junction genes were widely distributed in turbot tissues and significantly regulated after Vibrio anguillarum infection. Taken together, our research could provide valuable information for further exploration of the function of gap junction genes in teleost.

CXC chemokines and their receptors in black rockfish (Sebastes schlegelii): Characterization, evolution analyses, and expression pattern after Aeromonas salmonicida infection.

Chemokines are crucial regulators of cell mobilization for development, homeostasis, and immunity. Chemokines signal through binding to chemokine receptors, a superfamily of seven-transmembrane domain G-coupled receptors. In the present study, seventeen CXC chemokine ligands (SsCXCLs) and nine CXC chemokine receptors (SsCXCRs) were systematically identified from Sebastes schlegelii genome. Phylogeny, synteny, and evolutionary analyses were performed to annotate these genes, indicating that the tandem duplications (CXCL8, CXCL11, CXCL32, CXCR2, and CXCR3), the whole genome duplications (CXCL8, CXCL12, CXCL18, and CXCR4), and the teleost-specific members (CXCL18, CXCL19, and CXCL32) led to the expansion of SsCXCLs and SsCXCRs. In addition, SsCXCLs and SsCXCRs were ubiquitously expressed in nine examined healthy tissues, with high expression levels observed in head kidney, liver, gill and spleen. Moreover, most SsCXCLs and SsCXCRs were significantly differentially expressed in head kidney, liver, and gill after Aeromonas salmonicida infection, and exhibited tissue-specific and time-dependent manner. Finally, protein-protein interaction network (PPI) analysis indicated that SsCXCLs and SsCXCRs interacted with a few immune-related genes such as interleukins, cathepsins, CD genes, and TLRs, etc. These results should be valuable for comparative immunological studies and provide insights for further functional characterization of chemokines and receptors in teleost.

Comparative analysis of the miRNA-mRNA regulation networks in turbot (Scophthalmus maximus L.) following Vibrio anguillarum infection.

MicroRNAs could not only regulate posttranscriptional silencing of target genes in eukaryotic organisms, but also have positive effect on their target genes as well. These microRNAs have been reported to be involved in mucosal immune responses to pathogen infection in teleost. Therefore, we constructed the immune-related miRNA-mRNA networks in turbot intestine following Vibrio anguillarum infection. In our results, 1550 differentially expressed (DE) genes and 167 DE miRNAs were identified. 113 DE miRNAs targeting 89 DE mRNAs related to immune response were used to construct miRNA-mRNA interaction networks. Functional analysis showed that target genes were associated with synthesis and degradation of ketone bodies, mucin type O-Glycan biosynthesis, homologous recombination, biotin metabolism, and intestinal immune network for IgA production that were equivalent to the function of IgT and IgM in fish intestine. Finally, 10 DE miRNAs and 7 DE mRNAs were selected for validating the accuracy of high-throughput sequencing results by qRT-PCR. The results of this study will provide valuable information for the elucidation of the regulation mechanisms of miRNA-mRNA interactions involved in disease resistance in teleost mucosal immune system.

Galectins in turbot (Scophthalmus maximus L.): Characterization and expression profiling in mucosal tissues.

Galectins, a family of evolutionary conserved ß-galactoside-binding proteins, have been characterized in a wide range of species. Many reports have indicated vital roles of galectins in innate immunity, especially in the mucosal tissues against infection. However, the systematic identification of galectin gene family is still lacking in teleost. Here, we characterized the galectin gene family and investigated their expression profiles post bacterial challenge in turbot (Scophthalmus maximus L.). In this study, a total of 13 galectin genes were characterized in turbot, phylogenetic analyses revealed their strong relationships to half smooth tongue sole and puffer fish, and syntenic analyses confirmed the orthology suggested by the phylogenetic analysis. In addition, the copy number of galectin genes is similar across a broad spectrum of species from fish to amphibians, birds, and mammals, ranging from 8 to 16 genes. Furthermore, the galectin genes were widely expressed in all the examined turbot tissues, and most of the galectin genes were strongly expressed in mucosal tissues (skin, gill and intestine). Moreover, majority of the galectin genes were significantly regulated after Vibrio anguillarum infection in the intestine, gill and skin, suggesting that galectins were involved in the mucosal immune response to V. anguillarum infection in turbot. In addition, subcellular localization analysis showed lgals3a was distributed in the cytoplasm and nucleus. However, the knowledge of galectins are still limited in teleost species, further studies should be carried out to better characterize its detailed roles in teleost mucosal immunity.

Characterization and expression analysis of mitochondrial localization molecule: NOD-like receptor X1 (Nlrx1) in mucosal tissues of turbot (Scophthalmus maximus) following bacterial challenge.

The NOD-like receptor X1 (NLRX1) is a member of highly conserved nucleotide-binding domain (NBD)- and leucine-rich-repeat (LRR)-containing family (known as NLR), that localizes to the mitochondrial outer membrane and regulate the innate immunity by interacting with mitochondrial antiviral-signaling protein (MAVS). As one of cytoplasmic PRRs, NLRX1 plays key roles for pathogen recognition, autophagy and regulating of subsequent immune signaling pathways. In this study, we identified the nlrx1 in turbot as well as its expression profiles in mucosal surfaces following bacterial infection. In our results, the full-length nlrx1 transcript consists of an open reading frame (ORF) of 4,886 bp encoding the putative peptide of 966 amino acids. The phylogenetic analysis revealed the SmNlrx1 showed the closest relationship to Cynoglossus semilaevis. In addition, the Nlrx1 mRNA expression could be detected in all the examined tissues, with the most abundant expression level in head kidney, and the lowest expression level in liver. Moreover, Nlrx1 showed similar expression patterns following Vibrio anguillarum and Streptococcus iniae infection, that were both significantly up-regulated following challenge, especially post S. iniae challenge. Finally, fluorescence microscopy unveiled that the SmNlrx1 localized to mitochondria in HEK293T by N-terminal mitochondrial targeting sequence. Characterization of Nlrx1 might have an important implication in bioenergetic adaptation during metabolic stress, oncogenic transformation and innate immunity and will probably contribute to the development of novel intervention strategies for farming turbot.

Dynamics of MiRNA Transcriptome in Turbot (Scophthalmus maximus L.) Intestine Following Vibrio anguillarum Infection.

MicroRNAs (miRNAs) are a group of small non-coding RNAs, which could bind to the 3'-untranslated regions of their target mRNAs to regulate gene expression in various biological processes, including immune-regulated signaling pathways. Turbot (Scophthalmus maximus L.), an important commercial fish species in China, has been suffering with Vibrio anguillarum infection resulted in dramatic economic loss. Therefore, we investigated the expression profiles of miRNAs, as well as the immune-related miRNA-mRNA pairs in turbot intestine at 1 h, 4 h, and 12 h following V. anguillarum infection. As a result, 266 predicted novel miRNAs and 283 conserved miRNAs belonging to 92 miRNA families were detected. A total of 44 miRNAs were differentially expressed in the intestine following V. anguillarum infection. Following prediction, the potential target genes of differentially expressed miRNAs were grouped into a wide range of functional categories, including immune defense/evasion, inflammatory responses, RIG-I signaling pathway, and Toll-like receptor signaling pathway. Moreover, we selected 15 differentially expressed immune genes and their related differentially expressed miRNAs to construct an interaction network for V. anguillarum infection in turbot. These results suggested that in teleost, as in higher vertebrates, miRNAs prominently contribute to immune responses, protecting the host against infection. In addition, this is the first report of comprehensive identification of turbot miRNAs being differentially regulated in the intestine related to V. anguillarum infection. Our results provided an opportunity for further understanding of the molecular mechanisms of miRNA regulation in turbot host-pathogen interactions.

Identification and expression profiling analysis of microRNAs in Nile tilapia (Oreochromis niloticus) in response to Streptococcus agalactiae infection.

MicroRNAs (miRNAs) play vital regulatory roles in various biological processes, including in immune responses. Nile tilapia (Oreochromis niloticus) is an important commercial fish species in China. To identify immune-related miRNAs of O. niloticus, 4 libraries from liver during S. agalactiae infection (0 h, 5 h, 50 h, and 7 d) were sequenced by high-throughput sequencing technology in tilapia. We obtained 10,703,531, 11,507,163, 11,180,179 and 13,408,414 clean reads per library, respectively. In our results, a total of 482 miRNAs were identified through bioinformatic analysis, including 220 conserved miRNAs and 262 putative novel miRNAs. Moreover, 21 (4.36%), 50 (10.37%), and 46 (9.54%) miRNAs were significantly differentially expressed at 5 h, 50 h and 7 d, respectively. In addition, 6939 target genes regulated by these differentially expressed miRNAs were predicted, and their functional annotations were predicted by Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, which revealed that a majority of differentially expressed miRNAs were involved in apoptotic process, metabolic process, and immune responses. Finally, Real-time quantitative PCR experiments were performed for 7 miRNAs by stem-loop RT-PCR, and a general agreement was confirmed between the sequencing and RT-qPCR data. To our understanding, this is the first report of comprehensive identification of O. niloticus miRNAs being differentially regulated in liver related to S. agalactiae infection. This work provides an opportunity for further understanding of the molecular mechanisms of miRNA regulation in O. niloticus host-pathogen interactions, and genetic resources for molecular assistant selection for disease resistant breeding program.