Developmental Transcriptomics of the Hawaiian Anchialine Shrimp Halocaridina rubra Holthuis, 1963 (Crustacea: Atyidae)

Multiple transcriptome mining coupled with tissue specific molecular cloning and mass spectrometry provide insights into agatoxin-like peptide conservation in decapod crustaceans

Over the past decade, in silico genome and transcriptome mining has led to the identification of many new crustacean peptide families, including the agatoxin-like peptides (ALPs), a group named for their structural similarity to agatoxin, a spider venom component. Here, analysis of publicly accessible transcriptomes was used to expand our understanding of crustacean ALPs. Specifically, transcriptome mining was used to investigate the phylogenetic/structural conservation, tissue localization, and putative functions of ALPs in decapod species. Transcripts encoding putative ALP precursors were identified from one or more members of the Penaeoidea (penaeid shrimp), Sergestoidea (sergestid shrimps), Caridea (caridean shrimp), Astacidea (clawed lobsters and freshwater crayfish), Achelata (spiny/slipper lobsters), and Brachyura (true crabs), suggesting a broad, and perhaps ubiquitous, conservation of ALPs in decapods. Comparison of the predicted mature structures of decapod ALPs revealed high levels of amino acid conservation, including eight identically conserved cysteine residues that presumably allow for the formation of four identically positioned disulfide bridges. All decapod ALPs are predicted to have amidated carboxyl-terminals. Two isoforms of ALP appear to be present in most decapod species, one 44 amino acids long and the other 42 amino acids in length, both likely generated by alternative splicing of a single gene. In carideans, a gene or terminal exon duplication appears to have occurred, with alternative splicing producing four ALPs, two 44 and two 42 amino acid isoforms. The identification of ALP precursor-encoding transcripts in nervous system-specific transcriptomes (e.g., Homarus americanus brain, eyestalk ganglia, and cardiac ganglion assemblies, finding confirmed using RT-PCR) suggests that members of this peptide family may serve as locally-released and/or hormonally-delivered neuromodulators in decapods. Their detection in testis- and hepatopancreas-specific transcriptomes suggests that members of the ALP family may also play roles in male reproduction and innate immunity/detoxification.

Proteolytic profile of larval developmental stages of Penaeus vannamei: An activity and mRNA expression approach

In arthropods, the cleavage of specific proteins by peptidases has pivotal roles in multiple physiological processes including oogenesis, immunity, nutrition, and parasitic infection. These enzymes are also key players in the larval development, and well-described triggers of molting and metamorphosis. In this work the peptidase complement throughout the larvae development of Penaeus vannamei was quantified at the transcript and activity level using qPCR and fluorogenic substrates designed to be hydrolyzed by class-specific peptidases respectively, providing a detailed identification of the proteolytic repertoire in P. vannamei larvae. Significant changes in the peptidase activity profile were observed. During the lecithotrophic naupliar instars, the dominant peptidase activity and expression derive from cysteine peptidases, suggesting that enzymes of this class hydrolyze the protein components of yolk as the primary amino acid source. At the first feeding instar, zoea, dominant serine peptidase activity was found where trypsin activity is particularly high, supporting previous observations that during zoea the breakdown of food protein is primarily enzymatic. At decapodid stages the peptidase expression and activity is more diverse indicating that a multienzyme network achieves food digestion. Our results suggest that proteolytic enzymes fulfill specific functions during P. vannamei larval development.

Red Coloration in an Anchialine Shrimp: Carotenoids, Genetic Variation, and Candidate Genes

Red coloration is a widely distributed phenotype among animals, yet the pigmentary and genetic bases for this phenotype have been described in relatively few taxa. Here we show that the Hawaiian endemic anchialine shrimp Halocaridina rubra is red because of the accumulation of astaxanthin. Laboratory colonies of phylogenetically distinct lineages of H. rubra have colony-specific amounts of astaxanthin that are developmentally, and likely genetically, fixed. Carotenoid supplementation and restriction experiments failed to change astaxanthin content from the within-colony baseline levels, suggesting that dietary limitation is not a major factor driving coloration differences. A possible candidate gene product predicted to be responsible for the production of astaxanthin in H. rubra and other crustaceans is closely related to the bifunctional cytochrome P450 family 3 enzyme CrtS found in fungi. However, homologs to the enzyme thought to catalyze ketolation reactions in birds and turtles, CYP2J19, were not found. This work is one of the first steps in linking phenotypic variation in red coloration of H. rubra to genotypic variation. Future work should focus on (1) pinpointing the genes that function in the bioconversion of dietary carotenoids to astaxanthin, (2) examining what genomic variants might drive variation in coloration among discrete lineages, and (3) testing more explicitly for condition-dependent carotenoid coloration in crustaceans.

Disparate responses to salinity across species and organizational levels in anchialine shrimps

Environmentally induced plasticity in gene expression is one of the underlying mechanisms of adaptation to habitats with variable environments. For example, euryhaline crustaceans show predictable changes in the expression of ion-transporter genes during salinity transfers, although studies have typically been limited to specific genes, taxa and ecosystems of interest. Here, we investigated responses to salinity change at multiple organizational levels in five species of shrimp representing at least three independent invasions of the anchialine ecosystem, defined as habitats with marine and freshwater influences with spatial and temporal fluctuations in salinity. Although all five species were generally strong osmoregulators, salinity-induced changes in gill physiology and gene expression were highly species specific. While some species exhibited patterns similar to those of previously studied euryhaline crustaceans, instances of distinct and atypical patterns were recovered from closely related species. Species-specific patterns were found when examining: (1) numbers and identities of differentially expressed genes, (2) salinity-induced expression of genes predicted a priori to play a role in osmoregulation, and (3) salinity-induced expression of orthologs shared among all species. Notably, ion transport genes were unchanged in the atyid Halocaridina rubra while genes normally associated with vision and light perception were among those most highly upregulated. Potential reasons for species-specific patterns are discussed, including variation among anchialine habitats in salinity regimes and divergent evolution in anchialine taxa. Underexplored mechanisms of osmoregulation in crustaceans revealed here by the application of transcriptomic approaches to ecologically and taxonomically understudied systems are also explored.

