The Anadara trapezia transcriptome: a resource for molluscan physiological genomics

Molecular Insights into the Low Complexity Secreted Venom of Calliactis polypus

Sea anemones are venomous animals that rely on their venom for prey capture, defense against predators, and intraspecific competition. Currently, comprehensive molecular and evolutionary analyses of the toxin repertoire for sea anemones are limited by a lack of proteomic data for most species. In this study, proteo-transcriptomic analysis was used to expand our knowledge of the proteinaceous components of sea anemone venom by determining the secreted venom proteome of Calliactis polypus. Electromechanical stimulation was used to obtain the secreted venom of C. polypus. We identified a low complexity proteome that was dominated by toxins with similarity to known neurotoxins, as well as six novel toxin candidates. The novel putative toxin candidates were found to be taxonomically restricted to species from the superfamily Metridioidea. Furthermore, the secreted venom of C. polypus had only three putative toxins in common with the venom of acontia from the same species and little similarity with the secreted venom of closely related species. Overall, this demonstrates that regionalized and lineage-specific variability in toxin abundance is common among sea anemone species. Moreover, the limited complexity of the toxin repertoire found in C. polypus supports the idea that peptide neurotoxins make up the dominant toxin arsenal found in the venom of sea anemones.

An historical "wreck": A transcriptome assembly of the naval shipworm, Teredo navalis Linnaeus, 1978

Historically famous for their negative impact on human-built marine wood structures, mollusc shipworms play a central ecological role in marine ecosystems. Their association with bacterial symbionts, providing cellulolytic and nitrogen-fixing activities, underscores their exceptional wood-eating and wood-boring behaviours, improving energy transfer and the recycling of essential nutrients locked in the wood cellulose. Importantly, from a molecular standpoint, a minute of omic resources are available from this lineage of Bivalvia. Here, we produced and assembled a transcriptome from the globally distributed naval shipworm, Teredo navalis (family Teredinidae). The transcriptome was obtained by sequencing the total RNA from five equidistant segments of the whole body of a T. navalis specimen. The quality of the produced assembly was accessed with several statistics, revealing a highly contiguous (1194 N50) and complete (over 90% BUSCO scores for Eukaryote and Metazoan databases) transcriptome, with nearly 38,000 predicted ORF, more than half being functionally annotated. Our findings pave the way to investigate the unique evolutionary biology of these highly modified bivalves and lay the foundation for an adequate gene annotation of a full genome sequence of the species.

Venoms for all occasions: The functional toxin profiles of different anatomical regions in sea anemones are related to their ecological function

The phylum Cnidaria is the oldest extant venomous group and is defined by the presence of nematocysts, specialized organelles responsible for venom production and delivery. Although toxin peptides and the cells housing nematocysts are distributed across the entire animal, nematocyte and venom profiles have been shown to differ across morphological structures in actiniarians. In this study, we explore the relationship between patterns of toxin expression and the ecological roles of discrete anatomical structures in Telmatactis stephensoni. Specifically, using a combination of proteomic and transcriptomic approaches, we examined whether there is a direct correlation between the functional similarity of regions and the similarity of their associated toxin expression profiles. We report that the regionalization of toxin production is consistent with the partitioning of the ecological roles of venom across envenomating structures, and that three major functional regions are present in T. stephensoni: tentacles, epidermis and gastrodermis. Additionally, we find that most structures that serve similar functions not only have comparable putative toxin profiles but also similar nematocyst types. There was no overlap in the putative toxins identified using proteomics and transcriptomics, but the expression patterns of specific milked venom peptides were conserved across RNA-sequencing and mass spectrometry imaging data sets. Furthermore, based on our data, it appears that acontia of T. stephensoni may be transcriptionally inactive and only mature nematocysts are present in the distal portions of the threads. Overall, we find that the venom profile of different anatomical regions in sea anemones varies according to its ecological functions.

Tentacle Morphological Variation Coincides with Differential Expression of Toxins in Sea Anemones

