Comparative transcriptome analysis of red swamp crayfish (Procambarus clarkia) hepatopancreas in response to WSSV and Aeromonas hydrophila infection

The involvement of VEGF and VEGFR in bacterial recognition and regulation of antimicrobial peptides in Eriocheir sinensis

Vascular endothelial growth factor (VEGF) and VEGF reporter (VEGFR) are essential molecules in VEGF signalling pathway. Although the functions of VEGF and VEGFR have been well reported in vertebrates, their functions are still poorly understood in invertebrates. In this study, the open reading frame sequences of EsVEGF1 and EsVEGFR4 were cloned from Eriocheir sinensis, and their corresponding proteins shared typical structure characteristics with their counterparts in other species. EsVEGF1 were predominantly expressed in hepatopancreas and muscle while EsVEGFR4 mainly expressed in hemocytes and intestine. The expression levels of EsVEGF1 in hemocytes were rapidly induced by Staphylococcus aureus and Vibrio parahaemolyticus, and it also increased rapidly in hepatopancreas after being challenged with V. parahaemolyticus. The expression levels of EsVEGFR4 only increased in hepatopancreas of crabs injected with S. aureus. The extracellular immunoglobulin domain of EsVEGFR4 could bind with Gram-negative and Gram-positive bacteria as well as lipopolysaccharide and peptidoglycan. EsVEGF1 could act as the ligand for EsVEGFR4 and Toll-like receptor and regulate the expression of crustins and lysozyme with a tissue-specific manner, while have no regulatory function on that of anti-lipopolysaccharide factors. This study will provide new insights into the immune defense mechanisms mediated by VEGF and VEGFR in crustaceans.

Diapause-Linked Gene Expression Pattern and Related Candidate Duplicated Genes of the Mountain Butterfly Parnassius glacialis (Lepidoptera: Papilionidae) Revealed by Comprehensive Transcriptome Profiling

The mountain butterfly Parnassius glacialis is a representative species of the genus Parnassius, which probably originated in the high-altitude Qinhai–Tibet Plateau in the Miocene and later dispersed eastward into relatively low-altitude regions of central to eastern China. However, little is known about the molecular mechanisms underlying the long-term evolutionary adaptation to heterogeneous environmental conditions of this butterfly species. In this study, we obtained the high-throughput RNA-Seq data from twenty-four adult individuals in eight localities, covering nearly all known distributional areas in China, and firstly identified the diapause-linked gene expression pattern that is likely to correlate with local adaptation in adult P. glacialis populations. Secondly, we found a series of pathways responsible for hormone biosynthesis, energy metabolism and immune defense that also exhibited unique enrichment patterns in each group that are probably related to habitat-specific adaptability. Furthermore, we also identified a suite of duplicated genes (including two transposable elements) that are mostly co-expressed to promote the plastic responses to different environmental conditions. Together, these findings can help us to better understand this species’ successful colonization to distinct geographic areas from the western to eastern areas of China, and also provide us with some insights into the evolution of diapause in mountain Parnassius butterfly species.

Applying genetic technologies to combat infectious diseases in aquaculture

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

Transcriptome analysis of the immunomodulatory effects of Salvia miltiorrhiza polysaccharide on hemocyte immune response in Procambarus clarkii

Natural plant polysaccharide as immune modulator is considered an effective strategy for healthy aquaculture to reduce medicine treatment. Salvia miltiorrhiza polysaccharides (SMP) had applications to regulate immune activity and enhance antioxidant in vertebrates, but the potential function has been rarely reported in crustaceans. In this study, the immunological effects of SMP on hemocytes of Procambarus clarkii were analyzed. Results showed that total superoxide dismutase (T-SOD), phenoloxidase (PO) activity and respiratory burst were up-regulated after SMP treatment. After high-throughput sequencing, 2170 differentially expressed genes (DEGs) including 1294 up-regulated and 876 down-regulated genes were identified. KEGG function enrichment analysis indicated that DEGs are involved in crustaceans cellular immune-related signaling pathways, including lysosome, phagosome and endocytosis. Transcriptome mining and qRT-PCR showed that SMP up-regulated humoral immunity factors gene expression. Diets supplemented with 0.8% SMP significantly up-regulated the total number of hemocytes (THC), T-SOD and PO activity, improved the survival of crayfish after Citrobacter freundii infection. This study suggested that SMP could improve the cellular and humoral immunity of P. clarkii. Furthermore, this finding supplied a molecular foundation for further comprehending the immunopotentiator effects of plant polysaccharides in crustaceans.

