Transcriptome-wide identification and characterization of the Procambarus clarkii microRNAs potentially related to immunity against Spiroplasma eriocheiris infection

Comprehensive analysis of the mRNA and miRNA transcriptome implicated in the immune response of Procambarus clarkii to Spiroplasma eriocheiris

In recent years, the red swamp crayfish (Procambarus clarkii, P. clarkii) farming industry has suffered huge economic losses due to the pathogenic bacterium Spiroplasma eriocheiris (S. eriocheiris). To elucidate the immune response mechanism and identify hub immune genes as well as their associated microRNAs that regulate the host response of P. clarkii against S. eriocheiris infection, we conducted a comprehensive analysis on P. clarkii hemocyte mRNA and microRNA (miRNA) transcriptomes at different infection stages using third- and second-generation sequencing technologies. In full-length transcriptome functional annotation, 8155 unigenes were annotated, and 1168 potential new transcripts were predicted. In the mRNA transcriptome, a total of 3168 differentially expressed genes were identified at different infection stages, including 1492 upregulated and 1676 downregulated genes (duplicate genes excluded). Transcriptome analysis revealed 880 differentially expressed genes involved in multiple pathways and processes such as endocytosis, autophagy, lysosome, mTOR signaling, phagosome, and the Fanconi anemia pathway. Mfuzz analysis was employed to integrate and cluster the differential expression trends of genes across the three infection stages. In the miRNA transcriptome, 234 miRNAs and 966 predicted target genes were identified, with 86 differentially expressed miRNAs identified across the three time periods. A significant difference (P < 0.05) was observed for miRNAs including pcl-miR-146-3p, pcl-miR-74-3p, pcl-miR-225-5p, and pcl-miR-68-5p. These miRNAs are involved in multiple immune and autophagy-related pathways and have regulatory effects on immune genes including Vps26, lqf, and ERK-A. Based on the differentially expressed immune-related genes, we constructed a protein-protein interaction (PPI) network, which revealed the interactions among hub genes including Rac1, Akt1, Rho1, and Egfr. We also constructed a miRNA-gene interaction network in immune and autophagy-related processes, highlighting the potential regulatory effects of miRNAs including pcl-miR-183-5p, pcl-miR-146-3p, pcl-miR-176-5p, and pcl-miR-225-5p on proteins including LST8, SNAP29, Rab-7A, and ERK-A. To conclude, this study has identified hub immune genes and corresponding regulatory miRNAs in P. clarkii hemocytes in response to S. eriocheiris infection and explored the roles of these genes in selected pathways and processes. These findings are expected to provide further insights into the molecular mechanisms that confer resistance to S. eriocheiris infection in P. clarkii.

MiR-214_L-1R+4 regulate gossypol-induced immune response through MyD88-dependent signaling pathway in Cyprinus carpio

MicroRNAs (miRNAs) have been demonstrated to act as crucial modulators with considerable impacts on the immune system. Cottonseed meal is often used as a protein source in aqua feed, cottonseed meal contains gossypol, which is harmful to animals. However, there is a lack of research on the role of miRNAs in fish exposed to gossypol stress. To determine the regulatory effects of miRNAs on gossypol toxicity, Cyprinus carpio were given to oral administration of 20 mg/kg gossypol for 7 days, and the gossypol concentration in the tissues was tested. Then, we detected spleen index, histology, immune enzyme activities of fish induced by gossypol. The results of miRNA sequencing revealed 8 differentially expressed miRNAs in gossypol group, and miR-214_L-1R+4 was found involved in immune response induced by gossypol. The potential targets of miR-214_L-1R+4 were predicted, and found a putative miR-214_L-1R+4 binding site in the 3'UTR of MyD88a. Furthermore, dual-luciferase reporter assays displayed miR-214_L-1R+4 decreased MyD88a expression through binding to the 3'UTR of MyD88a. Moreover, miR-214_L-1R+4 antagomir were intraperitoneally administered to C. carpio, down-regulated miR-214_L-1R+4 could increase MyD88a expression, as well as inflammatory cytokines and anti-inflammatory cytokines expression. These findings revealed that miR-214_L-1R+4 via the MyD88-dependent signaling pathway modulate the immune response to gossypol in C. carpio spleen.

Effects of Dietary Melatonin on Antioxidant Capacity, Immune Defense, and Intestinal Microbiota in Red Swamp Crayfish (Procambarus clarkii)

The aim of this study was to investigate the effects of melatonin (MT) feed supplementation on the antioxidant capacity, immune defense, and intestinal flora in Procambarus clarkii (P. clarkii). Six groups of P. clarkii were fed test feeds containing different levels of MT: 0 mg/kg (control), 22.5, 41.2, 82.7, 165.1, and 329.2 mg/kg for a duration of 2 months. The specific growth rate, hepatosomatic index, and condition factor were recorded highest in the test group of shrimp fed an MT concentration of 165.1 mg/kg. Compared to the control group, the rate of apoptosis was lower in hepatopancreas cells of P. clarkii supplemented with high concentrations of MT. Analyses of antioxidant capacity and immune-response-related enzymes in the hepatopancreas indicated that dietary supplementation of MT significantly augmented both the antioxidant system and immune responses. Dietary MT supplementation significantly increased the expression levels of antioxidant-immunity-related genes and decreased the expression levels of genes linked to apoptosis. Dietary MT was associated with an elevation in the abundance of the Firmicutes and a reduction in the abundance of the Proteobacteria in the intestines; besides, resulting in an increase in the abundance of beneficial bacteria, such as Lactobacilli. The broken-line model indicated that the suitable MT concentration was 154.09-157.09 mg/kg. MT supplementation enhanced the growth performance of P. clarkii, exerting a positive influence on the intestinal microbiota, and bolstered both immune response and disease resistance. Thus, this study offered novel perspectives regarding the application of dietary MT supplementation within the aquaculture field.

MicroRNA sequencing analysis reveals immune responses in hepatopancreas of Fenneropenaeus penicillatus under white spot syndrome virus infection

