Search

The light-dark cycle significantly impacts the growth and development of animals. Mantis shrimps (Oratosquilla oratoria) receive light through their complex photoreceptors. To reveal the adaptive expression mechanism of the mantis shrimp induced in a dark environment, we performed comparative transcriptome analysis with O. oratoria cultured in a light environment (Oo-L) as the control group and O. oratoria cultured in a dark environment (Oo-D) as the experimental group. In the screening of differentially expressed genes (DEGs) between the Oo-L and Oo-D groups, a total of 88 DEGs with |log2FC| > 1 and FDR < 0.05 were identified, of which 78 were upregulated and 10 were downregulated. Then, FBP1 and Pepck were downregulated in the gluconeogenesis pathway, and MKNK2 was upregulated in the MAPK classical pathway, which promoted cell proliferation and differentiation, indicating that the activity of mantis shrimp was slowed and the metabolic rate decreases in the dark environment. As a result, the energy was saved for its growth and development. At the same time, we performed gene set enrichment analysis (GSEA) on all DEGs. In the KEGG pathway analysis, each metabolic pathway in the dark environment showed a slowing trend. GO was enriched in biological processes such as eye development, sensory perception and sensory organ development. The study showed that mantis shrimp slowed down metabolism in the dark, while the role of sensory organs prominent. It provides important information for further understanding the energy metabolism and has great significance to study the physiology of mantis shrimp in dark environment., (© 2023. The Author(s), under exclusive licence to Springer Nature Switzerland AG.)

Transcriptomes from nontraditional model organisms often harbor a wealth of unexplored data. Examining these data sets can lead to clarity and novel insights in traditional systems, as well as to discoveries across a multitude of fields. Despite significant advances in DNA sequencing technologies and in their adoption, access to genomic and transcriptomic resources for nontraditional model organisms remains limited. Crustaceans, for example, being among the most numerous, diverse, and widely distributed taxa on the planet, often serve as excellent systems to address ecological, evolutionary, and organismal questions. While they are ubiquitously present across environments, and of economic and food security importance, they remain severely underrepresented in publicly available sequence databases. Here, we present CrusTome, a multispecies, multitissue, transcriptome database of 201 assembled mRNA transcriptomes (189 crustaceans, 30 of which were previously unpublished, and 12 ecdysozoans for phylogenetic context) as an evolving and publicly available resource. This database is suitable for evolutionary, ecological, and functional studies that employ genomic/transcriptomic techniques and data sets. CrusTome is presented in BLAST and DIAMOND formats, providing robust data sets for sequence similarity searches, orthology assignments, phylogenetic inference, etc. and thus allowing for straightforward incorporation into existing custom pipelines for high-throughput analyses. In addition, to illustrate the use and potential of CrusTome, we conducted phylogenetic analyses elucidating the identity and evolution of the cryptochrome/photolyase family of proteins across crustaceans., Competing Interests: Conflicts of interest The author(s) declare no conflict of interest., (© The Author(s) 2023. Published by Oxford University Press on behalf of The Genetics Society of America.)

Copepods are zooplanktonic crustaceans ubiquitously widespread in aquatic systems. Although they are not the target, copepods are exposed to a wide variety of pollutants such as insect growth regulators (IGRs). The aim of this study was to investigate the molecular response of a non-targeted organism, the copepod Eurytemora affinis , to an IGR. Adult males and females were exposed to two sub-lethal concentrations of tebufenozide (TEB). Our results indicate a sex-specific response with a higher sensitivity in males, potentially due to a differential activation of stress response pathways. In both sexes, exposure to TEB triggered similar pathways to those found in targeted species by modulating the transcription of early and late ecdysone responsive genes. Among them were genes involved in cuticle metabolism, muscle contraction, neurotransmission, and gametogenesis, whose mis-regulation could lead to moult, locomotor, and reproductive impairments. Furthermore, genes involved in epigenetic processes were found in both sexes, which highlights the potential impact of exposure to TEB on future generations. This work allows identification of (i) potential biomarkers of ecdysone agonists and (ii) further assessment of putative physiological responses to characterize the effects of TEB at higher biological levels. The present study reinforces the suitability of using E. affinis as an ecotoxicological model.

Little is known about the mechanisms underlying the relationship between genetic variation and the adaptation of Oratosquilla oratoria populations to different habitat temperature. Here, the genome-wide genetic information of three O. oratoria populations were obtained by IIB restriction site-associated DNA (2b-RAD) sequencing and 2403 single-nucleotide polymorphisms (SNPs) were identified. Based on the 2403 SNPs, we found a remarkable genetic differentiation between the Yellow Sea and the East China Sea groups of O. oratoria. Furthermore, 63 SNPs are thought to be associated with different sea temperatures. Based on the 63 SNPs, it is hypothesised that the long-term temperature differences may contribute to the variation of genes associated with multiple biological functions, such as material metabolism, cytoskeleton, cellular processes, inflammatory response and hormonal regulation. This study provides new information for elucidating the molecular mechanisms underlying the relationship between genetic variation and the adaptation of Oratosquilla oratoria populations to different temperature., (Copyright © 2021 Elsevier Inc. All rights reserved.)

