Pathogenicity of non-O1/ O139 Vibrio cholerae and its induced immune response in Macrobrachium rosenbergii

Functional role of rpoN in regulating the virulence of non-O1/O139 Vibrio cholerae

Non-O1/O139 Vibrio cholerae is widely distributed in brackish and estuarine ecosystems, which can infect many aquatic animals. RpoN, an alternative sigma factor, plays a critical role in regulating cell functions such as motility, quorum sensing, and virulence. However, the function of rpoN in non-O1/O139 V. cholerae has rarely been reported. In the present study, we constructed the deletion mutant ΔrpoN of non-O1/O139 V. cholerae GXFL1-4 using recombination technology and investigated the function of rpoN through transcriptomic and phenotypic analyses. RNA-seq results showed that many major virulence-related genes were down-regulated in the ΔrpoN mutant, including the type VI secretion system (tssJ, tssA, tagO, tssG), type IV pilus assembly proteins (pilM, pilB), biofilm formation genes (vpsC, cheC), and hemolysin-related genes (hlyD, hlyD-PA). Additionally, phenotypic assays showed that the growth and motility of the ΔrpoN had no apparent change. The deletion of rpoN in non-O1/O139 V. cholerae led to decreased biofilm formation and reduced hemolytic activity. Furthermore, artificial infection tests showed that the virulence of the ΔrpoN mutant toward Macrobrachium rosenbergii was decreased. Our study provides essential insights into the regulatory function of rpoN, revealing that rpoN is a key determinant of virulence regulation in non-O1/O139 V. cholerae.

Pathogenicity of Citrobacter freundii Causing Mass Mortalities of Macrobrachium rosenbergii and Its Induced Host Immune Response

Citrobacter freundii is an opportunistic pathogen of freshwater aquatic animals, which severely restricts the sustainable development of the aquaculture industry. In this study, a dominant strain, named FSNM-1, was isolated from the hepatopancreas of diseased Macrobrachium rosenbergii. This strain was identified as C. freundii based on a comprehensive analysis of its morphological, physiological, and biochemical features and molecular identification. Challenge experiments were conducted to assess the pathogenicity of C. freundii to M. rosenbergii. The results showed that the FSNM-1 strain had high virulence to M. rosenbergii with a median lethal dose (LD50) of 1.1 × 106 CFU/mL. Histopathological analysis revealed that C. freundii infection caused different degrees of inflammation in the hepatopancreas, gills, and intestines of M. rosenbergii. The detection of virulence-related genes revealed that the FSNM-1 strain carried colonization factor antigen (cfa1, cfa2), ureases (ureG, ureF, ureD, ureE), and outer membrane protein (ompX), and virulence factor detection showed that the FSNM-1 strain had lecithinase, amylase, lipase, gelatinase, and hemolysin activities but did not produce protease and DNase activities. To investigate the immune response of M. rosenbergii to C. freundii, the expression levels of ALF3, MyD88, SOD, proPO, TRAF6, and TNF immune-related genes were monitored at different points of time in the hepatopancreas, gills, intestines, and hemocytes of M. rosenbergii after infection. The results demonstrated a significant upregulation in the expression levels of the ALF3, MyD88, SOD, proPO, TRAF6, and TNF genes in M. rosenbergii at the early stage of C. freundii infection. This study highlights C. freundii as a major pathogen causing mass mortality in M. rosenbergii and provides valuable insights into its virulence mechanisms and the host's immune response.

Complete genome sequence of Vibrio cholerae LK-18 isolated from tail-rotted Procambarus clarkii

Isolated from tail-rotted Procambarus clarkii, the pathogenic bacterium Vibrio cholerae LK-18 features two circular chromosomes: chromosome I (2,895,335 bp) and chromosome II (1,175,190 bp). The genome includes 3,522 open reading frames, 100 tRNA genes, and 31 rRNA genes, and it harbors the Vibrio cholera cytolysin and chitinase genes.