White Spot Disease is one of the most harmful diseases of the red tail shrimp, which can cause devastating economic losses due to the highest mortality up to 100% within a few days. MicroRNAs (miRNAs) are large class of small noncoding RNAs with the ability to post-transcriptionally repress the translation of target mRNAs. MiRNAs are considered to have a significant role in the innate immune response of crustaceans, particularly in relation to antiviral defense mechanisms. Numerous crustacean miRNAs have been verified to be required in host immune defense against viral infection, however, till present, the miRNAs functions of F. penicillatus defense WSSV infection have not been studied yet. Here in this study, for the first time, miRNAs involved in the F. penicillatus immune defense against WSSV infection were identified using high-throughput sequencing platform. A total of 432 miRNAs were obtained including 402 conserved miRNAs and 30 novel predicted miRNAs. Comparative analysis between the WSSV-challenged group and the control group revealed differential expression of 159 microRNAs in response to WSSV infection. Among these, 48 were up-regulated and 111 were down-regulated. Ten candidate MicroRNAs associated with immune activities were randomly selected for qRT-PCR analysis, which confirming the expression profiling observed in the MicroRNA sequencing data. As a result, most differentially expressed miRNAs were down-regulated lead to increase the expression of various target genes that mediated immune reaction defense WSSV infection, including genes related to signal transduction, Complement and coagulation cascade, Phagocytosis, and Apoptosis. Furthermore, the genes expression of the key members in Toll and Imd signaling pathways and apoptotic signaling were mediated by microRNAs to activate host immune responses including apoptosis against WSSV infection. These results will help to understand molecular defense mechanism against WSSV infection in F. penicillatus and to develop an effective WSSV defensive strategy in shrimp farming.

Analysis of mRNA and MicroRNA Expression Profiles of Nervous Tissues and Reproductive Tissues in Male Procambarus clarkii After Silencing IAG

The insulin-like androgenic gland hormone gene (IAG), primarily expressed in the androgenic gland (AG), plays a crucial role in controlling male sex differentiation and maintaining male secondary sexual characteristics in decapods. In this study, we investigated the mRNA and microRNA expression profiles of male Procambarus clarkii to understand the transcriptomic regulatory mechanism of IAG after the injection of an efficient siRNA (GsiRNA) designed based on IAG. The results revealed that several differentially expressed genes were enriched in reproduction-related pathways, such as the wnt signaling pathway, MAPK signaling pathway, and GnRH signaling pathway. In the testis (Te), the injection of GsiRNA led to the up-regulation of many ovary-related genes and down-regulation of testis-related genes. Moreover, the brain (Br) and abdominal nerve cord (AN) appeared to be involved in the regulation of IAG, with numerous differentially expressed genes found in Br and AN. Notably, the expression of five neuropeptide genes, Crustacean hyperglycemic hormone, pigment-dispersing hormone, red pigment concentrating hormone precursor, corazonin, and gonadotropin-releasing hormone II receptor isoform X1 in Br/AN, was significantly changed. Additionally, three ovary-related miRNAs (miR-263a, miR-263b, miR-133) highly expressed in Te/AG showed significant up-regulation after GsiRNA injection. Furthermore, the long-term interference of GsiRNA was found to inhibit the development of male external sexual characteristics during the juvenile stage and delay it during the adult stage. This research provides valuable insights into the molecular regulatory mechanism and function of IAG in P. clarkii.

MicroRNA sequencing analysis reveals injury-induced immune responses of Scylla paramamosain against cheliped autotomy

During pond culture or intensive culture system of crabs (mainly Eriocheir sinensis, Portunus trituberculatus and Scylla paramamosain), high-density farming has typically contributed to a higher limb autotomy level in juvenile animals, especially in S. paramamosain which has a high level of cannibalism. Due to the high limb autotomy level, the survival and growth rates in S. paramamosain farming are restricted, which limit the growth of the mud crab farming industry. MicroRNAs (miRNAs) are small noncoding RNAs that regulate a series of biological processes including innate immune responses by post-transcriptional suppression of their target genes. MiRNAs are believed to be crucial for innate immune process of host wound healing. Many miRNAs have been verified to be required in host immune responses to repair wound and to defense pathogen after tissue damage. However, to our best knowledge, the miRNAs functions of crustacean innate immune reactions against injury induced by limb autotomy have not been studied yet. Here in this study, for the first time, miRNAs involved in the S. paramamosain immune reactions against injury induced by cheliped autotomy were obtained by high-throughput sequencing. A total of 575 miRNAs (518 known miRNAs and 57 novel predicted miRNAs) were obtained, of which 141 differentially expressed microRNAs (93 up-regulated microRNAs and 48 down-regulated microRNAs) were revealed to be modified against cheliped autotomy, and the qPCR results of randomly selected miRNAs confirmed the expression patterns in the miRNAs sequencing data. Numerous immune-related target genes associated with innate immune system were mediated by miRNAs to induce host humoral immune and cellular immune defense to minimize acute physical damage. Furthermore, the genes expression in hemolymph coagulation and melanization pathways, as well as Toll and Imd signaling pathways were mediated by miRNAs to activate host immune responses including melanization and antimicrobial peptides for rapid wound healing and killing invaded pathogens. These results will help to understand injury-induced immune responses in crabs and to develop an effective control strategy of autotomy rate in crabs farming.

Identification and characterization of immune-related microRNAs in hybrid snakehead(Channa maculata♀ × Channa argus♂)after treated by Echinacea purpurea (Linn.) Moench

Echinacea purpurea (Linn.) Moench (EP) is a globally popular herbal medicine, which showed effects on growth promotion, antioxidant and immunomodulatory activities in fish culture world widely. However, there are few studies about the effects on miRNAs by EP in fish. The hybrid snakehead fish (Channa maculate♀ × Channa argus ♂) was new important economic specie of freshwater aquaculture in China with high market value and demand while there were only a few reports about its miRNAs. To overview immune-related miRNAs of the hybrid snakehead fish and to further understand the immune regulating mechanism of EP, we herein constructed and analyzed three small RNA libraries of immune tissues including liver, spleen and head kidney of the fish with or without EP treatment via Illumina high-throughput sequencing technology. Results showed that EP can affect the immune activities of fish by the miRNA-regulated ways. Totally, 67 (47 up and 20 down) miRNAs in liver, 138 (55 up and 83 down) miRNAs in spleen, and 251 (15 up and 236 down) miRNAs in spleen were detected, as well as 30, 60, 139 kinds of immune-related miRNAs belonging to 22, 35 and 66 families of the three tissues respectively. The expressions of 8 immune-related miRNA family members were found in all the three tissues, including miR-10, miR-133, miR-22 and etc. Some miRNAs have been identified involved in the innate and adaptive immune responses, such as the miR-125, miR-138, and miR-181 family. Ten miRNA families with antioxidant target genes were also discovered, including miR-125, miR-1306, and miR-138, etc. Results from Gene Ontology (GO) and KEGG pathway analysis further confirmed there are a majority immune response targets of the miRNAs involved in the EP treatment process. Our study deepened understanding roles of miRNAs in fish immune system and provides new ideas for the study of immune mechanism of EP.

Acute toxicity, biochemical and transcriptomic analysis of Procambarus clarkii exposed to avermectin

BACKGROUND

