A CCL25 chemokine functions as a chemoattractant and an immunomodulator in black rockfish, Sebastes schlegelii

Hypoxia stress induces hepatic antioxidant activity and apoptosis, but stimulates immune response and immune-related gene expression in black rockfish Sebastes schlegelii

Dissolved oxygen concentrations both in the open ocean and coast have been declining due to the interaction of global climate change and human activity. Fish have evolved different adaptative strategies to cope with possibly damage induced by hypoxic environments. Black rockfish as important economic fish widely reared in the offshore sea cage, whereas related physiological response subject to hypoxia stress remained unclear. In this study, hepatic anti-oxidant enzymes (superoxide dismutase [SOD], catalase [CAT], glutathione peroxidase [GSH-Px]), aminotransferase (AST) and alanine aminotransferase (ALT) activities, lipid peroxidation (LPO), malondialdehyde (MDA) and glutathione (GSH) content, immunological parameters and the expression of apoptosis (bax, bcl2, p53, caspase3, xiap) and immune-related genes (c3, il-1β, ccl25, saa, hap, isg15) of black rockfish were determined during hypoxia and reoxygenation to illustrate the underlying defense response mechanisms. Results showed that hypoxia stress remarkably increased hepatic LPO and MDA content, AST and ALT activity and proportion of pyknotic nucleus. Hepatic SOD, CAT and GSH-Px activity manifested similar results, whereas GSH levels significantly decreased under hypoxia stress. The apoptosis rate of hepatocyte increased during hypoxia stress and reoxygenation. Meanwhile, p53, caspase3, bax and xiap mRNAs and bax/bcl2 rations were significantly up-regulated under hypoxia stress. However, bcl2 mRNA was significantly down-regulated. Interestingly, hypoxia stress significantly increased NBT-positive cell percent, phagocytic index, respiratory burst and ACH50 activity, and lysozyme activity. The mRNA levels of c3, ilβ, ccl25, saa, hap and isg15 were significantly up-regulated in the liver, spleen and head-kidney under hypoxia stress. The above parameters recovered to normal status after reoxygenation for 24 h Thus, hypoxia stress impairs hepatic antioxidant capacity, induces oxidative damage and apoptosis via the xiap-p53-bax-bcl2 and the caspase-dependent pathways, but enhances host immunity by regulating nonspecific immune indices and related genes expression to maintain homeostasis in black rockfish. These findings will help fully understand the hypoxia tolerance mechanisms of black rockfish and provide more data for offshore open ocean farming.

Genome-wide identification, evolutionary analysis, and antimicrobial activity prediction of CC chemokines in allotetraploid common carp, Cyprinus carpio

Chemokines are a group of secreted small molecules which are essential for cell migration in physiological and pathological conditions by binding to specific chemokine receptors. They are structurally classified into five groups, namely CXC, CC, CX3C, XC and CX. CC chemokine group is the largest one among them. In this study, we identified and characterized 61 CC chemokines from allotetraploid common carp (Cyprinus carpio). The sequence analyses showed that the majority of CC chemokines had an N-terminal signal peptide, and an SCY domain, and all CC chemokines were located in the extracellular region. Phylogenetic, evolutionary and syntenic analyses confirmed that CC chemokines were annotated as 11 different types (CCL19, CCL20, CCL25, CCL27, CCL32, CCL33, CCL34, CCL35, CCL36, CCL39, and CCL44), which exhibited unique gene arrangement pattern and chromosomal location respectively. Furthermore, genome synteny analyses between common carp and four representative teleost species indicated expansion of common carp CC chemokines resulted from the whole genome duplication (WGD) event. Additionally, the continuous evolution of gene CCL25s in teleost afforded a novel viewpoint to explain the WGD event in teleost. Then, we predicted the three-dimensional structures and probable function regions of common carp CC chemokines. All the CC chemokines core structures were constituted of an N-loop, a three-stranded β-sheet, and a C-terminal helix. Finally, 43 CC chemokines were predicted to have probable general antimicrobial activity. Their tertiary structures, cationic and amphiphilic physicochemical property supported the viewpoint. To verify the prediction, six recombinant CCL19s proteins were prepared and the antibacterial activity against Escherichia coli and Aeromonas hydrophila were verified. The results supported our prediction that rCCL19a.1s (rCCL19a.1_a, rCCL19a.1_b) and rCCL19bs (rCCL19b_a, rCCL19b_b), especially rCCL19bs, exhibited extremely significant inhibition to the growth of both E. coli and A. hydrophila. On the contrary, two rCCL19a.2s had no significant inhibitory effect. These studies suggested that CC chemokines were essential in immune system evolution and not monofunctional during pathogen infection.