Isolation and the pathogenicity characterization of Decapod iridescent virus 1 (DIV1) from diseased Macrobrachium nipponense and its activation on host immune response

The high morbidity and mortality of Macrobrachium nipponense occurred in several farms in China, with cardinal symptoms of slow swimming, loss of appetite, empty of intestine, reddening of the hepatopancreas and gills. The pathogen has been confirmed as Decapod Iridescent Virus 1 (DIV1), namely DIV1-mn, by molecular epidemiology, histopathological examination, TEM observation, challenge experiment, and viral load detection. Histopathological analysis showed severe damage in hepatopancreas and gills of diseased prawns, exhibited few eosinophilic inclusions and pyknosis, and TEM of diseased prawns revealed that icosahedral virus particles existed in hepatopancreas and gill, which confirmed the disease of the farmed prawns caused by the DIV1 infection. Besides, challenge tests showed LD50 of DIV1 to M. nipponense was determined to be 2.14 × 104 copies/mL, and real-time PCR revealed that M. nipponense had a very high DIV1 load in the hemocytes, gills and hepatopancreas after infection. Furthermore, qRT-PCR was undertaken to investigated the expression of six immune-related genes in DIV1-infected M. nipponense after different time points, and the results revealed UCHL3, Relish, Gly-Cru2, CTL, MyD88 and Hemocyanin were significantly up-regulated in hemocytes, gills and hepatopancreas, which revealed various expression patterns in response to DIV1 infection. This study revealed that DIV1 infection is responsible for the mass mortality of M. nipponense, one of the important crustacean species, indicating its high susceptibility to DIV1. Moreover, this study will contribute to exploring the interaction between the host and DIV1 infection, specifically in terms of understanding how M. nipponense recognizes and eliminates the invading of DIV1.

Whole Genome Analysis of a Non-O1, Non-O139 Vibrio cholerae Detected from Human Blood in China

Non-O1, non-O139 Vibrio cholerae (NOVC) can cause cholera-like diarrhea, but it rarely causes extraintestinal infection, so it is easily overlooked. In this report, we present a case of NOVC detected through blood culture in a 58-year-old male patient with cirrhosis, resulting in severe infection. The patient had been diagnosed with cirrhosis seven years prior and was admitted to the hospital due to abdominal distension and gastrointestinal bleeding. Gram-negative bacilli were isolated from blood cultures and identified as V. cholerae using matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF MS) and average nucleotide identity (ANI). Moreover, the serum agglutination test showed that the strain was non-O1/non-O139. Further whole genome sequencing and analysis of the strain showed that the strain mainly carried virulence genes tox R, RTX, hly A, T3SS/T6SS, but no resistant genes such as sulII, dfrA1, strB were detected. It provides information for the study of the pathogenic mechanism and drug resistance mechanism of V. cholerae. The patient had severe symptoms and a poor prognosis, indicating that although the NOVC strain infected in this patient had few virulence genes, it was not weak in pathogenicity. It may be caused by the effect of some virulence genes, which should be paid attention to.

Transcriptome analysis of Macrobrachium rosenbergii hemocytes in response to Staphylococcus aureus infection

The aquaculture industry has suffered significant financial losses as a result of disease outbreaks. In particular, disease outbreaks have become a major problem that can seriously affect the sustainable development of the Macrobrachium rosenbergii aquaculture industry. It is crucial to determine the defense mechanism of the host after pathogenic invasion in order to provide effective defense measures after disease outbreaks. Shrimp, like other invertebrates, primarily depend on their innate immune systems to defend against pathogens, and recognize and resist pathogens through humoral and cellular immune responses. In this investigation, we used RNA-seq technology to investigate the transcriptome of hemocytes from M. rosenbergii induced by Staphylococcus aureus. Our main targets were immune pathways and genes related to innate immunity. RNA-seq identified 209,069 and 204,775 unigenes in the control and experimental groups, respectively. In addition, we identified 547 and 1734 differentially expressed genes (DEGs) following S. aureus challenge after 6 and 12 h (h), respectively. GO and KEGG enrichment analysis revealed that the DEGs were significantly enriched in several biological signalling pathways, including NOD-like receptor, PI3K-Akt, Toll and Imd, IL-17, TGF-beta, RIG-I-like receptor, cAMP, apoptosis, and C-type lectin receptor. Sixteen DEGs were chosen at random for qPCR verification; these results concurred with those from sequencing. Our findings revealed that immune-related genes play an important role in antibacterial activities and have specific functions for gram-positive bacteria. These results provide more data for the prevention of M. rosenbergii diseases and offer a basis for the better prevention of diseases.

