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.

Construction and Characterization of Two Novel Transcriptome Assemblies in the Congeneric Porcelain Crabs Petrolisthes cinctipes and P. manimaculis

Crustaceans have commonly been used as non-model systems in basic biological research, especially physiological regulation. With the recent and rapid adoption of functional genomic tools, crustaceans are increasingly becoming model systems for ecological investigations of development and evolution and for mechanistic examinations of genotype-phenotype interactions and molecular pathways of response to environmental stressors. Comparative transcriptomic approaches, however, remain constrained by a lack of sequence data in closely related crustacean taxa. We identify challenges in the use of functional genomics tools in comparative analysis among decapod crustacean in light of recent advances. We present RNA-seq data from two congeneric species of porcelain crabs (Petrolisthes cinctipes and P. manimaculis) used to construct two de novo transcriptome assemblies with ∼194K and ∼278K contigs, respectively. We characterize and contrast these assemblies and compare them to a previously generated EST sequence library for P. cinctipes We also discuss the potential use of these data as a case-study system in the broader context of crustacean comparative transcriptomics.

Accessing transcriptomic data for ecologically important genes in the goose barnacle (Pollicipes pollicipes), with particular focus on cement proteins

In this study 4310 expressed sequence tags (ESTs) were used to identify potentially useful transcripts for future studies in the gooseneck barnacle Pollicipes pollicipes (Gmelin, 1789). 119 ESTs were obtained in this work and 4191 were taken from Meusemann et al. (2010). The gooseneck barnacle is a sessile pedunculate cirripede of great economic importance that occurs in dense aggregations, and is harvested for human consumption. The assembly of these ESTs yielded 1805 unigenes (461 contigs and 1344 singlets). The identification of cement proteins in our data is particularly interesting for cirripedes. Only a small part of the assembled unigenes could be functionally annotated. However, our results greatly improve our understanding of the biological features of P. pollicipes. In addition to this, a large number of potentially interesting genes were identified in order to serve as the base for future evolutionary studies in P. pollicipes.

Characterization and analysis of a transcriptome from the boreal spider crab Hyas araneus

Research investigating the genetic basis of physiological responses has significantly broadened our understanding of the mechanisms underlying organismic response to environmental change. However, genomic data are currently available for few taxa only, thus excluding physiological model species from this approach. In this study we report the transcriptome of the model organism Hyas araneus from Spitsbergen (Arctic). We generated 20,479 transcripts, using the 454 GS FLX sequencing technology in combination with an Illumina HiSeq sequencing approach. Annotation by Blastx revealed 7159 blast hits in the NCBI non-redundant protein database. The comparison between the spider crab H. araneus transcriptome and EST libraries of the European lobster Homarus americanus and the porcelain crab Petrolisthes cinctipes yielded 3229/2581 sequences with a significant hit, respectively. The clustering by the Markov Clustering Algorithm (MCL) revealed a common core of 1710 clusters present in all three species and 5903 unique clusters for H. araneus. The combined sequencing approaches generated transcripts that will greatly expand the limited genomic data available for crustaceans. We introduce the MCL clustering for transcriptome comparisons as a simple approach to estimate similarities between transcriptomic libraries of different size and quality and to analyze homologies within the selected group of species. In particular, we identified a large variety of reverse transcriptase (RT) sequences not only in the H. araneus transcriptome and other decapod crustaceans, but also sea urchin, supporting the hypothesis of a heritable, anti-viral immunity and the proposed viral fragment integration by host-derived RTs in marine invertebrates.

Changes in extreme cold tolerance, membrane composition and cardiac transcriptome during the first day of thermal acclimation in the porcelain crab Petrolisthes cinctipes

Intertidal zone organisms can experience transient freezing temperatures during winter low tides, but their extreme cold tolerance mechanisms are not known. Petrolisthes cinctipes is a temperate mid-high intertidal zone crab species that can experience wintertime habitat temperatures below the freezing point of seawater. We examined how cold tolerance changed during the initial phase of thermal acclimation to cold and warm temperatures, as well as the persistence of cold tolerance during long-term thermal acclimation. Thermal acclimation for as little as 6 h at 8°C enhanced cold tolerance during a 1 h exposure to -2°C relative to crabs acclimated to 18°C. Potential mechanisms for this enhanced tolerance were elucidated using cDNA microarrays to probe for differences in gene expression in cardiac tissue of warm- and cold-acclimated crabs during the first day of thermal acclimation. No changes in gene expression were detected until 12 h of thermal acclimation. Genes strongly upregulated in warm-acclimated crabs represented immune response and extracellular/intercellular processes, suggesting that warm-acclimated crabs had a generalized stress response and may have been remodelling tissues or altering intercellular processes. Genes strongly upregulated in cold-acclimated crabs included many that are involved in glucose production, suggesting that cold acclimation involves increasing intracellular glucose as a cryoprotectant. Structural cytoskeletal proteins were also strongly represented among the genes upregulated in only cold-acclimated crabs. There were no consistent changes in composition or the level of unsaturation of membrane phospholipid fatty acids with cold acclimation, which suggests that neither short- nor long-term changes in cold tolerance are mediated by changes in membrane fatty acid composition. Overall, our study demonstrates that initial changes in cold tolerance are likely not regulated by transcriptomic responses, but that gene-expression-related changes in homeostasis begin within 12 h, the length of a tidal cycle.