Effects of air exposure stress on crustaceans: Histopathological changes, antioxidant and immunity of the red swamp crayfish Procambarus clarkii

Air exposure stress may result in oxidative damage and ultimately disease or death in crustaceans. Using the Procambarus clarkia, one of the main commercial aquaculture species in China, as a study model, the molecular mechanism including histopathological changes, antioxidant capacity and immunity response under the air exposure stress were firstly evaluated. Results showed that the surfaces of gill were wrinkled while the morphologies of the nuclei and mitochondria in the hepatopancreas were altered after 48 h of air exposure stress, and the damage of mitochondria was more serious after additional bacterial infection. Moreover, the activity of antioxidant enzymes increased at first and then decreased along with increasement of air exposure time. The concentration of malondialdehyde (MDA) in hepatopancreas was significantly increased under the air exposure stress, while the bacterial infection further aggravated such oxidative damage. The transcriptome analysis exhibited that the stress- and immunity-related genes in hepatopancreas altered when response to the air exposure stress. This study could help uncover the mechanisms of aerial exposure stress responses in Procambarus clarkii.

Transcriptome profiles of red swamp crayfish Procambarus clarkii hematopoietic tissue in response to WSSV challenge

The crayfish Procambarus clarkii could achieve a high cumulative mortality after WSSV infections. To better understand the immune response to WSSV in hematopoietic tissue, the present study investigated the immunological response of P. clarkii and analyzed the expression of some hematopoietic cytokines. After assembly, there was an average of 47,712,411 clean reads were obtained in control and treatment groups. A total of 35,945 unigenes were discovered with N50 length of 1554 bp. Under functional classification, enrichment, and pathway analysis using different database, there were about 257 differentially expressed genes (DEGs) identified, of which 139 were up-regulated and 118 were down-regulated. The GO function analysis of these DEGs were mostly participated in activation of immune response, complement activation, complement binding, negative regulation of humoral immune response and secretory granule membrane. Under KEGG analysis, these DEGs were involved in ECM-receptor interaction, HIF-1 signaling pathway, Glycolysis/Gluconeogenesis, Thyroid hormone signaling pathway and Glucagon signaling pathway. The real-time quantitative PCR (RT-qPCR) analysis of 9 selected genes confirmed the reliability of RNA-Seq results. The present research provide for the first time the transcriptomic profile of P. clarkii hematopoietic tissue in response to WSSV infection and reveals the astakines may play important roles in antiviral immune response. The results of the present study will further enrich the theoretical basis of the crayfish immune system and provide new ideas for disease prevention and control.

The immune response of fairy shrimp Streptocephalus sirindhornae against bacterial black disease by de novo transcriptome analysis

To enhance genomic resources and to understand the molecular immune mechanisms underlying the response of fairy shrimp (Streptocephalus sirindhornae) to pathogens, we first performed a comparative gene transcription analysis from Aeromonas hydrophila-immunized shrimp and from a control group through RNA sequencing. Meanwhile, the differentially expressed genes (DEGs) were investigated, and a total of 46,958,894 clean reads were obtained and then assembled into 73,297 unigenes with an average length of 993 bp and an N50 of 1,458 bp. Unigenes were annotated by comparison with the NR/NT/KO/SwissProt/PFAM/GO and KOG databases, and 28,198 unigenes (38.47%) were annotated in at least one database. After a bacterial challenge, 143 and 287 genes were identified as markedly up- or downregulated, respectively, and 345 were associated with 142 pathways, including the classic immune-related apoptosis, toll-like receptor and MAPK signaling pathways. Moreover, ten differently expressed immune-related genes were confirmed by using quantitative real-time PCR. This study characterized a gene expression pattern for normal and Aeromonas hydrophila-immunized S. sirindhornae for the first time and shed new light on its molecular mechanisms, thus enabling the future efforts of disease control programs for this valuable aquaculture species.

iTRAQ-based quantitative proteomic analysis of differentially expressed proteins in the hepatopancreas of Litopenaeus vannamei after WSSV infection

