Identification and comparative analysis of the Eriocheir sinensis microRNA transcriptome response to Spiroplasma eriocheiris infection using a deep sequencing approach

Current Disease Threats for Cultivated Crab Eriocheir sinensis in China

The Chinese mitten crab, Eriocheir sinensis, is a commercially important crustacean in China due to its great commercial value and compatibility in a variety of aquaculture systems. However, increases in its production have been accompanied by the emergence of various diseases affecting yield, profit, and trading potential. In this study, we review the pathogenic agents associated with E. sinensis since the start of its commercial culture. The history of crab cultivation implies that increased pathogen transfer can occur as E. sinensis aquaculture grows because polyculture of E. sinensis with other aquaculture species is a prevalent practice. With this in mind, a special focus of this review is placed on pathogens that were initially discovered in other crustacean species but have since been demonstrated to infect and cause disease in E. sinensis. We expect that this review will not only offer recommendations for disease management in the E. sinensis aquaculture sector but will also advance other crustacean cultivation.

Deep sequencing of microRNAs reveals circadian-dependent microRNA expression in the eyestalks of the Chinese mitten crab Eriocheir sinensis

MicroRNAs (miRNAs) are small endogenous non-coding RNAs. In crustaceans, miRNAs might be involved in the regulation of circadian rhythms. Many physiological functions of crustaceans including immunity and hormone secretion exhibit circadian rhythms, but it remains unclear whether specific miRNAs contribute to the alteration of crustacean physiological processes under circadian rhythms. This study investigated the mechanisms of miRNA regulation of circadian rhythms in the Chinese mitten crab (Eriocheir sinensis), one of China's most important aquaculture species. We obtained eyestalks from crab specimens at four time points (6:00; 12:00; 18:00; 24:00) during a 24-h period. We identified 725 mature miRNAs, with 23 known miRNAs differentially expressed depending on the time of day. The Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses revealed that the putative target genes for differentially expressed miRNAs were significantly enriched in the immune response and endocrine-related pathways. Numerous putative target genes are involved in the circadian-related pathways and enriched on circadian-control genes. These results suggest that the expression of miRNAs regulates some specific physiological functions in E. sinensis under circadian cycles. We also profiled various putative target genes enriched under the circadian-related pathway. This study performed miRNA expression in the eyestalks of E. sinensis during a 24-h daily cycle, providing insights into the molecular mechanism underlying crustacean circadian rhythms and suggesting miRNAs' role in studying crustacean physiology should not be overlooked.

Comparative transcriptome, digital gene expression and proteome profiling analyses provide insights into the brachyurization from the megalopa to the first juvenile in Eriocheir sinensis

Eriocheir sinensis larva normally experiences 11 stages. The reduced abdomen folded beneath the thorax is the most prominent characteristic of morphological alteration from megalopa to juvenile crab. Up to date, the molecular mechanisms of brachyurization remain a mystery. Here, transcriptome library, digital gene expression (DGE) libraries and proteome libraries at two developmental stages [the megalopa stage of E. sinensis (stage M) and the first stage of juvenile crab (stage J1)] of the Chinese mitten crab larva were constructed for RNA sequencing and iTRAQ approaches followed by bioinformatics analysis, respectively. In total, 1106 genes and 871 proteins were differentially expressed between the stage M and stage J1. Moreover, several important pathways were identified, including biosynthesis of secondary metabolites, metabolic pathways, focal adhesion, and some disease pathways. Besides, muscle contraction, oxidative phosphorylation, calcium signaling, PI3K-Akt, DNA replication pathway, and integrin signaling pathway also had important functions in brachyurization process. Furthermore, the components, actin, actin-related protein, collagens, filamin-A/B, laminin, integrins, paxillin, and fibronectin had up-regulated expression levels in M stage compared to J1 stage.

Full-length transcriptome sequencing and comparative transcriptome analysis of Eriocheir sinensis in response to infection by the microsporidian Hepatospora eriocheir

As a new generation of high-throughput sequencing technology, PacBio Iso-Seq technology (Iso-Seq) provides a better alternative sequencing method for the acquisition of full-length unigenes. In this study, a total of 22.27 gigabyte (Gb) subread bases and 128,614 non-redundant unigenes (mean length: 2,324 bp) were obtained from six main tissues of Eriocheir sinensis including the heart, nerve, intestine, muscle, gills and hepatopancreas. In addition, 74,732 unigenes were mapped to at least one of the following databases: Non-Redundant Protein Sequence Database (NR), Gene Ontology (GO), Kyoto Encyclopaedia of Genes and Genomes (KEGG), KEGG Orthology (KO) and Protein family (Pfam). In addition, 6696 transcription factors (TFs), 28,458 long non-coding RNAs (lncRNAs) and 94,230 mRNA-miRNA pairs were identified. Hepatospora eriocheir is the primary pathogen of E. sinensis and can cause hepatopancreatic necrosis disease (HPND); the intestine is the main target tissue. Here, we attempted to identify the key genes related to H. eriocheir infection in the intestines of E. sinensis. By combining Iso-Seq and Illumina RNA-seq analysis, we identified a total of 12,708 differentially expressed unigenes (DEUs; 6,696 upregulated and 6,012 downregulated) in the crab intestine following infection with H. eriocheir. Based on the biological analysis of these DEUs, several key processes were identified, including energy metabolism-related pathways, cell apoptosis and innate immune-related pathways. Twelve selected genes from these DEUs were subsequently verified by quantitative real-time PCR (qRT-PCR) analysis. Our findings enhance our understanding of the E. sinensis transcriptome and the specific association between E. sinensis and H. eriocheir 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.

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.

The miRNAs profiling revealed by high-throughput sequencing upon WSSV infection in mud crab Scylla paramamosain

microRNAs (miRNAs) are known to regulate various immune functions by silencing the target genes in both vertebrates and invertebrates. However, in mud crab Scylla paramamosain, the role of miRNAs during the response to virus invasion remains unclear. To investigate the roles of miRNAs in S. paramamosain during virus infection, the mud crab was challenged with white spot syndrome virus (WSSV) and then subjected to the transcriptional analysis at different conditions. The results of high-throughput sequencing revealed that 940,379 and 1,306,023 high-quality mappable reads were detected in the hemocyte of normal and WSSV-infected mud crabs, respectively. Besides, the total number of 261 unique miRNAs were identified. Among them, 131 miRNAs were specifically expressed in the hemocytes of normal mud crabs, 46 miRNAs were specifically transcribed in those of WSSV-infected individuals, the other 84 miRNAs were expressed in both normal and WSSV-infected individuals. Furthermore, a number of 152 (89 down-regulated and 63 up-regulated) miRNAs were found to be differentially expressed in the WSSV-infected hemocytes, normalized to the controls. The identified miRNAs were subjected to GO analysis and target gene prediction and the results suggested that the differentially regulated miRNAs were mainly correlated with the changes of the immune responses of the hemocytes, including phagocytosis, melanism, and apoptosis as well. Taken together, the results demonstrated that the expressed miRNAs during the virus infection were mainly involved in the regulation of immunological pathways in mud crabs. Our findings not only enrich the understanding of the functions of miRNAs in the innate immune system but also provide some novel potential targets for the prevention of WSSV infection in crustaceans.

