Yi Zhang
Mengqi Ni
Yunhui Bai
Qiao Shi
Jinbin Zheng
Zhaoxia Cui- 1School of Marine Sciences, Ningbo University, Ningbo, China
- 2Laboratory for Marine Biology and Biotechnology, Pilot Qingdao National Laboratory for Marine Science and Technology (Qingdao), Qingdao, China
Generally, invertebrates were thought to solely rely on their non-specific innate immune system to fight against invading microorganisms. However, increasing studies have implied that the innate immune response of invertebrates displayed diversity and specificity owing to the hyper-variable immune molecules in organisms. In order to get an insight into the diversity of immune-related genes in Portunus trituberculatus, a full-length transcriptome analysis of several immune-related tissues (hemocytes, hepatopancreas and gills) in P. trituberculatus was performed and the diversity of several immune-related genes was analyzed. The full-length transcriptome analysis of P. trituberculatus was conducted using a combination of SMRT long-read sequencing and Illumina short-read sequencing. A total of 17,433 nonredundant full-length transcripts with average length of 2,271 bp and N50 length of 2,841 bp were obtained, among which 13,978 (80.18%) transcripts were annotated. Moreover, numerous transcript variants of various immune-related genes were identified, including pattern recognition receptors, antimicrobial peptides, heat shock proteins (HSPs), antioxidant enzymes and vital molecules in prophenoloxidase (proPO)-activating system. Based on the full-length transcriptome analysis, open reading frames (ORFs) of four C-type lectins (CTLs) were cloned, and tissue distributions showed that the four CTLs were ubiquitously expressed in all the tested tissues, and mainly expressed in hepatopancreas and gills. The transcription of the four CTLs significantly increased in several immune-related tissues (hemocytes, hepatopancreas and gills) of P. trituberculatus challenged with Vibrio alginolyticus and displayed different profiles. Moreover, the four CTLs displayed distinct bacterial binding and antibacterial activities. The recombinant protein PtCTL-1 (rPtCTL-1) and rPtCTL-3 displayed bacterial binding and antibacterial activities against all tested bacteria. rPtCTL-2 only showed bacterial binding and antibacterial activities against V. alginolyticus. No obvious bacterial binding or antibacterial activities for PtCTL-4 was observed against the tested bacteria. This study enriches the transcriptomic information on P. trituberculatus and provides new insights into the innate immune system of crustaceans. Additionally, our study provided candidates of antibiotic agents for the prevention and treatment of bacteriosis.
Introduction
The swimming crab, Portunus trituberculatus, which is extensively distributed in East and Southeast Asia, is one of the most important commercial crabs in China (1). In the recent decade, the farming industry of P. trituberculatus has developed rapidly in China, and the annual production of P. trituberculatus in China is remain at approximately 100,000 tons (1). However, frequent outbreaks of diseases have resulted in severe economic losses and impede the sustainable development of the P. trituberculatus aquaculture industry. Mass mortalities of crustaceans attributed to outbreaks of vibriosis have been recorded from many regions (2–4). The vibriosis caused by Vibrio alginolyticus and V. parahaemolyticus has been responsible for the costliest disease (4–6). In recent years, in addition to Vibrio species, pathogenic bacteria such as Staphylococcus aureus has also become a ubiquitous pathogen in crab culture system (7). Coinfection with various pathogen seem to become a problematic for the crab farming industry. Since the abuse of antibiotics in aquaculture have resulted in serious antibiotics contamination and antibiotic resistance (8, 9), it is urgent to get a better understanding of the host immune system, which would be benefit for developing effective methods to control diseases.
Identification and functional characterization of immune molecules is basic for investigating immune characteristics of organisms. Moreover, several immune effector molecules are promising candidates of therapeutic agents (10–15). Recently, with the development of third-generation sequencing (TGS) technique, full-length transcriptome analysis has become an attractive high-throughput approach to excavate transcripts at isoform-level (16, 17). To date, an increasing number of full-length transcriptome analysis have been conducted in a variety of commercially important aquatic organisms using the TGS technique (18–22), which greatly facilitate the transcriptome research of aquatic animals lacking high-quality reference genome and provide technology for the development and utilization of antibacterial substances. However, the investigation of the full-length transcriptome has not yet been performed in P. trituberculatus.
In this context, with the purpose to get new insights into the innate immune system of P. trituberculatus and excavate candidate agents to control and cure the bacteriosis, the full-length transcriptome analysis of several immune-related tissues (hemocytes, hepatopancreas and gills) in P. trituberculatus was conducted. The tremendous transcript variants of immune-related molecules identified in this study provide new insight into the immune characteristics of P. trituberculatus. The data of full-length transcripts obtained in this study provided valuable genomic resources for further research into the new gene discovery and genomic research on P. trituberculatus. Moreover, this study provides candidates of antibiotic agents for the prevention and treatment of bacteriosis.
Materials and Methods
Animals and Samples
P. trituberculatus were purchased from a local seafood market in Ningbo (Zhejiang Province, China). The crabs were anesthetized by chilling on ice before dissection. For the full-length transcriptome analysis, three immune-related tissues (hemocytes, hepatopancreas and gills) were collected from P. trituberculatus with a body weight of 262.5 g. Crabs with a body weight of 87.72 ± 8.96 g were used for tissue distributions and immune response analyses. Crabs were acclimated in rectangular tanks (60 cm × 40 cm × 25 cm) with aerated seawater (25‰), each tank containing 8 individuals, and the seawater was renewed daily. Hemocytes, hepatopancreas, gills, stomach, intestine, heart, muscle and eyestalk were harvested from five crabs used for tissue distributions. For the immune challenge experiment, crabs were injected with 100 μL V. alginolyticus (2 × 108 CFU/mL) resuspended in 0.1 M phosphate buffered saline (PBS) into the arthrodial membrane of the last walking leg. All of the bacterial strains used in this study were provided by Xiamen University (Xiamen, China). Crabs injected with 0.1 M PBS served as controls. At 3 h and 6 h post injection (hpi), hemocytes, hepatopancreas and gills were collected from five randomly selected crabs from each group. All samples were immediately frozen in liquid nitrogen and stored at −80°C until the extraction of total RNA for analysis.
RNA Extraction and Qualification
Total RNA was extracted using RNAiso Plus (TaKaRa, Japan) according to the manufacturer’s instructions. RNA degradation was assessed on 1% agarose gels. The NanoPhotometer® spectrophotometer (IMPLEN, CA, USA), Qubit®2.0 Fluorometer (Life Technologies, CA, USA) and Agilent Bioanalyzer 2100 system (Agilent Technologies, CA, USA) were used to detect the purity, concentration and integrity of the RNA samples.
PacBio and Illumina Library Construction and Sequencing
Equal amounts of total RNA extracted from the three tissues were mixed to construct the library for PacBio and Illumina sequencing. The Isoform Sequencing (Iso-Seq) library was prepared according to the Iso-Seq protocol using the SMARTer PCR cDNA Synthesis Kit (Takara, Japan) and the BluePippin Size Selection System protocol as described by Pacific Biosciences (Sage science, USA). The qualified libraries were sequenced on a Pacbio Sequel platform.
Illumina sequencing library was constructed by pooling equal amounts of total RNA extracted from the three tissues. Sequencing libraries were generated using NEBNext® Ultra™ RNA Library Prep Kit for Illumina® (NEB, USA) following manufacturer’s recommendations and library quality was assessed on the Agilent Bioanalyzer 2100 system. The qualified libraries were sequenced on an Illumina Novaseq platform and 150 bp paired-end reads were generated.
Data Processing
Raw polymerase reads were processed using the SMRTlink 5.0 software to obtain sub-reads. Circular consensus (CCS) reads were generated from the sub-reads via self-correction, and CCS reads contained the 5’ primer, 3’ primer and polyA tail were considered as full-length transcript sequence. And then, full-length non-chimeric (FLNC) reads were subjected to cluster and non-redundant treatment using hierarchical n*log(n) algorithm and Arrow polishing, and polished consensus reads were obtained. Finally, LoRDEC was used to correct these polished consensus reads with short reads yielded by Illumina sequencing, and CD-HIT program was further used to remove redundancy. Through above data processing procedure, non-redundant and accuracy full-length transcripts were produced for following analysis.
Functional Annotation
For functional annotation, all transcripts were searched against the databases NCBI non-redundant nucleotide sequences (NT), NCBI non-redundant protein sequences (NR), Protein family (Pfam), Clusters of Orthologous Groups of proteins (COG), Swiss-Prot, Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) and KEGG Ortholog database (KO). The BLAST programs with an E-value threshold of 10−10 were used in the NT database analysis. The Diamond BLASTX programs with an E-value threshold of 10−10 was used in the NR, KOG, Swiss-Prot and KEGG databases analysis. The Hmmscan programs was used in the Pfam database analysis.
Prediction of Alternative Splicing (AS) Events
Alternative splicing analysis of non-redundant transcripts were processed with the Coding GENome reconstruction Tool (Cogent v3.1, https://github.com/Magdoll/Cogent).
Sequence Analysis
Identities of deduced aa sequences were determined using BioEdit software (v7.1.3). Functional domains were predicted by SMART program (http://smart.embl-heidelberg.de/).
Cloning of the Open Reading Frames (ORFs) Sequences of CTLs
Four CTLs, named PtCTL-1, PtCTL-2, PtCTL-3 and PtCTL-4, were selected for the immune function analysis. Extraction and qualification of total RNA was conducted as above mentioned. Reverse transcription to synthesize first-strand cDNA was carried out using PrimeScript™ RT reagent Kit with gDNA Eraser (TaKaRa, Japan) following the manufacturer’s protocols. The primers (Table 1) used for cloning were designed based on the sequence of the CTLs identified from the full-length transcriptome of P. trituberculatus. The amplification products were ligated into the pMD19-T vector (Takara, Japan) and sequenced (Sangon, China).
Table 1 Primers used in this study.
Quantitative Real-Time RT-PCR
Quantitative real-time PCR (qRT-PCR) was carried out on an ABI 7500 real-time PCR detection system (Applied Biosystems, USA) using TB Green Premix Dimer Eraser kit (Takara, Japan) according to the manufacturer’s instructions. The specific primers used for qRT-PCR were designed based on the ORF sequence of the CTLs and listed in Table 1. The β-actin gene of P. trituberculatus were served as the internal standardization. The melting curve was recorded to confirm that only a single PCR product was amplified. Three technical replicates were conducted. The relative gene expression levels were determined by the 2-△△Ct method (23).
Expression and Purification of Recombinant Proteins
The cDNA fragments encoding the mature peptides of the four PtCTLs were amplified using primers with the restriction sites (Table 1). The target cDNA fragments were ligated into the expression vector pET-32a and confirmed by sequencing. The recombinant plasmids were transformed into Escherichia coli BL21 (DE3) pLysS Chemically Competent Cells (TransGen Biotech, China) and induced at 16°C for 12 h with isopropyl-β-thiogalactopyranoside (IPTG) at a final concentration of 0.5 mM. The induced bacterial pellets were collected after centrifugation (10,000 ×g) for 5 min at 4°C and then suspended in Tris-HCl buffer (50 mM Tris-HCl, 500 mM NaCl, 20 mM imidazole, pH 7.8) and sonicated on ice for 30 min with burst duration of 3 s. The supernatants were collected after centrifugation (10,000 ×g) for 25 min at 4°C. The purification of recombinant proteins was conducted using the high affinity Ni-NTA resin (GenScript, China) according to the manufacturer’s protocols. The purity of the recombinant proteins was analyzed by sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) and Coomassie Blue staining, and the recombinant proteins were further confirmed by western blot analysis as previously described (24). The concentrations of recombinant proteins were examined using BCA Protein Assay Kit (Beyotimie, China).
Antibacterial and Bacterial Binding Activities Assay
Gram-negative bacteria (V. alginolyticus and V. parahaemolyticus) and Gram-positive bacteria (S. aureus) were tested in the antibacterial and bacterial binding activities assay. V. alginolyticus, V. parahaemolyticus and S. aureus were suspended in TBS (50 mM Tris-HCl, 50 mM NaCl, 10 mM CaCl2, pH 7.5) at a concentration of 6.5 × 107 CFU/mL, 6.8 × 107 CFU/mL and 6.0 × 107 CFU/mL, respectively. 50 μL of recombinant PtCTL proteins (rPtCTLs) (500 μg/mL) were mixed with the same volume of bacteria suspension and incubated at 37°C (V. alginolyticus and S. aureus) or 30°C (V. parahaemolyticus) for 2 h under slight rotation. The mixtures were diluted 104-fold and spread on Luria-Bertani (S. aureus) or 2216E (V. alginolyticus and V. parahaemolyticus) agar plates, and the CFU were calculated after culturing for 8 h at 37°C/30°C. In the binding assay, the rPtCTLs and bacteria suspension were mixed and incubated as above mentioned. The bacterial pellets were harvested by centrifugation at 5,000 ×g for 3 min and washed five times with TBS. Bindings between rPtCTLs and bacteria was confirmed by western blot analysis.
Statistical Analysis
All data were shown as the mean ± standard deviation (SD) and analyzed by the one-way analysis of variance (ANOVA) using the SPSS 17.0 software. 0.01 < p < 0.05 was considered statistically significant; p < 0.01 was considered to be extremely significant.
Results
Overview of the Sequencing Data
In total, 1,172,590 polymerase reads were obtained. The average length of polymerase reads was 75,948 bp, and the N50 length of polymerase reads was 141,113 bp (Supplementary Table 1). A total of 64,159,133 sub-reads with the average length of 1,314 bp and N50 length of 1,900 bp were generated from polymerase reads (Supplementary Table 1). Based on self-correction among sub-reads, 1,061,265 CCS reads were obtained, and the average length and N50 length of CCS reads were 2,006 bp and 2,471 bp, respectively (Supplementary Table 1). Subsequently, among the CCS reads, 668,918 full-length non-chimeric (FLNC) reads with average length of 1,760 bp and N50 length of 2,266 bp were identified (Supplementary Table 1). And then, 43,748 polished consensus reads with N50 length of 2,522 bp were produced from the FLNC reads (Supplementary Table 1). Finally, the polished consensus reads were subjected to correction using short reads produced by Illumina sequencing and removing redundancy. In total, 17,433 nonredundant full-length transcripts with average length of 2,271 bp and N50 length of 2,841 bp were obtained (Supplementary Table 1). Raw PacBio and Illumina sequencing reads are available at NCBI GenBank under accession numbers SRR12442560 and SRR12442561.
Functional Annotation of Transcripts
In total, 13,978 (80.18%) transcripts were annotated in at least one database. A total of 12,972 (74.41%), 6,371 (36.55%), 10,172 (58.35%), 11,296 (64.80%), 10,172 (58.35%), 10,126 (58.09%), 12,343 (70.80%) transcripts were annotated in the NR, NT, Pfam, SwissProt, GO, KOG and KEGG database, respectively (Supplementary Table 2). For main species distribution matched against NR database, most of the matched transcripts were similar with Hyalella azteca (4,612), followed by Zootermopsis nevadensis (649), Portunus trituberculatus (506), Scylla paramamosain (473), Limulus polyphemus (344), Daphnia magna (332), Eriocheir sinensis (276) and Litopenaeus vannamei (209) (Supplementary File 1).
The distribution of GO second-level terms showed that cellular process (4,907), metabolic process (4,506) and single-organism process (3,678) were the most represented terms in the biological process category (Figure 1 and Supplementary File 2). In the molecular function category, binding (6,075) and catalytic activity (4,173) comprised the large proportion (Figure 1 and Supplementary File 2). As to cellular component category, most transcripts were associated with cell (2,344), cell parts (2,344) and organelle (1,739) (Figure 1 and Supplementary File 2). The results of COG annotation showed that the cluster for general function prediction only (1,837), signal transduction mechanisms (1,368) and post-translation modification, protein turnover, chaperones (1,031) represented the top three class (Supplementary File 3). For the KEGG annotation, most transcripts were mapped to signal transduction (1,281) and transport and catabolism (911) related pathways (Figure 2 and Supplementary File 4). Meanwhile, a total of 561 transcripts were mapped to immune system related pathways, such as chemokine signaling pathway, complement and coagulation cascades, Toll-like receptor signaling pathway, NOD-like receptor signaling pathway, RIG-I-like receptor signaling pathway and IMD signaling pathway (Figure 2 and Supplementary File 4).
Figure 1 GO classification of the transcripts.
Figure 2 KEGG pathway classification of the transcripts. The number above the bars indicates the number of transcripts annotated in different KEGG pathways.
Alternative Splicing Events Analysis
In total, 3,960 unique transcript models (UniTransModels) were reconstructed by 8,563 non-redundant full-length transcripts. 1,684 UniTransModels that have more than one isoform were used for alternative event analysis. As a result, a total of 315 AS events were predicted based on 244 UniTransModels and generate 483 full-length transcripts. The retained intron (RI) (167) is the most abundant type of AS events, followed by alternative 3’ splice sites (A3) (72), alternative 5’ splice sites (A5) (61), skipping exon (SE) (9), alternative first exons (AF) (4), alternative last exons (AL) (1) and mutually exclusive exons (MX) (1).
Transcript Variants of Immune-Related Genes
In this study, transcript variants of various immune-related genes such as pattern recognition receptors (PRRs), antimicrobial peptides (AMPs), vital molecules in prophenoloxidase (proPO)-activating system, antioxidant enzymes, and heat shock proteins (HSPs) were identified (Table 2). The ORFs and deduced amino acid (aa) sequences of these transcripts are listed in Supplementary File 5. A total of twenty-five C-type lectins (CTLs) were identified in this study. The functional domain analysis showed that the deduced aa sequences of the CTLs all contain a carbohydrate recognition domain (CRD), moreover, a signal peptide and a low-density lipoprotein receptor class A domain were also observed in some CTLs (Figure 3A). Homology comparisons showed that the deduced aa sequences of these CTLs shared 2.2% to 99.6% identity with each other (Supplementary File 6). In total, four Toll-like receptors (TLRs) were identified, and these TLRs showed 19.8% to 94.2% identity with each other (Supplementary File 6). These TLRs exhibited a typical domain architecture of their counterparts, including extracellular leucine-rich repeats (LRRs), a transmembrane domain and an intracellular Toll/interleukin-1 receptor (TIR) domain (Figure 3B). Herein, a total of twelve crustins, seven antilipopolysaccharide factors (ALFs) and fourteen hyastatin were identified. Homology comparisons showed that the deduced aa sequences of the crustins, ALFs and hyastatin shared 7.2% to 99%, 4% to 95.7% and 24% to 93.1% identity with each other, respectively (Supplementary File 6). In the current study, twenty-four proPO-activating factors (PPAFs) and twenty-one serine proteinases transcript variants were identified. Homology comparisons showed that the deduced aa sequences of the PPAFs and serine proteinases shared 10.1% to 98.6% and 7.1% to 96% identity with each other, respectively (Supplementary File 6). Moreover, several antioxidant enzymes, including twelve copper/zinc SODs (Cu/ZnSODs), two cytosolic manganese SODs (CytMnSODs), three CATs, three Gpxs and fifteen GSTs were identified. Additionally, we also identified several HSPs, including a HSP10, two HSP40, two HSP60, six HSP70 and two HSP90.
Table 2 Statistics of transcript variants of immune-related genes.
Figure 3 Domain architecture of CTLs (A) and TLRs (B) predicted by SMART. The sequences and transcript IDs are listed in Supplementary File 5 and Supplementary File 6.
Cloning of PtCTLs cDNA
The cDNA sequences of PtCTL-1, PtCTL-2, PtCTL-3, and PtCTL-4 were deposited in GenBank with the accession number OL848444, OL848445, OL848446 and OL848447, respectively. The four PtCTLs encoded a polypeptide of 161 aa (18.08 kDa), 265 aa (29.68 kDa), 175 aa (20.62 kDa) and 164 aa (18.45 kDa), respectively. All of the four PtCTLs contained a single CRD with a characteristic carbohydrate binding motif (a QPN motif in PtCTL1, an EPS motif in PtCTL2, a FPR motif in PtCTL3 and an EPD motif in PtCTL4) and six cysteine residues (Cys23, Cys40, Cys57, Cys121, Cys153 and Cys156 in PtCTL1; Cys136, Cys147, Cys164, Cys231, Cys247 and Cys255 in PtCTL2; Cys33, Cys44, Cys61, Cys134, Cys150 and Cys158 in PtCTL3; Cys20, Cys35, Cys52, Cys119, Cys135 and Cys156 in PtCTL4) (Figure 4). Additionally, the PtCTL1, PtCTL2 and PtCTL3 also contained a signal peptide (Met1-Ala21 in PtCTL1, Met1-Ala16 in PtCTL2 and Met1-Gly22 in PtCTL3) in the N-terminus (Figures 4A–C).
Figure 4 cDNA and deduced amino acid sequence of PtCTLs. Red letters indicate the start sites and blue letters represent the stop codons. The asterisk indicates the stop codon. The signal peptide is shaded in yellow and the CRD is highlighted in green. Highly conserved cysteine residues are shown by black boxes. The carbohydrate binding motif is underlined in red. (A) PtCTL-1; (B) PtCTL-2; (C) PtCTL-3; (D) PtCTL-4.
Tissue Distribution and Expression Profiles of the Four CTLs in P. trituberculatus Challenged With V. alginolyticus
Tissue distribution analysis showed that the four PtCTLs were expressed in all tested tissues and showed relatively higher expression levels in the hepatopancreas and gills (Figure 5). Significant induction in the transcription of PtCTL-1, PtCTL-2, and PtCTL-4 were observed in the hemocytes, hepatopancreas and gills of P. trituberculatus challenged with V. alginolyticus, and significant expression of PtCTL-3 was detected in hemocytes and gills when V. alginolyticus challenged the P. trituberculatus (Figure 6). The transcriptional level of PtCTL-1 in the hemocytes was up-regulated 4.59-fold (P<0.01) at 3 hpi (Figure 6A), and it was up-regulated 4.50-fold (P<0.01) at 6 hpi (Figure 6B) in the hepatopancreas. The transcription of PtCTL-1 in the gills was up-regulated 2.07-fold and 2.65-fold (P<0.05) at 3 and 6 hpi (Figure 6C), respectively. The transcripts of PtCTL-2 in hemocytes increased 19.81-fold and 5.57-fold (P<0.05) at 3 hpi and 6 hpi (Figure 6D), respectively. The transcripts of PtCTL-2 in hepatopancreas increased 5.73-fold (P<0.05) at 3 hpi (Figure 6E), and it was up-regulated 1.98-fold (P<0.05) at 6 hpi in gills (Figure 6F). The mRNA expression of PtCTL-3 in hemocytes was significantly up-regulated 7.18-fold and 2.83-fold (P<0.01) at 3 hpi and 6 hpi, respectively (Figure 6G). The expression level of PtCTL-3 achieved a 1.84-fold (P<0.05) increment in gills (Figure 6I), while it was not changed in hepatopancreas (Figure 6H). The transcription of PtCTL-4 in hemocytes (Figure 6J) and hepatopancreas (Figure 6K) increased 3.23-fold and 4.01-fold, respectively. The transcriptional level of PtCTL-4 in gills was up-regulated 2.10-fold (P<0.05) and 2.62-fold (P<0.01) at 3 hpi and 6 hpi, respectively (Figure 6L).
Figure 5 Spatial expression analysis of PtCTLs in different tissues. (A) PtCTL-1; (B) PtCTL-2; (C) PtCTL-3; (D) PtCTL-4. Expression levels in hepatopancreas, gills, stomach, intestine, heart, muscle and eyestalk are normalized to those in the hemocytes, and β-actin is used as a reference gene. Vertical bars represent the means ± SD, n = 5. Results were analyzed by one-way ANOVA and different letters (a-d) above the bars presented significant differences between groups (p < 0.05).
Figure 6 Expression profiles of PtCTLs in V. alginolyticus-challenged P. trituberculatus. The relative expression of PtCTL-1 in hemocytes (A), hepatopancreas (B) and gills (C) after V. alginolyticus injection; The relative expression of PtCTL-2 in hemocytes (D), hepatopancreas (E) and gills (F) after V. alginolyticus injection; The relative expression of PtCTL-3 in hemocytes (G), hepatopancreas (H) and gills (I) after V. alginolyticus injection; The relative expression of PtCTL-4 in hemocytes (J), hepatopancreas (K) and gills (L) after V. alginolyticus injection. Each bar represents the mean ± SD (n = 5). Asterisks indicate statistically significant differences: *p < 0.05, **p < 0.01.
Bacterial Binding and Antibacterial Activities of rPtCTLs
Recombinant proteins of the four PtCTLs were successfully produced and purified. Coomassie Blue staining (Figure 7A) and western blot (Figure 7B) showed that the recombinant proteins exhibited the expected sizes of approximately 23 kDa (His-Trx), 35 kDa (rPtCTL-1), 34 kDa (rPtCTL-2), 35 kDa (rPtCTL-3) and 36 kDa (rPtCTL-4). Bacterial binding activity assays showed that rPtCTL-1 and rPtCTL-3 could bind to V. alginolyticus, V. parahaemolyticus and S. aureus (Figures 7C–E). rPtCTL-2 only exhibited binding activity to V. alginolyticus (Figures 7C–E), while rPtCTL-4 could not bind to any tested bacteria (Figures 7C–E). The antibacterial activity assay showed that the four rPtCTLs displayed various antibacterial spectra against the tested bacteria. The rPtCTL-1 and rPtCTL-3 exhibited antibacterial activity against V. alginolyticus, V. parahaemolyticus and S. aureus (Figures 7F–H). The CFU of V. alginolyticus treated with rPtCTL-1, rPtCTL-2 and rPtCTL-3 significantly decreased 2.41-fold, 5.35-fold and 18.88-fold, respectively. Compared with the blank control group, the CFU of V. parahaemolyticus treated with rPtCTL-1 and rPtCTL-3 were significantly decreased by 3.29-fold and 2.72-fold, respectively. The CFU of S. aureus treated with rPtCTL-1 and rPtCTL-3 significantly decreased 3.32-fold and 3.79-fold, respectively, compared with the blank control group. No obvious antimicrobial activity of rPtCTL-4 was observed (Figures 7F–H).
Figure 7 The bacteria binding activity and the antibacterial activities of recombinant CTLs. M, protein molecular standard (kDa); lane 1, rTrx-His; lane 2, PtCTL-1; lane 3, PtCTL-2; lane 4, PtCTL-3; lane 5, PtCTL-4. (A) SDS-PAGE analyses of purified recombinant proteins. (B) Western blot analyses of purified recombinant proteins. (C–E) Bacterial binding activity of recombinant PtCTLs proteins against V. alginolyticus (C), V. parahaemolyticus (D) and S. aureus (E). (F–H) Antibacterial activity of recombinant PtCTLs proteins against V. alginolyticus (F), V. parahaemolyticus (G) and S. aureus (H). Results were analyzed by one-way ANOVA and different letters (a-d) above the bars presented significant differences between groups at p < 0.05 (n = 3).
Discussion
Full-length cDNA sequence is basic for functional studies of genes. In the present study, to get an insight into the diversity of immune-related molecules in the P. trituberculatus, PacBio SMRT sequencing technology was used for identifying full-length transcripts in immune-related tissues of P. trituberculatus. In this study, a total of 17,433 nonredundant full-length transcripts with an average length of 2,271 bp and N50 length of 2,841 bp were obtained, which was longer than previously reported transcriptomes of P. trituberculatus obtained by NGS technologies (25–30). To obtain comprehensive information on gene function, transcripts were subjected to annotation analysis by searching against databases including NR, NT, Pfam, SwissProt, KOG, GO and KEGG. The results showed that the annotated rate (80.18%) of the full-length transcripts was much higher than that of transcripts obtained by next-generation sequencing (NGS) (25, 26, 29, 30).
Traditionally, be absent of adaptive immune system, invertebrates were thought to solely rely on their non-specific innate immune system to combat against invading microorganisms (31, 32). Recently, increasing studies revealed that the innate response of invertebrates displayed diversity and specificity (33–35). Unlike the antigen specificity and immunological memory of adaptive immunity, the hyper-variable immune molecules might be an underlying molecular mechanism of pathogen-specific innate immune responses in invertebrates (32, 36, 37). In this study, transcript variants of various immune-related genes were identified.
Pattern recognition is the primary step to implement an efficient immune response, in which PRRs sense the invading microorganisms and activate immune responses (38–40). To date, numerous PRRs have been identified and functionally studied in invertebrates, and PRRs have been found to display diverse and specific immune functions (41, 42). Herein, transcript variants of several PRRs, including CTL, TLR, galectin and scavenger receptor B (SRB), were identified. CTLs are a group of highly diverse PRRs that recognize various pathogens (41, 43). A total of 25 CTLs were identified in this study, most (20) of them contain a signal peptide and a CRD domain, and two CTLs with a low-density lipoprotein receptor class A domain (LDLa) were identified. The TLR-1 showed differences in the numbers and arrangement of LRRs compared with other three TLRs. Since the LRRs involved in the binding of microorganisms (44, 45), these TLRs may play diverse roles in recognizing invading pathogens. These highly abundant PRRs might contribute to the immune specificity of P. trituberculatus.
AMPs are effectors of the innate immune system to fight against invading microorganisms. Several AMPs in crab, such as crustin, ALF and hyastatin, have been identified and characterized in crab, and different types and isoforms of AMPs exhibit diverse and specific antimicrobial activities (25, 46–60). The massive variants of AMPs might be crucial for P. trituberculatus to eliminate diverse invading microorganisms.
ProPO-activating system is a crucial constituent of the immune system in invertebrates (51). To date, several genes associated with proPO-activating system, such as prophenoloxidase, PPAF, serine proteinase and its inhibitor, have been identified and characterized in crustacean species (52–57). In this study, numerous transcript variants of PPAFs and serine proteinases were identified. Considering the fact that proPO activation is mediated by PPAF and serine proteinase, the abundant variants of PPAFs and serine proteinases in P. trituberculatus may provide a molecular basis for the specific activation of proPO.
Reactive oxygen species (ROS) are considered as an important part of the innate immunity that participate in eliminating invading microorganisms (58). ROS can kill invading microorganisms; however, the overproduction of ROS will cause serious damage to cellular macromolecules (59). To maintain a balance of ROS, organisms have developed an antioxidant system to eliminate excessive ROS, which consist of various antioxidant enzymes such as superoxide dismutase (SOD), catalase (CAT), glutathione peroxidase (Gpx) and glutathione-S-transferase (GST) (60). In addition to antioxidant system, HSPs also play important roles in the maintenance of immune homeostasis by acting as molecular chaperones to maintain protein homeostasis (61). In the present study, various transcripts variants of antioxidant enzymes and HSPs were identified, which might equip P. trituberculatus to efficiently cope with adverse effects resulted from immune reaction.
Since abuse of antibiotics pose serious threats to the ecological environment and human health (8, 9), mining bioactive molecules to develop novel antimicrobial agents has attracted increasing attention of researchers. CTLs, which show high diversity in invertebrates, are promising candidates for antimicrobial agents for their antibacterial and antiviral activities (14, 15, 24, 35, 62). In this study, based on the full-length transcriptome analysis, the bacterial binding and antibacterial activities of four PtCTLs were investigated. The four PtCTLs were distinguished from previously identified CTLs in P. trituberculatus and shared less than 20.20% identity to CTLs identified in previous studies (GenBank: AGH68927.1, AHK59786.1, QEM39053.1, ATE51171.1, ATE51203.1, ACC86854.1, QIE05465.1, QGV11036.1). The four PtCTLs exhibited typical structures of CTLs, such as the CRD, a mannose-binding motif or a galactose-binding motif, four cysteine residues that were crucial for the formation of two disulfide bonds to stabilize the CRD structure. The tissue distribution showed that the four PtCTLs mainly expressed in the hepatopancreas and gills, which is consistent with the fact that the hepatopancreas and gills of crustaceans, two important tissues involved in immune response of crustaceans, are two main tissues that synthesizes CTLs (43, 63–66).
As multifunctional immune molecules in crustaceans, CTLs have been proven to be involved in diverse immune responses (43). In this study, to get a preliminary knowledge on the specific immune response of the four PtCTLs, the expression profiles of these PtCTLs in the hemocytes, hepatopancreas and gills of P. trituberculatus post-injection with V. alginolyticus were investigated. The transcriptional levels of the four PtCTLs showed different expression profiles after challenged with V. alginolyticus. The tissue-specific transcriptional profiles of the four PtCTLs implying the complicated and diverse innate immune response of P. trituberculatus.
T.he primary function of CTLs is act as PRRs that recognize and bind to the microorganisms (11, 15, 67–69). To further ascertain the roles of these PtCTLs, recombinant proteins of the PtCTLs were produced, and the bindings between the rPtCTLs and bacteria were analyzed. Bacterial binding assay showed that the four rPtCTLs displayed different bacterial binding activity. The rPtCTL-1 and rPtCTL-3 possessed the broadest spectrum of bacterial binding activity against the tested bacteria, including V. alginolyticus, V. parahaemolyticus and S. aureus, which was similar to previous studies (70). In addition to act as PRRs, CTLs also exhibit antibacterial activities (14, 24, 62, 64). Corresponding to the bacterial binding activities, antibacterial assays showed that rPtCTL-1 and rPtCTL-3 possessed the widest spectrum of antibacterial activity against the tested bacteria, which was similar with a CTL isolated from P. trituberculatus in previous study (70). The rPtCTL-2 only displayed bacterial binding and antibacterial activities against V. alginolyticus against the tested bacteria, whereas no bacterial binding and antibacterial activities of rPtCTL-4 was observed against the tested bacteria, which was different to previous studies (70–72). Differences in bacterial binding and antibacterial activities of the four rPtCTLs might due to the mutation of the carbohydrate binding motif (24, 73). Additionally, previous reports showed that the CRD recognized subtle differences in the arrangement and branching of the carbohydrate residues, which may result in specific sets of carbohydrate recognition profiles for CRD (74). These findings suggested that these diverse PtCTLs may play specific roles in the defense against bacteria, and the rPtCTL-1 and rPtCTL-3 might be a promising antibiotic agent.
Conclusion
In conclusion, a full-length transcriptome analysis of several immune-related tissues in P. trituberculatus was performed. In total, 17,433 nonredundant full-length transcripts were obtained. Massive transcript variants of various immune-related genes, including CTLs, TLRs, crustins, ALFs, hyastatins, SODs, CATs, Gpxs, GSTs, HSP40, HSP60, HSP70, HSP90, PPAFs and serine proteinase, were identified. Moreover, the bacterial binding and antibacterial activities of four PtCTLs were investigated. The four PtCTLs exhibited diverse transcriptional profiles when challenged with V. alginolyticus, and displayed different bacterial binding and antibacterial activities against bacteria. These results will facilitate the excavation of immune molecules in P. trituberculatus and provide clues for further investigating the diversity and specificity of innate immune response in crustaceans. Furthermore, the widest antibacterial spectrum exhibited by rPtCTL-1 and rPtCTL-3 implied their promising application prospects in prevention and treatment of bacteriosis.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.
Ethics Statement
Ethical review and approval was not required for the animal study because Portunus trituberculatus is not an endangered or protected species, and there is no requirement for ethical approval to perform experiments involving this species in China.
Author Contributions
Conceived and designed the experiments, ZC and JZ. Performed the experiments, YZ, MN, YB, and QS. Analyzed the data and wrote the paper, JZ and YZ. All authors contributed to the article and approved the submitted version.
Funding
This research was supported by the National Key R&D Program of China (No. 2018YFD0900303), the National Natural Science Foundation of China (No. 32072964), the Ten Thousand Talents Program and the K. C. Wong Magna Fund of Ningbo University.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2022.843347/full#supplementary-material
References
1. MOA (Ministry of Agriculture). China Fishery Statistical Yearbook. Beijing, China: China Agriculture Press (2021). (in Chinese with English abstract).
2. Jayasree L, Janakiram P, Madhavi R. Characterization of Vibrio Spp. Associated With Diseased Shrimp From Culture Ponds of Andhra Pradesh (India). J World Aquacult Soc (2006) 37:523–32. doi: 10.1111/j.1749-7345.2006.00066.x
3. Heenatigala PPM, Fernando MUL. Occurrence of Bacteria Species Responsible for Vibriosis in Shrimp Pond Culture Systems in Sri Lanka and Assessment of the Suitable Control Measures, Sri Lanka. J Aquat Sci (2016) 21:1–17. doi: 10.4038/sljas.v21i1.7481
4. Kumar BK, Deekshit VK, Raj JRM, Rai P, Shivanagowda BM, Karunasagar I, et al. Diversity of Vibrio Parahaemolyticus Associated With Disease Outbreak Among Cultured Litopenaeus Vannamei (Pacific White Shrimp) in India. Aquaculture (2014) 433:247–51. doi: 10.1016/j.aquaculture.2014.06.016
5. Chen LL, Lo CF, Chiu YL, Chang CF, Kou GH. Natural and Experimental Infection of White Spot Syndrome Virus (WSSV) in Benthic Larvae of Mud Crab Scylla Serrata. Dis Aquat Organ (2000) 40:157–61. doi: 10.3354/dao040157
6. Wang W. Bacterial Diseases of Crabs: A Review. J Invertebr Pathol (2011) 106:0–26. doi: 10.1016/j.jip.2010.09.018
7. Lin Q, Zhang BY, Liu L, Wang XP, Lei CL, Lin YJ, et al. The Role of P38 MAPK, JNK, and ERK in Antibacterial Responses of Chilo Suppressalis (Lepidoptera: Crambidae). J Econ Entomol (2017) 110:1460–4. doi: 10.1093/jee/tox126
8. Carvalho E, David GS, Silva RJ. Health and Environment in Aquaculture . Available at: https://www.intechopen.com/.
9. Watts JEM, Schreier HJ, Lanska L, Hale MS. The Rising Tide of Antimicrobial Resistance in Aquaculture: Sources, Sinks and Solutions. Mar Drugs (2017) 15(6):158. doi: 10.3390/md15060158
10. Pham LN, Dionne MS, Shirasu-Hiza M, Schneider DS. A Specific Primed Immune Response in Drosophila Is Dependent on Phagocytes. PloS Pathog (2007) 3(3):e26. doi: 10.1371/journal.ppat.0030026
11. Wang XW, Zhao XF, Wang JX. C-Type Lectin Binds to Beta-Integrin to Promote Hemocytic Phagocytosis in an Invertebrate. J Biol Chem (2014) 289(4):2405–14. doi: 10.1074/jbc.M113.528885
12. Zhang XW, Wang Y, Wang XW, Wang L, Mu Y, Wang JX. A C-Type Lectin With an Immunoglobulin-Like Domain Promotes Phagocytosis of Hemocytes in Crayfish Procambarus Clarkii. Sci Rep (2016) 6:29924. doi: 10.1038/srep29924
13. Yu XQ, Kanost MR. Immulectin-2, a Pattern Recognition Receptor That Stimulates Haemocyte Encapsulation and Melanization in the Tobacco Hornworm, Manduca Sexta. Dev Comp Immunol (2004) 28:891e900. doi: 10.1016/j.dci.2004.02.005
14. Guo XN, Jin XK, Li S, Yu AQ, Wu MH, Tan SJ, et al. A Novel C-Type Lectin From Eriocheir Sinensis Functions as a Pattern Recognition Receptor With Antibacterial Activity. Fish Shellfish Immunol (2013) 35(5):1554–65. doi: 10.1016/j.fsi.2013.08.021
15. Fang ZY, Li D, Li XJ, Zhang X, Zhu YT, Li WW, et al. A Single CRD C-Type Lectin From Eriocheir Sinensis (EsLecB) With Microbial-Binding, Antibacterial Prophenoloxidase Activation and Hem-Encapsulation Activities. Fish Shellfish Immunol (2016) 50:175–90. doi: 10.1016/j.fsi.2016.01.031
16. Rhoads A, Au KF. PacBio Sequencing and its Applications. Genom Proteom Bioinf (2015) 13:278–89. doi: 10.1016/j.gpb.2015.08.002
17. van Dijk EL, Jaszczyszyn Y, Naquin D, Thermes C. The Third Revolution in Sequencing Technology. Trends Genet (2018) 34:666–81. doi: 10.1016/j.tig.2018.05.008
18. Zhang XJ, Li GY, Jiang HY, Li LM, Ma JG, Li HM, et al. Full-Length Transcriptome Analysis of Litopenaeus Vannamei Reveals Transcript Variants Involved in the Innate Immune System. Fish Shellfish Immunol (2019) 87:346–59. doi: 10.1016/j.fsi.2019.01.023
19. Zhang JY, Liu CL, He MC, Xiang ZL, Yin YN, Liu SF, et al. A Full-Length Transcriptome of Sepia Esculenta Using a Combination of Single-Molecule Long-Read (SMRT) and Illumina Sequencing. Mar Genom (2018) 43:54–7. doi: 10.1016/j.margen.2018.08.008
20. Lin JL, Shi X, Fang SB, Zhang Y, You CH, Ma HY, et al. Comparative Transcriptome Analysis Combining SMRT and NGS Sequencing Provides Novel Insights Into Sex Differentiation and Development in Mud Crab (Scylla Paramamosain). Aquaculture (2019) 513:734447. doi: 10.1016/j.aquaculture.2019.734447
21. Pootakham W, Uengwetwanit T, Sonthirod C, Sittikankaew K, Karoonuthaisiri N. A Novel Full-Length Transcriptome Resource for Black Tiger Shrimp (Penaeus Monodon) Developed Using Isoform Sequencing (Iso-Seq). Front Mar Sci (2020) 7:172. doi: 10.3389/fmars.2020.00172
22. Yi SK, Zhou XY, Li J, Zhang MM, Luo SS. Full-Length Transcriptome of Misgurnus Anguillicaudatus Provides Insights Into Evolution of Genus Misgurnus. Sci Rep (2018) 8:11699. doi: 10.1038/s41598-018-29991-6
23. Livak KJ, Schmittgen TD. Analysis of Relative Gene Expression Data Using Real-Time Quantitative PCR and the 2-ΔΔct Method. Methods (2001) 25:402–8. doi: 10.1006/meth.2001.1262
24. Zheng JB, Mao Y, Su YQ, Wang J. Identification of Two Novel C-Type Lectins Involved in Immune Defense Against White Spot Syndrome Virus and Vibrio Parahaemolyticus From Marsupenaeus Japonicus. Aquaculture (2020) 519:734797. doi: 10.1016/j.aquaculture.2019.734797
25. Chen FY, Wang KJ. Characterization of the Innate Immunity in the Mud Crab Scylla Paramamosain. Fish Shellfish Immunol (2019) 93:436–48. doi: 10.1016/j.fsi.2019.07.076
26. Wang ZF, Sun LX, Guan WB, Zhou CL, Tang BP, Cheng YX, et al. De Novo Transcriptome Sequencing and Analysis of Male and Female Swimming Crab (Portunus Trituberculatus) Reproductive Systems During Mating Embrace (Stage II). BMC Genet (2018) 19:3. doi: 10.1186/s12863-017-0592-5
27. Li YQ, Lai SM, Wang RJ, Zhao YC, Qin H, Jiang LX, et al. RNA-Seq Analysis of the Antioxidant Status and Immune Response of Portunus Trituberculatus Following Aerial Exposure. Mar Biotechnol (2017) 19:89–101. doi: 10.1007/s10126-017-9731-2
28. Yang YX, Wang JT, Han T, Liu T, Wang CL, Xiao J, et al. Ovarian Transcriptome Analysis of Portunus Trituberculatus Provides Insights Into Genes Expressed During Phase III and IV Development. PloS One (2015) 10:e0138862. doi: 10.1371/journal.pone.0138862
29. Meng XL, Liu P, Jia FL, Li J, Gao BQ. De Novo Transcriptome Analysis of Portunus Trituberculatus Ovary and Testis by RNA-Seq: Identification of Genes Involved in Gonadal Development. PloS One (2015) 10:e0128659. doi: 10.1371/journal.pone.0128659
30. Lv JJ, Liu P, Gao BQ, Wang Y, Wang Z, Chen P, et al. Transcriptome Analysis of the Portunus Trituberculatus: De Novo Assembly, Growth-Related Gene Identification and Marker Discovery. PloS One (2014) 9:e94055. doi: 10.1371/journal.pone.0094055
31. Iwanaga S, Lee BL. Recent Advances in the Innate Immunity of Invertebrate Animals. J Biochem Mol Biol (2005) 38:128–50. doi: 10.5483/BMBRep.2005.38.2.128
32. Schulenburg H, Boehnisch C, Michiels NK. How do Invertebrates Generate a Highly Specific Innate Immune Response. Mol Immunol (2007) 44:3338–44. doi: 10.1016/j.molimm.2007.02.019
33. Tetreau G, Pinaud S, Portet A, Galinier R, Gourbal B, David D. Specific Pathogen Recognition by Multiple Innate Immune Sensors in an Invertebrate. Front Immunol (2017) 8:1249. doi: 10.3389/fimmu.2017.01249
34. Schulenburg H, Hoeppner MP, Iii JW, Bornberg-Bauer E. Specificity of the Innate Immune System and Diversity of C-Type Lectin Domain (CTLD) Proteins in the Nematode Caenorhabditis Elegans. Immunobiology (2008) 213:237–50. doi: 10.1016/j.imbio.2007.12.004
35. Pees B, Yang W, Zárate-Potes A, Schulenburg H, Dierking K. High Innate Immune Specificity Through Diversified C-Type Lectin-Like Domain Proteins in Invertebrates. J Innate Immun (2015) 8:129–42. doi: 10.1159/000441475
36. Cerenius L, Lee BL, Soderhall K. The proPO-System: Pros and Cons for its Role in Invertebrate Immunity. Trends Immunol (2008) 29:263–71. doi: 10.1016/j.it.2008.02.009
37. Ghosh J, Lun CM, Majeske AJ, Sacchi S, Schrankel CS, Smith LC. Invertebrate Immune Diversity. Dev Comp Immunol (2011) 35:959–74. doi: 10.1016/j.dci.2010.12.009
38. Brubaker SW, Bonham KS, Zanoni I, Kagan JC. Innate Immune Pattern Recognition: A Cell Biological Perspective. Annu Rev Immunol (2015) 33:257–90. doi: 10.1146/annurev-immunol-032414-112240
39. Medzhitov R, Janeway CJ. Innate Immune Recognition: Mechanisms and Pathways. Immunol Rev (2010) 173:89–97. doi: 10.1034/j.1600-065x.2000.917309.x
40. Thompson MR, Kaminski JJ, Kurtjones EA, Fitzgerald KA. Pattern Recognition Receptors and the Innate Immune Response to Viral Infection. Viruses (2011) 3:920–40. doi: 10.3390/v3060920
41. Wang XW, Wang JX. Pattern Recognition Receptors Acting in Innate Immune System of Shrimp Against Pathogen Infections. Fish Shellfish Immunol (2013) 34:981–9. doi: 10.1016/j.fsi.2012.08.008
42. Zhang LL, Li L, Guo XM, Litman GW, Dishaw LJ, Zhang GF. Massive Expansion and Functional Divergence of Innate Immune Genes in a Protostome. Sci Rep (2015) 5:8693–703. doi: 10.1038/srep08693
43. Wang XW, Wang JX. Diversity and Multiple Functions of Lectins in Shrimp Immunity. Dev Comp Immunol (2013) 39:27e38. doi: 10.1016/j.dci.2012.04.009
44. Botos I, Segal DM, Davies DR. The Structural Biology of Toll-Like Receptors. Structure (2011) 19:447–59. doi: 10.1016/j.str.2011.02.004
45. Hughes AL, Piontkivska H. Functional Diversification of the Toll-Like Receptor Gene Family. Immunogenetics (2008) 60:249–56. doi: 10.1007/s00251-008-0283-5
46. Liu Y, Cui ZX, Luan WS, Song CW, Nie Q, Wang SY, et al. Three Isoforms of Anti-Lipopolysaccharide Factor Identified From Eyestalk cDNA Library of Swimming Crab Portunus Trituberculatus. Fish Shellfish Immunol (2011) 30:583–91. doi: 10.1016/j.fsi.2010.12.005
47. Liu Y, Cui ZX, Li XH, Song CW, Li QQ, Wang SY. Molecular Cloning, Expression Pattern and Antimicrobial Activity of a New Isoform of Anti-Lipopolysaccharide Factor From the Swimming Crab Portunus Trituberculatus. Fish Shellfish Immunol (2012) 33:85–91. doi: 10.1016/j.fsi.2012.04.004
48. Liu Y, Cui ZX, Li XH, Song CW, Shi GH. A Newly Identified Anti-Lipopolysaccharide Factor From the Swimming Crab Portunus Trituberculatus With Broad Spectrum Antimicrobial Activity. Fish Shellfish Immunol (2013) 34:463–70. doi: 10.1016/j.fsi.2012.11.050
49. Zhang Y, Wang LL, Wang LL, Yang JL, Gai YC, Qiu LM, et al. The Second Anti-Lipopolysaccharide Factor (EsALF-2) With Antimicrobial Activity From Eriocheir Sinensis. Dev Comp Immunol (2010) 34:945–52. doi: 10.1016/j.dci.2010.04.002
50. Mu CK, Zheng PL, Zhao JM, Wang LL, Qiu LM, Zhang H, et al. A Novel Type III Crustin (CrusEs2) Identified From Chinese Mitten Crab Eriocheir Sinensis. Fish Shellfish Immunol (2011) 31:142–7. doi: 10.1016/j.fsi.2011.04.013
51. Cerenius L, Soderhall K. Variable Immune Molecules in Invertebrates. J Exp Biol (2013) 216:4313–9. doi: 10.1242/jeb.085191
52. Cui ZX, Liu Y, Wu DH, Luan WS, Wang SY, Li QQ, et al. Molecular Cloning and Characterization of a Serine Proteinase Homolog Prophenoloxidase-Activating Factor in the Swimming Crab Portunus Trituberculatus. Fish Shellfish Immunol (2010) 29:679–86. doi: 10.1016/j.fsi.2010.07.015
53. Liu HR, Liu Y, Song CW, Ning JH, Cui ZX. Functional Characterization of Two Clip-Domain Serine Proteases in the Swimming Crab Portunus Trituberculatus. Fish Shellfish Immunol (2019) 89:98–107. doi: 10.1016/j.fsi.2018.12.047
54. Chen P, Li JT, Li J, Liu P, Gao BQ, Wang QY. Molecular Cloning and Characterization of Prophenoloxidase Gene in Swimming Crab Portunus Trituberculatus. Fish Shellfish Immunol (2010) 28:106–12. doi: 10.1016/j.fsi.2009.10.002
55. Wang J, Jiang KJ, Zhang FY, Song W, Zhao M, Wei HQ, et al. Characterization and Expression Analysis of the Prophenoloxidase Activating Factor From the Mud Crab Scylla Paramamosain. Genet Mol Res (2015) 14:8847–60. doi: 10.4238/2015.August.3.8
56. Yang Y, Bao C, Liu A, Ye H, Li S. Immune Responses of Prophenoloxidase in the Mud Crab Scylla Paramamosain Against Vibrio Alginolyticus Infection: In Vivo and In Vitro Gene Silencing Evidence. Fish Shellfish Immunol (2014) 39:237–44. doi: 10.1016/j.fsi.2014.05.014
57. Amparyup P, Charoensapsri W, Tassanakajon A. Prophenoloxidase System and its Role in Shrimp Immune Responses Against Major Pathogens. Fish Shellfish Immunol (2013) 34:990–1001. doi: 10.1016/j.fsi.2012.08.019
58. Kohchi C, Inagawa H, Nishizawa T, Soma G. ROS and Innate Immunity. Anticancer Res (2009) 29:817–21. doi: 10.1111/j.1469-8749.2009.03380.x
59. Sies H. Oxidative Stress: Oxidants and Antioxidants. Exp Physiol (1997) 82:291–5. doi: 10.1113/expphysiol.1997.sp004024
60. Davies KJA. Oxidative Stress, Antioxidant Defenses, and Damage Removal, Repair, and Replacement Systems. IUBMB Life (2000) 50:279–89. doi: 10.1080/713803728
61. Li FH, Xiang JH. Recent Advances in Researches on the Innate Immunity of Shrimp in China. Dev Comp Immunol (2013) 39:11–26. doi: 10.1016/j.dci.2012.03.016
62. Li M, Li C, Ma C, Li H, Zuo H, Weng S, et al. Identification of a C-Type Lectin With Antiviral and Antibacterial Activity From Pacific White Shrimp Litopenaeus Vannamei. Dev Comp Immunol (2014) 46(2):231–40. doi: 10.1016/j.dci.2014.04.014
63. Ma THT, Benzie JAH, He JG, Chan SM. PmLT, a C-Type Lectin Specific to Hepatopancreas Is Involved in the Innate Defense of the Shrimp Penaeus Monodon. J Invertebr Pathol (2008) 99(3):332–41. doi: 10.1016/j.jip.2008.08.003
64. Sun YD, Fu LD, Jia YP, Du XJ, Wang Q, Wang YH, et al. A Hepatopancreas-Specific C-Type Lectin From the Chinese Shrimp Fenneropenaeus Chinensis Exhibits Antimicrobial Activity. Mol Immunol (2008) 45(2):348–61. doi: 10.1016/j.molimm.2007.06.355
65. Alenton RRR, Koiwai K, Nakamura R, Thawonsuwan J, Kondo H, Hirono I. A Hint of Primitive Mucosal Immunity in Shrimp Through Marsupenaeus Japonicus Gill C-Type Lectin. J Immunol (2019) 203(8):2310–8. doi: 10.4049/jimmunol.1900156
66. Mistry AC, Honda S, Hirose S. Structure, Properties and Enhanced Expression of Galactose-Binding C-Type Lectins in Mucous Cells of Gills From Freshwater Japanese Eels (Anguilla Japonica). Biochem J (2001) 360(1):107–15. doi: 10.1042/bj3600107
67. Dodd RB, Drickamer K. Lectin-Like Proteins in Model Organisms: Implications for Evolution of Carbohydrate-Binding Activity. Glycobiology (2001) 11(5):71R–9R. doi: 10.1093/glycob/11.5.71R
68. Lee RT, Hsu TL, Huang SK, Hsieh SL, Wong CH, Lee YC. Survey of Immune-Related, Mannose/Fucose-Binding C-Type Lectin Receptors Reveals Widely Divergent Sugar-Binding Specificities. Glycobiology (2011) 21:512e520. doi: 10.1093/glycob/cwq193
69. Wongpanya R, Sengprasert P, Amparyup P, Tassansksjon A. A Novel C-Type Lectin in the Black Tiger Shrimp Penaeus Monodon Functions as a Pattern Recognition Receptor by Binding and Causing Bacterial Agglutination. Fish Shellfish Immunol (2017) 60:103–13. doi: 10.1016/j.fsi.2016.11.042
70. Su Y, Liu Y, Gao FT, Cui ZX. A Novel C-Type Lectin With a YPD Motif From Portunus Trituberculatus (PtCLec1) Mediating Pathogen Recognition and Opsonization. Dev Comp Immunol (2020) 106:103609. doi: 10.1016/j.dci.2020.103609
71. Huang MM, Mu CK, Wu YH, Ye F, Wang D, Sun C, et al. The Functional Characterization and Comparison of Two Single CRD Containing C-Type Lectins With Novel and Typical Key Motifs From Portunus Trituberculatus. Fish Shellfish Immunol (2017) 70:398–407. doi: 10.1016/j.fsi.2017.09.029
72. Liu Y, Su Y, Zhang A, Cui ZX. A C-Type Lectin Highly Expressed in Portunus Trituberculatus Intestine Functions in AMP Regulation and Prophenoloxidase Activation. Antibiotics (2021) 10(5):541. doi: 10.3390/antibiotics10050541
73. Zhou Z, Sun L. CsCTL1, a Teleost C-Type Lectin That Promotes Antibacterial and Antiviral Immune Defense in a Manner That Depends on the Conserved EPN Motif. Dev Comp Immunol (2015) 50(2):69–77. doi: 10.1016/j.dci.2015.01.007
Keywords: Portunus trituberculatus, full-length transcriptome, innate immune, transcript variants, c-type lectin
Citation: Zhang Y, Ni M, Bai Y, Shi Q, Zheng J and Cui Z (2022) Full-Length Transcriptome Analysis Provides New Insights Into the Diversity of Immune-Related Genes in Portunus trituberculatus. Front. Immunol. 13:843347. doi: 10.3389/fimmu.2022.843347
Received: 31 December 2021; Accepted: 16 March 2022;
Published: 07 April 2022.
Edited by:
Chris Hauton, University of Southampton, United KingdomCopyright © 2022 Zhang, Ni, Bai, Shi, Zheng and Cui. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jinbin Zheng, zhengjinbin@nbu.edu.cn