White spot syndrome virus (WSSV) is the most destructive virus among invertebrates. In this study, we analyzed the immune response after WSSV infection in Pacific white shrimp Litopenaeus vannamei using isobaric tags for relative and absolute quantitation (iTRAQ). We identified 325 differentially expressed proteins (DEPs) in the hepatopancreas of L. vannamei. Among them, 212 were up-regulated proteins, and several of them might be related to immunity (e.g. arginine kinase and peroxiredoxin). Of the 113 down-regulated proteins, some were related to immunity (e.g. cathepsin C and cathepsin L) and others to the antioxidant defense process (e.g. glutathione peroxidase and catalase). One down-regulated DEP (C7M84_014268) and 3 up-regulated DEPs (C7M84_003456, C7M84_020702, and C7M84_007135) were randomly selected and analyzed using parallel reaction monitoring. This study is an important step for a comprehensive understanding of the immune relationship between L. vannamei and WSSV and provides valuable information for the prevention of viral diseases in the crustacean aquaculture industry.

Structure, gene expression, and putative functions of crustacean heat shock proteins in innate immunity

Heat shock proteins (HSPs) are molecular chaperones with critical roles in the maintenance of cellular proteostasis. HSPs, which regulate protein folding and refolding, assembly, translocation, and degradation, are induced in response to physiological and environmental stressors. In recent years, HSPs have been recognized for their potential role in immunity; in particular, these proteins elicit a variety of immune responses to infection and modulate inflammation. This review focuses on delineating the structural and functional roles of crustacean HSPs in the innate immune response. Members of crustacean HSPs include high molecular weight HSPs (HSP90, HSP70, and HSP60) and small molecular weight HSPs (HSP21 and HSP10). The sequences and structures of these HSPs are highly conserved across various crustacean species, indicating strong evolutionary links among this group of organisms. The expression of HSP-encoding genes across different crustacean species is significantly upregulated upon exposure to a wide range of pathogens, emphasizing the important role of HSPs in the immune response. Functional studies of crustacean HSPs, particularly HSP70s, have demonstrated their involvement in the activation of several immune pathways, including those mediating anti-bacterial resistance and combating viral infections, upon heat exposure. The immunomodulatory role of HSPs indicates their potential use as an immunostimulant to enhance shrimp health for control of disease in aquaculture.

Transcriptomic analysis of Procambarus clarkii affected by "Black May" disease

Each year from April to May, high mortality rates are reported in red swamp crayfish (Procambarus clarkii) cultured in Jiangsu and other regions, in China, and this phenomenon has come to be known as “Black May” disease (BMD). Therefore, in order to investigate the possible causes of this disease, this study gathered BMD-affected P. clarkii samples and performed transcriptome analysis on hepatopancreas, gill, and muscle tissues. A total of 19,995,164, 149,212,804, and 222,053,848 clean reads were respectively obtained from the gills, muscle, and hepatopancreas of BMD-affected P. clarkii, and 114,024 unigenes were identified. The number of differentially expressed genes (DEGs) in gill, muscle, and hepatopancreas was 1703, 964, and 476, respectively. GO and KEGG enrichment analyses of the DEGs were then conducted. Based on KEGG pathway enrichment analysis, the most significantly differentially expressed pathways were mainly those involved with metabolism, human disease, and cellular processes. Further analysis of the significantly DEGs revealed that they were mainly related to the mitochondrial-mediated apoptosis pathway and that the expression of these DEGs was mostly down-regulated. Moreover, the expression of genes related to immune and metabolism-related pathways was also significantly down-regulated, and these significantly-inhibited pathways were the likely causes of P. clarkii death. Therefore, our results provide a basis for the identification of BMD causes.

Comparative transcriptomics reveals eyestalk ablation induced responses of the neuroendocrine-immune system in the Pacific white shrimp Litopenaeus vannamei

In decapod crustaceans, eyestalk ablation is widely used to expedite ovarian maturation and spawning because of the removal of a gonad inhibiting hormone produced by the X-organ sinus gland. However, eyestalk ablation also results in negative impacts on the immunocompetence of the eyestalk-ablated females. In the current study, we investigated the impact of eyestalk ablation on the transcriptomic responses of three major nervous organs of shrimp, including the eyestalk ganglion, brain and thoracic ganglion, using the Illumina Hiseq™ 4000 platform. A total of 48,249 unigenes with an average length of 1253 bp and a N50 value of 2482 bp were obtained. Following eyestalk ablation treatment, a total of 2,983, 6325 and 6575 unigenes were detected as differentially expressed (log₂Ratio >1 and FDR <0.05) from the eyestalk, brain and thoracic ganglia, respectively. Functional GO and KEGG analysis of these differential expression genes (DEGs) showed that these DEGs were associated with a wide variety of biological processes and pathways. The distribution of DEGs among three comparison groups was similar, and many DEGs were mapped to the phagosome pathway, indicating that eyestalk ablation triggers activation of the neuroendocrine-immune (NEI) system. Interestingly, several important pathways were uniquely enriched in the brain tissue, suggesting that the brain may play a crucial role in the NEI system in response to eyestalk ablation. This is the first report on the transcriptomic regulation of the nervous system in response to eyestalk ablation in L. vannamei. The genes and pathways identified in this study will help to elucidate the molecular mechanisms of neuroendocrine-immune responses to eyestalk ablation in penaeid shrimp.