Transcriptome analysis of hepatopancreas of Eriocheir sinensis with hepatopancreatic necrosis disease (HPND)

Hepatopancreatic necrosis disease (HPND) is a newly emerging disease in the Chinese mitten crab, Eriocheir sinensis, which has resulted in large economic losses. However, the underlying cause of this disease remains unclear. To better understand the pathogenesis and pathogenic mechanism of HPND, we compared the transcriptome differences of the hepatopancreas of E. sinensis with and without HPND. The analysis yielded > 30 million reads for each sample of three test (with HPND) and three control groups (without HPND). We observed 978 downregulated genes and 644 upregulated genes. Among the gene ontology categories "biological process," "cellular component," and "molecular function", the subcategories cellular process, single-organism process, biological regulation, metabolic process, cell part, organelle, organelle part, binding, and catalytic were enriched. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis showed that "metabolism of xenobiotics by cytochrome P450," "drug metabolism-cytochrome P450," "chemical carcinogenesis," and "material metabolism" were the "five" most significantly enriched pathways in the hepatopancreas of E. sinensis with HPND. The results revealed that material metabolic abnormalities and drug effects from the external environment might be associated with HPND in the Chinese mitten crab. Considering the wide use of pyrethroids for pond cleaning in Xinghua city, we speculated that pyrethroids might cause HPND in the Chinese mitten crab. Our study provided useful information about the cause and pathogenetic mechanisms of HPND and could help to prevent this disease in production practice.

Identification and profiling of microRNAs responsive to cadmium toxicity in hepatopancreas of the freshwater crab Sinopotamon henanense

Background

Cadmium (Cd) is a ubiquitous environmental toxicant for aquatic animals. The freshwater crab, Sinopotamon henanense (S. henanense), is a useful model for monitoring Cd exposure since it is widely distributed in sediments whereby it tends to accumulate several toxicants, including Cd. In the recent years, the toxic effects of Cd in the hepatopancreas of S. henanense have been demonstrated by a series of biochemical analysis and ultrastructural observations as well as the deep sequencing approaches and gene expression profile analysis. However, the post-transcriptional regulatory network underlying Cd toxicity in S.henanense is still largely unknown.

Results

The miRNA transcriptional profile of the hepatopancreas of S. henanense was used to investigate the expression levels of miRNAs in response to Cd toxicity. In total, 464 known miRNAs and 191 novel miRNAs were identified. Among these 656 miRNAs, 126 known miRNAs could be matched with the miRNAs of Portunus trituberculatus, Eriocheir sinensis and Scylla paramamosain. Furthermore, a total of 24 conserved miRNAs were detected in these four crab species. Fifty-one differentially expressed miRNAs were identified in the Cd-exposed group, with 31 up-regulated and 20 down-regulated. Eight of the differentially expressed miRNAs were randomly selected and verified by the quantitative real-time PCR (qRT-PCR), and there was a general consistency (87.25%) between the qRT-PCR and miRNA transcriptome data. A total of 5258 target genes were screened by bioinformatics prediction. GO term analysis showed that, 17 GO terms were significantly enriched, which were mainly related to the regulation of oxidoreductase activity. KEGG pathway analysis showed that 18 pathways were significantly enriched, which were mainly associated with the biosynthesis, modification and degradation of proteins.

Conclusion

In response to Cd toxicity, in the hepatopancreas of S. henanense, the expressions of significant amount of miRNAs were altered, which may be an adaptation to resist the oxidative stress induced by Cd. These results provide a basis for further studies of miRNA-mediated functional adaptation of the animal to combat Cd toxicity.

The Single-molecule long-read sequencing of Scylla paramamosain

Scylla paramamosain is an important aquaculture crab, which has great economical and nutritional value. To the best of our knowledge, few full-length crab transcriptomes are available. In this study, a library composed of 12 different tissues including gill, hepatopancreas, muscle, cerebral ganglion, eyestalk, thoracic ganglia, intestine, heart, testis, ovary, sperm reservoir, and hemocyte was constructed and sequenced using Pacific Biosciences single-molecule real-time (SMRT) long-read sequencing technology. A total of 284803 full-length non-chimeric reads were obtained, from which 79005 high-quality unique transcripts were obtained after error correction and sequence clustering and redundant. Additionally, a total of 52544 transcripts were annotated against protein database (NCBI nonredundant, Swiss-Prot, KOG, and KEGG database). A total of 23644 long non-coding RNAs (lncRNAs) and 131561 simple sequence repeats (SSRs) were identified. Meanwhile, the isoforms of many genes were also identified in this study. Our study provides a rich set of full-length cDNA sequences for S. paramamosain, which will greatly facilitate S. paramamosain research.

Identification and expression profiling analysis of microRNAs in Nile tilapia (Oreochromis niloticus) in response to Streptococcus agalactiae infection

MicroRNAs (miRNAs) play vital regulatory roles in various biological processes, including in immune responses. Nile tilapia (Oreochromis niloticus) is an important commercial fish species in China. To identify immune-related miRNAs of O. niloticus, 4 libraries from liver during S. agalactiae infection (0 h, 5 h, 50 h, and 7 d) were sequenced by high-throughput sequencing technology in tilapia. We obtained 10,703,531, 11,507,163, 11,180,179 and 13,408,414 clean reads per library, respectively. In our results, a total of 482 miRNAs were identified through bioinformatic analysis, including 220 conserved miRNAs and 262 putative novel miRNAs. Moreover, 21 (4.36%), 50 (10.37%), and 46 (9.54%) miRNAs were significantly differentially expressed at 5 h, 50 h and 7 d, respectively. In addition, 6939 target genes regulated by these differentially expressed miRNAs were predicted, and their functional annotations were predicted by Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, which revealed that a majority of differentially expressed miRNAs were involved in apoptotic process, metabolic process, and immune responses. Finally, Real-time quantitative PCR experiments were performed for 7 miRNAs by stem-loop RT-PCR, and a general agreement was confirmed between the sequencing and RT-qPCR data. To our understanding, this is the first report of comprehensive identification of O. niloticus miRNAs being differentially regulated in liver related to S. agalactiae infection. This work provides an opportunity for further understanding of the molecular mechanisms of miRNA regulation in O. niloticus host-pathogen interactions, and genetic resources for molecular assistant selection for disease resistant breeding program.

Changes of microRNAs expression profiles from red swamp crayfish (Procambarus clarkia) hemolymph exosomes in response to WSSV infection

