De novo assembly and transcriptome analysis of osmoregulation in Litopenaeus vannamei under three cultivated conditions with different salinities

Effects of myo-inositol on growth and biomarkers of environmental stress and metabolic regulation in Pacific white shrimp (Litopenaeus vannamei) reared at low salinity

This study explored the role of myo-inositol in alleviating the low salinity stress of White Shrimp (Litopenaeus vannamei). Juvenile shrimp (0.4 ± 0.02 g) in low salinity (salinity 3) water were fed diets with myo-inositol levels of 0, 272, 518, 1020 and 1950 mg/kg (crude protein is 36.82 %, crude lipid is 7.58 %), fed shrimp in seawater at a salinity of 25 were fed a 0 mg/kg myo-inositol diet as a control (Ctrl). The experiment was carried out in tanks (50 L) with satiety feeding, and the experiment lasted for 6 weeks. After sampling, the serum was used to measure immune function, the hepatopancreas homogenate was used to measure the antioxidant capacity and hepatopancreas damage state, the hepatopancreas was used for transcriptomics analysis, and the gills were used for qPCR to measure osmotic pressure regulation. The results showed that the final weight and survival of the shrimp in the 1020 mg/kg group increased significantly compared with those in the other low salinity groups, but the final weight and biomass increase were significantly lower than those in the Ctrl group. Dietary myo-inositol improved the antioxidant capacity of shrimp under low salinity. B-cell hyperplasia and hepatic duct damage were observed in the hepatopancreas in the 0 mg/kg group. The results of transcriptome analysis showed that myo-inositol could participate in the osmotic pressure regulation of shrimp by regulating carbohydrate metabolism, amino acid metabolism, lipid metabolism and other related genes. Myo-inositol significantly affected the expression of related genes in ion transporter and G protein-coupled receptor-mediated pathways. This study demonstrated that myo-inositol can not only act as an osmotic pressure effector and participate in the osmolar regulation of shrimp through the phosphatidylinositol signaling pathway mediated by G protein-coupled receptors but also relieve low salinity stress by mediating physiological pathways such as immunity, antioxidation, and metabolism in shrimp. The binomial regression analysis of biomass increases and survival showed that the appropriate amount of myo-inositol in the L. vannamei diet was 862.50-1275.00 mg/kg under low salinity.

Analysis of adaptive molecular mechanisms in response to low salinity in antennal gland of mud crab, Scylla paramamosain

As an important marine aquaculture species, the mud crab (Scylla paramamosain) is a good candidate for studying the osmoregulatory mechanism of crustaceans. While previous studies have focused on the osmoregulatory function of the gills, this study aims to explore the osmoregulatory function of the antennal glands. By the comparative transcriptomic analysis, we found the pathways of ion regulation including "proximal tubule bicarbonate reclamation" and "mineral absorption" were activated in the antennal glands of the crabs long-term dwelling in low salinity. The enhanced ionic reabsorption was associated with up-regulated ion transport genes such as NKA, CA-c, VPA, and NHE, and with energy metabolism genes such as MDH, SLC25, and PEPCK. The upregulation of NKA and CA-c was also verified by the increased enzyme activity. The lowered osmolality and ion concentration of the hemolymph and the enlarged labyrinth lumen and hemolymph capillary inside the antennal glands indicated the infiltration of external water and the responsively increase of urine excretion, which explained the requirement of enhanced ionic reabsorption. To further confirm these findings, we examined the change of gene expression, enzyme activity, internal ion concentration, and external ion concentration during a 96 h low salinity challenge with seven intervals. The results were basically consistent with the results as shown in the long-term low salinity adaptation. The present study provides valuable information on the osmoregulatory function of the antennal glands of S. paramamosain. The implication of this study in marine aquaculture is that it provides valuable information on the osmoregulatory mechanism of mud crabs, which can be used to improve their culture conditions and enhance their tolerance to salinity stress. The identified genes and pathways involved in osmoregulation can also be potential targets for genetic selection and breeding programs to develop more resilient mud crab strains for aquaculture.

Growth, immunity and transcriptome response to different stocking densities in Litopenaeus vannamei