New insights into the immune regulation and tissue repair of Litopenaeus vannamei during temperature fluctuation using TMT-based proteomics

To investigate shrimp immunoregulation and tissue self-repair mechanism during temperature fluctuation stage, Litopenaeus vannamei (L. vannamei) was treated under conditions of gradual cooling from an acclimation temperature (28 °C, C group) to 13 °C (T group) in 2 days with a cooling rate of 7.5 °C/d and then rewarmed to 28 °C (R group) with the same rate. Tandem mass tags (TMT) -based proteomics technology was used to investigate the protein abundance changes of intestine in L. vannamei during temperature fluctuation. The results showed that a total of 5796 proteins with function annotation were identified. Of which, the abundances of 1978 proteins (34%) decreased after cooling and then increased after rewarming, 1498 proteins (26%) increased during the whole stage, 1263 proteins (22%) increased after cooling and then decreased after rewarming and 1057 proteins (18%) decreased during the whole stage. Differentially expressed proteins such as C-lectin, NFκBIA and Caspase may contributed to the regulation of immunity and tissue repair of shrimp intestine during the temperature fluctuation stage. These findings contribute to the better understanding of shrimp' regulatory mechanism against adverse environment.

Acute BPA exposure-induced oxidative stress, depressed immune genes expression and damage of hepatopancreas in red swamp crayfish Procambarus clarkii

Bisphenol A is a typical endocrine disrupting chemicals (EDCs) and produce various toxic effects on animals due to its potential endocrine disruption, oxidative damage effect, mutagenic effect and hypomethylation. To study its effect on the immune system of crustaceans, the Procambarus clarkii were utilized to detect the immune related indicators after 225 μg/L BPA exposure for 1 week. Hepatopancreatic histology and ultrastructure analysis showed that the brush border disappeared, the lumen increased, and the connection between the hepatic tubules fade away in BPA treated group. BPA could significantly increase the level of ROS, inhibit the activities of antioxidant-related enzymes (SOD, POD, and CAT), and thereby cause the oxidative stress. The enzyme activities of AKP, ACP and lysozyme in hepatopancreas after BPA exposure were also depressed even after Aeromonas hydrophila infections. The relative expression profiles of immune-related genes after BPA exposure and bacterial infection showed suppressed trends of most selected genes. Under A. hydrophila infections, the cumulative mortality of 225 μg/L BPA-treated crayfish was significantly higher than other groups. All these results indicated that BPA exposure had adverse effects on the immune ability of P. clarkii. The present study will provide an important foundation for further understanding the effects of EDCs on crustacean immune functions.

Pathogenicity of Aeromonas hydrophila causing mass mortalities of Procambarus clarkia and its induced host immune response

Outbreaks of mass mortalities among cultured Procambarus clarkia occurred in a commercial hatchery during the spring of 2019 in Jiangsu province of China. Here, we exploit the pathogenicity and immune response of Aeromonas hydrophila (GPC1-2), which was isolated from diseased P. clarkia. Crayfish challenged showed similar pathological signs to the naturally diseased P. clarkia, lethal dose 50% (LD₅₀) of the strain GPC1-2 to P. clarkia was 3.8 × 10⁶ CFU/mL. Detection of virulence-associated genes by PCR indicated that the strain GPC1-2 carried hlyA, aerA, alt, ast, act, aha, ahp, ahpA, and ahpB. Histopathological analysis of hepatopancreas revealed that the hepatic tubule lumen and the gap between the hepatic tubules became larger, and the brush border disappeared in the P. clarkia infected by GPC1-2. Quantitive real-time PCR (qRT-PCR) was undertaken to measure mRNA expression levels for six immune-related genes in P. clarkia after A. hydrophila infection. The expression level of proPO, NOS, ALF1, TLR2, PX, and AST were detected in hemolymph, hepatopancreas, gill and intestine tissues, and clear transcriptional activation of these genes were observed in the infected individuals. These results revealed pathogenicity of A. hydrophila and its activation of host immune response, which will provide a scientific reference for the breeding and disease prevention in P. clarkia culture.