MicroRNAs (miRNAs) as short noncoding RNAs play important regulatory roles in diverse biological processes by degrading the target mRNAs, and could be delivered by exosomes. WSSV is a highly pathogenic and prevalent virus, and has brought high mortality of P. clarkia. Till present, no studies focus on the miRNAs changes in exosomes during WSSV infection. To understand the different virulence of WSSV on miRNAs expression in P. clarkia hemolymph exosome, the deep sequencing was performed to compare the small RNA libraries from the hemolymph exosome of P. clarkia individuals with or without WSSV infections. From the TEM observations, NTA and Western Blot analysis, the extracted exosomes were well identified with classic characteristics. The 209 conserved miRNAs and 250 novel miRNAs were identified from the small RNA libraries. In response to WSSV infection, there were about 98 miRNAs significantly up-regulated and 59 miRNAs significantly down-regulated. The target genes prediction, GO and KEGG enrichment analysis revealed that some target genes of P. clarkia miRNAs were grouped mainly into the categories of biological regulation, immune system process, signal pathway and other more functions. This is the first report of comprehensive identification of P. clarkia hemolymph exosome miRNAs being differentially regulated in response to WSSV infection. These results will help to understand the hemolymph exosome miRNAs response to different virulence WSSV infection.

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.

Differential expression of microRNAs in Portunus trituberculatus in response to Hematodinium parasites

Latest studies have indicated that microRNAs (miRNAs) play important roles in defending against bacterial and viral infections in marine crustacean, whereas little is known regarding the immunological roles of crustacean miRNAs in response to parasitic infection. To further reveal the host-parasite interactions between the parasitic dinoflagellate Hematodinium and its crustacean hosts, we applied the high-throughput sequencing technology to identify and characterize miRNAs in the Chinese swimming crab Portunus trituberculatus challenged with the Hematodinium parasite at a timescale of 16 days (d). A total of 168 miRNAs were identified and 51 miRNAs were differentially expressed in the hepatopancreas tissues of affected hosts. Eleven of the differentially expressed miRNAs were selected and verified by the quantitative real-time RT-PCR (qRT-PCR), manifesting the consistency between the high throughout sequencing and qRT-PCR assays. Further analysis of the putative target genes indicated that various immune-related pathways (e.g. endocytosis, Fc gamma R-mediated phagocytosis, lysosome, ECM-receptor interaction, complement and coagulation cascades, antigen processing and presentation, focal adhesion, etc.) and signal transduction pathways (e.g. JAK-STAT signaling pathway, MAPK signaling pathway, p53 signaling pathway, etc.) were mediated by the differentially expressed miRNAs. The results presented fundamental knowledge on the immunological roles of crustacean miRNAs and contributed to the better understanding of hosts' miRNAs-mediated immunity against the parasitic infection.

MicroRNA repertoire and comparative analysis of Andrias davidianus infected with ranavirus using deep sequencing

Andrias davidianus is a large and economically important amphibian in China. Ranavirus infection causes serious losses in A. davidianus farming industry. MicroRNA mediated host-pathogen interactions are important in antiviral defense. In this study, five small-RNA libraries from ranavirus infected and non-infected A. davidianus spleens were sequenced using high throughput sequencing. The miRNA expression pattern, potential functions, and target genes were investigated. In total, 1356 known and 431 novel miRNAs were discovered. GO and KEGG analysis revealed that certain miRNA target genes are associated with apoptotic, signal pathway, and immune response categories. Analysis identified 82 downregulated and 9 upregulated differentially expressed miRNAs, whose putative target genes are involved in pattern-recognition receptor signaling pathways and immune response. These findings suggested miRNAs play key roles in A. davidianus's response to ranavirus and could provide a reference for further miRNA functional identification, leading to novel approaches to improve A. davidianus ranavirus resistance.

Analysis of the microRNA transcriptome of Daphnia pulex during aging

Daphnia pulex is an important food organism that exhibits a particular mode of reproduction known as cyclical parthenogenesis (asexual) and sexual reproduction. Regulation of the aging process by microRNAs (miRNAs) is a research hotspot in miRNA studies. To investigate a possible role of miRNAs in regulating aging and senescence, we used Illumina HiSeq to sequence two miRNA libraries from 1-day-old (1d) and 25-day-old (25d) D. pulex specimens. In total, we obtained 11,218,097 clean reads and 28,569 unique miRNAs from 1d specimens and 11,819,106 clean reads and 44,709 unique miRNAs from 25d specimens. Bioinformatic analyses was used to identify 1335 differentially expressed miRNAs from known miRNAs, including 127 miRNAs that exhibited statistically significant differences (P < 0.01); 92 miRNAs were upregulated and 35 were downregulated. Quantitative real-time (qRT)-PCR experiments were performed for nine miRNAs from five samples (1d, 5d, 10d, 15d, 20d and 25d) during the aging process, and the sequencing and qRT-PCR data were found to be consistent. Ninety-four miRNAs were predicted to correspond to 2014 target genes in known miRNAs with 4032 target gene sites. Sixteen pathways changed significantly (P < 0.05) at different developmental stages, revealing many important principles of the miRNA regulatory aging network of D. pulex. Overall, the difference in miRNA expression profile during aging of D. pulex forms a basis for further studies aimed at understanding the role of miRNAs in regulating aging, reproductive transformation, senescence, and longevity.

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.

Differential expression of microRNAs in hemocytes from white shrimp Litopenaeus vannamei under copper stress

MicroRNAs (miRNAs) are small noncoding RNAs that regulate diverse cellular processes, including organismal stress response, through posttranscriptional repression of gene transcripts. They are known to have antiviral functions in aquatic crustacean species, but little is known about the role of miRNAs against environmental stress caused by Cu, a common chemical contaminant in aquatic environment. We performed small RNA sequencing to characterize the differentially expressed microRNAs in Cu exposed shrimp. A total of 4524 known miRNAs and 73 novel miRNAs were significantly (P < .05) differentially expressed after Cu exposure. The peak size of miRNAs was 22 nt. Among them, 218 miRNAs were conserved across 115 species. The validation of 12 miRNAs by stem-loop quantitative RT-PCR were found to be coherent with the expression profile of deep sequencing data as evaluated with Pearson's correlation coefficient (r = 0.707). Target genes of these differentially expressed miRNAs related to immune defense, apoptosis, and xenobiotics metabolism also showed significant changes in expression under Cu stress. The present study provides the first characterization of L. vannamei miRNAs and some target genes expression in response to Cu stress, and the findings support the hypothesis that certain miRNAs along with their target genes might be essential in the intricate adaptive response regulation networks. Our current study will provide valuable information to take an insight into molecular mechanism of L. vannamei against environmental stress.

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.

Modulation of Atlantic salmon miRNome response to sea louse infestation