Pesticides are extensively applied globally. Pesticide residues induce calamitous effects on the environment and untargeted organisms. Public concerns for the safety of freshwater organisms and the challenges posed by aquatic contaminants remain high. In the present study, the acute toxicity of avermectins (AVMs) to the crayfish, Procambarus clarkii was evaluated. We also evaluated the potential effects of AVM on the biochemical and transcriptomic status of the hepatopancreas and gastrointestinal tract in P. clarkii.

RESULTS

The 24, 48, 72, 96 h median lethal concentrations (LC₅₀ ) of AVM on crayfish were 2.626, 1.162, 0.723, 0.566 mg L^-1 , respectively. The crayfish were then exposed to 0.65 mg L^-1 of AVM for 96 h. AVM significantly altered biochemical parameters including AChE and CAT activities in the hepatopancreas, and AChE, SOD and Na + -K + -ATPase activities in the gastrointestinal tract at several time points. Furthermore, transcriptomic analysis identified 953 and 1851 differentially-expressed genes (DEGs) in the hepatopancreas and gastrointestinal tract, respectively. KEGG enrichment showed that the gene expression profiles of the hepatopancreas and gastrointestinal tract were distinct from each other. The DEGs in the hepatopancreas were mostly enriched with stress-response pathways, while the majority of the DEGs in the gastrointestinal tract belonged to metabolism-related pathways.

CONCLUSION

We demonstrated that the AVM induced acute toxicity, oxidative stress, osmoregulation disturbance, neurotoxicity and transcriptome imbalance in crayfish. These findings unraveled the detrimental effects of AVMs exposure on crayfish. © 2022 Society of Chemical Industry.

Differential expression of miRNAs in the body wall of the sea cucumber Apostichopus japonicus under heat stress

MicroRNAs, as one of the post-transcriptional regulation of genes, play an important role in the development process, cell differentiation and immune defense. The sea cucumber Apostichopus japonicus is an important cold-water species, known for its excellent nutritional and economic value, which usually encounters heat stress that affects its growth and leads to significant economic losses. However, there are few studies about the effect of miRNAs on heat stress in sea cucumbers. In this study, high-throughput sequencing was used to analyze miRNA expression in the body wall of sea cucumber between the control group (CS) and the heat stress group (HS). A total of 403 known miRNAs and 75 novel miRNAs were identified, of which 13 miRNAs were identified as significantly differentially expressed miRNAs (DEMs) in response to heat stress. A total of 16,563 target genes of DEMs were predicted, and 101 inversely correlated target genes that were potentially regulated by miRNAs in response to heat stress of sea cucumbers were obtained. Based on these results, miRNA-mRNA regulatory networks were constructed. The expression results of high-throughput sequencing were validated in nine DEMs and four differentially expressed genes (DEGs) by quantitative real-time polymerase chain reaction (RT-qPCR). Moreover, pathway enrichment of target genes suggested that several important regulatory pathways may play an important role in the heat stress process of sea cucumber, including ubiquitin-mediated proteolysis, notch single pathway and endocytosis. These results will provide basic data for future studies in miRNA regulation and molecular adaptive mechanisms of sea cucumbers under heat stress.

Lipid metabolism, immune and apoptosis transcriptomic responses of the hepatopancreas of Chinese mitten crab to the exposure to microcystin-LR

Global warming is favouring the incidence, intensity and duration of harmful cyanobacterial blooms. Microcystin-LR (MC-LR), a hepatotoxic agent, is produced during cyanobacterial blooms. To understand the molecular mechanisms of acute hepatotoxic effect of low doses of MC-LR in crab, we examined differentially expressed genes in samples of the hepatopancreas of Chinese mitten crab (Eriocheir sinensis) collected in 48 h after injections of MC-LR at doses of 0, 25, 50, and 75 µg/kg. The results revealed that MC-LR induced changes in corresponding gene led to the accumulation of triglycerides. MC-LR exposure affected sterol metabolism. Apoptosis-related genes such as Fas-L, Bcl-XL, Cytc, AiF, p53, PERK, calpain, CASP2, CASP7, α-tubulin, PARP, GF, G12, and PKC were upregulated. Conversely, expression levels of CASP10 and ASK1 were downregulated. Genes related to the regulation of actin cytoskeleton (Rho, ROCK, MLCP, MLC, PAK, and PFN) were upregulated. Further, expression levels of genes encoding fatty acid elongation-related enzymes were upregulated, but the expression of genes related to fatty acid synthesis was slightly down regulated. Taken together, these results demonstrated the hepatic toxicity and molecular mechanisms of changes in lipid metabolism, immune and apoptosis in Chinese mitten crab under the MC-LR-induced stress, which is the first report on crabs and performs a comprehensive analysis and a new insight of the molecular toxicological responses in crabs.

Integrated transcriptome analysis of immune-related mRNAs and microRNAs in Macrobrachium rosenbergii infected with Spiroplasma eriocheiris

Macrobrachium rosenbergii (M. rosenbergii), is a major aquaculture species in China and Southeast Asia. However, infection with Spiroplasma eriocheiris (S. eriocheiris) has caused huge economic losses to the cultivation of M. rosenbergii. Currently, there are few reports on the immune response mechanism of M. rosenbergii that are infected with S. eriocheiris. To clarify the immune response mechanism of M. rosenbergii infected with S. eriocheiris, the key immune genes which respond to the infection with the pathogen and the regulation of related microRNAs (miRNAs) on them were identified. In this study, the mRNA and miRNA transcriptome of hepatopancreas of M. rosenbergii at different infection stages were analyzed using high-throughput sequencing and qRT-PCR. In the mRNA transcriptome, 27,703 and 33,402 genes were expressed in healthy and susceptible M. rosenbergii, respectively. By digital gene-expression profiling analysis, 23,929 and 24,325 genes were expressed, and 223 and 373 genes were significantly up-regulated and down-regulated, respectively. A total of 145 key genes related to Toll, IMD, JAK/STAT and MAPK were excavated from the transcriptome. In the miRNA transcriptome, 549 miRNAs (Conserved: 41, PN-type: 83, PC-type: 425) were sequenced, of which 87 were significantly up-regulated and 23 were significantly down-regulated. Among the related immune pathways, there are 259 miRNAs involved in the regulation of target genes in the Toll and IMD pathways, 231 JAK/STAT pathways and 122 MAPK pathways. qRT-PCR differential detection of immune-related miRNAs and mRNAs showed that 22 miRNAs with significant differences (P < 0.05) such as mro-miR-100, PC-mro-3p-27 and PN-mro-miR-316 had corresponding regulatory relationships with 22 important immune genes such as TLR2, TLR3, TLR4, TLR5, MyD88, Pelle and Relish in different stages after infection. In this study, the immune genes and related regulatory miRNAs of M. rosenbergii in response to S. eriocheiris infection were obtained. The results can provide basic data to further reveal the immune defense mechanism of M. rosenbergii against S. eriocheiris infection.