One fundamental character common to pancrustaceans (Crustacea and Hexapoda) is a mineralized rigid exoskeleton whose principal organic components are chitin and proteins. In contrast to traditional research in the field that has been devoted to the structural and physicochemical aspects of biomineralization, the present study explores transcriptomic aspects of biomineralization as a first step towards adding a complementary molecular layer to this field. The rigidity of the exoskeleton in pancrustaceans dictates essential molt cycles enabling morphological changes and growth. Thus, formation and mineralization of the exoskeleton are concomitant to the timeline of the molt cycle. Skeletal proteinaceous toolkit elements have been discovered in previous studies using innovative molt-related binary gene expression patterns derived from transcriptomic libraries representing the major stages comprising the molt cycle of the decapod crustacean Cherax quadricarinatus. Here, we revisited some prominent exoskeleton-related structural proteins encoding and, using the above molt-related binary pattern methodology, enlarged the transcriptomic database of C. quadricarinatus. The latter was done by establishing a new transcriptomic library of the cuticle forming epithelium and molar tooth at four different molt stages (i.e., inter-molt, early pre-molt, late pre-molt and post-molt) and incorporating it to a previous transcriptome derived from the gastroliths and mandible. The wider multigenic approach facilitated by the newly expanded transcriptomic database not only revisited single genes of the molecular toolkit, but also provided both scattered and specific information that broaden the overview of proteins and gene clusters which are involved in the construction and biomineralization of the exoskeleton in decapod crustaceans., (Copyright © 2020 The Authors. Published by Elsevier Inc. All rights reserved.)

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.

Chinese mitten crab (Eriocheir sinensis) is an important aquaculture species in Crustacea. Functional analysis, although essential, has been hindered due to the lack of sufficient genomic or transcriptomic resources. In this study, transcriptome sequencing was conducted on 59 samples representing diverse developmental stages (fertilized eggs, zoea, megalopa, three sub-stages of larvae, juvenile crabs, and adult crabs) and different tissues (eyestalk, hepatopancreas, and muscle from juvenile crabs, and eyestalk, hepatopancreas, muscle, heart, stomach, gill, thoracic ganglia, intestine, ovary, and testis from adult crabs) of E. sinensis. A comprehensive reference transcriptome was assembled, including 19,023 protein-coding genes. Hierarchical clustering based on 128 differentially expressed cuticle-related genes revealed two distinct expression patterns during the early larval developmental stages, demonstrating the distinct roles of these genes in "crab-like" cuticle formation during metamorphosis and cuticle calcification after molting. Phylogenetic analysis of 1406 one-to-one orthologous gene families identified from seven arthropod species and Caenorhabditis elegans strongly supported the hypothesis that Malacostraca and Branchiopoda do not form a monophyletic group. Furthermore, Branchiopoda is more phylogenetically closely related to Hexapoda, and the clade of Hexapoda and Branchiopoda and the clade of Malacostraca belong to the Pancrustacea. This study offers a high-quality transcriptome resource for E. sinensis and demonstrates the evolutionary relationships of major arthropod groups. The differentially expressed genes identified in this study facilitate further investigation of the cuticle-related gene expression networks which are likely associated with "crab-like" cuticle formation during metamorphosis and cuticle calcification after molting., (Copyright © 2020 The Authors. Published by Elsevier B.V. All rights reserved.)

Salinity fluctuation may detrimentally affect the composition and biological processes of crustaceans. As a euryhaline crustacean, Oratosquilla oratoria can survive at salinities ranging from 20 psu to 40 psu. Therefore, we designed five salinity gradients (20, 25, 30, 35, and 40 psu) and 66.39 Gb clean transcriptome data were obtained after O. oratorias were exposed to each gradient for 24 h. All clean data were spliced into 50,482 unigenes and 17,035 unigenes were annotated in at least one database. Compared with 30 psu, 1010, 851, 1733 and 2188 differentially expressed genes were obtained at 20, 25, 35 and 40 psu, respectively. Results also showed that the osmoregulation of O. oratoria is primarily regulated by lipid and amino acid metabolism, amongst others. No significant up-regulated pathways were enriched at 25 psu and 35 psu, although more significant down-regulated pathways were obtained at 35 psu. Therefore, we assumed that the optimum survival salinity of O. oratoria may range from 25 psu to 35 psu. However, 35 psu may be more suitable for O. oratoria. In addition, 55 unigenes that encode putative inorganic ion exchanges were identified. This study aims to provide fundamental information for understanding the osmoregulation mechanisms of crustaceans., (Copyright © 2019 Elsevier B.V. All rights reserved.)

As the representative species of Stomatopoda, Oratosquilla oratoria can tolerate the complex benthic environment. In this study, our goal was to develop the transcriptomic resource of O. oratoria that would support adaptation mechanism studies. We generated the whole transcriptome of O. oratoria from combined tissues (eyestalk, muscle, sexual and viscus) using Hiseq technology. A total of 51,305,284 high-quality clean reads were assembled to produce 59,054 non-redundant transcripts with a mean length of 987nt using Trinity and CORSET software. Among the predictable unigenes, a total of 32,451 unigenes were annotated based on protein databases. Finally, we predicted the coding sequences of 31,822 unigenes and obtained 19,057 SSRs in the present study. The present study will provide an important resource and foundational understanding for future genomic research of O. oratoria., (Copyright © 2017 Elsevier B.V. All rights reserved.)