MicroRNAs are non-coding RNA that plays a crucial role in post-transcriptional regulation and immune system regulation. On other hand, sea lice are prevalent parasites that affect salmon farming, generating different degrees of immune suppression depending on the salmon and sea louse species. Caligus rogercresseyi for example, which affects the salmon industry in Chile, decreases Th1 response, macrophage activation, TLR-mediated response and iron regulation in infected fish. In this study, we explore Atlantic salmon miRNome during infestation by C. rogercresseyi. Using small RNA sequencing, we annotated 1718 miRNAs for skin and head kidney from infected Atlantic salmon. The most abundant families identified were mir-10, mir-21, mir-30, mir-181 and let7. Significant differences were found between tissue, with 1404 annotated miRNA in head kidney and 529 in skin. Differential analysis of transcript expression indicated that at an early stage of infestation miRNA expression was higher in head kidney than in skin tissue, revealing tissue-specific expression patterns. In parallel, miRNA target prediction using 3'UTRs from highly regulated immune-related genes and iron metabolism showed that mir-140-4 and mir-181a-2-5 modulate the expression of TLR22 and Aminolevulinic acid synthase, respectively. This study contributes knowledge about the immune response of Atlantic salmon during infestation with sea lice.

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.

Transcriptional responses in the hepatopancreas of Eriocheir sinensis exposed to deltamethrin

Deltamethrin is an important pesticide widely used against ectoparasites. Deltamethrin contamination has resulted in a threat to the healthy breeding of the Chinese mitten crab, Eriocheir sinensis. In this study, we investigated transcriptional responses in the hepatopancreas of E. sinensis exposed to deltamethrin. We obtained 99,087,448, 89,086,478, and 100,117,958 raw sequence reads from control 1, control 2, and control 3 groups, and 92,094,972, 92,883,894, and 92,500,828 raw sequence reads from test 1, test 2, and test 3 groups, respectively. After filtering and quality checking of the raw sequence reads, our analysis yielded 79,228,354, 72,336,470, 81,859,826, 77,649,400, 77,194,276, and 75,697,016 clean reads with a mean length of 150 bp from the control and test groups. After deltamethrin treatment, a total of 160 and 167 genes were significantly upregulated and downregulated, respectively. Gene ontology terms "biological process," "cellular component," and "molecular function" were enriched with respect to cell killing, cellular process, other organism part, cell part, binding, and catalytic. Pathway analysis using the Kyoto Encyclopedia of Genes and Genomes showed that the metabolic pathways were significantly enriched. We found that the CYP450 enzyme system, carboxylesterase, glutathione-S-transferase, and material (including carbohydrate, lipid, protein, and other substances) metabolism played important roles in the metabolism of deltamethrin in the hepatopancreas of E. sinensis. This study revealed differentially expressed genes related to insecticide metabolism and detoxification in E. sinensis for the first time and will help in understanding the toxicity and molecular metabolic mechanisms of deltamethrin in E. sinensis.

Integrated analysis of mRNA-seq in the haemocytes of Eriocheir sinensis in response to Spiroplasma eriocheiris infection

The Chinese mitten crab Eriocheir sinensis is an important economic crustacean that has been exposed to various diseases. Spiroplasma eriocheiris, isolated from tremor-diseased E. sinensis, was first identified as a lethal pathogen of freshwater crustaceans. To understand the pathogenesis of S. eriocheiris to E. sinensis, the transcriptomic profiles of haemocytes in the experimental and control groups at 1 d and 7 d post-injection were obtained using Illumina HiSeq 2500. These results showed that 40,358,724, 44,462,112, 45,516,576 and 37,713,728 paired-end clean reads were obtained from the cDNA libraries of DZ1 (the control group at 1 d), DZ7 (the control group at 7 d), SY1 (the experimental group at 1 d) and SY7 (the experimental group at 7 d), respectively. In total, 106,641 unique transcript fragments (unigenes) were assembled, with an average length of 710 bp. On the first day of stimulation, 33,084 up-regulated transcripts and 19,208 down-regulated transcripts were found in the experimental group compared with those in the control group. On the seventh day of stimulation, 40,198 up-regulated transcripts and 12,032 down-regulated transcripts were found in the experimental group compared with those in the control group. Some canonical immune-related pathways were identified via KEGG pathway analysis, including complement and coagulation cascades, the VEGF signalling pathway, the Wnt signalling pathway, natural killer cell-mediated cytotoxicity, the MAPK signalling pathway, neuroactive ligand-receptor interactions, and the Lysosome pathway. We found important immune-related genes (GNPTAB, MASP2, F7, F5, NFATC, TRAF6, MAP3K5, and TRa) in the KEGG pathway, and those genes were confirmed by qRT-PCR analysis. In addition, the significantly enriched neuroactive ligand-receptor interaction pathway was associated with intense paroxysmal tremors of infected crabs. Our results provide valuable information for the further analysis of the mechanisms of E. sinensis defence against S. eriocheiris invasion.

Interaction of heat shock protein 60 (HSP60) with microRNA in Chinese mitten crab during Spiroplasma eriocheiris infection

Heat shock protein 60 from the Chinese mitten crab Eriocheir sinensis (EsHSP60) was previously identified in relation to Spiroplasma eriocheiris infection by isobaric tags for relative and absolute quantitation labelling followed by liquid chromatography-tandem mass spectrometry. In the present study, to validate the immune function of this protein, the cDNA of the EsHSP60 gene was cloned. Various crab tissues were assessed using real-time PCR, which showed that EsHSP60 transcription occurred in all tissues examined. The expression profiles of EsHSP60 in haemolymph at transcription and protein levels when infected with S. eriocheiris were investigated by real-time PCR and Western blot analysis, respectively. A significant increase of EsHSP60 transcription and protein expression appeared post-injection in response to S. eriocheiris infection when compared to the control group. The double-luciferase reporter gene assay showed that the microRNA PC-533-3p interacted with the 3'-untranslated region of EsHSP60 and inhibited the translation of EsHSP60. The expression profiles of PC-533-3p during S. eriocheiris infection were also investigated by real-time PCR. However, the change tendency of PC-533-3p was opposite to that of the EsHSP60 after S. eriocheiris challenge. These data indicate that the EsHSP60 proteins may play an important role in mediating the immune responses of E. sinensis to an S. eriocheiris challenge.

Differential expression of microRNAs of Litopenaeus vannamei in response to different virulence WSSV infection