The effects of different stocking densities on Litopenaeus vannamei were investigated from the aspects of growth performance, immune response and transcriptome in this experiment. L. vannamei (initial body weight: 0.30 ± 0.02 g) were reared for 8 weeks at three stocking densities of 100 (LSD), 200 (MSD) and 300 (HSD) shrimp/m³, respectively. The results showed that the survival rate (SR), final body weight (FBW), weight gain rate (WGR), specific growth ratio (SGR) and protein efficiency ratio (PER) of L. vannamei significantly decreased, while the feed factor (FCR) significantly increased with the increase of stocking density. After Vibrio parahemolyticus infection, the SR of L. vannamei in the HSD group was significantly lower than that in the LSD and MSD groups. Increasing stocking density significantly increased the activities of aspartate aminotransferase (AST), alanine aminotransferase (ALT) and lysozyme (LYS) while significantly decreased the activities of catalase (CAT) and phenol oxidase (PO) in the serum of L. vannamei. Similar changes of the gene expression as the activities of immune enzymes were found in the hemocytes. Pairwise comparison between the LSD, MSD and HSD group in the transcriptome analysis identified that there were 304, 1376 and 2083 differentially expressed genes (DEGs) in LSD vs MSD, MSD vs HSD and LSD vs HSD, respectively. Among them, most of the immune-related DEGs were down-regulated and metabolism-related DEGs were up-regulated with the increasing stocking density. In addition, KEGG enrichment pathway analysis revealed that several immune and metabolic related pathways including PI3K-Akt signaling pathway and AMPK signaling pathway were significantly enriched. Of these, the PI3K-Akt signaling pathway had the most DEGs and was also the most significantly enriched pathway. Furthermore, 16 DEGs (such as FOXO, PCK2 and CTSC, etc.) and partial immune enzyme activity (such as AST, CAT and PO, etc.) changes were closely correlated with the increase of stocking density when partial immune-related DEGs and immune-related enzymes were analyzed jointly. All these results indicated that changes in stocking density had a significant effect on the growth performance, immunity and transcriptome of L. vannamei.

RNA-Seq comparative study reveals molecular effectors linked to the resistance of Pinna nobilis to Haplosporidium pinnae parasite

With the intensification of maritime traffic, recently emerged infectious diseases have become major drivers in the decline and extinction of species. Since 2016, mass mortality events have decimated the endemic Mediterranean Sea bivalve Pinna nobilis, affecting ca. 100% of individuals. These events have largely been driven by Haplosporidium pinnae's infection, an invasive species which was likely introduced by shipping. While monitoring wild populations of P. nobilis, we observed individuals that survived such a mass mortality event during the summer of 2018 (France). We considered these individuals resistant, as they did not show any symptoms of the disease, while the rest of the population in the area was devastated. Furthermore, the parasite was not detected when we conducted a PCR amplification of a species-specific fragment of the small subunit ribosomal DNA. In parallel, the transcriptomic analysis showed evidence of some parasite RNA indicating that the resistant individuals had been exposed to the parasite without proliferating. To understand the underlying mechanisms of resistance in these individuals, we compared their gene expression with that of susceptible individuals. We performed de novo transcriptome assembly and annotated the expressed genes. A comparison of the transcriptomes in resistant and susceptible individuals highlighted a gene expression signature of the resistant phenotype. We found significant differential expressions of genes involved in immunity and cell architecture. This data provides the first insights into how individuals escape the pathogenicity associated with infection.

Transcriptome Analysis to Study the Molecular Response in the Gill and Hepatopancreas Tissues of Macrobrachium nipponense to Salinity Acclimation

Macrobrachium nipponense is an economically important prawn species and common in Chinese inland capture fisheries. During aquaculture, M. nipponense can survive under freshwater and low salinity conditions. The molecular mechanism underlying the response to salinity acclimation remains unclear in this species; thus, in this study, we used the Illumina RNA sequencing platform for transcriptome analyses of the gill and hepatopancreas tissues of M. nipponense exposed to salinity stress [0.4‰ (S0, control group), 6‰ (S6, low salinity group), and 12‰ (S12, high salinity group)]. Differentially expressed genes were identified, and several important salinity adaptation-related terms and signaling pathways were found to be enriched, such as “ion transport,” “oxidative phosphorylation,” and “glycometabolism.” Quantitative real-time PCR demonstrated the participation of 12 key genes in osmotic pressure regulation in M. nipponense under acute salinity stress. Further, the role of carbonic anhydrase in response to salinity acclimation was investigated by subjecting the gill tissues of M. nipponense to in situ hybridization. Collectively, the results reported herein enhance our understanding of the mechanisms via which M. nipponense adapts to changes in salinity.

Antibacterial effects of extracted corn zein with garlic extract-based nanoemulsion on the shelf life of Vannamei prawn (Litopenaeus vannamei) at refrigerated temperature

BACKGROUND

Litopenaeus vannamei is one of the most perishable foods due to microorganism growth. Using essential oil-based nanoemulsion as a biodegradable and edible coating can enhance the shelf life of shrimp at refrigerated temperature through retarding microbial growth compared to synthetic coatings.

METHODS

Zein was extracted from 50 g dry milled corn with ethanol. Garlic essential oil was prepared by mixing with Tween 80. Nanoemulsion was prepared in an ultrasonic bath. DLS analysis, turbidity, and stability of nanoemulsions were performed. Radical scavenging activity, and total phenol content were done for evaluation of garlic essential oil and nanoemulsion. A 10% zein solution was prepared using ethanol 95% and glycerol plasticizer (2.5%). In the end, various microbial analysis, peroxide value, the thiobarbituric acid reactive substance (TBARS) value, the total volatile basic nitrogen (TVB-N) values, and sensory evaluations of different shrimp samples were determined.

RESULTS