Small RNA analysis of Perna viridis after exposure to Prorocentrum lima, a DSP toxins-producing dinoflagellate

Diarrheic shellfish poisoning toxins (DSP toxins) are a set of the most important phycotoxins produced by some dinoflagellates. Studies have shown that DSP toxins have various toxicities such as genotoxicity, cytotoxicity, and immunotoxicity to bivalve mollusks. However, these toxicities appear decreasing with exposure time and concentration of DSP toxins. The underlying mechanism involved remains unclear. In this study, small RNA sequencing was performed in the digestive gland of the mussel Perna viridis after exposure to DSP toxins-producing dinoflagellate Prorocentrum lima for different time periods. The potential roles of miRNAs in response and detoxification to DSP toxins in the mussel were analyzed. Small RNA sequencing of 12 samples from 72 individuals was conducted by BGISEQ-500. A total of 123 mature miRNAs were identified, including 90 conserved miRNAs and 33 potential novel miRNAs. After exposure to P. lima, multiple important miRNAs displayed some alterations. Further miRNA target prediction revealed some important genes involved in cytoskeleton, apoptosis, complement system and immune stress. qPCR demonstrated that miR-71_5, miR-750_1 and novel_mir4 were significantly up-regulated at 6 h after exposure to P. lima, while miR-100_2 was significantly down-regulated after 96 h of exposure. Accordingly, putative target genes of these differentially expressed miRNAs experienced some changes. After 6 h of DSP toxins exposure, NHLRC2 and C1q-like were significantly down-regulated. After 96 h of DSP toxins exposure, NHLRC2 was significantly up-regulated. It is reasonable to speculate that the mussel P. viridis might respond to DSP toxins through miR-750_1, novel_mir4 and miR-71_5 regulating the expression of relevant target genes involved in apoptosis, cytoskeleton, and immune response, etc. This study might provide new clues to uncover the toxic response of bivalve to DSP toxins and lay a foundation for revealing the roles of miRNAs in the environmental adaptation in shellfish.

Comparative transcriptomic analysis of the different developmental stages of ovary in red swamp crayfish Procambarus clarkii

Background

The red swamp crayfish Procambarus clarkii is a freshwater species that possesses high adaptability, environmental tolerance, and fecundity. P. clarkii is artificially farmed on a large scale in China. However, the molecular mechanisms of ovarian development in P. clarkii remain largely unknown. In this study, we identified four stages of P. clarkii ovary development, the previtellogenic stage (stage I), early vitellogenic stage (stage II), middle vitellogenic stage (stage III), and mature stage (stage IV) and compared the transcriptomics among these four stages through next-generation sequencing (NGS).

Results

The total numbers of clean reads of the four stages ranged from 42,013,648 to 62,220,956. A total of 216,444 unigenes were obtained, and the GC content of most unigenes was slightly less than the AT content. Principal Component Analysis (PCA) and Anosim analysis demonstrated that the grouping of these four stages was feasible, and each stage could be distinguished from the others. In the expression pattern analysis, 2301 genes were continuously increase from stage I to stage IV, and 2660 genes were sharply decrease at stage IV compared to stages I-III. By comparing each of the stages at the same time, four clusters of differentially expressed genes (DEGs) were found to be uniquely highly expressed in stage I (136 genes), stage II (43 genes), stage III-IV (49 genes), and stage IV (22 genes), thus exhibiting developmental stage specificity. Moreover, in comparisons between adjacent stages, the number of DEGs between stage III and IV was the highest. GO enrichment analysis demonstrated that nutrient reservoir activity was highest at stage II and that this played a foreshadowing role in ovarian development, and the GO terms of cell, intracellular and organelle participated in the ovary maturation during later stages. In addition, KEGG pathway analysis revealed that the early development of the ovary was mainly associated with the PI3K-Akt signaling pathway and focal adhesion; the middle developmental period was related to apoptosis, lysine biosynthesis, and the NF-kappa B signaling pathway; the late developmental period was involved with the cell cycle and the p53 signaling pathway.

Conclusion

These transcriptomic data provide insights into the molecular mechanisms of ovarian development in P. clarkii. The results will be helpful for improving the reproduction and development of this aquatic species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07537-x.

Integrated analysis of miRNAome and transcriptome reveals miRNA-mRNA network regulation in Vibrio alginolyticus infected thick shell mussel Mytilus coruscus

The thick shell mussel Mytilus coruscus has developed into a model species for studying the interaction between molluscs and environmental stimuli. Herein, integrated analysis of miRNAome and transcriptome was performed to reveal miRNA-mRNA network regulation in Vibrio alginolyticus infected M. coruscus. There have detected some histological abnormalities in digestive gland and gills of V. alginolyticus challenged mussels, ascertaining the effective irritation by the present bacterial strain. A total of 265 novel miRNAs were finally predicted, of which 26 were differentially expressed miRNAs (DEMs). Additionally, 667 differentially expressed genes (DEGs) were detected, which may be potentially associated with innate immune response to V. alginolyticus infection. A regulatory network linked to 22 important pathways and 16 DEMs and 34 OGs was constructed. Some traditional immune-related signaling pathways such as toll-like receptor signaling pathway (TLR) signaling pathway, transforming growth factor-beta (TGF-beta) signaling pathway, peroxisome, phagosome, lysosome, mammalian target of rapamyoin (mTOR) signaling pathway were linked to specific miRNAs and genes in this network. Further, interactional relationship between certain miRNAs and TLR pathway was dissected, which the results predicted that a number of TLRs and TLR-associated signaling genes including TLR1, TLR2, TLR4, TLR6, IRAK1, TRAF6, MAPK, and IL-17 were negatively regulated by novel_miR_11, novel_miR_145, novel_miR_196, novel_miR_5, novel_miR_163 and novel_miR_217 in the TLR pathway. Additionally, interactional relationship between novel_miR_145 and TLR2 was validated by laboratory experiment. The integrated analysis of mRNA and microRNA deep sequencing data exhibited a sophisticated miRNA-mRNA regulation network in M. coruscus in response to V. alginolyticus challenge, which shed a new light on the underlying mechanism of molluscan confronting bacterial infection.

Effects of hesperidin on the growth performance, antioxidant capacity, immune responses and disease resistance of red swamp crayfish (Procambarus clarkii)