Phylum Cnidaria is an ancient venomous group defined by the presence of cnidae, specialised organelles that serve as venom delivery systems. The distribution of cnidae across the body plan is linked to regionalisation of venom production, with tissue-specific venom composition observed in multiple actiniarian species. In this study, we assess whether morphological variants of tentacles are associated with distinct toxin expression profiles and investigate the functional significance of specialised tentacular structures. Using five sea anemone species, we analysed differential expression of toxin-like transcripts and found that expression levels differ significantly across tentacular structures when substantial morphological variation is present. Therefore, the differential expression of toxin genes is associated with morphological variation of tentacular structures in a tissue-specific manner. Furthermore, the unique toxin profile of spherical tentacular structures in families Aliciidae and Thalassianthidae indicate that vesicles and nematospheres may function to protect branched structures that host a large number of photosynthetic symbionts. Thus, hosting zooxanthellae may account for the tentacle-specific toxin expression profiles observed in the current study. Overall, specialised tentacular structures serve unique ecological roles and, in order to fulfil their functions, they possess distinct venom cocktails.

Comparative molecular analyses of select pH- and osmoregulatory genes in three freshwater crayfish Cherax quadricarinatus, C. destructor and C. cainii

Systemic acid-base balance and osmotic/ionic regulation in decapod crustaceans are in part maintained by a set of transport-related enzymes such as carbonic anhydrase (CA), Na⁺/K⁺-ATPase (NKA), H⁺-ATPase (HAT), Na⁺/K⁺/2Cl⁻ cotransporter (NKCC), Na⁺/Cl⁻/HCO\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}${}_{3}^{-}$\end{document}3− cotransporter (NBC), Na⁺/H⁺ exchanger (NHE), Arginine kinase (AK), Sarcoplasmic Ca⁺²-ATPase (SERCA) and Calreticulin (CRT). We carried out a comparative molecular analysis of these genes in three commercially important yet eco-physiologically distinct freshwater crayfish, Cherax quadricarinatus, C. destructor and C. cainii, with the aim to identify mutations in these genes and determine if observed patterns of mutations were consistent with the action of natural selection. We also conducted a tissue-specific expression analysis of these genes across seven different organs, including gills, hepatopancreas, heart, kidney, liver, nerve and testes using NGS transcriptome data. The molecular analysis of the candidate genes revealed a high level of sequence conservation across the three Cherax sp. Hyphy analysis revealed that all candidate genes showed patterns of molecular variation consistent with neutral evolution. The tissue-specific expression analysis showed that 46% of candidate genes were expressed in all tissue types examined, while approximately 10% of candidate genes were only expressed in a single tissue type. The largest number of genes was observed in nerve (84%) and gills (78%) and the lowest in testes (66%). The tissue-specific expression analysis also revealed that most of the master genes regulating pH and osmoregulation (CA, NKA, HAT, NKCC, NBC, NHE) were expressed in all tissue types indicating an important physiological role for these genes outside of osmoregulation in other tissue types. The high level of sequence conservation observed in the candidate genes may be explained by the important role of these genes as well as potentially having a number of other basic physiological functions in different tissue types.

Expression patterns of two carbonic anhydrase genes, Na+/K+-ATPase and V-type H+-ATPase, in the freshwater crayfish, Cherax quadricarinatus, exposed to low pH and high pH

Carbonic anhydrase (CA), Na+/K+-ATPase (NKA) and Vacuolar-type H+-ATPase (HAT) play vital roles in osmoregulation and pH balance in decapod crustaceans. As variable pH levels have a significant impact on the physiology of crustaceans, it is crucial to understand the mechanisms by which an animal maintains its internal pH. We examined expression patterns of cytoplasmic (CAc) and membrane-associated form (CAg) of CA, NKA α subunit and HAT subunit a in gills of freshwater crayfish, Cherax quadricarinatus, at three pH levels – 6.2, 7.2 (control) and 8.2 – over 24 h. Expression levels of CAc were significantly increased at low pH and decreased at high pH conditions 24 h after transfer. Expression increased at low pH after 12 h, and reached its maximum level by 24 h. CAg showed a significant increase in expression at 6 h after transfer at low pH. Expression of NKA significantly increased at 6 h after transfer to pH 6.2 and remained elevated for up to 24 h. Expression for HAT and NKA showed similar patterns, where expression significantly increased 6 h after transfer to low pH and remained significantly elevated throughout the experiment. Overall, CAc, CAg, NKA and HAT gene expression is induced at low pH conditions in freshwater crayfish.