Transcriptome-wide identification and characterization of the Macrobrachium rosenbergii microRNAs potentially related to immunity against non-O1 Vibrio cholerae infection

Non-O1 Vibrio cholerae, a member of the Vibrio family, could cause gastrointestinal infection of Macrobrachium rosenbergii and result in significant economic losses. However, few studies on microRNA immunity related to non-O1 V. cholerae infection of M. rosenbergii. The aim of this study was to elucidate the mechanism of miRNA in the potential immune response of M. rosenbergii. to non-O1 V. cholerae MSVC-GY01 infection by transcriptome sequencing. Following quality screening, the control group received 10, 616, 712 clean reads, whereas the infected group received 9,727,616. The miRNA sequences in the two samples are extremely consistent and have a length of roughly 23 nt. In all, 871 known miRNAs were discovered, with 279 differentially expressed miRNAs (DEMs). Meanwhile, 62 novel miRNAs were predicted, including 43 DEMs. In order to understand the immune-related biological functions of DEMs, target genes were predicted. Pathway function annotation analysis showed that non-O1 V. cholerae affected the NOD-like receptor signaling pathway, RIG-I-like receptor signaling pathway, and Toll-like receptor signaling pathway, suggesting that miRNAs in the hepatopancreas play a key role in immune responses to pattern recognition receptors. Twelve DEMs were randomly selected for Quantitative Real-time PCR (qRT-PCR). Overall, the expression trends of qRT-PCR were consistent with the sequencing results. These findings corroborate the immunomodulatory function of miRNA in M. rosenbergii against non-O1 V. cholerae infection and provide guidance for the prevention and treatment of related illnesses.

Identification and genomic analysis of a Vibrio cholerae strain isolated from a patient with bloodstream infection

Vibrio cholerae is a bacterium ubiquitous in aquatic environments which can cause widespread infection worldwide. V. cholerae gradually became a rare species of bacteria in clinical microbiology laboratories with the control of the cholera epidemic. In this study, we isolated a V. cholerae strain, named VCHL017, from the blood of an elderly patient without gastrointestinal symptoms. The patient had a history of hookworm infection and multiple myeloma. Furthermore, she was immunocompromised, and received long-term chemotherapy and antimicrobial agents. VCHL017 was inoculated on blood agar and thiosulfate citrate bile salt sucrose plates (TCBS) to observe morphological characteristics. Then this isolate was identified by matrix-assisted laser desorption/ionization time-of-flight spectrometry (MALDI-TOF MS). The minimum inhibitory concentrations (MICs) for cefazolin, ceftazidime, cefepime, meropenem, tetracycline, ciprofloxacin, chloramphenicol, and gentamicin of VCHL017 were determined by the microbroth dilution method. PCR and serum agglutination tests were used to determine whether the serogroups of the isolate belonged to the O1/O139 and cholera toxin encoding genes. Finally, the genomic features and phylogeny of VCHL017 were analyzed by whole genome sequencing (WGS). VCHL017 was a non-O1/O139 V cholerae strain that did not carry the ctxA gene. Antimicrobial susceptibility tests revealed that VCHL017 was susceptive to chloramphenicol and tetracycline. Although it did not carry the genes encoding the cholera toxin, WGS indicated that VCHL017 carried a variety of other virulence factors. By calculating the average nucleotide identity (ANI), we precisely identified the species of VCHL017 as V. cholerae. There are also A171S and A202S missense mutations in gyrA of VCHL017. The phylogenetic analysis indicated that VCHL017 was closely related to V. cholerae strains isolated from aquatic environments. Our results suggest that continuous monitoring is necessary for non-O1/O139 V cholerae strains isolated from outside the digestive tract, which could be pathogenic through multiple virulence factors.