We evaluated the effects of hesperidin on the nonspecific immunity, antioxidant capacity and growth performance of red swamp crayfish (Procambarus clarkii). A total of 900 healthy crayfish were randomly divided into six groups: the control group (fed the basal diet) and the HES25, HES50, HES75, HES100 and HES150 groups, which were fed the basal diet supplemented with 25, 50, 75, 100 and 150 mg kg^-1 hesperidin, respectively. The feeding experiment lasted 8 weeks. The results indicated that compared with the control group, the crayfish groups supplemented with 50-150 mg kg^-1 hesperidin had a decreased feed conversion ratio (FCR) and increased final body weight (FBW), specific growth rate (SGR) and weight gain (WG) (P < 0.05). The protein carbonyl content (PCC), reactive oxygen species (ROS) level and malondialdehyde (MDA) level in the hepatopancreas and hemocytes were significantly lower, while the total antioxidant capacity (T-AOC), glutathione peroxidase (GPx) activity, and superoxide dismutase (SOD) activity were significantly higher in the crayfish groups supplemented with 50-150 mg kg^-1 hesperidin than in the control group. Supplementation with 50-150 mg kg^-1 hesperidin significantly increased the activities of acid phosphatase (ACP), alkaline phosphatase (AKP), lysozyme (LZM), and phenoloxidase (PO) compared with the control group (P < 0.05); upregulated the mRNA expression of cyclophilin A (CypA), extracellular copper-zinc superoxide dismutase (ecCuZnSOD), GPxs, crustin, astacidin, Toll3 and heat shock protein 70 (HSP70) (P < 0.05); and decreased crayfish mortality following white spot syndrome virus (WSSV) infection. These findings indicate that dietary hesperidin supplementation at an optimum dose of 50-150 mg kg^-1 may effectively improve nonspecific immunity, antioxidant capacity and growth performance in crayfish.

The immune function of prophenoloxidase from red swamp crayfish (Procambarus clarkii) in response to bacterial infection

Prophenoloxidase (proPO) is the zymogen form of phenoloxidase (PO), a key enzyme in melanization cascade that has been co-opted in invertebrate immune reactions. There have been reported that proPO plays many essential roles in the crustacean immune system. However, little is known about the function of proPO from red swamp crayfish (Procambarus clarkii) which is an important cultured species worldwide. Here, we cloned and expressed proPO gene from red swamp crayfish (PcproPO). Subsequently, specific antibody against PcproPO was generated. The immune function of PcproPO was further characterized in vitro and in vivo. The results showed that the expression of PcproPO mRNA could be significantly up-regulated during the challenge of Gram-positive-negative (Vibrio parahaemolyticus) and Gram-positive-positive bacterial (Staphylococcus aureus). Furthermore, the purified recombinant PcproPO protein had a strong affinity binding to both bacteria and polysaccharides. In vivo knockdown of PcproPO could significantly reduce the crayfish bacterial clearance ability, resulting in the higher mortality of the crayfish during V. parahaemolyticus infection. In addition, in vitro knockdown of PcproPO in the hemocytes significantly reduced the phenoloxidase (PO) activity and the bacterial clearance ability, indicating that PcproPO might involve in hemocyte-mediated melanization. Our results will shed a new light on the immune function of PcproPO in the crayfish.

Integrated analysis of differentially expressed microRNAs and mRNAs to screen miRNAs and genes related to reproduction in Macrobrachium nipponense

Female Macrobrachium nipponense has the characteristic of short sexual maturity during the breeding season, which can increase breeding risk and lead to prevalent female individual miniaturization. In this study, we characterized micro (mi)RNA-seq data of the eyestalk (E) and cerebral ganglia (B) of female M. nipponense during breeding and non-breeding seasons. A total of 393 and 189 differentially expressed miRNAs (DE miRNAs) were identified in BSE vs. NBSE and BSB vs. NBSB, respectively. The most abundant up- and down-regulated DE miRNAs were miR-124, miR-14, and miR-7. Enrichment analysis showed that DE miRNA target genes were mainly involved in 'metabolic process' and 'binding', and were associated with 'neurohormonal regulation' and 'photoreceptor activity' signaling pathways. Integrated analysis of miRNA-mRNA expression showed that the most abundant DE miRNAs were miR-14 and miR-278 in BSE vs. NBSE and BSB vs. NBSB, respectively. Four pairs of DE miRNAs and their corresponding target annotated genes were selected from the DE miRNA-mRNA interaction network (bmo-miR-316-5p/opsin protein, ame-miR-125/skeletal muscle actin 8, dmo-miR-278/sugar transporter, and tca-miR-3885-5p/5-HT1 receptor). Gene expression analysis of these four pairs in different ovary development stages showed their potential regulatory roles in ovary maturation.

microRNAs in Macrobrachium olfersii embryos: Identification, their biogenesis components and potential targets

In embryonic development, microRNAs (miRNAs) regulate the complex gene expression associated with the complexity of embryogenesis. Today, few studies have been conducted on the identification of miRNAs and components of miRNA biogenesis on embryonic development in crustaceans, especially in prawns. In this context, the aim of this study was to identify in silico components of miRNA biogenesis, and miRNAs and potential target genes during embryonic development in the prawn Macrobrachium olfersii through small RNAs and transcriptome analyses. Using the miRDeep2 program, we identified 17 miRNA precursors in M. olfersii, which seven (miR-9, miR-10, miR-92, miR-125, miR-305, miR-1175, and miR-2788) were reported in the miRBase database, indicating high evolutionary conservation of these sequences among animals. The other 10 miRNAs of M. olfersii were novel miRNAs and only similar to Macrobrachium niponnense miRNAs, indicating genus-specific miRNAs. In addition, eight key components of miRNA biogenesis (DROSHA, PASHA/DGCR8, XPO5, RAN, DICER, TRBP2, AGO, and PIWI) were identified in M. olfersii embryos unigenes. In the annotation of miRNA targets, 516 genes were similar to known sequences in the GenBank database. Regarding the conserved miRNAs, we verified that they were differentially expressed during embryonic development in M. olfersii. In conclusion, this is the first study that identifies conserved and novel miRNAs in the prawn M. olfersii with some miRNA target genes involved in embryonic development. Our results will allow further studies on the function of these miRNAs and miRNA biogenesis components during embryonic development in M. olfersii and other prawns of commercial interest.

High-throughput sequencing analysis of microRNAs in gills of red swamp crayfish, Procambarus clarkii infected with white spot syndrome virus

MicroRNAs (miRNAs) are important posttranscriptional regulators. They play an important role in the antiviral innate immunity of invertebrates. In the present study, high-throughput small RNAs Illumina sequencing systems were carried out to identify differentially expressed miRNAs (DEMs) in the gills of Procambarus clarkii, which was challenged with white spot syndrome virus (WSSV). Our results identified 11,617 known and 6 novel miRNAs in normal group (NG) and WSSV-challenged group (WG) small RNA libraries. Additionally, 27 DEMs were shown to participate in the antiviral innate immunity of P. clarkii and were significantly upregulated or downregulated. In addition, the results of the KEGG pathway prediction of the DEMs target genes showed that putative target genes of these 27 DEMs were related mainly to the RNA transport pathway, tight junction pathway, mRNA surveillance pathway, regulation actin cytoskeleton pathway, focal adhesion pathway, and MAPK signaling pathway. These results provide important information for future studies about the antiviral innate immunity of crustaceans.