WSSV is one of the most harmful pathogeny in the pacific white shrimp, and genetic variations caused the strains of different virulence. MicroRNAs (miRNAs) involved in the regulation of virus defense. To understand the different virulence of WSSV on miRNA expression in Litopeneaus vannamei, the deep sequencing was performed to compare two small RNA libraries prepared from hepatopancreas of Litopeneaus vannamei infected with normal-virulence or low-virulence WSSV. Approximately 29,398,623 raw reads from normal-virulence library and 35,291,803 raw reads from low-virulence library were obtained. There were about 37 miRNAs homologs identified. Sixteen miRNAs were significantly up-regulated and twenty-one miRNAs were significantly down-regulated in normal-virulence infection library compared with low-virulence infection library. Of these, Igi-miR-1175-3p was the most significant different miRNA, followed by bmo-miR-1175-3p and ipu-miR-26b, respectively. The putative target genes for differentially expressed miRNAs were concerned with biological processes, signal meditated, cell differentiation and apoptosis, immune recognition and other more functions. The results will help to understand the miRNAs response to different virulence WSSV infection.

Comprehensive identification and profiling of Nile tilapia (Oreochromis niloticus) microRNAs response to Streptococcus agalactiae infection through high-throughput sequencing

MicroRNAs are a kind of small non-coding RNAs that participate in various biological processes. Deregulated microRNA expression is associated with several types of diseases. Tilapia (Oreochromis niloticus) is an important commercial fish species in China. To identify miRNAs and investigate immune-related miRNAs of O. niloticus, we applied high-throughput sequencing technology to identify and analyze miRNAs from tilapia infected with Streptococcus agalactiae at a timescale of 72 h divided into six different time points. The results showed that a total of 3009 tilapia miRNAs were identified, including in 1121 miRNAs which have homologues in the currently available databases and 1878 novel miRNAs. The expression levels of 218 tilapia miRNAs were significantly altered at 6 h-72 h post-bacterial infection (pi), and these miRNAs were therefore classified as differentially expressed tilapia miRNAs. For the 1121 differentially expressed tilapia miRNAs target 41961 genes. GO and KEGG enrichment analysis revealed that some target genes of tilapia miRNAs were grouped mainly into the categories of apoptotic process, signal pathway, and immune response. This is the first report of comprehensive identification of O. niloticus miRNAs being differentially regulated in spleen in normal conditions relating to S. agalactiae infection. This work provides an opportunity for further understanding of the molecular mechanisms of miRNA regulation in O. niloticus host-pathogen interactions.

MicroRNAs of the mesothorax in Qinlingacris elaeodes, an alpine grasshopper showing a wing polymorphism with unilateral wing form

MicroRNAs (miRNAs) are now recognized as key post-transcriptional regulators in regulation of phenotypic diversity. Qinlingacris elaeodes is a species of the alpine grasshopper, which is endemic to China. Adult individuals have three wing forms: wingless, unilateral-winged and short-winged. This is an ideal species to investigate the phenotypic plasticity, development and evolution of insect wings because of its case of unilateral wing form in both the sexes. We sequenced a small RNA library prepared from mesothoraxes of the adult grasshoppers using the Illumina deep sequencing technology. Approximately 12,792,458 raw reads were generated, of which the 854,580 high-quality reads were used only for miRNA identification. In this study, we identified 49 conserved miRNAs belonging to 41 families and 69 species-specific miRNAs. Moreover, seven miRNA*s were detected both for conserved miRNAs and species-specific miRNAs, which were supported by hairpin forming precursors based on polymerase chain reaction. This is the first description of miRNAs in alpine grasshoppers. The results provide a useful resource for further studies on molecular regulation and evolution of miRNAs in grasshoppers. These findings not only enrich the miRNAs for insects but also lay the groundwork for the study of post-transcriptional regulation of wing forms.

High throughput sequencing of small RNAs transcriptomes in two Crassostrea oysters identifies microRNAs involved in osmotic stress response

Increasing evidence suggests that microRNAs post-transcriptionally regulate gene expression and are involved in responses to biotic and abiotic stress. However, the role of miRNAs involved in osmotic plasticity remains largely unknown in marine bivalves. In the present study, we performed low salinity challenge with two Crassostrea species (C. gigas and C. hongkongensis), and conducted high-throughput sequencing of four small RNA libraries constructed from the gill tissues. A total of 202 and 87 miRNAs were identified from C. gigas and C. hongkongensis, respectively. Six miRNAs in C. gigas and two in C. hongkongensis were differentially expressed in response to osmotic stress. The expression profiles of these eight miRNAs were validated by qRT-PCR. Based on GO enrichment and KEGG pathway analysis, genes associated with microtubule-based process and cellular component movement were enriched in both species. In addition, five miRNA-mRNA interaction pairs that showed opposite expression patterns were identified in the C. hongkongensis, Differential expression analysis identified the miRNAs that play important regulatory roles in response to low salinity stress, providing insights into molecular mechanisms that are essential for salinity tolerance in marine bivalves.

In silico genome wide mining of conserved and novel miRNAs in the brain and pineal gland of Danio rerio using small RNA sequencing data

MicroRNAs (miRNAs) are small, non-coding RNA molecules that bind to the mRNA of the target genes and regulate the expression of the gene at the post-transcriptional level. Zebrafish is an economically important freshwater fish species globally considered as a good predictive model for studying human diseases and development. The present study focused on uncovering known as well as novel miRNAs, target prediction of the novel miRNAs and the differential expression of the known miRNA using the small RNA sequencing data of the brain and pineal gland (dark and light treatments) obtained from NCBI SRA. A total of 165, 151 and 145 known zebrafish miRNAs were found in the brain, pineal gland (dark treatment) and pineal gland (light treatment), respectively. Chromosomes 4 and 5 of zebrafish reference assembly GRCz10 were found to contain maximum number of miR genes. The miR-181a and miR-182 were found to be highly expressed in terms of number of reads in the brain and pineal gland, respectively. Other ncRNAs, such as tRNA, rRNA and snoRNA, were curated against Rfam. Using GRCz10 as reference, the subsequent bioinformatic analyses identified 25, 19 and 9 novel miRNAs from the brain, pineal gland (dark treatment) and pineal gland (light treatment), respectively. Targets of the novel miRNAs were identified, based on sequence complementarity between miRNAs and mRNA, by searching for antisense hits in the 3'-UTR of reference RNA sequences of the zebrafish. The discovery of novel miRNAs and their targets in the zebrafish genome can be a valuable scientific resource for further functional studies not only in zebrafish but also in other economically important fishes.

Identification and Characterization of MicroRNAs in Snakehead Fish Cell Line upon Snakehead Fish Vesiculovirus Infection