Construction and characterization of a normalized cDNA library from the river snail Bellamya aeruginosa after exposure to copper

The construction of a normalized cDNA library is a popular tool for identifying novel biomarkers for monitoring environmental pollution. In the present study, a normalized cDNA library was constructed from the river snail Bellamya aeruginosa after exposure to Cu(2+) by using the SMART technique. The titer of the cDNA library was 1.78 × 10(6) pfu/ml, with a recombinant efficiency of 95.8%. In addition, from 6,000 randomly selected and sequenced clones, 5,473 high-quality ESTs were identified. After processing the sequences, 3,961 unigenes representing 897 contigs and 3,064 singlets were obtained with 27.6% redundancy. Analysis of expressed sequenced tags using COG and GO annotation and KEGG pathway data showed that a large group of genes related to growth and development, signal transduction, and defense mechanisms were present in the cDNA library. Based on our findings, this normalized cDNA library will provide a valuable resource for further research on functional genes and ecotoxicology in B. aeruginosa.

EST analysis on the gonad development related organs and microarray screen for differentially expressed genes in mature ovary and testis of Scylla paramamosain

A total of 5160 high quality ESTs (expressed sequence tags) averaging 357 bp were collected from normalized cDNA libraries created from testis, ovary and mixed organs of mud crab Scylla paramamosain. Clustering and assembly of these ESTs resulted in a total of 3837 unique sequences with 576 overlapping contigs and 3261 singletons. Comparisons with the GenBank non-redundant (Nr) protein database (BLASTx, e-values <10(-5)) revealed putative functions or matched homologs from other organisms for 847 (22%) of the ESTs. Several gonad development related genes such as cathepsin C, thioredoxin peroxidase, vitellogenin receptor precursor, 50S ribosomal protein L24 and ubiquitin-conjugating enzyme E2 isoform 2 were identified from this EST project and demonstrated as gonad differential expression genes by rqRT-PCR. Sixty five different types of SSRs (simple sequence repeats) were identified from the total 411 EST-SSR motifs. A home-made cDNA microarray containing 5664 spots was developed and the hybridization results indicated that 39 unique transcripts were differentially expressed in testis and ovaries (P<0.05). The expression levels of eleven unique transcripts examined by rqRT-PCR were matched with microarray fairly. These results will provide a useful resource for functional genomic studies on the biology of reproduction of mud crab.

Using comparative genomic hybridization to survey genomic sequence divergence across species: a proof-of-concept from Drosophila

Background

Genome-wide analysis of sequence divergence among species offers profound insights into the evolutionary processes that shape lineages. When full-genome sequencing is not feasible for a broad comparative study, we propose the use of array-based comparative genomic hybridization (aCGH) in order to identify orthologous genes with high sequence divergence. Here we discuss experimental design, statistical power, success rate, sources of variation and potential confounding factors. We used a spotted PCR product microarray platform from Drosophila melanogaster to assess sequence divergence on a gene-by-gene basis in three fully sequenced heterologous species (D. sechellia, D. simulans, and D. yakuba). Because complete genome assemblies are available for these species this study presents a powerful test for the use of aCGH as a tool to measure sequence divergence.

Results

We found a consistent and linear relationship between hybridization ratio and sequence divergence of the sample to the platform species. At higher levels of sequence divergence (< 92% sequence identity to D. melanogaster) ~84% of features had significantly less hybridization to the array in the heterologous species than the platform species, and thus could be identified as "diverged". At lower levels of divergence (≥ 97% identity), only 13% of genes were identified as diverged. While ~40% of the variation in hybridization ratio can be accounted for by variation in sequence identity of the heterologous sample relative to D. melanogaster, other individual characteristics of the DNA sequences, such as GC content, also contribute to variation in hybridization ratio, as does technical variation.

Conclusions

Here we demonstrate that aCGH can accurately be used as a proxy to estimate genome-wide divergence, thus providing an efficient way to evaluate how evolutionary processes and genomic architecture can shape species diversity in non-model systems. Given the increased number of species for which microarray platforms are available, comparative studies can be conducted for many interesting lineages in order to identify highly diverged genes that may be the target of natural selection.