Corn zein along with garlic-based nanoemulsion, had a great impact on the TVB-N, TBARS, and peroxide value (which indicated that z+24% garlic nanoemulsion group was the lowest among other groups on days 3, 7, and 14 (p < 0.05) as well as microbial properties (garlic EO nanoemulsion had significantly better antibacterial effectiveness compared to other groups (p < 0.05)), and sensory evaluation (the z+24% garlic nanoemulsion sample received a significantly higher score than other groups (p < 0.05)) of Litopenaeus vannamei.

CONCLUSIONS

Corn zein nanoemulsion functioned as an antioxidant and antimicrobial agent, increasing the shelf life of Litopenaeus vannamei at refrigerated temperature.

PRACTICAL APPLICATION

Emerging nanotechnology-based approaches with no side-effects on immune system of consumers plays a vital role in bioactive packaging, and on reduction of food spoilage or food poisoning in the transportation, exporting, and distribution stages of food products, especially marine-based products with cold chain transportation. Additionally, it can reduce aquaculture and environmental risks due to the usage of chemical agents used in packaging. Our results showed that administration of a corn zein with nanoemulsion of garlic extract can prolong shelf life of Litopenaeus vannamei. The paper should be of interest to readers such as food microbiologists, aquaculture scientists, fisheries scientists, marine biologists, biomaterial scientists, food packaging industrialists, medical microbiologists, public health managers, and health system managers.

Comparative Physiological and Transcriptomic Profiling Offers Insight into the Sexual Dimorphism of Hepatic Metabolism in Size-Dimorphic Spotted Scat (Scatophagus argus)

The spotted scat (Scatophagus argus) is an economically important cultured marine fish that exhibits a typical sexual size dimorphism (SSD). SSD has captivated considerable curiosity for farmed fish production; however, up till now the exact underlying mechanism remains largely unclear. As an important digestive and metabolic organ, the liver plays key roles in the regulation of fish growth. It is necessary to elucidate its significance as a downstream component of the hypothalamic-pituitary-liver axis in the formation of SSD. In this study, the liver physiological differences between the sexes were evaluated in S. argus, and the activity of several digestive and metabolic enzymes were affected by sex. Females had higher amylase, protease, and glucose-6-phosphate dehydrogenase activities, while males exhibited markedly higher hepatic lipase and antioxidant enzymes activities. A comparative transcriptomics was then performed to characterize the responsive genes. Illumina sequencing generated 272.6 million clean reads, which were assembled into 79,115 unigenes. A total of 259 differentially expressed genes were identified and a few growth-controlling genes such as igf1 and igfbp1 exhibited female-biased expression. Further analyses showed that several GO terms and pathways associated with metabolic process, particularly lipid and energy metabolisms, were significantly enriched. The male liver showed a more active mitochondrial energy metabolism, implicating an increased energy expenditure associated with reproduction. Collectively, the female-biased growth dimorphism of S. argus may be partially attributed to sexually dimorphic metabolism in the liver. These findings would facilitate further understanding of the nature of SSD in teleost fish.

Comparative genomics of the coconut crab and other decapod crustaceans: exploring the molecular basis of terrestrial adaptation

Background

The complex life cycle of the coconut crab, Birgus latro, begins when an obligate terrestrial adult female visits the intertidal to hatch zoea larvae into the surf. After drifting for several weeks in the ocean, the post-larval glaucothoes settle in the shallow subtidal zone, undergo metamorphosis, and the early juveniles then subsequently make their way to land where they undergo further physiological changes that prevent them from ever entering the sea again. Here, we sequenced, assembled and analyzed the coconut crab genome to shed light on its adaptation to terrestrial life. For comparison, we also assembled the genomes of the long-tailed marine-living ornate spiny lobster, Panulirus ornatus, and the short-tailed marine-living red king crab, Paralithodes camtschaticus. Our selection of the latter two organisms furthermore allowed us to explore parallel evolution of the crab-like form in anomurans.

Results

All three assembled genomes are large, repeat-rich and AT-rich. Functional analysis reveals that the coconut crab has undergone proliferation of genes involved in the visual, respiratory, olfactory and cytoskeletal systems. Given that the coconut crab has atypical mitochondrial DNA compared to other anomurans, we argue that an abundance of kif22 and other significantly proliferated genes annotated with mitochondrial and microtubule functions, point to unique mechanisms involved in providing cellular energy via nuclear protein-coding genes supplementing mitochondrial and microtubule function. We furthermore detected in the coconut crab a significantly proliferated HOX gene, caudal, that has been associated with posterior development in Drosophila, but we could not definitively associate this gene with carcinization in the Anomura since it is also significantly proliferated in the ornate spiny lobster. However, a cuticle-associated coatomer gene, gammacop, that is significantly proliferated in the coconut crab, may play a role in hardening of the adult coconut crab abdomen in order to mitigate desiccation in terrestrial environments.

Conclusion

The abundance of genomic features in the three assembled genomes serve as a source of hypotheses for future studies of anomuran environmental adaptations such as shell-utilization, perception of visual and olfactory cues in terrestrial environments, and cuticle sclerotization. We hypothesize that the coconut crab exhibits gene proliferation in lieu of alternative splicing as a terrestrial adaptation mechanism and propose life-stage transcriptomic assays to test this hypothesis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07636-9.

Multi-omic approach provides insights into osmoregulation and osmoconformation of the crab Scylla paramamosain

Osmoregulation and osmoconformation are two mechanisms through which aquatic animals adapt to salinity fluctuations. The euryhaline crab Scylla paramamosain, being both an osmoconformer and osmoregulator, is an excellent model organism to investigate salinity adaptation mechanisms in brachyurans. In the present study, we used transcriptomic and proteomic approaches to investigate the response of S. paramamosain to salinity stress. Crabs were transferred from a salinity of 25 ppt to salinities of 5 ppt or 33 ppt for 6 h and 10 days. Data from both approaches revealed that exposure to 5 ppt resulted in upregulation of ion transport and energy metabolism associated genes. Notably, acclimation to low salinity was associated with early changes in gene expression for signal transduction and stress response. In contrast, exposure to 33 ppt resulted in upregulation of genes related to amino acid metabolism, and amino acid transport genes were upregulated only at the early stage of acclimation to this salinity. Our study reveals contrasting mechanisms underlying osmoregulation and osmoconformation within the salinity range of 5–33 ppt in the mud crab, and provides novel candidate genes for osmotic signal transduction, thereby providing insights on understanding the salinity adaptation mechanisms of brachyuran crabs.

CrustyBase: an interactive online database for crustacean transcriptomes

Transcriptome sequencing has opened the field of genomics to a wide variety of researchers, owing to its efficiency, applicability across species and ability to quantify gene expression. The resulting datasets are a rich source of information that can be mined for many years into the future, with each dataset providing a unique angle on a specific context in biology. Maintaining accessibility to this accumulation of data presents quite a challenge for researchers.The primary focus of conventional genomics databases is the storage, navigation and interpretation of sequence data, which is typically classified down to the level of a species or individual. The addition of expression data adds a new dimension to this paradigm - the sampling context. Does gene expression describe different tissues, a temporal distribution or an experimental treatment? These data not only describe an individual, but the biological context surrounding that individual. The structure and utility of a transcriptome database must therefore reflect these attributes. We present an online database which has been designed to maximise the accessibility of crustacean transcriptome data by providing intuitive navigation within and between datasets and instant visualization of gene expression and protein structure.The site is accessible at https://crustybase.org and currently holds 10 datasets from a range of crustacean species. It also allows for upload of novel transcriptome datasets through a simple web interface, allowing the research community to contribute their own data to a pool of shared knowledge.

Probiotic Role of Salt Pan Bacteria in Enhancing the Growth of Whiteleg Shrimp, Litopenaeus vannamei

Development of probiotics to improve the growth of cultured species is a key to sustainable aquaculture. The present study investigates the potential of salt pan bacteria as probiotics for Litopenaeus vannamei. Halotolerant bacteria (100) were screened for enzyme production and mucus adhesion in vitro. The bacteria (SK07, SK27, ABSK55, FSK444, TSK17, TSK71) exhibiting promising enzyme activity and adhesive property in vitro were selected to study their effect on the growth and metabolism of L. vannamei in vivo. When administered to shrimps individually as a water additive in experiment I, SK07, SK27 and TSK71 significantly (p < 0.05) increased shrimp weight as compared to the control. In experiment II, a lyophilized bacterial consortium (test) prepared with the four best isolates (SK07, SK27, ABSK55, TSK71), exhibited significantly higher weight gain of shrimps, better feed efficiency and final yield as compared to control. Total enzyme activity (amylase, protease, lipase) in the shrimp gut was significantly higher in the test than the control. The four isolates showed 99% nBLAST similarity with Bacillus subtilis, Bacillus amyloliquefaciens, Bacillus licheniformis and Pseudomonas sp. Presence of these bacteria in the shrimp gut was confirmed by using specific PCR-based molecular probes and 16S rDNA sequencing. Safety evaluation by antibiotic susceptibility test and hemolytic activity test indicated that the bacteria are safe as bioinoculants. The increased enzyme activity by colonisation of the isolates in the shrimp gut, along with improved growth and feed utilisation efficiency, strongly confirms that these salt pan bacteria are prospective probiotics in shrimp aquaculture.

A new insight into the intestine of Pacific white shrimp: Regulation of intestinal homeostasis and regeneration in Litopenaeus vannamei during temperature fluctuation

Litopenaeus vannamei (L. vannamei) is an essential aquaculture shrimp throughout the world, but its aquaculture industry is threatened by temperature fluctuation. In this study, our histological results indicated that the shrimp intestine has a self-repairing ability during temperature fluctuation; however the potential mechanisms were still unknown. Therefore, transcriptome profiles of the intestine were collected from shrimp at 28 °C (C28), 13 °C (T13) and 28 °C after their temperature rose back (R28) and were analyzed. A total of 2229 differentially expressed genes (DEGs) (986 up- and 1243 downregulated) were identified in the C28 group, and 1790 DEGs (933 up- and 857 downregulated) were identified in the R28 group when compared to their expression levels in the T13 group. According to the functional annotation using KEGG, we found that the immune system was the most enriched section of organismal systems and that the shrimp can mobilize the body's immune response to regulate organism homeostasis during temperature fluctuation, although cold stress decreased the immunity. Additionally, metabolic inhibition is a strategy to cope with cold stress, and the regulation of lipid metabolism was especially important for shrimp during temperature fluctuation. Remarkably, the Hippo signaling pathway might help the repair of intestinal structure. Our research provides the first histological analysis and transcriptome profiling for the L. vannamei intestine during the temperature fluctuation stage. These results enrich our understanding of the mechanism of intestinal self-repair and homeostasis and could provide guidance for shrimp farming.

Transcriptome Analysis Reveals the Tolerance Mechanism of Mantis Shrimp (Oratosquilla oratoria) under a Lipopolysaccharide Challenge

Lipopolysaccharide (LPS), a major cell wall component of Gram-negative bacteria, is considered to lead to some disease development in commercial crustaceans. However, mantis shrimps Oratosquilla oratoria (Crustacea: Stomatopoda) have a strong vitality and ability to resist disease. To study the tolerance mechanism of mantis shrimp, transcriptome analyses were conducted in hepatopancreas of O. oratoria under LPS challenge investigation. Totally, 84 547 044 clean reads were obtained from transcriptomes (43 159 230 in OP (control), 41 387 814 in OL (treatment), respectively). Unigenes, the longest transcript of each gene, with a total length of 68 318 880 bp and the total number of 100 978 were obtained. 8369 (8.28%) of unigenes were successfully annotated in all databases and 54 888 (54.35%) were annotated in at least one database. Finally, 1012 differentially expressed genes (DEGs) including 439 and 573 showed significantly upregulated and downregulated were determined between OL and OP, respectively. Moreover, those DEGs only expressed in OL or OP accounted for 8.99%. The functional classification based on GO and KEGG indicated that the common enrichment categories for the DEGs are "amino sugar metabolic" and "cellular homeostasis" and that the progress of nutrient metabolic and homeostasis in cells is important in facing variable environmental conditions. Protein-protein interaction analysis elucidated proteins, β-actin (ACTB_G1), T-complex protein subunits (TCPs), heat shock proteins (HSPs), hydroxysteroid dehydrogenase-like protein 2 (HSDL2), kinesin family member 5 (KIF5), methylglutaconyl-CoA hydratase (AUH), and myosin heavy chain (MYH) may play key roles in response to an LPS challenge. This study laid a foundation to further investigate the possible adaptation way that O. oratoria survives in a bacterial challenge.

Effects of a Sudden Drop in Salinity on Immune Response Mechanisms of Anadara kagoshimensis

In this experiment, the effects of a sudden drop of salinity on the immune response mechanisms of the ark shell Anadara kagoshimensis were examined by simulating the sudden drop of salinity that occurs in seawater after a rainstorm. Additionally, the differentially expressed genes (DEGs) were identified using transcriptome sequencing. When the salinity dropped from 30‱ (S30) to 14‱ (S14), the phagocytic activity of blood lymphocytes, the O₂^- levels produced from respiratory burst, the content of reactive oxygen species, and the activities of lysozymes and acid phosphatases increased significantly, whereas the total count of blood lymphocytes did not increase. Total count of blood lymphocytes in 22‱ salinity (S22) was significantly higher than that in any other group. The raw data obtained from sequencing were processed with Trimmomatic (Version 0.36). The expression levels of unigenes were calculated using transcripts per million (TPM) based on the effects of sequencing depth, gene length, and sample on reads. Differential expression analysis was performed using DESeq (Version 1.12.4). Transcriptome sequencing revealed 269 (101 up-regulated, 168 down-regulated), 326 (246 up-regulated, 80 down-regulated), and 185 (132 up-regulated, 53 down-regulated) significant DEGs from comparison of the S14 vs. S22, S22 vs. S30, and S14 vs. S30 groups, respectively. Gene Ontology enrichment analysis of the DEGs in these salinity comparison groups revealed that the cellular amino acid metabolic process, the regulation of protein processing, the regulation of response to stress, and other terms were significantly enriched. Kyoto Encyclopedia of Genes and Genomes enrichment analysis showed that nucleotide-binding, oligomerization domain (NOD)-like receptor signaling pathway (ko04621), apoptosis-multiple species (ko04215), Toll and Imd signaling pathway (ko04624), NF-κB signaling pathway (ko04064), apoptosis (ko04210), and focal adhesion (ko04510) were significantly enriched in all salinity comparison groups. qRT-PCR verification of 12 DEGs in the above six pathways was conducted, and the results were consistent with the transcriptome sequencing results in terms of up-regulation and down-regulation, which illustrates that the transcriptome sequencing data are credible. These results were used to preliminarily explore the effects of a sudden drop of salinity on blood physiological and biochemical indexes and immunoregulatory mechanisms of A. kagoshimensis. They also provide a theoretical basis for the selection of bottom areas optimal for release and proliferation of A. kagoshimensis required to restore the declining populations of this species.

Comparative transcriptome analysis reveals the different roles between hepatopancreas and intestine of Litopenaeus vannamei in immune response to aflatoxin B1 (AFB1) challenge

Aflatoxin B1 (AFB1) is a mycotoxin mainly produced by Aspergillus flavus and Aspergillus parasiticus contaminating food, feed ingredients and products of animal origin. In mammals, this toxin causes widespread organ-specific damage; it is immunotoxicity and could promote hepatotoxicity, alter intestinal functions and so on. In this study, we conducted transcriptome and histomorphology analyses of hepatopancreas and intestinal in Litopenaeus vannamei (L. vannamei) challenged with AFB1. Totally 12,014 and 1387 differentially expression genes (DEGs) were identified in the hepatopancreas and intestine, respectively. In hepatopancreas, a total of 1995 DEGs were mainly annotated and grouped into 18 processes or pathways related to animal immune system. With respect to intestine, a total of 152 DEGs were mainly annotated to 7 processes or pathways related to animal immune system. Meanwhile, we determined the relative mRNA expression of several crucial representative immune genes including Toll, immune deficiency (IMD), prophenoloxidase (proPO), Rab and glutathione S-transferase (GST) in the hepatopancreas and intestines of shrimp at 3-, 6-, 12-, 18-, 24- and 30-d after challenged by AFB1. Exposure to AFB1 increased mortality, decrease weight gain rate, severely destroyed the histomorphology of hepatopancreas and intestine, and resulted in the damaged of immune system of shrimp. The present data reveals the different roles between hepatopancreas and intestine of L. vannamei in immune response to AFB1 challenge, and provides insight into the molecular basis of the relationship between hepatopancreas and intestinal immunity during either homeostasis or inflammation.

The Effect of Disease and Season to Hepatopancreas and Intestinal Mycobiota of Litopenaeus vannamei

Increasing evidence has manifested that the gut bacterial microbiota of shrimps is closely related to the environmental factors, host developmental stage and health status like that of humans and animals does. These studies have provided an important guidance for improving shrimp culture benefits. In practice, aside from bacteria, eukaryotic microorganisms dominated by fungal microbiota (mycobiota), also play a key role in host growth, metabolism and homeostasis. However, little so far is known about the mycobiota in the digestive tract of shrimp. In this study, we used high-throughput sequencing of internal transcribed spacer 1 region to characterize the hepatopancreas and intestinal mycobiota of Pacific white shrimp and their connections with disease incidence and seasonal variation. The results showed that the hepatopancreas and intestinal mycobiota of Litopenaeus vannamei are dominated by the phyla Ascomycota and Basidiomycota, and the genera Alternaria, Tuber, Hortaea, Sarocladium, and Stagonospora. The fungal microbiota significantly varies under the influence of disease and seasonal variation. Sick shrimps had a higher level of potential pathogenic fungus, Candida in the intestine. Healthy shrimps had a higher abundance of the genera Didymella and Filobasidium in the gut, and Pyrenochaetopsis in the hepatopancreas. Of note, most of the fungi carried by Pacific white shrimps were pathogens to humans. This study has revealed the intestinal and hepatopancreas mycobiota of L. vannamei and the effects of diseases and seasonal variation to the mycobiota. Our study provides important guidance for Pacific white shrimp farming and sheds further insight on the fungal microbiota.

Time-series transcriptomic analysis of the kelp grouper Epinephelus moara in response to low salinity stress

The Kelp grouper Epinephelus moara is one of the most widely consumed and economically important marine fish in China. The species can tolerate a wide range of salinity, but genomic resources are not available, and the molecular mechanisms underlying adaptation to salinity at the transcriptomic level remain largely unclear. In this study, the transcriptomic responses of the liver of E. moara under low salinity were investigated using the Illumina digital gene expression system. After de novo assembly, 499,356 transcripts were generated and contributed 445,068 unigenes. A total of 14, 19, 33 and 3101 genes were differentially expressed following exposure to low salinity stress for 2, 6, 24 and 48 h, respectively. Only two genes were differentially expressed in all groups. Four genes related to metabolism and ambient salinity adaption were randomly selected to validate the differentially expressed genes (DEGs) by real-time PCR. Gene Ontology (GO) and Kyoto Encyclopaedia of Genes and Genomes (KEGG) pathway enrichment analysis were used to analyse the functional significance of DEGs, including those responding to salinity through diverse biological processes, cellular components, molecular functions, and pathways associated with metabolic and osmotic responses. This work provides new insight into the response to salinity challenges in E. moara, and the findings expand our knowledge of the molecular basis of metabolic regulation mechanisms in this species. Additionally, the transcriptional data provide a valuable resource for future molecular and genetic studies on E. moara.

De novo assembly, characterization, functional annotation and expression patterns of the black tiger shrimp (Penaeus monodon) transcriptome

The black tiger shrimp (Penaeus monodon) remains the second most widely cultured shrimp species globally; however, issues with disease and domestication have seen production levels stagnate over the past two decades. To help identify innovative solutions needed to resolve bottlenecks hampering the culture of this species, it is important to generate genetic and genomic resources. Towards this aim, we have produced the most complete publicly available P. monodon transcriptome database to date based on nine adult tissues and eight early life-history stages (BUSCO - Complete: 98.2% [Duplicated: 51.3%], Fragmented: 0.8%, Missing: 1.0%). The assembly resulted in 236,388 contigs, which were then further segregated into 99,203 adult tissue specific and 58,678 early life-history stage specific clusters. While annotation rates were low (approximately 30%), as is typical for a non-model organisms, annotated transcript clusters were successfully mapped to several hundred functional KEGG pathways. Transcripts were clustered into groups within tissues and early life-history stages, providing initial evidence for their roles in specific tissue functions, or developmental transitions. We expect the transcriptome to provide an essential resource to investigate the molecular basis of commercially relevant-significant traits in P. monodon and other shrimp species.

The State of "Omics" Research for Farmed Penaeids: Advances in Research and Impediments to Industry Utilization

Elucidating the underlying genetic drivers of production traits in agricultural and aquaculture species is critical to efforts to maximize farming efficiency. "Omics" based methods (i.e., transcriptomics, genomics, proteomics, and metabolomics) are increasingly being applied to gain unprecedented insight into the biology of many aquaculture species. While the culture of penaeid shrimp has increased markedly, the industry continues to be impeded in many regards by disease, reproductive dysfunction, and a poor understanding of production traits. Extensive effort has been, and continues to be, applied to develop critical genomic resources for many commercially important penaeids. However, the industry application of these genomic resources, and the translation of the knowledge derived from "omics" studies has not yet been completely realized. Integration between the multiple "omics" resources now available (i.e., genome assemblies, transcriptomes, linkage maps, optical maps, and proteomes) will prove critical to unlocking the full utility of these otherwise independently developed and isolated resources. Furthermore, emerging "omics" based techniques are now available to address longstanding issues with completing keystone genome assemblies (e.g., through long-read sequencing), and can provide cost-effective industrial scale genotyping tools (e.g., through low density SNP chips and genotype-by-sequencing) to undertake advanced selective breeding programs (i.e., genomic selection) and powerful genome-wide association studies. In particular, this review highlights the status, utility and suggested path forward for continued development, and improved use of "omics" resources in penaeid aquaculture.

Transcriptome analysis of Pacific white shrimp (Litopenaeus vannamei) challenged by Vibrio parahaemolyticus reveals unique immune-related genes

Pacific white shrimp (Litopenaeus vannamei) is an important cultural species worldwide. However, Vibrio spp. infections have caused a great economic loss in Pacific white shrimp culture industry. The immune responses of Pacific white shrimp to the Vibrio spp. is not fully characterized. In this study, the transcriptomic profiles of L. vannamei hemocytes were explored by injecting with or without Vibrio parahaemolyticus. Totally, 42,632 high-quality unigenes were obtained from RNAseq data. Comparative genome analysis showed 2258 differentially expressed genes (DEGs) following the Vibrio challenge, including 1017 up-regulated and 1241 down-regulated genes. Eight DEGs were randomly selected for further validation by quantitative real-time RT-PCR (qRT-PCR) and the results showed that are consistent with the RNA-seq data. Due to the lack of predictable adaptive immunity, shrimps rely on an innate immune system to defend themselves against invading microbes by recognizing and clearing them through humoral and cellular immune responses. Here we focused our studies on the humoral immunity, five genes (SR, MNK, CTL3, GILT, and ALFP) were selected from the transcriptomic data, which were significantly up-regulated by V. parahaemolyticus infection. These genes were widely expressed in six different tissues and were up-regulated by both Gram negative bacteria (V. parahaemolyticus) and Gram positive bacteria (Staphylococcus aureus). To further extend our studies, we knock-down those five genes by dsRNA in L. vannamei and analyzed the functions of specific genes against V. parahaemolyticus and S. aureus by bacterial clearance analysis. We found that the ability of L. vannamei was significantly reduced in bacterial clearance when treated with those specific dsRNA. These results indicate that those five genes play essential roles in antibacterial immunity and have its specific functions against different types of pathogens. The obtained data will shed a new light on the immunity of L. vannamei and pave a new way for fighting against the bacterial infection in Pacific white shrimp.

Composition, diversity and function of intestinal microbiota in pacific white shrimp (Litopenaeus vannamei) at different culture stages

Intestinal microbiota is an integral component of the host and plays important roles in host health. The pacific white shrimp is one of the most profitable aquaculture species commercialized in the world market with the largest production in shrimp consumption. Many studies revealed that the intestinal microbiota shifted significantly during host development in other aquaculture animals. In the present study, 22 shrimp samples were collected every 15 days from larval stage (15 day post-hatching, dph) to adult stage (75 dph) to investigate the intestinal microbiota at different culture stages by targeting the V4 region of 16S rRNA gene, and the microbial function prediction was conducted by PICRUSt. The operational taxonomic unit (OTU) was assigned at 97% sequence identity. A total of 2,496 OTUs were obtained, ranging from 585 to 1,239 in each sample. Forty-three phyla were identified due to the classifiable sequence. The most abundant phyla were Proteobacteria, Cyanobacteria, Tenericutes, Fusobacteria, Firmicutes, Verrucomicrobia, Bacteroidetes, Planctomycetes, Actinobacteria and Chloroflexi. OTUs belonged to 289 genera and the most abundant genera were Candidatus_Xiphinematobacter, Propionigenium, Synechococcus, Shewanella and Cetobacterium. Fifty-nine OTUs were detected in all samples, which were considered as the major microbes in intestine of shrimp. The intestinal microbiota was enriched with functional potentials that were related to transporters, ABC transporters, DNA repair and recombination proteins, two component system, secretion system, bacterial motility proteins, purine metabolism and ribosome. All the results showed that the intestinal microbial composition, diversity and functions varied significantly at different culture stages, which indicated that shrimp intestinal microbiota depended on culture stages. These findings provided new evidence on intestinal microorganism microecology and greatly enhanced our understanding of stage-specific community in the shrimp intestinal ecosystem.

Transcriptomic analysis of Crassostrea sikamea × Crassostrea angulata hybrids in response to low salinity stress

Hybrid oysters often show heterosis in growth rate, weight, survival and adaptability to extremes of salinity. Oysters have also been used as model organisms to study the evolution of host-defense system. To gain comprehensive knowledge about various physiological processes in hybrid oysters under low salinity stress, we performed transcriptomic analysis of gill tissue of Crassostrea sikamea ♀ × Crassostrea angulata♂ hybrid using the deep-sequencing platform Illumina HiSeq. We exploited the high-throughput technique to delineate differentially expressed genes (DEGs) in oysters maintained in hypotonic conditions. A total of 199,391 high quality unigenes, with average length of 644 bp, were generated. Of these 35 and 31 genes showed up- and down-regulation, respectively. Functional categorization and pathway analysis of these DEGs revealed enrichment for immune mechanism, apoptosis, energy metabolism and osmoregulation under low salinity stress. The expression patterns of 41 DEGs in hybrids and their parental species were further analyzed by quantitative real-time PCR (qRT-PCR). This study will serve as a platform for subsequent gene expression analysis regarding environmental stress. Our findings will also provide valuable information about gene expression to better understand the immune mechanism, apoptosis, energy metabolism and osmoregulation in hybrid oysters under low salinity stress.

Comparative transcriptomics enlarges the toolkit of known developmental genes in mollusks

BACKGROUND

Mollusks display a striking morphological disparity, including, among others, worm-like animals (the aplacophorans), snails and slugs, bivalves, and cephalopods. This phenotypic diversity renders them ideal for studies into animal evolution. Despite being one of the most species-rich phyla, molecular and in silico studies concerning specific key developmental gene families are still scarce, thus hampering deeper insights into the molecular machinery that governs the development and evolution of the various molluscan class-level taxa.

RESULTS

Next-generation sequencing was used to retrieve transcriptomes of representatives of seven out of the eight recent class-level taxa of mollusks. Similarity searches, phylogenetic inferences, and a detailed manual curation were used to identify and confirm the orthology of numerous molluscan Hox and ParaHox genes, which resulted in a comprehensive catalog that highlights the evolution of these genes in Mollusca and other metazoans. The identification of a specific molluscan motif in the Hox paralog group 5 and a lophotrochozoan ParaHox motif in the Gsx gene is described. Functional analyses using KEGG and GO tools enabled a detailed description of key developmental genes expressed in important pathways such as Hedgehog, Wnt, and Notch during development of the respective species. The KEGG analysis revealed Wnt8, Wnt11, and Wnt16 as Wnt genes hitherto not reported for mollusks, thereby enlarging the known Wnt complement of the phylum. In addition, novel Hedgehog (Hh)-related genes were identified in the gastropod Lottia cf. kogamogai, demonstrating a more complex gene content in this species than in other mollusks.

CONCLUSIONS

The use of de novo transcriptome assembly and well-designed in silico protocols proved to be a robust approach for surveying and mining large sequence data in a wide range of non-model mollusks. The data presented herein constitute only a small fraction of the information retrieved from the analysed molluscan transcriptomes, which can be promptly employed in the identification of novel genes and gene families, phylogenetic inferences, and other studies using molecular tools. As such, our study provides an important framework for understanding some of the underlying molecular mechanisms involved in molluscan body plan diversification and hints towards functions of key developmental genes in molluscan morphogenesis.