MicroRNAs (miRNAs) play important roles in mediating multiple biological processes in eukaryotes and are being increasingly studied to evaluate their roles associated with cellular changes following viral infection. Snakehead fish Vesiculovirus (SHVV) has caused mass mortality in snakehead fish during the past few years. To identify specific miRNAs involved in SHVV infection, we performed microRNA deep sequencing on a snakehead fish cell line (SSN-1) with or without SHVV infection. A total of 205 known miRNAs were identified when they were aligned with the known zebrafish miRNAs, and nine novel miRNAs were identified using MiRDeep2 software. Eighteen and 143 of the 205 known miRNAs were differentially expressed at three and 24 h post-infection (poi), respectively. From the differentially-expressed miRNAs, five were randomly selected to validate their expression profiles using quantitative reverse transcription polymerase chain reaction (qRT-PCR), and their expression profiles were consistent with the microRNA sequencing results. In addition, the target gene prediction of the SHVV genome was performed for the differentially-expressed host miRNAs, and a total of 10 and 58 differentially-expressed miRNAs were predicted to bind to the SHVV genome at three and 24 h poi, respectively. The effects of three selected miRNAs (miR-130-5p, miR-214 and miR-216b) on SHVV multiplication were evaluated using their mimics and inhibitors via qRT-PCR and Western blotting. The results showed that all three miRNAs were able to inhibit the multiplication of SHVV; whereas the mechanisms underlying the SHVV multiplication inhibited by the specific miRNAs need to be further characterized in the future.

The identification and characteristics of salinity-related microRNAs in gills of Portunus trituberculatus

MicroRNAs (miRNAs) are small noncoding RNAs that post-transcriptionally regulate gene expression in organisms. To understand the underlying mechanisms behind the molecular response of the crab to low salt-stress, high-throughput Illumina/Solexa deep sequencing technology was used to investigate the expression profiles of miRNAs under low salinity challenged. Two mixed RNA pool libraries of gill tissues from low salinity challenged (LC) and the control groups (NC) were sequenced on the Illumina platform. A total of 6,166,057 and 7,032,973 high-quality reads were obtained from the NC and LC libraries, respectively. Sixty-seven miRNAs consisting of 16 known and 51 novel ones were identified, among which, 12 miRNAs were differentially expressed in LC compared to NC. Thirty-four of the target genes predicted were differentially expressed in the opposite direction to the miRNAs, which were involved in crucial processes related to osmoregulation by gene ontology (GO) functional enrichment analysis, such as anion transport processes (GO:0006820) and chitin metabolic process (GO:0006030). These results provide a basis for further investigation of the miRNA-modulating networks in osmoregulation of Portunus trituberculatus.

Identification of microRNAs in the Toxigenic Dinoflagellate Alexandrium catenella by High-Throughput Illumina Sequencing and Bioinformatic Analysis

Micro-ribonucleic acids (miRNAs) are a large group of endogenous, tiny, non-coding RNAs consisting of 19–25 nucleotides that regulate gene expression at either the transcriptional or post-transcriptional level by mediating gene silencing in eukaryotes. They are considered to be important regulators that affect growth, development, and response to various stresses in plants. Alexandrium catenella is an important marine toxic phytoplankton species that can cause harmful algal blooms (HABs). To date, identification and function analysis of miRNAs in A. catenella remain largely unexamined. In this study, high-throughput sequencing was performed on A. catenella to identify and quantitatively profile the repertoire of small RNAs from two different growth phases. A total of 38,092,056 and 32,969,156 raw reads were obtained from the two small RNA libraries, respectively. In total, 88 mature miRNAs belonging to 32 miRNA families were identified. Significant differences were found in the member number, expression level of various families, and expression abundance of each member within a family. A total of 15 potentially novel miRNAs were identified. Comparative profiling showed that 12 known miRNAs exhibited differential expression between the lag phase and the logarithmic phase. Real-time quantitative RT-PCR (qPCR) was performed to confirm the expression of two differentially expressed miRNAs that were one up-regulated novel miRNA (aca-miR-3p-456915), and one down-regulated conserved miRNA (tae-miR159a). The expression trend of the qPCR assay was generally consistent with the deep sequencing result. Target predictions of the 12 differentially expressed miRNAs resulted in 1813target genes. Gene ontology (GO) analysis and the Kyoto Encyclopedia of Genes and Genomes pathway database (KEGG) annotations revealed that some miRNAs were associated with growth and developmental processes of the alga. These results provide insights into the roles that miRNAs play in the growth of A. catenella, and they provide the basis for further studies of the molecular mechanisms that underlie bloom growth in red tides species.

Hepatitis From Spiroplasma sp. in an Immunocompromised Patient

A 70-year-old lung transplant recipient patient was admitted with fever, nausea, abdominal pain, peripheral edema and pronounced weakness. An initial work-up for presumed infection revealed cholestatic hepatitis, leukocytosis and thrombocytopenia, but failed to detect a pathogen. An increased glucose uptake exclusively in the liver was demonstrated by positron emission tomography. Liver biopsy showed basophilic inclusions in the cytoplasm of hepatocytes. Broad- range 16S rRNA gene PCR followed by sequence analysis yielded Spiroplasma sp. in two independent blood samples and the liver biopsy, confirming Spiroplasma sp. as the causative agent. Antibiotic treatment with doxycycline and azithromycin led to complete recovery.

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.

Expression profiling analysis of the microRNA response of Cynoglossus semilaevis to Vibrio anguillarum and other stimuli

To investigate the roles of microRNAs (miRNA) of Cynoglossus semilaevis in response to Vibrio anguillarum that were previously identified using high-throughput sequencing, microarray analyses was performed on three small RNA libraries (CG, NOSG, and HOSG) prepared from C. semilaevis immune tissues. In total, of 1279 designed probes, 739 (57.78 %) were detectable. The expression levels of these miRNAs were analyzed using pairwise comparisons among the three libraries, and a total of 99 miRNAs were observed to be significantly differentially expressed. The expression patterns of 10 differentially expressed miRNAs were validated by real-time quantitative PCR (RT-qPCR). In addition, expression of miR-142-5p, miR-223, and miR-181a in response to V. anguillarum at numerous time-points in four tissues, as well as the responses to lipopolysaccharide (LPS), polyinosinic:polycytidylic acid (poly I:C), peptidoglycan (PGN), and red-spotted grouper nervous necrosis virus (RGNNV) in head kidney cells, were studied by qRT-PCR. Taken together, all of the expression profiles showed significant differences compared to the control group; both similarities and differences in the expression responses to the same pathogen were observed. Collectively, these findings highlighted the putative roles for miRNAs in the context of the innate immune response of C. semilaevis exposing to pathogens and that further studies are needed to understand the molecular mechanisms of miRNA regulation in C. semilaevis host-pathogen interactions.

Identification and analysis of differentially-expressed microRNAs in Japanese encephalitis virus-infected PK-15 cells with deep sequencing