Transcriptional Identification of Related Proteins in the Immune System of the Crayfish Procambarus clarkii

The freshwater crayfish Procambarus clarkii is an animal model employed for physiological and immunological studies and is also of great economic importance in aquaculture. Although it is a species of easy husbandry, a high percentage of its production is lost annually as a result of infectious diseases. Currently, genetic information about the immune system of crustaceans is limited. Therefore, we used the abdominal nerve cord from P. clarkii to obtain its transcriptome using Next Generation Sequencing (NGS) to identify proteins that participate in the immune system. The reads were assembled de novo and consensus sequences with more than 3000 nucleotides were selected for analysis. The transcripts of the sequences of RNA were edited for annotation and sent to the GenBank database of the National Center for Biotechnology Information (NCBI). We made a list of accession numbers of the sequences which were organized by the putative role of the immune system pathway in which they participate. In this work, we report on 80 proteins identified from the transcriptome of crayfish related to the immune system, 74 of them being the first reported for P. clarkii. We hope that the knowledge of these sequences will contribute significantly to the development of future studies of the immune system in crustaceans.

Acute Hepatopancreatic Necrosis Disease (AHPND) related microRNAs in Litopenaeus vannamei infected with AHPND-causing strain of Vibrio parahemolyticus

BACKGROUND

Acute hepatopancreatic necrosis disease (AHPND) has emerged as a major debilitating disease that causes massive shrimp death resulting in substantial economic losses in shrimp aquaculture. Given that several diseases and infections have been associated with microRNAs (miRNAs), we conducted a comparative transcriptomic analysis using the AHPND (VA) and non-AHPND (VN) strains of Vibrio parahemolyticus to identify miRNAs potentially involved in AHPND pathogenesis in Litopenaeus vannamei.

RESULTS

A total of 83 miRNAs (47 upregulated and 36 downregulated) were significantly differentially expressed between the VA and VN challenged groups, while 222 target genes of these miRNAs were predicted. Functional enrichment analysis revealed that the miRNAs target genes were involved in multiple biological processes including metabolic pathways, amoebiasis, Vibrio cholerae infection etc. Finally, interaction network and qPCR (Real-time Quantitative PCR) analysis of 12 potential key AHPND-related miRNAs and their predicted target genes, revealed their possible involvement in modulating several immune-related processes in the pathogenesis of AHPND.

CONCLUSIONS

We have shown using comparative transcriptomic analysis, miRNAs and their target genes that are responsive to AHPND V. parahemolyticus infection in shrimp, therefore suggesting their possible role in defense response to AHPND V. parahemolyticus infection.

A potential microRNA regulation of immune-related genes in invertebrate haemocytes

Bivalve mollusks have been employed as sentinel organisms in environmental health programs due to their sedentary lifestyle, filter-feeding behavior and their ability to accumulate pathogens or toxin molecules inside tissues. Endocrine disrupting chemicals (EDCs) can be up taken and bioaccumulated, and due to sensibility of mollusks to these EDCs, being able to cause immune alterations. Recently, microRNAs (miRNAs) were shown to be involved in modulation and buffering developmental processes against the effects of environmental alterations and pathogenic microorganisms. Moreover, it is suggested that this miRNAs are incorporated into the estrogen-controlled immune network, regulating mechanism of immune gene expression at the posttranscriptional level, modulating immune responses as phagocytosis, redox reaction and apoptosis in bivalve haemocytes. Thus, miRNAs can be used as biomarkers that specifically elucidate immunotoxic effects caused by exogenous biotic or abiotic factors, and can act as useful tools in integrated monitoring environmental health programs. In this review, we aim to describe the investigations that have been carried out on miRNAs in bivalve mollusks, especially those associated with immune responses against infectious agents and xenobiotic exposure.

Antibacterial activity of hemocyanin from red swamp crayfish (Procambarus clarkii)

Hemocyanins (HMC): the copper-containing respiratory proteins present in invertebrate hemolymph, which plays many essential roles in the immune system. Currently, little is known about the HMC domains of Procambarus clarkii (P. clarkii) and their function in antimicrobial immune response. In this present study, we comparatively studied the expression pattern of native PcHMC with the three recombinant proteins of variable domains of crayfish hemocyanin (PcHMC-N, N-terminal domain of hemocyanin; PcHMC-T, tyrosinase domain of hemocyanin; PcHMC-C, C-terminal domain of hemocyanin). The results showed that three purified recombinant proteins had a strong binding to various bacteria and lipopolysaccharides that further highly agglutinated. The HMCs recombinant proteins showed strong antibacterial activity against V. parahaemolyticus and S. aureus by bacterial growth inhibition, phenoloxidase (PO) and phagocytosis assays. Specifically, rPcHMC1-T and rPcHMC1-C inhibited both the bacteria efficiently, rPcHMC1-T was highly upregulated the PO activity than the other recombinant proteins. Whereas, recombinant proteins pretreated crayfish hemocytes participated in phagocytosis activity, rPcHMC1-N and rPcHMC1-C proteins had a profound effect than the rPcHMC1-T on S. aureus and V. parahaemolyticus phagocytosis. The crayfish hemocyanin domains clearly exhibited antibacterial and phagocytic activities against both the bacteria, suggesting that its variable domains of hemocyanin have the different function on specific pathogen during the assault of pathogens.

Novel Insights into Antiviral Gene Regulation of Red Swamp Crayfish, Procambarus clarkii, Infected with White Spot Syndrome Virus

White spot syndrome virus (WSSV), one of the major pathogens of Procambarus clarkii, has caused severe disruption to the aquaculture industry of P. clarkii in China. To reveal the gene regulatory mechanisms underlying WSSV infection, a comparative transcriptome analysis was performed among WSSV-infected susceptible individuals (GS), viral resistant individuals (GR), and a non-infected control group (GC). A total of 61,349 unigenes were assembled from nine libraries. Subsequently, 515 and 1033 unigenes exhibited significant differential expression in sensitive and resistant crayfish individuals compared to the control group (GC). Many differentially expressed genes (e.g., C-type lectin 4, Peroxinectin, Prophenoloxidase, and Serine/threonine-protein kinase) observed in GR and GS play critical roles in pathogen recognition and viral defense reactions after WSSV infection. Importantly, the glycosaminoglycan biosynthesis-chondroitin sulfate/dermatan sulfate pathway was identified to play critical roles in defense to WSSV infection for resistant crayfish individuals by upregulating the chondroitin sulfate related genes for the synthesis of WSSV-sensitive, functional chondroitin sulfate chains containing E units. Numerous genes and the key pathways identified between resistant and susceptible P. clarkii individuals provide valuable insights regarding antiviral response mechanisms of decapoda species and may help to improve the selective breeding of P. clarkii WSSV-resistance.