The porcelain crab transcriptome and PCAD, the porcelain crab microarray and sequence database

Background

With the emergence of a completed genome sequence of the freshwater crustacean Daphnia pulex, construction of genomic-scale sequence databases for additional crustacean sequences are important for comparative genomics and annotation. Porcelain crabs, genus Petrolisthes, have been powerful crustacean models for environmental and evolutionary physiology with respect to thermal adaptation and understanding responses of marine organisms to climate change. Here, we present a large-scale EST sequencing and cDNA microarray database project for the porcelain crab Petrolisthes cinctipes.

Methodology/Principal Findings

A set of ∼30K unique sequences (UniSeqs) representing ∼19K clusters were generated from ∼98K high quality ESTs from a set of tissue specific non-normalized and mixed-tissue normalized cDNA libraries from the porcelain crab Petrolisthes cinctipes. Homology for each UniSeq was assessed using BLAST, InterProScan, GO and KEGG database searches. Approximately 66% of the UniSeqs had homology in at least one of the databases. All EST and UniSeq sequences along with annotation results and coordinated cDNA microarray datasets have been made publicly accessible at the Porcelain Crab Array Database (PCAD), a feature-enriched version of the Stanford and Longhorn Array Databases.

Conclusions/Significance

The EST project presented here represents the third largest sequencing effort for any crustacean, and the largest effort for any crab species. Our assembly and clustering results suggest that our porcelain crab EST data set is equally diverse to the much larger EST set generated in the Daphnia pulex genome sequencing project, and thus will be an important resource to the Daphnia research community. Our homology results support the pancrustacea hypothesis and suggest that Malacostraca may be ancestral to Branchiopoda and Hexapoda. Our results also suggest that our cDNA microarrays cover as much of the transcriptome as can reasonably be captured in EST library sequencing approaches, and thus represent a rich resource for studies of environmental genomics.

Ecological genetics in the North Atlantic: environmental gradients and adaptation at specific loci

The North Atlantic intertidal community provides a rich set of organismal and environmental material for the study of ecological genetics. Clearly defined environmental gradients exist at multiple spatial scales: there are broad latitudinal trends in temperature, meso-scale changes in salinity along estuaries, and smaller scale gradients in desiccation and temperature spanning the intertidal range. The geology and geography of the American and European coasts provide natural replication of these gradients, allowing for population genetic analyses of parallel adaptation to environmental stress and heterogeneity. Statistical methods have been developed that provide genomic neutrality tests of population differentiation and aid in the process of candidate gene identification. In this paper, we review studies of marine organisms that illustrate associations between an environmental gradient and specific genetic markers. Such highly differentiated markers become candidate genes for adaptation to the environmental factors in question, but the functional significance of genetic variants must be comprehensively evaluated. We present a set of predictions about locus-specific selection across latitudinal, estuarine, and intertidal gradients that are likely to exist in the North Atlantic. We further present new data and analyses that support and contradict these simple selection models. Some taxa show pronounced clinal variation at certain loci against a background of mild clinal variation at many loci. These cases illustrate the procedures necessary for distinguishing selection driven by internal genomic vs. external environmental factors. We suggest that the North Atlantic intertidal community provides a model system for identifying genes that matter in ecology due to the clarity of the environmental stresses and an extensive experimental literature on ecological function. While these organisms are typically poor genetic and genomic models, advances in comparative genomics have provided access to molecular tools that can now be applied to taxa with well-defined ecologies. As many of the organisms we discuss have tight physiological limits driven by climatic factors, this synthesis of molecular population genetics with marine ecology could provide a sensitive means of assessing evolutionary responses to climate change.

Gene discovery from an ovary cDNA library of oriental river prawn Macrobrachium nipponense by ESTs annotation

The oriental river prawn, Macrobrachium nipponense, is an important crustacean species in aquaculture. However, early gonad maturity is a ubiquitous problem which devalues the product quality. While husbandry and nutritional management have achieved little success in tackling this issue, a molecular approach may discover the genes involved in reproduction and development, which will provide the basic knowledge on reproductive control. In this study, a high-quality cDNA library of prawn was constructed from the ovary tissue. A total of 3294 successful sequencing reactions yielded 3256 expressed sequence tags (ESTs) longer than 100 bp. The cluster and assembly analyses yielded 1514 unique sequences including 414 contigs and 1168 singletons. About 719 (47.49%) unique sequences were identified as orthologs of genes from other organisms. By sequence comparability analysis, 28 important genes including cathepsin B, chromobox protein, Cdc2, cyclin B, DEAD box protein and ADF/cofilin protein were expressed. These genes may be involved in reproductive and developmental functions in prawn. Peritrophin consisting of cortical rods was also found in this species. The identification of these EST sequences in M. nipponense would improve our understanding on the genes that regulate reproduction and development in prawn species. This study also lays the groundwork for development of molecular markers related to ovary development in other prawn species.