Japanese encephalitis virus (JEV), a mosquito-borne Flavivirus, causes acute viral encephalitis with high morbidity and mortality in humans and animals. MicroRNAs (miRNAs) are small noncoding RNAs that are important modulators of the intricate host-pathogen interaction networks. However, our knowledge of the changes that occur in miRNAs in host cells after JEV infection is still limited. To understand the molecular pathogenesis of JEV at the level of posttranscriptional regulation, we used Illumina deep sequencing to sequence two small RNA libraries prepared from PK-15 cells before and after JEV infection. We identified 522 and 427 miRNAs in the infected and uninfected cells, respectively. Overall, 132 miRNAs were expressed significantly differently after challenge with JEV: 78 were upregulated and 54 downregulated. The sequencing results for selected miRNAs were confirmed with RT-qPCR. GO analysis of the host target genes revealed that these dysregulated miRNAs are involved in complex cellular pathways, including the metabolic pathway, inflammatory response and immune response. To our knowledge, this is the first report of the comparative expression of miRNAs in PK-15 cells after JEV infection. Our findings will underpin further studies of miRNAs' roles in JEV replication and identify potential candidates for antiviral therapies against JEV.

MicroRNA profile analysis of Epithelioma papulosum cyprini cell line before and after SVCV infection

MicroRNAs (miRNAs) play significant roles in regulating almost all of the biological processes in eukaryotes. An accumulating body of evidence shows that miRNAs are associated with cellular changes following viral infection. Spring viremia of carp virus (SVCV) is the pathogen of Spring viremia of carp (SVC), which results in heavy losses in the cultured common carp (Cyprinus carpio) industry in many countries. To study the involvement of miRNAs during SVCV infection, we adopted the Solexa sequencing technology to sequence small RNA libraries from the Epithelioma papulosum cyprini (EPC) cell line before and after infection with SVCV. In this study, a total of 161 conserved and 26 novel miRNAs were identified. Subsequently, the expression patterns of these miRNAs were compared between the uninfected (control library, M) and SVCV-infected (infection library, E) libraries. In addition, to verify the Solexa sequencing results, the expression patterns of 14 randomly selected miRNAs were validated by qRT-PCR. The targets of the significantly differentially expressed miRNAs were then predicted, and the miRNAs that could directly target the SVCV genome were also predicted. No miRNA encoded by SVCV itself was detected. To the best of our knowledge, this study presents the first miRNA profiling assessment in association with fish rhabdovirus infection, and the data presented lay a foundation for further investigations to determine the roles of miRNAs in regulating the molecular mechanism during SVCV infection.

Transcriptome and expression profiling analysis of the hemocytes reveals a large number of immune-related genes in mud crab Scylla paramamosain during Vibrio parahaemolyticus infection

BACKGROUND

Mud crab Scylla paramamosain is an economically important marine species in China. However, frequent outbreaks of infectious diseases caused by marine bacteria, such as Vibrio parahaemolyticus, result in great economic losses.

METHODOLOGY/PRINCIPAL FINDINGS

Comparative de novo transcriptome analysis of S. paramamosain infected with V. parahaemolyticus was carried out to investigate the molecular mechanisms underlying the immune response to pathogenic bacteria by using the Illumina paired-end sequencing platform. A total of 52,934,042 clean reads from the hemocytes of V. parahaemolyticus-infected mud crabs and controls were obtained and assembled into 186,193 contigs. 59,120 unigenes were identified from 81,709 consensus sequences of mud crabs and 48,934 unigenes were matched proteins in the Nr or Swissprot databases. Among these, 10,566 unigenes belong to 3 categories of Gene Ontology, 25,349 to 30 categories of KEGG, and 15,191 to 25 categories of COG database, covering almost all functional categories. By using the Solexa/Illumina's DGE platform, 1213 differentially expressed genes (P<0.05), including 538 significantly up-regulated and 675 down-regulated, were detected in V. parahaemolyticus-infected crabs as compared to that in the controls. Transcript levels of randomly-chosen genes were further measured by quantitative real-time PCR to confirm the expression profiles. Many differentially expressed genes are involved in various immune processes, including stimulation of the Toll pathway, Immune Deficiency (IMD) pathway, Ras-regulated endocytosis, and proPO-activating system.

CONCLUSIONS/SIGNIFICANCE

Analysis of the expression profile of crabs under infection provides invaluable new data for biological research in S. paramamosain, such as the identification of novel genes in the hemocytes during V. parahaemolyticus infection. These results will facilitate our comprehensive understanding of the mechanisms involved in the immune response to bacterial infection and will be helpful for diseases prevention in crab aquaculture.

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.

Identification and characterization of Cynoglossus semilaevis microRNA response to Vibrio anguillarum infection through high-throughput sequencing

MicroRNAs (miRNA) play key regulatory roles in diverse biological processes. Cynoglossus semilaevis is an important commercial mariculture fish species in China. To identify miRNAs and investigate immune-related miRNAs of C. semilaevis, we performed high-throughput sequencing on three small RNA libraries prepared from C. semilaevis immune tissues (liver, head kidney, spleen, and intestine). One library was prepared under normal conditions (control, CG); two were prepared during Vibrio anguillarum infection, where vibriosis symptoms were obvious and non-obvious (HOSG and NOSG, respectively). We obtained 11,216,875, 12,313,404, and 11,398,695 clean reads per library, respectively. Bioinformatic analysis identified 452 miRNAs, including 24 putative novel miRNAs. We analyzed differentially expressed miRNAs between two libraries using pairwise comparison. For NOSG-CG, there was significant differential expression of 175 (38.72%) miRNAs. There was significant differential expression of 215 (47.57%) miRNAs between HOSG and CG. Compared with CG, The HOSG-NOSG comparison revealed significantly different expression of 122 (26.99%) miRNAs respectively. Real-time quantitative PCR (RT-qPCR) experiments were performed for 10 miRNAs of the three samples, and agreement was found between the sequencing and RT-qPCR data. For miRNAs that were significantly differentially expressed, functional annotation of target genes by Gene Ontology enrichment and Kyoto Encyclopedia of Genes and Genomes pathway analysis indicated that a set of miRNAs that were expressed highly abundantly and significantly differentially were might involved in immune system development and immune response. To our understanding, this is the first report of comprehensive identification of C. semilaevis miRNAs being differentially regulated in immune tissues (liver, head kidney, spleen, and intestine) in normal conditions relating to V. anguillarum infection. Many miRNAs were differentially regulated upon pathogen exposure. This work provides an opportunity for further understanding of the molecular mechanisms of miRNA regulation in C. semilaevis host-pathogen interactions.

Bioinformatics analysis of hemocyte miRNAs of scallop Chlamys farreri against acute viral necrobiotic virus (AVNV)