Insights into the innate immunome of actiniarians using a comparative genomic approach

Background

Innate immune genes tend to be highly conserved in metazoans, even in early divergent lineages such as Cnidaria (jellyfish, corals, hydroids and sea anemones) and Porifera (sponges). However, constant and diverse selection pressures on the immune system have driven the expansion and diversification of different immune gene families in a lineage-specific manner. To investigate how the innate immune system has evolved in a subset of sea anemone species (Order: Actiniaria), we performed a comprehensive and comparative study using 10 newly sequenced transcriptomes, as well as three publically available transcriptomes, to identify the origins, expansions and contractions of candidate and novel immune gene families.

Results

We characterised five conserved genes and gene families, as well as multiple novel innate immune genes, including the newly recognised putative pattern recognition receptor CniFL. Single copies of TLR, MyD88 and NF-κB were found in most species, and several copies of IL-1R-like, NLR and CniFL were found in almost all species. Multiple novel immune genes were identified with domain architectures including the Toll/interleukin-1 receptor (TIR) homology domain, which is well documented as functioning in protein-protein interactions and signal transduction in immune pathways. We hypothesise that these genes may interact as novel proteins in immune pathways of cnidarian species. Novelty in the actiniarian immunome is not restricted to only TIR-domain-containing proteins, as we identify a subset of NLRs which have undergone neofunctionalisation and contain 3–5 N-terminal transmembrane domains, which have so far only been identified in two anthozoan species.

Conclusions

This research has significance in understanding the evolution and origin of the core eumetazoan gene set, including how novel innate immune genes evolve. For example, the evolution of transmembrane domain containing NLRs indicates that these NLRs may be membrane-bound, while all other metazoan and plant NLRs are exclusively cytosolic receptors. This is one example of how species without an adaptive immune system may evolve innovative solutions to detect pathogens or interact with native microbiota. Overall, these results provide an insight into the evolution of the innate immune system, and show that early divergent lineages, such as actiniarians, have a diverse repertoire of conserved and novel innate immune genes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3204-2) contains supplementary material, which is available to authorized users.

Sequencing and de novo assembly of visceral mass transcriptome of the critically endangered land snail Satsuma myomphala: Annotation and SSR discovery

Satsuma myomphala is critically endangered through loss of natural habitats, predation by natural enemies, and indiscriminate collection. It is a protected species in Korea but lacks genomic resources for an understanding of varied functional processes attributable to evolutionary success under natural habitats. For assessing the genetic information of S. myomphala, we performed for the first time, de novo transcriptome sequencing and functional annotation of expressed sequences using Illumina Next-Generation Sequencing (NGS) platform and bioinformatics analysis. We identified 103,774 unigenes of which 37,959, 12,890, and 17,699 were annotated in the PANM (Protostome DB), Unigene, and COG (Clusters of Orthologous Groups) databases, respectively. In addition, 14,451 unigenes were predicted under Gene Ontology functional categories, with 4581 assigned to a single category. Furthermore, 3369 sequences with 646 having Enzyme Commission (EC) numbers were mapped to 122 pathways in the Kyoto Encyclopedia of Genes and Genomes Pathway database. The prominent protein domains included the Zinc finger (C2H2-like), Reverse Transcriptase, Thioredoxin-like fold, and RNA recognition motif domain. Many unigenes with homology to immunity, defense, and reproduction-related genes were screened in the transcriptome. We also detected 3120 putative simple sequence repeats (SSRs) encompassing dinucleotide to hexanucleotide repeat motifs from >1kb unigene sequences. A list of PCR primers of SSR loci have been identified to study the genetic polymorphisms. The transcriptome data represents a valuable resource for further investigations on the species genome structure and biology. The unigenes information and microsatellites would provide an indispensable tool for conservation of the species in natural and adaptive environments.