D- and L-lactate dehydrogenases during invertebrate evolution

BACKGROUND

The L-lactate and D-lactate dehydrogenases, which are involved in the reduction of pyruvate to L(-)-lactate and D(+)-lactate, belong to evolutionarily unrelated enzyme families. The genes encoding L-LDH have been used as a model for gene duplication due to the multiple paralogs found in eubacteria, archaebacteria, and eukaryotes. Phylogenetic studies have suggested that several gene duplication events led to the main isozymes of this gene family in chordates, but little is known about the evolution of L-Ldh in invertebrates. While most invertebrates preferentially oxidize L-lactic acid, several species of mollusks, a few arthropods and polychaetes were found to have exclusively D-LDH enzymatic activity. Therefore, it has been suggested that L-LDH and D-LDH are mutually exclusive. However, recent characterization of putative mammalian D-LDH with significant similarity to yeast proteins showing D-LDH activity suggests that at least mammals have the two naturally occurring forms of LDH specific to L- and D-lactate. This study describes the phylogenetic relationships of invertebrate L-LDH and D-LDH with special emphasis on crustaceans, and discusses gene duplication events during the evolution of L-Ldh.

RESULTS

Our phylogenetic analyses of L-LDH in vertebrates are consistent with the general view that the main isozymes (LDH-A, LDH-B and LDH-C) evolved through a series of gene duplications after the vertebrates diverged from tunicates. We report several gene duplication events in the crustacean, Daphnia pulex, and the leech, Helobdella robusta. Several amino acid sequences with strong similarity to putative mammalian D-LDH and to yeast DLD1 with D-LDH activity were found in both vertebrates and invertebrates.

CONCLUSION

The presence of both L-Ldh and D-Ldh genes in several chordates and invertebrates suggests that the two enzymatic forms are not necessarily mutually exclusive. Although, the evolution of L-Ldh has been punctuated by multiple events of gene duplication in both vertebrates and invertebrates, a shared evolutionary history of this gene in the two groups is apparent. Moreover, the high degree of sequence similarity among D-LDH amino acid sequences suggests that they share a common evolutionary history.

Expression analysis of light-regulated genes isolated from a full-length-enriched cDNA library of Onosma paniculatum cell cultures

Shikonin and its derivatives are formed in large amounts in dark-cultured Onosma paniculatum cells. In order to isolate and identify the genes regulating shikonin biosynthesis, we constructed and characterized a full-length-enriched cDNA library of dark-cultured cells by using the SMART (Switching Mechanism At 5'-end of RNA Transcript) cDNA synthesis and LD-PCR (long-distance PCR) strategies. The titer of the primary cDNA library was 1.04 x 10(6)pfu/mL with a recombination rate of 99.60%. Most of the cDNA inserts ranged from 1.0 to 2.5 kb, and 78.33% of the 76 randomly selected clones contained full-length coding regions. Expression analysis of randomly selected genes by small scale microarray revealed that 23 genes were down-regulated, including 17 genes with known functions, 2 genes with putative functions, and 4 novel genes, and that 3 genes were up-regulated (two-fold) in cells cultured under white light as compared with those cultured in the dark. Interestingly, two of the down-regulated genes, encoding aci-reductone dioxygenase (ARD)-like protein and ethylene responsive factor (ERF), are involved in ethylene biosynthesis and signal transduction, implying that ethylene might play an important role as a signal molecule in light-regulated shikonin formation. These data contribute to a better understanding of light-involvement in regulating the formation of plant secondary metabolites.

Hepatopancreatic multi-transcript expression patterns in the crayfish Cherax quadricarinatus during the moult cycle

Alterations of hepatopancreatic multi-transcript expression patterns, related to induced moult cycle, were identified in male Cherax quadricarinatus through cDNA microarray hybridizations of hepatopancreatic transcript populations. Moult was induced by X-organ sinus gland extirpation or by repeated injections of 20-hydroxyecdysone. Manipulated males were sacrificed at premoult or early postmoult, and a reference population was sacrificed at intermoult. Differentially expressed genes among the four combinations of two induction methods and two moult stages were identified. Biologically interesting clusters revealing concurrently changing transcript expressions across treatments were selected, characterized by a general shift of expression throughout premoult and early postmoult vs. intermoult, or by different premoult vs. postmoult expressions. A number of genes were differentially expressed in 20-hydroxyecdysone-injected crayfish vs. X-organ sinus gland extirpated males.