Effects of high fat in the diet on growth, antioxidant, immunity and fat deposition of Macrobrachium rosenbergii post-larvae

Lipids are essential nutrients for organisms, and high-fat feeds for shrimp may cause oxidative stress. This study evaluated the effects of feeding high fat in the diet on the growth, antioxidant, immunity, and liver fat accumulation of Macrobrachium rosenbergii post-larvae. Five groups with an initial body weight of 0.0084 ± 0.001 g were fed five isonitrogenous and isoenergetic diets (47.01% crude protein and 18.40 kJ/g gross energy) containing 8%, 10%, 12%, 14% and 16% (named L8, L10, L12, L14 and L16) lipid for 8 weeks, respectively. The results showed that the weight gain rate (WGR) and specific growth rate (SGR) of L8 group were significantly higher than those of L10, L12, L14 and L16 group (P < 0.05), and the feed coefficient (FCR) of L8 group was significantly lower than that of other groups (P < 0.05). With the increase of dietary fat level, the content of MDA and the activity of SOD increased significantly, and the activities of T-AOC and CAT decreased significantly (P < 0.05). H&E staining clearly revealed the occurrence of hepatocyte swelling, hepatocyte vacuolization and nucleus displacement to the peripheral cell vacuolization in the L16 group, and hepatic lipid accumulation was further observed in the L14 and L16 group by Oil red O staining. In addition, high-fat diet significantly upregulated the expression of Dorsal, Relish and IκBα mRNA, and also upregulated the expression of fat synthesis-related genes FAS, ACC, DGAT and fat transport-related gene FABP (P < 0.05), and significantly downregulated the expression of fat metabolism-related genes AMPK and CPT-1 (P < 0.05) compared to that of the L8 group. In conclusion, this study showed that feeding a high-fat diet could induce oxidative stress, inhibit growth performance, alter antioxidant capacity, cause hepatic fat deposition and affect the immune system of M. rosenbergii post-larvae.

A Novel Disease (Water Bubble Disease) of the Giant Freshwater Prawn Macrobrachium rosenbergii Caused by Citrobacter freundii: Antibiotic Treatment and Effects on the Antioxidant Enzyme Activity and Immune Responses

The giant freshwater prawn, Macrobrachium rosenbergii, is an important and economical aquaculture species widely farmed in tropical and subtropical areas of the world. A new disease, "water bubble disease (WBD)", has emerged and resulted in a large loss of M. rosenbergii cultured in China. A water bubble with a diameter of about 7 mm under the carapace represents the main clinical sign of diseased prawns. In the present study, Citrobacter freundii was isolated and identified from the water bubble. The optimum temperature, pH, and salinity of the C. freundii were 32 °C, 6, and 1%, respectively. A challenging experiment showed that C. freundii caused the same typical signs of WBD in prawns. Median lethal dose of the C. freundii to prawn was 104.94 CFU/g. According to the antibiogram tests of C. freundii, florfenicol and ofloxacin were selected to evaluate their therapeutic effects against C. freundii in prawn. After the challenge with C. freundii, 86.67% and 72.22% survival of protective effects against C. freundii were evaluated in the oral florfenicol pellets and oral ofloxacin pellets feding prawns, respectively, whereas the mortality of prawns without fed antibiotics was 93%. After antibiotic treatment and C. freundii infection, the activities of superoxide dismutase (SOD), catalase (CAT), glutathione peroxidase (GPx), glutathione S-transferase (GST), malondialdehyde (MDA), acid phosphatase (ACP), alkaline phosphatase (ALP), and lysozyme (LZM) in the hemolymph and hepatopancreas of the prawns and the immune-related gene expression levels of Cu/Zn-SOD, CAT, GPx, GST, LZM, ACP, anti-lipopolysaccharide factor, crustin, cyclophilin A, and C-type lectin in hepatopancreas were all significantly changed, indicating that innate immune responses were induced by C. freundii. These results can be beneficial for the prevention and control of C. freundii in prawns.