Sodium butyrate can improve intestinal integrity and immunity in juvenile Chinese mitten crab (Eriocheir sinensis) fed glycinin

Butyrate is a fermentation byproduct of gut microbiota and is susceptible to chronic oxidative stress. This study investigates the mitigative effects of sodium butyrate (SBT) on growth inhibition and intestinal damage induced by glycinin in juvenile Chinese mitten crab (Eriocheir sinensis). All four experimental diets containing 80 g/kg glycinin were formulated with 0, 10, 20 and 40 g/kg SBT respectively. There was no glycinin or SBT in the control diet. Juvenile crabs (0.33 ± 0.01g) were respectively fed with these five diets for eight weeks. The diets with 10 and 20 g/kg SBT significantly improved the survival and weight gain of the crabs compared with those in the 0 g/kg SBT group, and showed no difference with the control group. The crabs fed diets containing glycinin without SBT had lower glutathione and glutathione peroxidase activities but higher malondialdehyde in the intestine than those in the control group. Moreover, dietary glycinin decreased the lysozyme and phenoloxidase activities and improved the level of histamine in the intestine compared with the control group, while the supplementation of SBT counteracted these negative effects. The addition of SBT could also restore the impaired immunity and morphological structure of the intestine. Dietary SBT could increase the mRNA expression of antimicrobial peptides genes (anti-lipopolysaccharide factor 1 and 2) and decrease the content of pro-inflammatory factor TNF-α. The SBT could restore the intestinal microbial community disorganized by glycinin. The abundance of pathogenic bacteria (Aeromonas, Vibrio and Pseudomonas) decreased significantly and the potential probiotic bacteria (Bacillus, Lactobacillus, Chitinibacter and Dysgonomonas) increased significantly in the 10 g/kg SBT group. This study suggests that sodium butyrate supplementation can mitigate the negative effects induced by glycinin such as growth inhibition, intestinal inflammation and reduction of beneficial flora in the gut.

Transcriptome analysis of Procambarus clarkii infected with infectious hypodermal and haematopoietic necrosis virus

Infectious hypodermal and haematopoietic necrosis virus (IHHNV) is a major viral pathogen in cultured penaeid shrimp. IHHNV has many hosts, mainly including crustaceans. It has recently been reported that Procambarus clarkii can be infected by IHHNV. In the present study, we studied the hepatopancreas of P. clarkii by transcriptome high-throughput sequencing to analyze the response of P. clarkii to IHHNV infection. After de novo assembly, there were 400,340,760 clean reads. A total of 237 differentially expressed genes (DEGs) were obtained, including 77 significantly up-regulated unigenes and 160 significantly down-regulated ones. The expression levels of 12 immune-related DEGs were validated by qRT-PCR, substantiating the reliability of RNA-Seq results. The enrichment analysis of DEGs showed that the immune-related pathways were closely related to apoptosis and phagocytosis. Moreover, a large number of pathways related to metabolic function were down-regulated, suggesting that IHHNV infection might affect the growth of P. clarkii.

Characterization and expression analysis of immune-related genes in the red swamp crayfish, Procambarus clarkii in response to lipopolysaccharide challenge

To learn more about red swamp crayfish related genes in response to bacterial infections, we investigated immune-related genes induced by lipopolysaccharide (LPS) in the hepatopancreas using high-throughput sequencing method. In present the study, a total of 55,107 unigenes were identified, with an average length of 678 bp. A total of 2215 differentially expressed genes (DEGs) were found, including 669 up-regulated genes and 1546 down-regulated genes. The result of Gene ontology (GO) analysis revealed that 3017 DEGs were enriched in 19 biological process subcategories, 17 cellular component subcategories and 15 molecular function subcategories. The top 20 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways showed that "ribosome" was the most abundant group, which had 34 DEGs. KEGG enrichment analysis identified several immune response pathways. Real-time quantitative reverse transcription-PCR (qRT-PCR) results exhibited that several immune responsive genes were greatly up-regulated following LPS stimulation as observed in the results of high-throughput sequencing. Overall, this study provides new insight into the immune defense mechanisms of P. clarkii against LPS infection.

Characterization of a classical 2-cysteine peroxiredoxin1 gene from Chinese soft-shelled turtle Pelodiscus sinensis with its potent antioxidant properties and putative immune response