A novel full-length transcriptome resource from multiple immune-related tissues in turbot (Scophthalmus maximus) using Pacbio SMART sequencing

Turbot (Scophthalmus maximus) is an important cold-water economic fish. However, the production and development of turbot industry has been constantly hindered by the frequent occurrence of some diseases. Lacking full-length transcriptome for turbot limits immune gene discoveries and gene structures analysis. Therefore, we generated a full-length transcriptome using mixed immune-related tissues of turbot with PacBio Sequel platform. In this study, a total of 31.7 Gb high quality data were generated with the average subreads length of 2618 bp. According to the presence of 5' and 3' primers as well as poly (A) tails, FL (Full-length) and NFL (Non-full-length) isoforms were obtained. Meanwhile, we identified 32,003 non-redundant transcripts, 76.02% of which was novel isoforms of known genes. In addition, 12,176 alternative splicing (AS) events, 6614 polyadenylation (APA) events, 1905 transcription factors, and 2703 lncRNAs were identified. This work is a comprehensive report on the full-length transcriptome of immune-related tissues of turbot, and it also provides valuable molecular resources for future research on the adaptation mechanisms and functional genomics of turbot.

DNA vaccine dual-expressing viral hemorrhagic septicemia virus glycoprotein and C-C motif chemokine ligand 19 induces the expression of immune-related genes in zebrafish (Danio rerio)

Glycoprotein (G protein)-based DNA vaccines are effective in protecting aquaculture fish from rhabdoviruses but the degree of immune response they elicit depends on plasmid concentration and antigen cassette. Here, we developed a DNA vaccine using the viral hemorrhagic septicemia virus G (VG) gene and chemokine (C-C motif) ligand 19 (CCL19)a.2 regulated by the CMV promoter as the molecular adjuvant. After transfection of the prepared plasmid (pVG + CCL19) into epithelioma papulosum cyprini cells, mRNA expression was confirmed through quantitative real-time polymerase chain reaction. The vaccine was intramuscularly injected into zebrafish (Danio rerio), and 28 days after immunization, viral hemorrhagic septicemia virus (10⁵ TCID₅₀/10 µl/fish) was intraperitoneally injected. A survival rate of 68% was observed in the pVG + CCL19 group but this was not significantly different from the survival rate of fish treated with pVG alone, that is, without the adjuvant. However, the expression of interferon- and cytokine-related genes in the spleen and kidney tissues of zebrafish was significantly increased (p < 0.05) on days 1, 3, 7, and 14 after immunization. Thus, CCL19a.2 induced an initial immune response as a molecular adjuvant, which may provide initial protection against virus infection before vaccination-induced antibody formation. This study provides insights on the functions of CCL19a.2 adjuvant in DNA vaccines.

CC and CXC chemokines in turbot (Scophthalmus maximus L.): Identification, evolutionary analyses, and expression profiling after Aeromonas salmonicida infection

Chemokines are a superfamily of structurally related cytokines, which exert essential roles in guiding cell migration in development, homeostasis, and immunity. CC and CXC chemokines are the two major subfamilies in teleost species. In this study, a total of seventeen CC and CXC chemokines, with inclusion of twelve CC and five CXC chemokines, were systematically identified from the turbot genome, making turbot the teleost harboring the least number of CC and CXC chemokines among all teleost species ever reported. Phylogeny, synteny, and genomic organization analyses were performed to annotate these genes, and multiple chemokine genes were identified in the turbot genome, due to the tandem duplications (CCL19 and CCL20), the whole genome duplications (CCL20, CCL25, and CXCL12), and the teleost-specific members (CCL34-36, CCL44, and CXCL18). In addition, chemokines were ubiquitously expressed in nine examined healthy tissues, with high expression levels observed in liver, gill, and spleen. Moreover, most chemokines were significantly differentially expressed in gill and spleen after Aeromonas salmonicida infection, and exhibited tissue-specific and time-dependent manner. Finally, protein-protein interaction network (PPI) analysis indicated that turbot chemokines interacted with a few immune-related genes such as interleukins, cathepsins, stats, and TLRs. These results should be valuable for comparative immunological studies and provide insights for further functional characterization of chemokines in teleost.

The dynamic immune responses of Mandarin fish (Siniperca chuatsi) to ISKNV in early infection based on full-length transcriptome analysis and weighted gene co-expression network analysis