Complete Genome Analysis of Highly Pathogenic Non-O1/O139 Vibrio cholerae Isolated From Macrobrachium rosenbergii Reveals Pathogenicity and Antibiotic Resistance-Related Genes

Non-O1/O139 Vibrio cholerae is a highly virulent pathogen that causes mass mortalities of various aquatic animals. In the present study, we sequenced the whole genome of non-O1/O139 V. cholerae GXFL1-4, isolated from Macrobrachium rosenbergii, to reveal the pathogenicity and antibiotic resistance. The result showed its genome contained two circular chromosomes and one plasmid with a total size of 4,282,243 bp, which harbored 3,869 coding genes. Among them, 3,047, 2,659, and 3,661 genes were annotated in the Clusters of Orthologous Genes (COG), Gene Ontology (GO), and Kyoto Encyclopedia of Genes and Genomes (KEGG), respectively. In addition, 372 potential virulence genes were predicted based on the Virulence Factor Database (VFDB) database, such as type II, III, IV, and VI secretion systems related genes, flagella genes, and pilus formation or motility-related genes. Blast results in the Comprehensive Antibiotic Resistance Database (CARD) database showed that the strain contained 148 antibiotic resistance-related genes belonging to 27 categories, such as efflux pump complex antibiotic resistance genes and antibiotic resistance gene cluster genes. The Pathogen-Host Interaction (PHI) database annotated 320 genes related to pathogen-host interaction, such as T3SS, virulence regulatory factors, transcriptional regulators, and two-component response regulator related genes. The whole-genome analysis suggested that the pathogenic non-O1/O139 V. cholerae strain GXFL1-4 might have a complex molecular mechanism of pathogenicity and antibiotic resistance. This study provides a wealth of information about non-O1/O139 V. cholerae genes related to its pathogenicity and drug resistance and will facilitate the understanding of its pathogenesis as well as the development of prevention and treatment strategies for the pathogen.

Survival Characteristics and Transcriptomic Analyses Reveal the Adaptive Response of the Aquatic Pathogen Non-O1/O139 Vibrio cholerae to Starvation Stress

Non-O1/O139 Vibrio cholerae is a pathogen of various aquatic organisms but requires major self-regulation to overcome environmental stress in the aquatic environment. However, its survival strategies under environmental stress are not well understood. The objective of this study was to describe the survival characteristics and changes in expression of stress resistance-related genes of non-O1/O139 V. cholerae after 6 months of starvation at room temperature. The results demonstrated that starved cells were still viable, exhibited shortened rods and shrinking surface, and maintained virulence to Macrobrachium rosenbergii. To investigate the changes in gene expression in non-O1/O139 V. cholerae under starvation stress, especially those involved in stress resistance, transcriptome profiles of starved and wild-type cells were determined. The differentially expressed genes (DEGs) in starved cells were identified, including 191 upregulated genes and 180 downregulated genes. Among these DEGs, the well-known stress resistance-related genes were upregulated significantly, including rpoS, rpoD, rpoN, rpoE, uspA, uspC, cspD, hslJ, etc. Gene Ontology (GO) analysis of the DEGs demonstrated that environmental adaptation-related categories, such as response to stimulus and signal transduction, were upregulated significantly in the starved cells, while cell motility was downregulated significantly. These DEGs were also enriched into 54 KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways, including biofilm formation, two-component system, quorum sensing, flagellar assembly, bacterial chemotaxis stress resistance-related pathways, etc. The potential existence of long-starved non-O1/O139 V. cholerae bacteria in the aquatic environment may raise new concerns about this devastating pathogen in aquaculture.