The sustainable development of the scallop Chlamys farreri industry in China is hindered by mass mortality mainly caused by a novel pathogen known as acute viral necrosis virus (AVNV). A better understanding of host-virus interactions, especially those at the molecular level, may facilitate the prevention and cure of AVNV infections. MicroRNAs (miRNAs) represent a class of small RNA molecules involved in several biological processes, including mediating host-pathogen responses. In this study, two hemocyte small RNA libraries were constructed from control (control library, CL) and AVNV-infected (infection library, IL) C. farreri for high throughput sequencing using Solexa technology. Acquired data were further used to identify conserved and novel miRNAs, screen differentially expressed miRNAs, and predict their target genes through bioinformatics analysis. Solexa sequencing produced 19,485,719 and 20,594,513 clean reads representing 2,248,814 and 1,510,256 unique small RNAs from CL and IL, respectively. A total of 57 conserved miRNAs were identified in both libraries, among which only two were unique to IL. Novel miRNA prediction using the Crassostrea gigas genome as a reference revealed 11 candidate miRNAs, 10 of which were validated by RT-PCR. Differential expression (p < 0.001) between libraries was observed in 37 miRNAs, among which 30 and 7 miRNAs were up- and downregulated, respectively. Expression differences were further confirmed by qRT-PCR. A sequence homology search against available C. farreri ESTs showed that these differentially expressed miRNAs may target 177 genes involved in a broad range of biological processes including immune defense and stress response. This study is the first to characterize C. farreri miRNAs and miRNA expression profiles in response to AVNV infection by deep sequencing. The results presented here will deepen our understanding of the pathogenesis of AVNV at the molecular level and provide new insights into the molecular basis of host-pathogen interactions in C. farreri.

The identification and characteristics of immune-related microRNAs in haemocytes of oyster Crassostrea gigas

Background

MicroRNAs (miRNAs) represent a class of small ncRNAs that repress gene expression on the post-transcriptional level by the degradation or translation inhibition of target mRNA.

Methodology

Three small RNA libraries from oyster haemocytes were sequenced on the Illumina platform to investigate the latent immunomodulation of miRNAs after bacteria challenge and heat stress. Totally, 10,498,663, 8,588,606 and 9,679,663 high-quality reads were obtained in the control, bacteria and bacteria+heat library, respectively, from which 199 oyster miRNAs including 71 known and 128 novel ones were identified. Among these miRNAs, 6 known and 23 novel ones were predicted to possess more than one precursor-coding region, and cgi-miR-10a, cgi-miR-184b, cgi-miR-100, cgi-miR-1984 and cgi-miR-67a were observed to be the most abundant miRNAs in the control library. The expression levels of 22 miRNAs in the bacteria library were significantly higher than those in the control library, while there were another 33 miRNAs whose expression levels were significantly lower than that in the control library. Meanwhile, the expression levels of 65 miRNAs in the bacteria+heat library changed significantly compared to those in the bacteria library. The target genes of these differentially expressed miRNAs were annotated, and they fell in immune and stress-related GO terms including antioxidant, cell killing, death, immune system process, and response to stimulus. Furthermore, there were 42 differentially expressed miRNAs detected in both control/bacteria and bacteria/bacteria+heat comparisons, among which 9 miRNAs displayed the identical pattern in the two comparisons, and the expression alterations of 8 miRNAs were confirmed using quantitative RT-PCR.

Conclusions

These results indicated collectively that immune challenge could induce the expression of immune-related miRNAs, which might modulate the immune response such as redox reaction, phagocytosis and apoptosis, and the expression of some immune-related miRNAs could be also regulated by heat stress to improve the environmental adaption of oyster.

Transcriptome profiling analysis on whole bodies of microbial challenged Eriocheir sinensis larvae for immune gene identification and SNP development

To study crab immunogenetics of individuals, newly hatched Eriocheir sinensis larvae were stimulated with a mixture of three pathogen strains (Gram-positive bacteria Micrococcus luteus, Gram-negative bacteria Vibrio alginolyticus and fungi Pichia pastoris; 10⁸ cfu·mL^-1). A total of 44,767,566 Illumina clean reads corresponding to 4.52 Gb nucleotides were generated and assembled into 100,252 unigenes (average length: 1,042 bp; range: 201-19,357 bp). 17,097 (26.09%) of 65,535 non-redundant unigenes were annotated in NCBI non-redundant protein (Nr) database. Moreover, 23,188 (35.38%) unigenes were assigned to three Gene Ontology (GO) categories, 15,071 (23.00%) to twenty-six Clusters of orthologous Groups (COG) and 8,574 (13.08%) to six Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, respectively. Numerous genes were further identified to be associated with multiple immune pathways, including Toll, immune deficiency (IMD), janus kinase (JAK)-signal transducers and activators of transcription (STAT) and mitogen-activated protein kinase (MAPK) pathways. Some of them, such as tumor necrosis factor receptor associated factor 6 (TRAF6), fibroblast growth factor (FGF), protein-tyrosine phosphatase (PTP), JNK-interacting protein 1 (JIP1), were first identified in E. sinensis. TRAF6 was even first discovered in crabs. Additionally, 49,555 single nucleotide polymorphisms (SNPs) were developed from over 13,309 unigenes. This is the first transcriptome report of whole bodies of E. sinensis larvae after immune challenge. Data generated here not only provide detail information to identify novel genes in genome reference-free E. sinensis, but also facilitate our understanding on host immunity and defense mechanism of the crab at whole transcriptome level.

Characterization of microRNAs in mud crab Scylla paramamosain under Vibrio parahaemolyticus infection

BACKGROUND

Infection of bacterial Vibrio parahaemolyticus is common in mud crab farms. However, the mechanisms of the crab's response to pathogenic V. parahaemolyticus infection are not fully understood. MicroRNAs (miRNAs) are a class of small noncoding RNAs that function as regulators of gene expression and play essential roles in various biological processes. To understand the underlying mechanisms of the molecular immune response of the crab to the pathogens, high-throughput Illumina/Solexa deep sequencing technology was used to investigate the expression profiles of miRNAs in S. paramamosain under V. parahaemolyticus infection.

METHODOLOGY/PRINCIPAL FINDINGS

Two mixed RNA pools of 7 tissues (intestine, heart, liver, gill, brain, muscle and blood) were obtained from V. parahaemolyticus infected crabs and the control groups, respectively. By aligning the sequencing data with known miRNAs, we characterized 421 miRNA families, and 133 conserved miRNA families in mud crab S. paramamosain were either identical or very similar to existing miRNAs in miRBase. Stem-loop qRT-PCRs were used to scan the expression levels of four randomly chosen differentially expressed miRNAs and tissue distribution. Eight novel potential miRNAs were confirmed by qRT-PCR analysis and the precursors of these novel miRNAs were verified by PCR amplification, cloning and sequencing in S. paramamosain. 161 miRNAs (106 of which up-regulated and 55 down-regulated) were significantly differentially expressed during the challenge and the potential targets of these differentially expressed miRNAs were predicted. Furthermore, we demonstrated evolutionary conservation of mud crab miRNAs in the animal evolution process.

CONCLUSIONS/SIGNIFICANCE

In this study, a large number of miRNAs were identified in S. paramamosain when challenged with V. parahaemolyticus, some of which were differentially expressed. The results show that miRNAs might play some important roles in regulating gene expression in mud crab under V. parahaemolyticus infection, providing a basis for further investigation of miRNA-modulating networks in innate immunity of mud crab.