Identification and characterization of intestine microRNAs and targets in red swamp crayfish, Procambarus clarkii infected with white spot syndrome virus

MicroRNAs (miRNAs) are small non-coding endogenous RNA molecules that play important roles in the innate immunity system of invertebrates, especially in the aspect of antivirus. In the present study, high-throughput small RNA Illumina sequencing systems were used to identify differentially expressed miRNAs (DEMs) from the intestines of Procambarus clarkii that were infected with white spot syndrome virus (WSSV). As a result, 39 known and 12 novel miRNAs were identified in both NG and WG small RNA libraries. Seven DEMs were determined to be involved in the antiviral innate immunity in the intestines of P. clarkii. The results of the target gene predictions of the DEMs showed that the putative target genes of these 7 DEMs are related to tight junctions, vascular smooth muscle contraction regulation of the actin cytoskeleton, focal adhesion, RNA transport, mRNA surveillance, viral carcinogenesis, and Salmonella infection. These results provide theoretical insights for future studies on the antiviral immunity of crustaceans.

Identification and characterization of lymph organ microRNAs in red swamp crayfish, Procambarus clarkii infected with white spot syndrome virus

MicroRNAs (miRNAs) were important post-transcriptional regulators and played vital roles in innate immunity system of invertebrates, especially in the aspect of antivirus. In this study, using high-throughput small RNAs Illumina sequencing system, differentially expressed miRNAs (DEMs) from lymph organs in red swamp crayfish, Procambarus clarkii, infected with white spot syndrome virus, were identified. As a result, 32 known miRNAs and 7 novel miRNAs were identified in crayfish lymph organ small RNAs library of NG and WG. Among them, 7 differentially expressed miRNAs (DEMs) were predicted to be involved in the lymph organ antiviral innate immunity of P. clarkii. Besides, the results showed that putative target genes of these DEMs were related with tight junction, RNA transport, regulation of actin cytoskeleton, focal adhesion, vascular smooth muscle contraction, mRNA surveillance pathway, NOD-like receptor signaling pathway, leukocyte transendothelial migration, and protein processing in endoplasmic reticulum. These results might provide the guiding theoretical foundation for future studies about crustaceans' antiviral innate immunity.

Sex and tissue specific gene expression patterns identified following de novo transcriptomic analysis of the Norway lobster, Nephrops norvegicus

BACKGROUND

The Norway lobster, Nephrops norvegicus, is economically important in European fisheries and is a key organism in local marine ecosystems. Despite multi-faceted scientific interest in this species, our current knowledge of genetic resources in this species remains very limited. Here, we generated a reference de novo transcriptome for N. norvegicus from multiple tissues in both sexes. Bioinformatic analyses were conducted to detect transcripts that were expressed exclusively in either males or females. Patterns were validated via RT-PCR.

RESULTS

Sixteen N. norvegicus libraries were sequenced from immature and mature ovary, testis and vas deferens (including the masculinizing androgenic gland). In addition, eyestalk, brain, thoracic ganglia and hepatopancreas tissues were screened in males and both immature and mature females. RNA-Sequencing resulted in >600 million reads. De novo assembly that combined the current dataset with two previously published libraries from eyestalk tissue, yielded a reference transcriptome of 333,225 transcripts with an average size of 708 base pairs (bp), with an N50 of 1272 bp. Sex-specific transcripts were detected primarily in gonads followed by hepatopancreas, brain, thoracic ganglia, and eyestalk, respectively. Candidate transcripts that were expressed exclusively either in males or females were highlighted and the 10 most abundant ones were validated via RT-PCR. Among the most highly expressed genes were Serine threonine protein kinase in testis and Vitellogenin in female hepatopancreas. These results align closely with gene annotation results. Moreover, a differential expression heatmap showed that the majority of differentially expressed transcripts were identified in gonad and eyestalk tissues. Results indicate that sex-specific gene expression patterns in Norway lobster are controlled by differences in gene regulation pattern between males and females in somatic tissues.

CONCLUSIONS

The current study presents the first multi-tissue reference transcriptome for the Norway lobster that can be applied to future biological, wild restocking and fisheries studies. Sex-specific markers were mainly expressed in males implying that males may experience stronger selection than females. It is apparent that differential expression is due to sex-specific gene regulatory pathways that are present in somatic tissues and not from effects of genes located on heterogametic sex chromosomes. The N. norvegicus data provide a foundation for future gene-based reproductive studies.

The miiuy croaker microRNA transcriptome and microRNA regulation of RIG-I like receptor signaling pathway after poly(I:C) stimulation

MicroRNAs (miRNAs) as endogenous small non-coding RNAs play key regulatory roles in diverse biological processes via degrading the target mRNAs or inhibiting protein translation. Previously many researchers have reported the identification, characteristic of miRNAs and the interaction with its target gene. But, the study on the regulation of miRNAs to biological processes via regulatory the key signaling pathway was still limited. In order to comprehend the regulatory mechanism of miRNAs, two small RNA libraries from the spleen of miiuy croaker individuals with or without poly(I:C) infection were constructed. The 197 conserved miRNAs and 75 novel miRNAs were identified, and 14 conserved and 8 novel miRNAs appeared significant variations. Those differently expressed miRNAs relate to immune regulation of miiuy croaker. Furthermore, expressions of four differently expressed miRNAs were validated by qRT-PCR, and the result was consistent with sequencing data. The target genes of the differently expressed miRNAs in the two libraries were predicted, and some candidate target genes were involved in the RIG-I-like receptor (RLR) signaling pathway. The negative regulation of miRNAs to target genes were confirmed by comparing the expression pattern of miRNAs and their target genes. The results of regulating target genes were that firstly directly or indirectly activating the downstream signaling cascades and subsequent inducting the type I interferon, inflammatory cytokines and apoptosis. These studies could help us to deeper understand the roles of miRNAs played in the fish immune system, and provide a new way to investigate the defense mechanism of fish.

Identification and profiling of growth-related microRNAs of the swimming crab Portunus trituberculatus by using Solexa deep sequencing