IMPORTANCE Non-O1/O139 V. cholerae is a causal agent of vibriosis that can be subject to nutrient insufficiency and cause high rates of mortality in aquatic animals. However, its molecular mechanisms of survival in response to starvation stress have been investigated only partially. Here, we demonstrate that under starvation stress, non-O1/O139 V. cholerae can survive over the long term and cause disease by dwarfing of the cell structure, upregulation of a series of stress resistance-related genes, and downregulation of flagellum assembly-related genes. This knowledge can help the development of intervention strategies to control non-O1/O139 V. cholerae infection in aquaculture.

Role of peroxinectin in the antibacterial immune response of the Chinese mitten crab, Eriocheir sinensis

To elucidate the antibacterial role of peroxinectin (referred to as PXN) and its molecular mechanism in Chinese mitten crab Eriocheir sinensis, we analyzed the bacterial binding and removal of the peroxinectin recombinant protein in vitro and the interaction of peroxinectin with integrin and CuZn-SOD through GST-pulldown and bimolecular fluorescence complementation methods. Concurrently, the effect of peroxinectin interference on the expression of other immune-related genes was studied using RNA interference. The results showed that the recombinant peroxinectin protein could bind to Bacillus subtilis, Staphylococcus aureus, Aeromonas hydrophila, and Vibrio parahaemolyticus with different affinities in vitro and could eliminate Vibrio parahaemolyticus in vivo. The findings also indicated that peroxinectin could establish interactions with integrin and CuZn-SOD in vitro. Furthermore, 48 h after the injection of the peroxinectin gene siRNA in vivo, the expression of peroxinectin mRNA decreased significantly (P < 0.05), integrin mRNA expression decreased by 16.8%, and CuZn-SOD mRNA expression decreased by 62.84% (P < 0.01). The expression levels of Dorsal, GPx, GST, PPAF, and Relish (P < 0.01), as well as that of lectin (P < 0.001) were significantly decreased. When peroxinectin siRNA was injected in vivo for 48 h and Aeromonas hydrophila was injected into mitten crabs, the expression of immune-related genes significantly increased. All data indicate that the recombinant peroxinectin protein in Chinese mitten crabs can recognize and bind different bacteria and promote the elimination of Vibrio parahaemolyticus from the body. Furthermore, peroxinectin may establish interactions with integrin and CuZn-SOD to activate the expression of related immune genes to elicit responses to bacterial infections and achieve immune protection.

One health pathogen surveillance demonstrated the dissemination of gut pathogens within the two coastal regions associated with intensive farming

BACKGROUND

Intensive aquaculture farming has caused significant degradation of coastal wetlands and has been proposed as a reservoir for pathogenic Vibrio spp.

RESULTS

Gut pathogens including Vibrio spp., Salmonella spp., and Klebsiella spp. were isolated from bird feces, shrimp and wetland water in two typical coastal regions of China in 2015 and 2017 and were subsequently subjected to whole-genome sequencing. Meanwhile, local patient isolates were also selected to confirm the epidemiological links. Bacterial community composition analyses of the sediments that were sampled in 2015 and 2017 were conducted by the hypervariable region 4 of the 16S rRNA gene. Together with the local clinical isolates, we observed highly related Vibrio isolates from waterbirds, wetlands and shrimp. Phylogenetic genome comparisons also demonstrated that sequence types ST3 and ST2414 Vibrio parahaemolyticus isolates obtained from aquatic animals were clonally related to patient isolates. Likewise, three Salmonella typhimurium isolates were also genomically related to one clinical strain. The results showed that farming activities significantly altered the community composition and resulted in the emergence of several pathogens, including Acinetobacter, Mycobacterium and Legionella.

CONCLUSIONS

In conclusion, our results demonstrated that intensive shrimp farming in wetlands has two devastating impacts: pathogen dissemination from aquatic animals into migratory birds and transmission of foodborne pathogens into local communities.