Transcriptome sequencing and de novo characterization of Korean endemic land snail, Koreanohadra kurodana for functional transcripts and SSR markers

The Korean endemic land snail Koreanohadra kurodana (Gastropoda: Bradybaenidae) found in humid areas of broadleaf forests and shrubs have been considered vulnerable as the number of individuals are declining in recent years. The species is poorly characterized at the genomic level that limits the understanding of functions at the molecular and genetics level. In the present study, we performed de novo transcriptome sequencing to produce a comprehensive transcript dataset of visceral mass tissue of K. kurodana by the Illumina paired-end sequencing technology. Over 234 million quality reads were assembled to a total of 315,924 contigs and 191,071 unigenes, with an average and N50 length of 585.6 and 715 bp and 678 and 927 bp, respectively. Overall, 36.32 % of the unigenes found matches to known protein/nucleotide sequences in the public databases. The direction of the unigenes to functional categories was determined using COG, GO, KEGG, and InterProScan protein domain search. The GO analysis search resulted in 22,967 unigenes (12.02 %) being categorized into 40 functional groups. The KEGG annotation revealed that metabolism pathway genes were enriched. The most prominent protein motifs include the zinc finger, ribonuclease H, reverse transcriptase, and ankyrin repeat domains. The simple sequence repeats (SSRs) identified from >1 kb length of unigenes show a dominancy of dinucleotide repeat motifs followed with tri- and tetranucleotide motifs. A number of unigenes were putatively assessed to belong to adaptation and defense mechanisms including heat shock proteins 70, Toll-like receptor 4, AMP-activated protein kinase, aquaporin-2, etc. Our data provide a rich source for the identification and functional characterization of new genes and candidate polymorphic SSR markers in K. kurodana. The availability of transcriptome information ( http://bioinfo.sch.ac.kr/submission/ ) would promote the utilization of the resources for phylogenetics study and genetic diversity assessment.

Transcriptomic Analysis of the Endangered Neritid Species Clithon retropictus: De Novo Assembly, Functional Annotation, and Marker Discovery

An aquatic gastropod belonging to the family Neritidae, Clithon retropictus is listed as an endangered class II species in South Korea. The lack of information on its genomic background limits the ability to obtain functional data resources and inhibits informed conservation planning for this species. In the present study, the transcriptomic sequencing and de novo assembly of C. retropictus generated a total of 241,696,750 high-quality reads. These assembled to 282,838 unigenes with mean and N50 lengths of 736.9 and 1201 base pairs, respectively. Of these, 125,616 unigenes were subjected to annotation analysis with known proteins in Protostome DB, COG, GO, and KEGG protein databases (BLASTX; E ≤ 0.00001) and with known nucleotides in the Unigene database (BLASTN; E ≤ 0.00001). The GO analysis indicated that cellular process, cell, and catalytic activity are the predominant GO terms in the biological process, cellular component, and molecular function categories, respectively. In addition, 2093 unigenes were distributed in 107 different KEGG pathways. Furthermore, 49,280 simple sequence repeats were identified in the unigenes (>1 kilobase sequences). This is the first report on the identification of transcriptomic and microsatellite resources for C. retropictus, which opens up the possibility of exploring traits related to the adaptation and acclimatization of this species.

Transcriptome of the Australian Mollusc Dicathais orbita Provides Insights into the Biosynthesis of Indoles and Choline Esters

Dicathais orbita is a mollusc of the Muricidae family and is well known for the production of the expensive dye Tyrian purple and its brominated precursors that have anticancer properties, in addition to choline esters with muscle-relaxing properties. However, the biosynthetic pathways that produce these secondary metabolites in D. orbita are not known. Illumina HiSeq 2000 transcriptome sequencing of hypobranchial glands, prostate glands, albumen glands, capsule glands, and mantle and foot tissues of D. orbita generated over 201 million high quality reads that were de novo assembled into 219,437 contigs. Annotation with reference to the Nr, Swiss-Prot and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases identified candidate-coding regions in 76,152 of these contigs, with transcripts for many enzymes in various metabolic pathways associated with secondary metabolite biosynthesis represented. This study revealed that D. orbita expresses a number of genes associated with indole, sulfur and histidine metabolism pathways that are relevant to Tyrian purple precursor biosynthesis, and many of which were not found in the fully annotated genomes of three other molluscs in the KEGG database. However, there were no matches to known bromoperoxidase enzymes within the D. orbita transcripts. These transcriptome data provide a significant molecular resource for gastropod research in general and Tyrian purple producing Muricidae in particular.

Transcriptome Analysis of the Sydney Rock Oyster, Saccostrea glomerata: Insights into Molluscan Immunity

BACKGROUND

Oysters have important ecological functions in their natural environment, acting as global carbon sinks and improving water quality by removing excess nutrients from the water column. During their life-time oysters are exposed to a variety of pathogens that can cause severe mortality in a range of oyster species. Environmental stressors encountered in their habitat can increase the susceptibility of oysters to these pathogens and in general have been shown to impact on oyster immunity, making immune parameters expressed in these marine animals an important research topic.

RESULTS

Paired-end Illumina high throughput sequencing of six S. glomerata tissues exposed to different environmental stressors resulted in a total of 484,121,702 paired-end reads. When reads and assembled transcripts were compared to the C. gigas genome, an overall low level of similarity at the nucleotide level, but a relatively high similarity at the protein level was observed. Examination of the tissue expression pattern showed that some transcripts coding for cathepsins, heat shock proteins and antioxidant proteins were exclusively expressed in the haemolymph of S. glomerata, suggesting a role in innate immunity. Furthermore, analysis of the S. glomerata ORFs showed a wide range of genes potentially involved in innate immunity, from pattern recognition receptors, components of the Toll-like signalling and apoptosis pathways to a complex antioxidant defence mechanism.

CONCLUSIONS

This is the first large scale RNA-Seq study carried out in S. glomerata, showing the complex network of innate immune components that exist in this species. The results confirmed that many of the innate immune system components observed in mammals are also conserved in oysters; however, some, such as the TLR adaptors MAL, TRIF and TRAM are either missing or have been modified significantly. The components identified in this study could help explain the oysters' natural resilience against pathogenic microorganisms encountered in their natural environment.

De novo Transcriptome Generation and Annotation for Two Korean Endemic Land Snails, Aegista chejuensis and Aegista quelpartensis, Using Illumina Paired-End Sequencing Technology

Aegista chejuensis and Aegista quelpartensis (Family-Bradybaenidae) are endemic to Korea, and are considered vulnerable due to declines in their population. The limited genetic resources for these species restricts the ability to prioritize conservation efforts. We sequenced the transcriptomes of these species using Illumina paired-end technology. Approximately 257 and 240 million reads were obtained and assembled into 198,531 and 230,497 unigenes for A. chejuensis and A. quelpartensis, respectively. The average and N50 unigene lengths were 735.4 and 1073 bp, respectively, for A. chejuensis, and 705.6 and 1001 bp, respectively, for A. quelpartensis. In total, 68,484 (34.5%) and 77,745 (33.73%) unigenes for A. chejuensis and A. quelpartensis, respectively, were annotated to databases. Gene Ontology terms were assigned to 23,778 (11.98%) and 26,396 (11.45) unigenes, for A. chejuensis and A. quelpartensis, respectively, while 5050 and 5838 unigenes were mapped to 117 and 124 pathways in the Kyoto Encyclopedia of Genes and Genomes database. In addition, we identified and annotated 9542 and 10,395 putative simple sequence repeats (SSRs) in unigenes from A. chejuensis and A. quelpartensis, respectively. We designed a list of PCR primers flanking the putative SSR regions. These microsatellites may be utilized for future phylogenetics and conservation initiatives.

Sequencing, De Novo Assembly, and Annotation of the Transcriptome of the Endangered Freshwater Pearl Bivalve, Cristaria plicata, Provides Novel Insights into Functional Genes and Marker Discovery

BACKGROUND

The freshwater mussel Cristaria plicata (Bivalvia: Eulamellibranchia: Unionidae), is an economically important species in molluscan aquaculture due to its use in pearl farming. The species have been listed as endangered in South Korea due to the loss of natural habitats caused by anthropogenic activities. The decreasing population and a lack of genomic information on the species is concerning for environmentalists and conservationists. In this study, we conducted a de novo transcriptome sequencing and annotation analysis of C. plicata using Illumina HiSeq 2500 next-generation sequencing (NGS) technology, the Trinity assembler, and bioinformatics databases to prepare a sustainable resource for the identification of candidate genes involved in immunity, defense, and reproduction.

RESULTS

The C. plicata transcriptome analysis included a total of 286,152,584 raw reads and 281,322,837 clean reads. The de novo assembly identified a total of 453,931 contigs and 374,794 non-redundant unigenes with average lengths of 731.2 and 737.1 bp, respectively. Furthermore, 100% coverage of C. plicata mitochondrial genes within two unigenes supported the quality of the assembler. In total, 84,274 unigenes showed homology to entries in at least one database, and 23,246 unigenes were allocated to one or more Gene Ontology (GO) terms. The most prominent GO biological process, cellular component, and molecular function categories (level 2) were cellular process, membrane, and binding, respectively. A total of 4,776 unigenes were mapped to 123 biological pathways in the KEGG database. Based on the GO terms and KEGG annotation, the unigenes were suggested to be involved in immunity, stress responses, sex-determination, and reproduction. A total of 17,251 cDNA simple sequence repeats (cSSRs) were identified from 61,141 unigenes (size of >1 kb) with the most abundant being dinucleotide repeats.

CONCLUSIONS

This dataset represents the first transcriptome analysis of the endangered mollusc, C. plicata. The transcriptome provides a comprehensive sequence resource for the conservation of genetic information in this species and enrichment of the genetic database. The development of molecular markers will assist in the genetic improvement of C. plicata.

Analysis of synonymous codon usage patterns in sixty-four different bivalve species

Synonymous codon usage bias (CUB) is a defined as the non-random usage of codons encoding the same amino acid across different genomes. This phenomenon is common to all organisms and the real weight of the many factors involved in its shaping still remains to be fully determined. So far, relatively little attention has been put in the analysis of CUB in bivalve mollusks due to the limited genomic data available. Taking advantage of the massive sequence data generated from next generation sequencing projects, we explored codon preferences in 64 different species pertaining to the six major evolutionary lineages in Bivalvia. We detected remarkable differences across species, which are only partially dependent on phylogeny. While the intensity of CUB is mild in most organisms, a heterogeneous group of species (including Arcida and Mytilida, among the others) display higher bias and a strong preference for AT-ending codons. We show that the relative strength and direction of mutational bias, selection for translational efficiency and for translational accuracy contribute to the establishment of synonymous codon usage in bivalves. Although many aspects underlying bivalve CUB still remain obscure, we provide for the first time an overview of this phenomenon in this large, commercially and environmentally important, class of marine invertebrates.

Transcriptome analysis and characterisation of gill-expressed carbonic anhydrase and other key osmoregulatory genes in freshwater crayfish Cherax quadricarinatus

The pH and salinity balance mechanisms of crayfish are controlled by a set of transport-related genes. We identified a set of the genes from the gill transcriptome from a freshwater crayfish Cherax quadricarinatus using the Illumina NGS-sequencing technology. We identified and characterized carbonic anhydrase (CA) genes and some other key genes involved in systematic acid-base balance and osmotic/ionic regulation. We also examined expression patterns of some of these genes across different sublethal pH levels [1]. A total of 72,382,710 paired-end Illumina reads were assembled into 36,128 contigs with an average length of 800 bp. About 37% of the contigs received significant BLAST hits and 22% were assigned gene ontology terms. These data will assist in further physiological-genomic studies in crayfish.

Comparative Analysis and Distribution of Omega-3 lcPUFA Biosynthesis Genes in Marine Molluscs

Recent research has identified marine molluscs as an excellent source of omega-3 long-chain polyunsaturated fatty acids (lcPUFAs), based on their potential for endogenous synthesis of lcPUFAs. In this study we generated a representative list of fatty acyl desaturase (Fad) and elongation of very long-chain fatty acid (Elovl) genes from major orders of Phylum Mollusca, through the interrogation of transcriptome and genome sequences, and various publicly available databases. We have identified novel and uncharacterised Fad and Elovl sequences in the following species: Anadara trapezia, Nerita albicilla, Nerita melanotragus, Crassostrea gigas, Lottia gigantea, Aplysia californica, Loligo pealeii and Chlamys farreri. Based on alignments of translated protein sequences of Fad and Elovl genes, the haeme binding motif and histidine boxes of Fad proteins, and the histidine box and seventeen important amino acids in Elovl proteins, were highly conserved. Phylogenetic analysis of aligned reference sequences was used to reconstruct the evolutionary relationships for Fad and Elovl genes separately. Multiple, well resolved clades for both the Fad and Elovl sequences were observed, suggesting that repeated rounds of gene duplication best explain the distribution of Fad and Elovl proteins across the major orders of molluscs. For Elovl sequences, one clade contained the functionally characterised Elovl5 proteins, while another clade contained proteins hypothesised to have Elovl4 function. Additional well resolved clades consisted only of uncharacterised Elovl sequences. One clade from the Fad phylogeny contained only uncharacterised proteins, while the other clade contained functionally characterised delta-5 desaturase proteins. The discovery of an uncharacterised Fad clade is particularly interesting as these divergent proteins may have novel functions. Overall, this paper presents a number of novel Fad and Elovl genes suggesting that many mollusc groups possess most of the required enzymes for the synthesis of lcPUFAs.

Analysis, characterisation and expression of gill-expressed carbonic anhydrase genes in the freshwater crayfish Cherax quadricarinatus

Changes in water quality parameters such as pH and salinity can have a significant effect on productivity of aquaculture species. Similarly, relative osmotic pressure influences various physiological processes and regulates expression of a number of osmoregulatory genes. Among those, carbonic anhydrase (CA) plays a key role in systemic acid-base balance and ion regulation. Redclaw crayfish (Cherax quadricarinatus) are unique in their ability to thrive in environments with naturally varied pH levels, suggesting unique adaptation to pH stress. To date, however, no studies have focused on identification and characterisation of CA or other osmoregulatory genes in C. quadricarinatus. Here, we analysed the redclaw gill transcriptome and characterized CA genes along with a number of other key osmoregulatory genes that were identified in the transcriptome. We also examined patterns of gene expression of these CA genes when exposed to three pH treatments. In total, 72,382,710 paired end Illumina reads were assembled into 36,128 contigs with an average length of 800bp. Approximately 37% of contigs received significant BLAST hits and 22% were assigned gene ontology terms. Three full length CA isoforms; cytoplasmic CA (ChqCAc), glycosyl-phosphatidylinositol-linked CA (ChqCAg), and β-CA (ChqCA-beta) as well as two partial CA gene sequences were identified. Both partial CA genes showed high similarity to ChqCAg and appeared to be duplicated from the ChqCAg. Full length coding sequences of Na(+)/K(+)-ATPase, V-type H(+)-ATPase, sarcoplasmic Ca(+)-ATPase, arginine kinase, calreticulin and Cl(-) channel protein 2 were also identified. Only the ChqCAc gene showed significant differences in expression across the three pH treatments. These data provide valuable information on the gill expressed CA genes and their expression patterns in freshwater crayfish. Overall our data suggest an important role for the ChqCAc gene in response to changes in pH and in systemic acid-base balance in freshwater crayfish.