MicroRNAs (miRNAs) are a class of endogenous small non-coding RNAs that regulate gene expression by post-transcriptional repression of mRNAs. The swimming crab Portunus trituberculatus is one of the most important crustacean species for aquaculture in China. However, to date no miRNAs have been reported to for modulating growth in P. trituberculatus. To investigate miRNAs involved in the growth of this species, we constructed six small RNA libraries for big individuals (BIs) and small individuals (SIs) from a highly inbred family. Six mixed RNA pools of five tissues (eyestalk, gill, heart, hepatopancreas, and muscle) were obtained. By aligning sequencing data with those for known miRNAs, a total of 404 miRNAs, including 339 known and 65 novel miRNAs, were identified from the six libraries. MiR-100 and miR-276a-3p were among the most prominent miRNA species. We identified seven differentially expressed miRNAs between the BIs and SIs, which were validated using real-time PCR. Preliminary analyzes of their putative target genes and GO and KEGG pathway analyzes showed that these differentially expressed miRNAs could play important roles in global transcriptional depression and cell differentiation of P. trituberculatus. This study reveals the first miRNA profile related to the body growth of P. trituberculatus, which would be particularly useful for crab breeding programs.

Comprehensive identification and profiling of host miRNAs in response to Singapore grouper iridovirus (SGIV) infection in grouper (Epinephelus coioides)

microRNAs (miRNAs) are an evolutionarily conserved class of non-coding RNA molecules that participate in various biological processes. Employment of high-throughput screening strategies greatly prompts the investigation and profiling of miRNAs in diverse species. In recent years, grouper (Epinephelus spp.) aquaculture was severely affected by iridoviral diseases. However, knowledge regarding the host immune responses to viral infection, especially the miRNA-mediated immune regulatory roles, is rather limited. In this study, by employing Solexa deep sequencing approach, we identified 116 grouper miRNAs from grouper spleen-derived cells (GS). As expected, these miRNAs shared high sequence similarity with miRNAs identified in zebrafish (Danio rerio), pufferfish (Fugu rubripes), and other higher vertebrates. In the process of Singapore grouper iridovirus (SGIV) infection, 45 and 43 miRNAs with altered expression (>1.5-fold) were identified by miRNA microarray assays in grouper spleen tissues and GS cells, respectively. Furthermore, target prediction revealed 189 putative targets of these grouper miRNAs.

Roles of small RNAs in the immune defense mechanisms of crustaceans

Small RNAs, 21-24 nucleotides in length, are non-coding RNAs found in most multicellular organisms, as well as in some viruses. There are three main types of small RNAs including microRNA (miRNA), small-interfering RNA (siRNA), and piwi-interacting RNA (piRNA). Small RNAs play key roles in the genetic regulation of eukaryotes; at least 50% of all eukaryote genes are the targets of small RNAs. In recent years, studies have shown that some unique small RNAs are involved in the immune response of crustaceans, leading to lower or higher immune responses to infections and diseases. SiRNAs could be used as therapy for virus infection. In this review, we provide an overview of the diverse roles of small RNAs in the immune defense mechanisms of crustaceans.

The environmental quality of Doñana surrounding areas affects the immune transcriptional profile of inhabitant crayfish Procambarus clarkii

This study aimed to identify differentially expressed genes in Procambarus clarkii crayfish collected from locations of different environmental qualities in the Doñana National Park surrounding areas. The pollution sustained by the crayfish was confirmed by their hepatopancreatic metal concentration. We generated forward and reverse libraries by suppression subtractive hybridization (SSH) to analyze the transcriptional profiles of crayfish from moderately and highly polluted zones in comparison with the control site within the Doñana Biological Reserve. Forty-three differentially expressed genes were detected, and most of them were identified as genes involved in a variety of biological functions, particularly in the innate immune response. To verify the SSH results and assess interindividual variability nine transcripts (ALP, AST, BTF3, CHIT, CTS, ferritin, HC, HC2, and SPINK4) were selected for absolute quantification by real-time qRT-PCR. The qRT-PCR data revealed substantial differences in the absolute amounts of the nine transcripts and confirmed their up- or down-regulation in the polluted sites. Additionally, a positive and significant linear correlation was found between the hepatopancreatic copper concentration and the levels of the transcripts encoding hemocyanins. Finally, the transcriptomic study was complemented with a detailed analysis of SNP profiles of the selected transcripts that revealed point mutations that might underlie adaptive response to environmental stress in P. clarkii. Overall, this work provides novel insights into the molecular pathways that could mediate the response to environmental pollutants in P. clarkii emphasizing the central role of the immune function and thus, should clearly benefit further immunotoxicological research in this organism.

Identification and comparative analysis of the Pseudosciaena crocea microRNA transcriptome response to poly(I:C) infection using a deep sequencing approach

Two sRNA libraries with or without poly(I:C) infection of large yellow croaker Pseudosciaena crocea were constructed and sequenced using the high-throughput Illumina/Solexa deep sequencing technology. The high-throughput sequencing pipeline yielded 163,79,272 and 217,07,070 raw reads corresponding to 132,27,594 and 206,86,409 clean reads for the normal and infected libraries, respectively. Bioinfromatic analysis identified 534 miRNAs, of which, 158 miRNAs were known in miRBase 20.0 and the remaining 376 were not found homology to any known metazoan miRNAs, suggesting a possible species-specificity. We analyzed the significance of differently expressed miRNAs between two libraries using pairwise comparison. There was significant differential expression of 112 miRNAs (p < 0.001) between two libraries. Thereinto, a number of known miRNAs were identified immune-related. Real-time quantitative PCR experiments (RT-qPCR) were preformed for 6 miRNAs of the two samples, and agreement was found between the sequencing and RT-qPCR data. To our knowledge, this is the first comprehensive study of miRNAs in P. crocea and of expression analysis of P. crocea miRNAs in response to poly(I:C) infection, and many miRNAs were differentially regulated under normal and infection conditions. These findings deepened our understanding of the role of miRNAs in the intricate host's immune system, and should be useful to develop new control strategies for host immune defense against various foreign infection in P. crocea.

Identification and differential expression of hepatopancreas microRNAs in red swamp crayfish fed with emodin diet

Using high-throughput Illumina Solexa system, the differential miRNA expressions from hepatopancreas in red swamp crayfish (Procambarus clarkii) fed with diets containing 0 (control) and 75 mg emodin kg(-1) (trial) were identified, respectively. As a result, 13,335,928 raw reads from the control sample and 14,938,951 raw reads from the trial sample were obtained while 13,053,344 (98.77%) and 14,517,522 (98.34%) small RNA were identified, respectively. 106 mature miRNAs (belonging to 68 miRNA gene families) were identified. 35 miRNAs displayed significantly differential expressions between two libraries. Of these, comparing to the control library, 6 miRNAs were significantly up-regulated and 29 miRNAs were significantly down-regulated. Moreover, 5 novel miRNAs (2 from control sample, 3 from trial sample) and target genes were predicted. GO analysis suggested that these miRNAs might be involved in innate immune response, growth, metabolism, cellular process, biological regulation and stimulus response. Our knowledge from this study could contribute to a better understanding of the miRNAs roles in regulating innate immune response and the study of miRNA function in crayfish.