Aeromonas hydrophila associated with red spot disease in Macrobrachium nipponense and host immune-related gene expression profiles

In September 2018, a serious disease causing high mortality with red spot syndrome occurred in a Macrobrachium nipponense aquaculture farm in Jintan County, Jiangsu Province, China. In this study, a pathogenic isolate 5-S3 was isolated from diseased M. nipponense and was identified as Aeromonas hydrophila by phenotypically and molecularly. The pathogenicity of the isolate 5-S3 to M. nipponense was determined by challenge experiments. Results of artificial challenge showed A. hydrophila was pathogenic to M. nipponense, the LD₅₀ was 9.58 × 10⁴ CFU/mL, and histopathological analysis revealed that the hepatopancreas of infected M. nipponense exhibited obvious inflammatory responses to A. hydrophila infection. The isolate showed significant phenotypical activities such as the lecithinase, esterase, caseinase and hemolysin which are indicative of their virulence potential. Besides, virulence genes such as aerA, act, fla, ahpβ, alt, lip, eprCAI, hlyA, acg and gcaT were detected in the isolate 5-S3. Subsequently, the immune-related genes expression in M. nipponense were evaluated by quantitative real-time PCR (qRT-PCR), and the results showed that the expression levels of dorsal, relish, crustin1, crustin2, anti-lipopolysaccharide factors 1 (ALF1), anti-lipopolysaccharide factors 2 (ALF2), hemocyanin, i-lysozyme and prophenoloxidase were significantly up-regulated in hepatopancreas of M. nipponense after A. hydrophila infection, the stat, p38, crustin3, anti-lipopolysaccharide factors 3 (ALF3) genes had no significant change during the infection. The present results reveal that A. hydrophila was an etiological agent causing red spot syndrome and mass mortality of M. nipponense and the influence of A. hydrophila infection on host immune genes.

Transcriptome analysis and immune-related genes expression reveals the immune responses of Macrobrachium rosenbergii infected by Enterobacter cloacae

Macrobrachium rosenbergii is an important cultural species in China and other Southeast Asian countries. However, Enterobacter cloacae infection has caused a great economic loss in M. rosenbergii culture industry. The immune responses of M. rosenbergii to the E. cloacae infection is not fully characterized. To investigate the immune response of M. rosenbergii against E. cloacae, we performed transcriptome analysis of the M. rosenbergii hepatopancreas with and without E. cloacae infection using RNA-seq. After assembly and annotation, 29,731 high quality unigenes were obtained from RNA-seq data. Differential expression analysis revealed the existence of 2498 significantly differently expressed genes (DEGs) at 12 h post infection, with 1365 up-regulated and 1133 down-regulated genes. Among these DEGs, some well-known immune-related genes were up-regulated significantly, including C-type lectin 1, lectin 3, anti-lipopolysaccharide factor 2, Cu/Zn superoxide dismutase and heat shock protein 70. GO analysis demonstrated 24 biological process subcategories, 14 cellular component subcategories, and 12 molecular function subcategories that were enriched among these DEGs, and some DEGs were clustered into immune related subcategories such as immune system process, response to stimulus, biological adhesion, and antioxidant activity. These DEGs were enriched into 216 KEGG pathways including a core set of immune correlated pathways notably in phagosome and lysosome. In addition, 5 up-regulated and 5 down-regulated immune-related DEGs were selected for further validation by quantitative real-time PCR and the results showed consistence with the RNA-seq data. Additionally, the expression level of six selected immune-related genes (ALF2, CLEC1, LEC3, hemocyanin1, HSP70 and SOD) based on the transcriptomic data were monitored at different point of time in hepatopancreas, gill, hemolymph and intestine. Results revealed these immune-related genes were significantly up-regulated in different tissues from 6 to 24 h after E. cloacae infection. Overall, these results provided valuable information for further studying the immune response of M. rosenbergii against E. cloacae infection.