Peroxiredoxin family members could function in host defense against oxidative stress, and modulate immune response. In the present study, a 2-cysteine peroxiredoxin gene named PsPrx1 was isolated from Chinese soft-shelled turtle Pelodiscus sinensis. The PsPrx1 cDNA was composed of 1130 bp, consisted of 199 amino acid residues and included a Redoxin and AphC-TSA domain. As detected by qPCR, PsPrx1 was ubiquitously expressed in the examined tissues with the higher levels in liver and spleen. Upon the immune challenge with A. jandaei bacteria and oxidative stress with ammonia pressure, both mRNA and protein expression level in liver could be significantly enhanced. The results of immunohistochemical examinations showed PsPrx1 was mainly distributed at the junction between the hepatic cells. The general functional properties of PsPrx1 were confirmed using purified rPsPrx1 protein. From the results, rPsPrx1 protein was confirmed to exhibit antioxidant activity and antibacterial properties. The potential for scavenging extracellular H₂O₂ was evidenced by the purified rPsPrx1 protein in vitro system. In the mixed-function oxidase assay, rPsPrx1 also exhibited a dose-dependent inhibition of DNA damage. These results suggest that rPsPrx1 was implicated defense against microbial pathogens and oxidants, and would provide important information to further understand the functional mechanism of Prx1 in P. sinensis immunity.

Full-length transcriptome analysis of Litopenaeus vannamei reveals transcript variants involved in the innate immune system

To better understand the immune system of shrimp, this study combined PacBio isoform sequencing (Iso-Seq) and Illumina paired-end short reads sequencing methods to discover full-length immune-related molecules of the Pacific white shrimp, Litopenaeus vannamei. A total of 72,648 nonredundant full-length transcripts (unigenes) were generated with an average length of 2545 bp from five main tissues, including the hepatopancreas, cardiac stomach, heart, muscle, and pyloric stomach. These unigenes exhibited a high annotation rate (62,164, 85.57%) when compared against NR, NT, Swiss-Prot, Pfam, GO, KEGG and COG databases. A total of 7544 putative long noncoding RNAs (lncRNAs) were detected and 1164 nonredundant full-length transcripts (449 UniTransModels) participated in the alternative splicing (AS) events. Importantly, a total of 5279 nonredundant full-length unigenes were successfully identified, which were involved in the innate immune system, including 9 immune-related processes, 19 immune-related pathways and 10 other immune-related systems. We also found wide transcript variants, which increased the number and function complexity of immune molecules; for example, toll-like receptors (TLRs) and interferon regulatory factors (IRFs). The 480 differentially expressed genes (DEGs) were significantly higher or tissue-specific expression patterns in the hepatopancreas compared with that in other four tested tissues (FDR <0.05). Furthermore, the expression levels of six selected immune-related DEGs and putative IRFs were validated using real-time PCR technology, substantiating the reliability of the PacBio Iso-seq results. In conclusion, our results provide new genetic resources of long-read full-length transcripts data and information for identifying immune-related genes, which are an invaluable transcriptomic resource as genomic reference, especially for further exploration of the innate immune and defense mechanisms of shrimp.

A novel ferritin gene from Procambarus clarkii involved in the immune defense against Aeromonas hydrophila infection and inhibits WSSV replication

Ferritin is a protein related to the storage of iron and widely distributed in animals. It participates in many biological process, including antioxidation, cell activation, angiogenesis, regulation of iron metabolic balance and immune defense. In the present study, a novel ferritin gene was identified from red swamp crayfish Procambarus clarkii, with a cDNA sequence encoding a predicted 221 amino-acid residues. The ferritin protein contains a 19-residue signal peptide and 145-residue classic ferritin domain. The NJ phylogenetic analysis showed PcFer clustered with other crustacean peptides. The recombinant PcFer protein was produced and purified in E. coli, and the anti-rabbit polyclonal antibody was obtained. The rPcFer exhibited iron binding activity at a dose-dependent effect. The qPCR and western blot analysis revealed that PcFer was highly expressed in hemocytes, hepatopancreas, and gills. After challenged with WSSV and Aeromonas hydrophila, the mRNA and protein expression patterns of PcFer were significantly up-regulated in hemocytes and hepatopancreas. dsRNA interfering technique was utilized to silence the expression of PcFer gene. The WSSV copy number in PcFer silenced shrimp was much higher than that in the control group. The present study indicated that PcFer was involved in the immune defense against WSSV and Aeromonas hydrophila, and might inhibit WSSV replication in P. clarkii. These results will provide important data support for further study of the functional role of the ferritin gene.