Mandarin fish (Siniperca chuatsi) been seriously harmed by infectious spleen and kidney necrosis virus (ISKNV) in recent years, but the early immune response mechanism of infection is still unknown. Here, we performed RNA sequencing on the spleens of mandarin fish infected with ISKNV at 0, 12, 24, 48, and 72 h post-infection (hpi) using short-read Illumina RNA sequencing and long-read Pacific Biosciences isoform sequencing to generate a full-length transcriptome. The immune responses of mandarin fish infected with ISKNV at the molecular level were characterized by RNA-seq analysis and weighted gene co-expression network analysis (WGCNA). A total of 26,528 full-length transcript sequences were obtained. There were 2,729 (1,680 up-regulated and 1,112 down-regulated), 1,874 (1,136 up-regulated and 738 down-regulated), 2,032 (1,158 up-regulated and 847 down-regulated), and 4,176 (2,233 up-regulated and 1,943 down-regulated) differentially expressed genes (DEGs) in mandarin fish at 12, 24, 48, and 72 hpi, compared with uninfected fish, respectively. A total of four modules of co-expressed DEGs identified by WGCNA were significantly positively correlated to the four time points after infection, respectively. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis revealed that the immune-related DEGs in all these modules were mainly enriched in Phagosome, Endocytosis, Herpes simplex infection, and Cytokine-cytokine receptor interaction pathways. Further analysis showed that oher signaling pathways, including CAMs, NOD-like receptor and ER protein processing, Intestinal immune network for IgA production, TLR pathway, and Apoptosis significantly enriched in four modules corresponding to 12, 24, 48, and 72 hpi respectively, had specifically participated in the immune response. Hub genes identified based on the high-degree nodes in the WGCN, including CAM3, IL-8, CCL21, STING, SNX1, PFR and TBK1, and some DEGs such as MHCI, MHCII, TfR, STING, TNF α, TBK1, IRF1, and NF-kB, BCR, IgA and Bcl-XL had involved in dynamic molecular response of mandarin fish to ISKNV infection. In sum, this study provides a set of full-length transcriptome of the spleen tissue of mandarin fish for the first time and revealed a group of immune genes and pathways involved in different temporal responses to ISKNV infection, which has implications for resource conservation and aiding the development of strategies to prevent virus early infection for mandarin fish.

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.

Molecular characterization and chemotaxis assay of a CC motif chemokine ligand 25 from Japanese sea bass (Lateolabrax japonicus)

CC motif chemokine ligand 25 (CCL25) is a key chemokine that attracts various types of leukocytes, such as activated peritoneal macrophages. However, information on CCL25 in fish is limited. Here, a CCL25 gene (LjCCL25) was identified from Japanese sea bass (Lateolabrax japonicus), showing upregulation in multiple tissues against Vibrio harveyi infection. The recombinant LjCCL25 (rLjCCL25) only significantly induced the migration of monocytes/macrophages (MO/MΦ) both in vitro and in vivo, but didn't induce that of neutrophils or lymphocytes. Additionally, rLjCCL25 only induced migration of the lipopolysaccharide-stimulated MO/MΦ (M1 type). Knockdown of Japanese sea bass CC chemokine receptor 9 (LjCCR9) expression in MO/MФ by RNA interference inhibited the LjCCL25-induced chemotaxis of resting and M1 type MO/MФ. Moreover, administration of 300 ng/g rLjCCL25 effectively increased the survival of V. harveyi-infected fish and decreased bacterial load. Our study demonstrates that LjCCL25 functions as an MO/MФ chemoattractant via LjCCR9 in Japanese sea bass against V. harveyi.

Identification of Antimicrobial Peptide Genes in Black Rockfish Sebastes schlegelii and Their Responsive Mechanisms to Edwardsiella tarda Infection