Resources and Recommendations for Using Transcriptomics to Address Grand Challenges in Comparative Biology

High-throughput RNA sequencing (RNA-seq) technology has become an important tool for studying physiological responses of organisms to changes in their environment. De novo assembly of RNA-seq data has allowed researchers to create a comprehensive catalog of genes expressed in a tissue and to quantify their expression without a complete genome sequence. The contributions from the "Tapping the Power of Crustacean Transcriptomics to Address Grand Challenges in Comparative Biology" symposium in this issue show the successes and limitations of using RNA-seq in the study of crustaceans. In conjunction with the symposium, the Animal Genome to Phenome Research Coordination Network collated comments from participants at the meeting regarding the challenges encountered when using transcriptomics in their research. Input came from novices and experts ranging from graduate students to principal investigators. Many were unaware of the bioinformatics analysis resources currently available on the CyVerse platform. Our analysis of community responses led to three recommendations for advancing the field: (1) integration of genomic and RNA-seq sequence assemblies for crustacean gene annotation and comparative expression; (2) development of methodologies for the functional analysis of genes; and (3) information and training exchange among laboratories for transmission of best practices. The field lacks the methods for manipulating tissue-specific gene expression. The decapod crustacean research community should consider the cherry shrimp, Neocaridina denticulata, as a decapod model for the application of transgenic tools for functional genomics. This would require a multi-investigator effort.

Tapping the Power of Crustacean Transcriptomics to Address Grand Challenges in Comparative Biology: An Introduction to the Symposium

Crustaceans, and decapods in particular (i.e., crabs, shrimp, and lobsters), are a diverse and ecologically and commercially important group of organisms. Understanding responses to abiotic and biotic factors is critical for developing best practices in aquaculture and assessing the effects of changing environments on the biology of these important animals. A relatively small number of decapod crustacean species have been intensively studied at the molecular level; the availability, experimental tractability, and economic relevance factor into the selection of a particular species as a model. Transcriptomics, using high-throughput next generation sequencing (NGS, coupled with RNA sequencing or RNA-seq) is revolutionizing crustacean biology. The 11 symposium papers in this volume illustrate how RNA-seq is being used to study stress response, molting and limb regeneration, immunity and disease, reproduction and development, neurobiology, and ecology and evolution. This symposium occurred on the 10th anniversary of the symposium, "Genomic and Proteomic Approaches to Crustacean Biology", held at the Society for Integrative and Comparative Biology 2006 meeting. Two participants in the 2006 symposium, the late Paul Gross and David Towle, were recognized as leaders who pioneered the use of molecular techniques that would ultimately foster the transcriptomics research reviewed in this volume. RNA-seq is a powerful tool for hypothesis-driven research, as well as an engine for discovery. It has eclipsed the technologies available in 2006, such as microarrays, expressed sequence tags, and subtractive hybridization screening, as the millions of "reads" from NGS enable researchers to de novo assemble a comprehensive transcriptome without a complete genome sequence. The symposium series concludes with a policy paper that gives an overview of the resources available and makes recommendations for developing better tools for functional annotation and pathway and network analysis in organisms in which the genome is not available or is incomplete.

Here We Are, But Where Do We Go? A Systematic Review of Crustacean Transcriptomic Studies from 2014-2015

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.

Salinity-induced changes in gene expression from anterior and posterior gills of Callinectes sapidus (Crustacea: Portunidae) with implications for crustacean ecological genomics

Decapods represent one of the most ecologically diverse taxonomic groups within crustaceans, making them ideal to study physiological processes like osmoregulation. However, prior studies have failed to consider the entire transcriptomic response of the gill - the primary organ responsible for ion transport - to changing salinity. Moreover, the molecular genetic differences between non-osmoregulatory and osmoregulatory gill types, as well as the hormonal basis of osmoregulation, remain underexplored. Here, we identified and characterized differentially expressed genes (DEGs) via RNA-Seq in anterior (non-osmoregulatory) and posterior (osmoregulatory) gills during high to low salinity transfer in the blue crab Callinectes sapidus, a well-studied model for crustacean osmoregulation. Overall, we confirmed previous expression patterns for individual ion transport genes and identified novel ones with salinity-mediated expression. Notable, novel DEGs among salinities and gill types for C. sapidus included anterior gills having higher expression of structural genes such as actin and cuticle proteins while posterior gills exhibit elevated expression of ion transport and energy-related genes, with the latter likely linked to ion transport. Potential targets among recovered DEGs for hormonal regulation of ion transport between salinities and gill types included neuropeptide Y and a KCTD16-like protein. Using publically available sequence data, constituents for a "core" gill transcriptome among decapods are presented, comprising genes involved in ion transport and energy conversion and consistent with salinity transfer experiments. Lastly, rarefication analyses lead us to recommend a modest number of sequence reads (~10-15M), but with increased biological replication, be utilized in future DEG analyses of crustaceans.