Cr (VI), the pathogenicity factor, is widely known to cause toxic effects in living organisms. Given the economic importance of the mantis shrimp ( Oratosquilla oratoria), the understanding of impacts by Cr (VI) is considered important. In this study, transcriptome of mantis shrimp was characterized by a comparison between control and Cr (VI)-treated samples using RNA-seq approach. Totally, 88 234 826 bp and 13.24G clean reads were obtained. The total length and number of unigenes were 68 411 206 bp and 100 918, respectively. The maximal and average length of unigenes was 24 906 bp and 678 bp, respectively (N50, 798 bp). 7115 of these unigenes accounted for 7.05% of the total that were annotated in all databases. After annotation of assembled unigenes, 35 619 of them were assigned into 3 functional categories and 56 subcategories using Gene Ontology; 18 580 of them were assigned into 26 functional categories using Clusters of Orthologous Groups of proteins; 16 864 of them were assigned into 5 major categories and 32 subclasses using KEGG. Finally, 1730 genes were differentially expressed (DGEs), 9 up-regulated pathways (protein digestion and absorption, neuroactive ligand-receptor interaction, pancreatic secretion, tyrosine metabolism, amoebiasis, ECM-receptor interaction, riboflavin metabolism, amino sugar and nucleotide sugar metabolism and AGE-RAGE signaling pathway in diabetic complications) were significantly enriched ( q < 0.05), and one down-regulated pathway ( Staphylococcus aureus infection) was significantly enriched ( q < 0.05). Up-regulation of genes in pathways of protein digestion/absorption ( PepT1/SLC15A and ATP1B) and environment information processing ( COL1AS, COL4A; LAMA3&#95;5, LAMB3; FN1 and TN) may imply the potentially positive toxicity resistance mechanisms.

Background: Crustacea, the second largest subphylum of Arthropoda, includes species of major ecological and economic importance, such as crabs, lobsters, crayfishes, shrimps, and barnacles. With the rapid development of crustacean aquaculture and biodiversity loss, understanding the gene regulatory mechanisms of growth, reproduction, and development of crustaceans is crucial to both aquaculture development and biodiversity conservation of this group of organisms. In these biological processes, transcription factors (TFs) play a vital role in regulating gene expression. However, crustacean transcription factors are still largely unknown, because the lack of complete genome sequences of most crustacean species hampers the studies on their transcriptional regulation on a system-wide scale. Thus, the current TF databases derived from genome sequences contain TF information for only a few crustacean species and are insufficient to elucidate the transcriptional diversity of such a large animal group., Results: Our database CrusTF ( http://qinlab.sls.cuhk.edu.hk/CrusTF ) provides comprehensive information for evolutionary and functional studies on the crustacean transcriptional regulatory system. CrusTF fills the knowledge gap of transcriptional regulation in crustaceans by exploring publicly available and newly sequenced transcriptomes of 170 crustacean species and identifying 131,941 TFs within 63 TF families. CrusTF features three categories of information: sequence, function, and evolution of crustacean TFs. The database enables searching, browsing and downloading of crustacean TF sequences. CrusTF infers DNA binding motifs of crustacean TFs, thus facilitating the users to predict potential downstream TF targets. The database also presents evolutionary analyses of crustacean TFs, which improve our understanding of the evolution of transcriptional regulatory systems in crustaceans., Conclusions: Given the importance of TF information in evolutionary and functional studies on transcriptional regulatory systems of crustaceans, this database will constitute a key resource for the research community of crustacean biology and evolutionary biology. Moreover, CrusTF serves as a model for the construction of TF database derived from transcriptome data. A similar approach could be applied to other groups of organisms, for which transcriptomes are more readily available than genomes.

In crustaceans, muscle growth and development is complicated, and to date substantial knowledge gaps exist. In this study, the claw muscle, hepatopancreas and nervous tissue of the mud crab (Scylla paramamosain) were collected at three fattening stages for sequence by the Illumina sequencing. A total of 127.87 Gb clean data with no less than 3.94 Gb generated for each sample and the cycleQ30 percentages were more than 86.13% for all samples. De Bruijn assembly of these clean data produced 94,853 unigenes, thereinto, 50,059 unigenes were found in claw muscle. A total of 121 differentially expressed genes (DEGs) were revealed in claw muscle from the three fattening stages with a Padj value < 0.01, including 63 genes with annotation. Functional annotation and enrichment analysis showed that the DEGs clusters represented the predominant gene catalog with roles in biochemical processes (glycolysis, phosphorylation and regulation of transcription), molecular function (ATP binding, 6-phosphofructokinase activity, and sequence-specific DNA binding) and cellular component (6-phosphofructokinase complex, plasma membrane, and integral component of membrane). qRT-PCR was employed to further validate certain DEGs. Single nucleotide polymorphism (SNP) analysis obtained 159,322, 125,963 and 166,279 potential SNPs from the muscle transcriptome at stage B, stage C and stage D, respectively. In addition, there were sixteen neuropeptide transcripts being predicted in the claw muscle. The present study provides a comprehensive transcriptome of claw muscle of S. paramamosain during fattening, providing a basis for screening the functional genes that may affect muscle growth of S. paramamosain.

Molecular markers are powerful tools for studying patterns of relatedness and parentage within populations and for making inferences about social evolution. However, the development of molecular markers for simultaneous study of multiple species presents challenges, particularly when species exhibit genome duplication or polyploidy. We developed microsatellite markers for Synalpheus shrimp, a genus in which species exhibit not only great variation in social organization, but also interspecific variation in genome size and partial genome duplication. From the four primary clades within Synalpheus, we identified microsatellites in the genomes of four species and in the consensus transcriptome of two species. Ultimately, we designed and tested primers for 143 microsatellite markers across 25 species. Although the majority of markers were disomic, many markers were polysomic for certain species. Surprisingly, we found no relationship between genome size and the number of polysomic markers. As expected, markers developed for a given species amplified better for closely related species than for more distant relatives. Finally, the markers developed from the transcriptome were more likely to work successfully and to be disomic than those developed from the genome, suggesting that consensus transcriptomes are likely to be conserved across species. Our findings suggest that the transcriptome, particularly consensus sequences from multiple species, can be a valuable source of molecular markers for taxa with complex, duplicated genomes., (© 2017 John Wiley & Sons Ltd.)

The Pacific white shrimp (Litopenaeus vannamei) is an important cultured crustacean species worldwide. However, little is known about the molecular mechanism of this species involved in the response to cold stress. In this study, four separate RNA-Seq libraries of L. vannamei were generated from 13°C stress and control temperature. Total 29,662 of Unigenes and overall of 19,619 annotated genes were obtained. Three comparisons were carried out among the four libraries, in which 72 of the top 20% of differentially-expressed genes were obtained, 15 GO and 5 KEGG temperature-sensitive pathways were fished out. Catalytic activity (GO: 0003824) and Metabolic pathways (ko01100) were the most annotated GO and KEGG pathways in response to cold stress, respectively. In addition, Calcium, MAPK cascade, Transcription factor and Serine/threonine-protein kinase signal pathway were picked out and clustered. Serine/threonine-protein kinase signal pathway might play more important roles in cold adaptation, while other three signal pathway were not widely transcribed. Our results had summarized the differentially-expressed genes and suggested the major important signaling pathways and related genes. These findings provide the first profile insight into the molecular basis of L. vannamei response to cold stress.

Elucidating the genetic mechanisms of adaptation to the hydrothermal vent in organisms at genomic level is significant for understanding the adaptive evolution process in the extreme environment. We performed RNA-seq on four different tissues of a vent crab species, Austinograea alayseae, producing 725,461 unigenes and 134,489 annotated genes. Genes related to sensory, circadian rhythm, hormone, hypoxia stress, metal detoxification and immunity were identified. It was noted that in the degenerated eyestalk, transcription of phototransduction related genes which are important for retinal function was greatly reduced; three crucial neuropeptide hormones, one molt-inhibiting and two crustacean hyperglycemic hormone precursors were characterized with conserved domains; hypoxia-inducible factor 1 and two novel isoforms of metallothioneins in the vent crabs were discovered. An analysis of 6,932 orthologs among three crabs A. alayseae, Portunus trituberculutus and Eriocheir sinensis revealed 19 positive selected genes (PSGs). Most of the PSGs were involved in immune responses, such as crustins and anti-lipopolysaccharide factor, suggesting their function in the adaptation to environment. The characterization of the first vent crab transcriptome provides abundant resources for genetic and evolutionary studies of this species, and paves the way for further investigation of vent adaptation process in crabs.

Bathynellacea (Crustacea, Syncarida, Parabathynellidae) are subterranean aquatic crustaceans that typically inhabit freshwater interstitial spaces (e.g., groundwater) and are occasionally found in caves and even hot springs. In this study, we sequenced the whole transcriptome of Allobathynella bangokensis using RNA-seq. De novo sequence assembly produced 74,866 contigs including 28,934 BLAST hits. Overall, the gene sequences were most similar to those of the waterflea Daphnia pulex. In the A. bangokensis transcriptome, no opsin or related sequences were identified, and no contig aligned to the crustacean visual opsins and non-visual opsins (i.e. arthropsins, peropsins, and melaopsins), suggesting potential regressive adaptation to the dark environment. However, A. bangokensis expressed conserved gene family sets, such as heat shock proteins and those related to key innate immunity pathways and antioxidant defense systems, at the transcriptional level, suggesting that this species has evolved adaptations involving molecular mechanisms of homeostasis. The transcriptomic information of A. bangokensis will be useful for investigating molecular adaptations and response mechanisms to subterranean environmental conditions., Competing Interests: The authors have declared that no competing interests exist.

Despite their economic, ecological, and experimental importance, genomic resources remain scarce for crustaceans. In lieu of genomes, many researchers have taken advantage of technological advancements to instead sequence and assemble crustacean transcriptomes de novo However, there is little consensus on what standard operating procedures are, or should be, for the field. Here, we systematically reviewed 53 studies published during 2014-2015 that utilized transcriptomic resources from this taxonomic group in an effort to identify commonalities as well as potential weaknesses that have applicability beyond just crustaceans. In general, these studies utilized RNA-Seq data, both novel and publicly available, to characterize transcriptomes and/or identify differentially expressed genes (DEGs) between treatments. Although the software suite Trinity was popular in assembly pipelines and other programs were also commonly employed, many studies failed to report crucial details regarding bioinformatic methodologies, including read mappers and the utilized parameters in identifying and characterizing DEGs. Annotation percentages for assembled transcriptomic contigs were low, averaging 32% overall. While other metrics, such as numbers of contigs and DEGs reported, correlated with the number of sequence reads utilized per sample, these did reach apparent saturation with increasing sequencing depth. Most disturbingly, a number of studies (55%) reported DEGs based on non-replicated experimental designs and single biological replicates for each treatment. Given this, we suggest future RNA-Seq experiments targeting transcriptome characterization conduct deeper (i.e., 50-100 M reads) sequencing while those examining differential expression instead focus more on increased biological replicates at shallower (i.e., ∼10-20 M reads/sample) sequencing depths. Moreover, the community must avoid submitting for review, or accepting for publication, non-replicated differential expression studies. Finally, mining the ever growing publicly available transcriptomic data from crustaceans will allow future studies to focus on hypothesis-driven research instead of continuing to simply characterize transcriptomes. As an example of this, we utilized neurotoxin sequences from the recently described remipede venom gland transcriptome in conjunction with publicly available crustacean transcriptomic data to derive preliminary results and hypotheses regarding the evolution of venom in crustaceans., (© The Author 2016. Published by Oxford University Press on behalf of the Society for Integrative and Comparative Biology. All rights reserved. For permissions please email: journals.permissions@oup.com.)

Background: The giant freshwater prawn, Macrobrachium rosenbergii, is a decapod crustacean that is commercially important as a food source. Farming of commercial crustaceans requires an efficient management strategy because the animals are easily subjected to stress and diseases during the culture. Autophagy, a stress response process, is well-documented and conserved in most animals, yet it is poorly studied in crustaceans., Results: In this study, we have performed an in silico search for transcripts encoding autophagy-related (Atg) proteins within various tissue transcriptomes of M. rosenbergii. Basic Local Alignment Search Tool (BLAST) search using previously known Atg proteins as queries revealed 41 transcripts encoding homologous M. rosenbergii Atg proteins. Among these Atg proteins, we selected commonly used autophagy markers, including Beclin 1, vacuolar protein sorting (Vps) 34, microtubule-associated proteins 1A/1B light chain 3B (MAP1LC3B), p62/sequestosome 1 (SQSTM1), and lysosomal-associated membrane protein 1 (Lamp-1) for further sequence analyses using comparative alignment and protein structural prediction. We found that crustacean autophagy marker proteins contain conserved motifs typical of other animal Atg proteins. Western blotting using commercial antibodies raised against human Atg marker proteins indicated their presence in various M. rosenbergii tissues, while immunohistochemistry localized Atg marker proteins within ovarian tissue, specifically late stage oocytes., Conclusions: This study demonstrates that the molecular components of autophagic process are conserved in crustaceans, which is comparable to autophagic process in mammals. Furthermore, it provides a foundation for further studies of autophagy in crustaceans that may lead to more understanding of the reproduction- and stress-related autophagy, which will enable the efficient aquaculture practices.

We sequenced the whole transcriptome of the marine mysid Neomysis awatschensis (Crustacea: Mysidacea) using Illumina RNA-seq. De novo assembly was performed with 14,018,702 raw reads using Trinity, resulting in 82,434 contigs. Transdecoder found 22,141 candidate coding contigs with homology to other species by BLAST analysis. Functional gene annotation was performed by Gene Ontology, InterProScan, and KEGG pathway analyses. We generated an expressed gene catalog for the mysid N. awatschensis to serve as a resource for marine environmental genomic and ecotoxicogenomic studies focused on uncovering the molecular mechanisms underlying the responses of N. awatschensis to environmental stressors and chemicals., (Copyright © 2016 Elsevier B.V. All rights reserved.)

Acclimation to low O2 in many organisms involves changes at the level of the transcriptome. Here we used high-throughput RNA sequencing (RNA-Seq) to explore the global transcriptomic response and specific involvement of a suite of hemocyanin (Hc) subunits to low O2 alone and in combination with high CO2, which naturally co-occurs with low O2. Hepatopancreas mRNA of juvenile L. vannamei exposed to air-saturated water, low O2, or low O2/high CO2 for 4 or 24 h was pooled, sequenced (HiSeq 2500) and assembled (Trinity: 52,190 contigs) to create a deep strand-specific reference transcriptome. Annotation of the assembly revealed sequences encoding the previously described small Hc subunit (HcS), and three full-length isoforms of the large subunit (HcL1-3). In addition to this, a previously unidentified full-length Hc subunit was discovered. Phylogenetic analysis demonstrated the subunit to be a β-type Hc subunit (denoted HcB), making this the first report of a β-type hemocyanin subunit in the Penaeoidea. RNAs of individual shrimp were sequenced; regulated genes identified from pairwise comparisons demonstrated a distinct pattern of regulation between prolonged low O2 and low O2/high CO2 treatments by GO term enrichment analysis (Roff-Bentzen, P < 0.0001), showcasing the stabilization of energetically costly translational machinery, mobilization of energy stores, and downregulation of the ubiquitin/proteasomal degradation machinery. Exposure to hypoxia for 24 h resulted in an increase in all of the full-length hemocyanin subunits (HcS, HcL1, HcL2, HcL3, and HcB). The addition of CO2 to hypoxia muted the transcriptomic response of all the Hc subunits to low O2, except for the β-type subunit., (Copyright © 2015 the American Physiological Society.)

Enhancing the production of aquatic animals is crucial for fishery management and aquaculture applications. Ovaries are specialized tissues that play critical roles in producing oocytes and hormones. Significant biochemical changes take place during the sexual maturation of Portunus trituberculatus, but the genetics of this process has not been extensively studied. Transcriptome sequencing can be used to determine gene expression changes within specific periods. In the current study, we used transcriptome sequencing to produce a comprehensive transcript dataset for the ovarian development of P. trituberculatus. Approximately 100 million sequencing reads were generated, and 126,075 transcripts were assembled. Functional annotation of the obtained transcripts revealed important pathways in ovarian development, such as those involving the vitellogenin gene. Also, we performed deep sequencing of ovaries in phases III and IV of sexual maturation in P. trituberculatus. Differential analysis of gene expression identified 506 significantly differentially expressed genes, which belong to 20 pathway, transporters, development, transcription factors, metabolism of other amino acids, carbohydrate and lipid, solute carrier family members, and enzymes. Taken together, our study provides the first comprehensive transcriptomic resource for P. trituberculatus ovaries, which will strengthen understanding of the molecular mechanisms underlying the sexual maturation process and advance molecular nutritional studies of P. trituberculatus.

Portunus trituberculatus is an important farmed crustacean in China. Growth traits are important economic traits that are controlled by many genes. The aim of this study was to screen growth-related genes using Illumina paired-end sequencing technology, and to identify growth-related SNPs in P. trituberculatus by association analysis. A total of 36,599,712 and 29,247,160 qualified Illumina read pairs for small size (SG) and large size crabs (LG) were obtained, giving rise to total clean bases of 7.14 and 5.72 G, respectively. A gene expression comparison showed that a total of 117 unigenes were differentially expressed between SG and LG, including 49 up-regulated and 68 down-regulated unigenes in SG and LG, respectively. GO functional enrichment analysis results showed up-regulated genes in LG were involved in three biological processes. These were cell communication, ovarian follicle cell development and cell division. Down-regulated genes in LG were involved in five biological processes, which were chitin metabolism, carbohydrate metabolism, ATP hydrolysis coupled proton transport, the steroid hormone mediated signaling pathway and fatty acid biosynthesis. Association analysis identified one growth-associated SNP, which was located in hemocyanin (comp58070-R31). This research has provided information on key genes, SNPs and biological processes related to growth, which will increase the efficiency of selective breeding programs for this important species., (Copyright © 2015 Elsevier B.V. All rights reserved.)

Background: ATP-binding cassette (ABC) protein family encode for membrane proteins involved in the transport of various biomolecules through the cellular membrane. These proteins have been identified in all taxa and present important physiological functions, including the process of insecticide detoxification in arthropods. For that reason the ectoparasite Caligus rogercresseyi represents a model species for understanding the molecular underpinnings involved in insecticide drug resistance., Methods: llumina sequencing was performed using sea lice exposed to 2 and 3 ppb of deltamethrin and azamethiphos. Contigs obtained from de novo assembly were annotated by Blastx. RNA-Seq analysis was performed and validated by qPCR analysis., Results: From the transcriptome database of C. rogercresseyi, 57 putative members of ABC protein sequences were identified and phylogenetically classified into the eight subfamilies described for ABC transporters in arthropods. Transcriptomic profiles for ABC proteins subfamilies were evaluated throughout C. rogercresseyi development. Moreover, RNA-Seq analysis was performed for adult male and female salmon lice exposed to the delousing drugs azamethiphos and deltamethrin. High transcript levels of the ABCB and ABCC subfamilies were evidenced. Furthermore, SNPs mining was carried out for the ABC proteins sequences, revealing pivotal genomic information., Conclusions: The present study gives a comprehensive transcriptome analysis of ABC proteins from C. rogercresseyi, providing relevant information about transporter roles during ontogeny and in relation to delousing drug responses in salmon lice. This genomic information represents a valuable tool for pest management in the Chilean salmon aquaculture industry.

Within the larval period of Eriocheir sinensis, there is pronounced morphological changes upon the molt from the fifth zoeae (Z5) to megalopae (M), and low survival rate exists during this transition, which is typical in crab species. RNA sequencing was applied to Z5 and M of E. sinensis, resulting in the discovery of 19,186 unigenes and 652 differentially expressed genes (DEGs, 3.40% of the unigenes). The important metabolic pathways that might play roles in the larval development of E. sinensis from Z5 to M were detected to be 'Xenobiotics Biodegradation and Metabolism (8.16%)', 'Metabolism of Cofactors and Vitamins (6.70%)', 'Lipid Metabolism (6.36%)', and 'Amino Acid Metabolism (6.28%)'. Further, 19 DEGs possibly contributing to the morphological and sensory capability changes of the larvae were identified, like multiple copies of cuticle protein genes, retinaldehyde-binding protein 1 (RLBP1), envelope protein (Envelope) and hormone-related gene ecdysteroid-regulated 16kDa protein (ESR16). Moreover, 62 DEGs were identified to be related to carbohydrate, lipid and protein digestion and metabolism, such as glucose dehydrogenases (GDHs), lipases (LIPs) and serine proteases (SPs). Among these DEGs, more genes related to the substance metabolism were found up-regulated in Z5 than M, suggesting that more energy might be essential to be released for Z5 to complete the transition into M. Characterization of the crucial DEGs by real-time quantitative PCR re-conformed their expression pattern. This study provides the first genome-wide transcriptomic analysis of E. sinensis Z5 and M for studying the molecular basis of the larvae metamorphosis and nutrition metabolism., (Copyright © 2014 Elsevier Inc. All rights reserved.)

MicroRNAs have been identified and analyzed in various model species, but an investigation of miRNAs in nonmodel species is required for a more complete understanding of miRNA evolution. In this study, we investigated the miRNAs of the nonmodel species Triops cancriformis (tadpole shrimp), a "living fossil," whose morphological form has not changed in almost 200 million years. Dramatic ontogenetic changes occur during its development. To clarify the evolution of miRNAs, we comparatively analyzed its miRNAs and the components of its RNAi machinery. We used deep sequencing to analyze small RNA libraries from the six different developmental stages of T. cancriformis (egg, first-fourth instars, and adult), and also analyzed its genomic DNA with deep sequencing. We identified 180 miRNAs (87 conserved miRNAs and 93 novel candidate miRNAs), and deduced the components of its RNAi machinery: the DICER1, AGO1-3, PIWI, and AUB proteins. A comparative miRNA analysis of T. cancriformis and Drosophila melanogaster showed inconsistencies in the expression patterns of four conserved miRNAs. This suggests that although the miRNA sequences of the two species are very similar, their roles differ across the species. An miRNA conservation analysis revealed that most of the conserved T. cancriformis miRNAs share sequence similarities with those of arthropods, although T. cancriformis is called a "living fossil." However, we found that let-7 and DICER1 of T. cancriformis are more similar to those of the vertebrates than to those of the arthropods. These results suggest that miRNA systems of T. cancriformis have evolved in a unique fashion., (© 2015 Ikeda et al.; Published by Cold Spring Harbor Laboratory Press for the RNA Society.)

We present a new transcriptome assembly of the Pacific whiteleg shrimp (Litopenaeus vannamei), the species most farmed for human consumption. Its functional annotation, a substantial improvement over previous ones, is provided freely. RNA-Seq with Illumina HiSeq technology was used to analyze samples extracted from shrimp abdominal muscle, hepatopancreas, gills and pleopods. We used the Trinity and Trinotate software suites for transcriptome assembly and annotation, respectively. The quality of this assembly and the affiliated targeted homology searches greatly enrich the curated transcripts currently available in public databases for this species. Comparison with the model arthropod Daphnia allows some insights into defining characteristics of decapod crustaceans. This large-scale gene discovery gives the broadest depth yet to the annotated transcriptome of this important species and should be of value to ongoing genomics and immunogenetic resistance studies in this shrimp of paramount global economic importance.

This article documents the public availability of transcriptome sequence data, assembled contigs and annotated contigs and ORFS for five anchialine shrimp species (Crustacea: Atyidae and Alpheidae)., (© 2014 John Wiley & Sons Ltd.)

Penaeid shrimp has a distinctive metamorphosis stage during early development. Although morphological and biochemical studies about this ontogeny have been developed for decades, researches on gene expression level are still scarce. In this study, we have investigated the transcriptomes of five continuous developmental stages in Pacific white shrimp (Litopenaeus vannamei) with high throughput Illumina sequencing technology. The reads were assembled and clustered into 66,815 unigenes, of which 32,398 have putative homologues in nr database, 14,981 have been classified into diverse functional categories by Gene Ontology (GO) annotation and 26,257 have been associated with 255 pathways by KEGG pathway mapping. Meanwhile, the differentially expressed genes (DEGs) between adjacent developmental stages were identified and gene expression patterns were clustered. By GO term enrichment analysis, KEGG pathway enrichment analysis and functional gene profiling, the physiological changes during shrimp metamorphosis could be better understood, especially histogenesis, diet transition, muscle development and exoskeleton reconstruction. In conclusion, this is the first study that characterized the integrated transcriptomic profiles during early development of penaeid shrimp, and these findings will serve as significant references for shrimp developmental biology and aquaculture research.

Animal venoms have evolved many times. Venomous species are especially common in three of the four main groups of arthropods (Chelicerata, Myriapoda, and Hexapoda), which together represent tens of thousands of species of venomous spiders, scorpions, centipedes, and hymenopterans. Surprisingly, despite their great diversity of body plans, there is no unambiguous evidence that any crustacean is venomous. We provide the first conclusive evidence that the aquatic, blind, and cave-dwelling remipede crustaceans are venomous and that venoms evolved in all four major arthropod groups. We produced a three-dimensional reconstruction of the venom delivery apparatus of the remipede Speleonectes tulumensis, showing that remipedes can inject venom in a controlled manner. A transcriptomic profile of its venom glands shows that they express a unique cocktail of transcripts coding for known venom toxins, including a diversity of enzymes and a probable paralytic neurotoxin very similar to one described from spider venom. We screened a transcriptomic library obtained from whole animals and identified a nontoxin paralog of the remipede neurotoxin that is not expressed in the venom glands. This allowed us to reconstruct its probable evolutionary origin and underlines the importance of incorporating data derived from nonvenom gland tissue to elucidate the evolution of candidate venom proteins. This first glimpse into the venom of a crustacean and primitively aquatic arthropod reveals conspicuous differences from the venoms of other predatory arthropods such as centipedes, scorpions, and spiders and contributes valuable information for ultimately disentangling the many factors shaping the biology and evolution of venoms and venomous species.

Pacific white shrimp (Litopenaeus vannamei) is the most extensively farmed crustacean species in the world. White spot syndrome virus (WSSV) is one of the major pathogens in the cultured shrimp. However, the molecular mechanisms of the host-virus interaction remain largely unknown. In this study, the impact of WSSV infection on host gene expression in the hepatopancreas of L. vannamei was investigated through the use of 454 pyrosequencing-based RNA-Seq of cDNA libraries developed from WSSV-challenged shrimp or normal controls. By comparing the two cDNA libraries, we show that 767 host genes are significantly up-regulated and 729 genes are significantly down-regulated by WSSV infection. KEGG analysis of the differentially expressed genes indicated that the distribution of gene pathways between the up- and down-regulated genes is quite different. Among the differentially expressed genes, several are found to be involved in various processes of animal defense against pathogens such as apoptosis, mitogen-activated protein kinase (MAPK) signaling, toll-like receptor (TLR) signaling, Wnt signaling and antigen processing and presentation pathways. The present study provides valuable information on differential expression of L. vannamei genes following WSSV infection and improves our current understanding of this host-virus interaction. In addition, the large number of transcripts obtained in this study provides a strong basis for future genomic research on shrimp.

Background: The Chinese mitten crab Eriocheir sinensis is an important economic crustacean and has been seriously attacked by various diseases, which requires more and more information for immune relevant genes on genome background. Recently, high-throughput RNA sequencing (RNA-seq) technology provides a powerful and efficient method for transcript analysis and immune gene discovery., Methods/principal Findings: A cDNA library from hepatopancreas of E. sinensis challenged by a mixture of three pathogen strains (Gram-positive bacteria Micrococcus luteus, Gram-negative bacteria Vibrio alginolyticus and fungi Pichia pastoris; 10(8) cfu·mL(-1)) was constructed and randomly sequenced using Illumina technique. Totally 39.76 million clean reads were assembled to 70,300 unigenes. After ruling out short-length and low-quality sequences, 52,074 non-redundant unigenes were compared to public databases for homology searching and 17,617 of them showed high similarity to sequences in NCBI non-redundant protein (Nr) database. For function classification and pathway assignment, 18,734 (36.00%) unigenes were categorized to three Gene Ontology (GO) categories, 12,243 (23.51%) were classified to 25 Clusters of Orthologous Groups (COG), and 8,983 (17.25%) were assigned to six Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. Potentially, 24, 14, 47 and 132 unigenes were characterized to be involved in Toll, IMD, JAK-STAT and MAPK pathways, respectively., Conclusions/significance: This is the first systematical transcriptome analysis of components relating to innate immune pathways in E. sinensis. Functional genes and putative pathways identified here will contribute to better understand immune system and prevent various diseases in crab.

An ambitious, yet fundamental goal for comparative biology is to understand the evolutionary relationships for all of life. However, many important taxonomic groups have remained recalcitrant to inclusion into broader scale studies. Here, we focus on collection of 9 new 454 transcriptome data sets from Ostracoda, an ancient and diverse group with a dense fossil record, which is often undersampled in broader studies. We combine the new transcriptomes with a new morphological matrix (including fossils) and existing expressed sequence tag, mitochondrial genome, nuclear genome, and ribosomal DNA data. Our analyses lead to new insights into ostracod and pancrustacean phylogeny. We obtained support for three epic pancrustacean clades that likely originated in the Cambrian: Oligostraca (Ostracoda, Mystacocarida, Branchiura, and Pentastomida); Multicrustacea (Copepoda, Malacostraca, and Thecostraca); and a clade we refer to as Allotriocarida (Hexapoda, Remipedia, Cephalocarida, and Branchiopoda). Within the Oligostraca clade, our results support the unresolved question of ostracod monophyly. Within Multicrustacea, we find support for Thecostraca plus Copepoda, for which we suggest the name Hexanauplia. Within Allotriocarida, some analyses support the hypothesis that Remipedia is the sister taxon to Hexapoda, but others support Branchiopoda + Cephalocarida as the sister group of hexapods. In multiple different analyses, we see better support for equivocal nodes using slow-evolving genes or when excluding distant outgroups, highlighting the increased importance of conditional data combination in this age of abundant, often anonymous data. However, when we analyze the same set of species and ignore rate of gene evolution, we find higher support when including all data, more in line with a "total evidence" philosophy. By concatenating molecular and morphological data, we place pancrustacean fossils in the phylogeny, which can be used for studies of divergence times in Pancrustacea, Arthropoda, or Metazoa. Our results and new data will allow for attributes of Ostracoda, such as its amazing fossil record and diverse biology, to be leveraged in broader scale comparative studies. Further, we illustrate how adding extensive next-generation sequence data from understudied groups can yield important new phylogenetic insights into long-standing questions, especially when carefully analyzed in combination with other data.

Background: The Pacific white shrimp, Litopenaeus vannamei, is a worldwide cultured crustacean species with important commercial value. Over the last two decades, Taura syndrome virus (TSV) has seriously threatened the shrimp aquaculture industry in the Western Hemisphere. To better understand the interaction between shrimp immune and TSV, we performed a transcriptome analysis in the hepatopancreas of L. vannamei challenged with TSV, using the 454 pyrosequencing (Roche) technology., Methodology/principal Findings: We obtained 126919 and 102181 high-quality reads from TSV-infected and non-infected (control) L. vannamei cDNA libraries, respectively. The overall de novo assembly of cDNA sequence data generated 15004 unigenes, with an average length of 507 bp. Based on BLASTX search (E-value <10-5) against NR, Swissprot, GO, COG and KEGG databases, 10425 unigenes (69.50% of all unigenes) were annotated with gene descriptions, gene ontology terms, or metabolic pathways. In addition, we identified 770 microsatellites and designed 497 sets of primers. Comparative genomic analysis revealed that 1311 genes differentially expressed in the infected shrimp compared to the controls, including 559 up- and 752 down- regulated genes. Among the differentially expressed genes, several are involved in various animal immune functions, such as antiviral, antimicrobial, proteases, protease inhibitors, signal transduction, transcriptional control, cell death and cell adhesion., Conclusions/significance: This study provides valuable information on shrimp gene activities against TSV infection. Results can contribute to the in-depth study of candidate genes in shrimp immunity, and improves our current understanding of this host-virus interaction. In addition, the large amount of transcripts reported in this study provide a rich source for identification of novel genes in shrimp.

Background: Arthropods are the most diverse animal phylum, but their genomic resources are relatively few. While the genome of the branchiopod Daphnia pulex is now available, no other large-scale crustacean genomic resources are available for comparison. In particular, genomic resources are lacking for the most tractable laboratory model of crustacean development, the amphipod Parhyale hawaiensis. Insight into shared and divergent characters of crustacean genomes will facilitate interpretation of future developmental, biomedical, and ecological research using crustacean models., Results: To generate a transcriptome enriched for maternally provided and zygotically transcribed developmental genes, we created cDNA from ovaries and embryos of P. hawaiensis. Using 454 pyrosequencing, we sequenced over 1.1 billion bases of this cDNA, and assembled them de novo to create, to our knowledge, the second largest crustacean genomic resource to date. We found an unusually high proportion of C2H2 zinc finger-containing transcripts, as has also been reported for the genome of the pea aphid Acyrthosiphon pisum. Consistent with previous reports, we detected trans-spliced transcripts, but found that they did not noticeably impact transcriptome assembly. Our assembly products yielded 19,067 unique BLAST hits against nr (E-value cutoff e-10). These included over 400 predicted transcripts with significant similarity to D. pulex sequences but not to sequences of any other animal. Annotation of several hundred genes revealed P. hawaiensis homologues of genes involved in development, gametogenesis, and a majority of the members of six major conserved metazoan signaling pathways., Conclusions: The amphipod P. hawaiensis has higher transcript complexity than known insect transcriptomes, and trans-splicing does not appear to be a major contributor to this complexity. We discuss the importance of a reliable comparative genomic framework within which to consider findings from new crustacean models such as D. pulex and P. hawaiensis, as well as the need for development of further substantial crustacean genomic resources.