The black rockfish, Sebastes schlegelii, is a typical viviparous teleost, which belongs to the family Scorpaenidae. Due to its high economic and ecological values, S. schlegelii has been widely cultured in East Asian countries. With the enlargement of cultivation scale, bacterial and viral diseases have become the main threats to the farming industry of S. schlegelii, which have resulted in significant economic losses. In this study, Illumina shotgun sequencing, single-molecule real-time (SMRT) sequencing, 10× genomics and high-throughput chromosome conformation capture (Hi-C) technologies were collectively applied to assemble the genome of S. schlegelii. Then, we identified the antimicrobial peptide genes (AMPs) in the S. schlegelii genome. In total, 214 AMPs were identified in the S. schlegelii genome, which can be divided into 33 classes according to the annotation and cataloging of the Antimicrobial Peptides Database (APD3). Among these AMPs, thrombin-derived C-terminal peptide (TCP) was the dominant type, followed by RegIIIgamma and chemokine. The amino acid sequences of the TCP, cgUbiquitin, RegIIIalpha, RegIIIgamma, chemokine shared 32.55%, 42.63%, 29.87%, 28.09%, and 32.15% similarities among the same type in S. schlegelii. Meanwhile, the expression patterns of these AMPs in nine healthy tissues and at different infection time points in intestine were investigated. The results showed that the numbers and types of AMPs that responded to Edwardsiella tarda infection gradually increased as the infection progressed. In addition, we analyzed the phylogenetic relationships of hepcidins in teleost. The identification of AMPs based on the whole genome could provide a comprehensive database of potential AMPs, and benefit for the understanding of the molecular mechanisms of immune responses to E. tarda infection in S. schlegelii. This would further offer insights into an accurate and effective design and development of AMP for aquaculture therapy in the future.

The CC and CXC chemokine receptors in turbot (Scophthalmus maximus L.) and their response to 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, eleven CC chemokine receptors (CCRs) and seven CXC chemokine receptors (CXCRs) were identified from turbot genome. Phylogenetic and syntenic analyses were performed to annotate these genes, indicating the closest relationship between the turbot chemokine receptors and their counterparts of Japanese flounders (Paralichthys olivaceus). Evolutionary analyses revealed that the tandem duplications of CCR8 and CXCR3, the whole genome duplications of CCR6, CCR9, CCR12, and CXCR4, and the teleost-specific CCR12 led to the expansion of turbot chemokine receptors. In addition, turbot chemokine receptors were ubiquitously expressed in nine examined healthy tissues, with high expression levels observed in spleen, gill, and head kidney. Moreover, most turbot chemokine receptors were significantly differentially expressed in spleen and gill after Aeromonas salmonicida infection, and exhibited general down-regulations at early time points and then gradually up-regulated. Finally, protein-protein interaction network (PPI) analyses indicated that chemokine receptors interacted with a few immune-related genes such as interleukins, Grk genes, CD genes, etc. These results should be valuable for comparative immunological studies and provide insights for further functional characterization of chemokine receptors in turbots.

Molecular cloning and characterization of chemokine C-C motif ligand 34 (CCL34) genes from olive flounder (Paralichthys olivaceus)

Chemokines are a superfamily of chemotactic cytokines that regulate the migration and immune responses of leukocytes. Depending on the arrangement of the first two cysteine residues, chemokines are divided into four groups: CXC (α), CC (β), C (γ), and CX3C (δ). Chemokine C-C motif ligand 34 (CCL34) is a member of the CC chemokine family and is known as a fish-specific CC chemokine. In this experiment, we analyzed the molecular cloning and characterization of the PoCCL34 gene in olive flounder (Paralichthys olivaceus), including CCL34a.3 (PoCCL34a.3) and CCL34b.3 (PoCCL34b.3). The amino acid sequence of PoCCL34 has four highly conserved cysteine residues and it has a C-C motif. Phylogenetic analysis revealed that PoCCL34 was phylogenetically clustered in the fish CCL34 subcluster. Recombinant PoCCL34 induced chemotaxis of head kidney leukocytes in a dose-dependent manner. Head kidney leukocytes stimulated with PoCCL34 also exhibited significant respiratory burst activity and increased expression of pro-inflammatory cytokines (IL-1β, IL-6, and CXCL8), but the overall expression of interferon-related genes (IFN-α/β, IFN-γ, Mx, and ISG15) did not increase. Olive flounder injected with recombinant PoCCL34 demonstrated increased expression of pro-inflammatory cytokines (IL-1β and IL-6) in the head kidney. However, there was no increase in the expression of interferon-related genes (IFN-α/β, IFN-γ, Mx, and ISG15). Additionally, recombinant PoCCL34 induced high lysozyme activity in the serum of the flounder. These results indicate that although PoCCL34 is not involved in the antiviral response, it may play a significant role in the overall immune response of the flounder, particularly in mediating the inflammatory response.

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.

Integrated Analysis of circRNA-miRNA-mRNA Regulatory Networks in the Intestine of Sebastes schlegelii Following Edwardsiella tarda Challenge

Sebastes schlegelii, an important aquaculture species, has been widely cultured in East Asian countries. With the increase in the cultivation scale, various diseases have become major threats to the industry. Evidence has shown that non-coding RNAs (ncRNAs) have remarkable functions in the interactions between pathogens and their hosts. However, little is known about the mechanisms of circular RNAs (circRNAs) and coding RNAs in the process of preventing pathogen infection in the intestine in teleosts. In this study, we aimed to uncover the global landscape of mRNAs, circRNAs, and microRNAs (miRNAs) in response to Edwardsiella tarda infection at different time points (0, 2, 6, 12, and 24 h) and to construct regulatory networks for exploring the immune regulatory mechanism in the intestine of S. schlegelii. In total, 1,794 mRNAs, 87 circRNAs, and 79 miRNAs were differentially expressed. The differentially expressed RNAs were quantitatively validated using qRT-PCR. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis showed that most of the differentially expressed mRNA genes and the target genes of ncRNAs were related to immune signaling pathways, such as the NF-κB signal pathway, pathogen recognition receptors related to signaling pathways (Toll-like receptors and Nod-like receptors), and the chemokine signaling pathway. Based on these differentially expressed genes, 624 circRNA-miRNA pairs and 2,694 miRNA-mRNA pairs were predicted using the miRanda software. Integrated analyses generated 25 circRNA-miRNA-mRNA interaction networks. In a novel_circ_0004195/novel-530/IκB interaction network, novel_530 was upregulated, while its two targets, novel_circ_0004195 and IκB, were downregulated after E. tarda infection. In addition, two circRNA-miRNA-mRNA networks related to apoptosis (novel_circ_0003210/novel_152/apoptosis-stimulating of p53 protein 1) and interleukin (novel_circ_0001907/novel_127/interleukin-1 receptor type 2) were also identified in our study. We thus speculated that the downstream NF-κB signaling pathway, p53 signaling pathway, and apoptosis pathway might play vital roles in the immune response in the intestine of S. schlegelii. This study revealed a landscape of RNAs in the intestine of S. schlegelii during E. tarda infection and provided clues for further study on the immune mechanisms and signaling networks based on the circRNA-miRNA-mRNA axis in S. schlegelii.

Full length transcriptome profiling reveals novel immune-related genes in black rockfish (Sebastes schlegelii)

Lacking full-length transcriptome for black rockfish (Sebastes schlegelii) limits novel gene discoveries and gene structures analysis. Therefore, we constructed the full-length transcriptome of black rockfish using Single-Molecule Real-Time Sequencing technology. Totally, we produced 21.73 Gb raw reads containing 298,904 circular consensus sequence (CCS) reads. Full-length (FL) and Non-full-length (NFL) isoforms were obtained based on the presence of 5' and 3' primers as well as poly (A) tails. The results showed 70.71% reads were identified as FL isoforms. Moreover, the average length of these PacBio isoforms is 2,632 bp, which is much longer than the length of the unigenes with the average length of 589 bp which generated from Illumina platform. Meanwhile, we identified 43,068 non-redundant transcripts, 12,485 alternative splicing (AS), 6,320 polyadenylation (APA) and 499 gene fusions as well as numerous long non-coding RNAs based on mapped FL isoforms. In addition, we identified 147 and 528 immune-related genes from novel genes and unmapped transcripts. The provided dataset can be utilized to discover novel genes and construct a comprehensive transcript dataset for black rockfish.

Genome-wide characterization of Toll-like receptors in black rockfish Sebastes schlegelii: Evolution and response mechanisms following Edwardsiella tarda infection

As one of the key components of pattern recognition receptors, Toll-like receptors (TLRs) play pivotal roles in the innate immune system. However, little information is available about the TLR genes in Sebastes schlegelii. In the present study, 17 TLR genes were identified and classified based on the whole genome database. Tandem duplication events in TLR1, TLR2, TLR5 and TLR13 played major role in the expansion of S. schlegelii TLR genes; both TLR2-3 and TLR2-4 had the same largest number of introns/exons, 11 exons and 10 introns. The syntenic analysis showed neighboring genes of TLR genes were most conserved in S. schlegelii and in L. crocea. Phylogenetic and evolutionary analysis showed that these TLR genes were divided into five subfamilies and exhibited different selection pressures. Meanwhile, the expression patterns of TLR genes in the intestine after E. tarda infection were investigated by qRT-PCR. Finally, protein and protein interaction (PPI) network analysis indicated that TLR genes interacted with IFN-inducible genes, inflammatory cytokines, and participated in MyD88-dependent pathway. In summary, this study provided valuable information for further functional characterization of TLR genes in S. schlegelii.