Next Article in Journal
New Cyclotetrapeptides and a New Diketopiperzine Derivative from the Marine Sponge-Associated Fungus Neosartorya glabra KUFA 0702
Next Article in Special Issue
The Role of Spongia sp. in the Discovery of Marine Lead Compounds
Previous Article in Journal
New Polyketides and New Benzoic Acid Derivatives from the Marine Sponge-Associated Fungus Neosartorya quadricincta KUFA 0081
Previous Article in Special Issue
Guanidine Alkaloids from the Marine Sponge Monanchora pulchra Show Cytotoxic Properties and Prevent EGF-Induced Neoplastic Transformation in Vitro
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

by
Abdul Baten
1,†,
Ajit Kumar Ngangbam
2,†,
Daniel L. E. Waters
1 and
Kirsten Benkendorff
2,*
1
Southern Cross Plant Science, Southern Cross University, Lismore NSW 2480, Australia
2
Marine Ecology Research Centre, School of Environment, Science and Engineering, Southern Cross University, Lismore NSW 2480, Australia
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Mar. Drugs 2016, 14(7), 135; https://doi.org/10.3390/md14070135
Submission received: 16 May 2016 / Revised: 8 July 2016 / Accepted: 12 July 2016 / Published: 20 July 2016
(This article belongs to the Collection Bioactive Compounds from Marine Invertebrates)

Abstract

:
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.

Graphical Abstract

1. Introduction

Dicathais orbita is well known for the production of the dye Tyrian purple, which is a historically important colourant exclusively produced by the Muricidae family of marine molluscs [1,2]. Tyrian purple is not produced within the mollusc but is formed from oxidative and photolytic reactions from a precursor tyrindoxyl sulfate (Figure 1A), which is stored as a salt of the choline ester murexine (Figure 1B) [3,4]. Once the salt is liberated by an aryl sulfatase enzyme, a range of intermediate brominated indole precursors are produced, which have anticancer and antimicrobial properties [5,6,7,8,9]. The major dye component of Tyrian Purple, 6,6′-dibromoindigo, was the first marine natural product to be structurally elucidated [10]; however, a century later, limited information is available on the biosynthesis or gene regulation of this secondary metabolite.
The “post-genomics era” has seen an expansion in the application of bioinformatics to the fields of transcriptomics, proteomics and metabolomics [11]. Transcriptomics provides information on the genes expressed by an organism under certain circumstances or stages of development, in a particular tissue or cell type [12,13]. Mollusc transcriptome studies have been instrumental in establishing the gene expression events associated with shell formation [14], host parasite interactions [15,16], nervous system function [17,18,19], immune defence [20,21], developmental processes [22,23] and cellular and physiological mechanisms [24,25]. Only a few transcriptome profiling studies have been undertaken on predatory marine neogastropods including the gonadal tissues of Reishia (Thais) clavigera [23], the mantle, foot, gills and gonadal tissues of Concholepas concholepas [26], the alimentary canal and salivary glands of Colubraia reticulata [27], the venom glands of Conus consors [28] and the venom ducts of C. tribblei, C. lenavati [29] and Lophiotoma olangoensis, a Turrid snail [30].
Transcriptomics can identify genes involved in the biosynthesis of secondary metabolites [31,32]. The Australasian mollusc D. orbita is an ideal model species for gene expression studies of biosynthetic pathways that may be involved in the synthesis of biologically active secondary metabolites [3]. A preliminary transcriptome study of D. orbita used suppressive subtractive hybridisation to identify genes that were upregulated in the hypobranchial gland, the biosynthetic organ where Tyrian purple is produced [33]. This study confirmed that the hypobranchial gland is a significant site for protein synthesis and regulation, but the only enzyme associated with Tyrian purple production identified was arylsulphatase [33]. However, the study was limited by short read lengths and a low total number of reads. Therefore, the aim of this study was to generate a more comprehensive transcriptome from the hypobranchial glands, prostate glands, capsule glands, albumen glands, and mantle and foot tissues of D. orbita (Muricidae Neogastropoda) and then search these transcriptomes for potential metabolic pathways that could contribute to indole and choline ester biosynthesis using the Kyoto Encyclopedia of Genes and Genomes (KEGG) PATHWAYS database.

2. Results and Discussion

2.1. De Novo Transcriptome Assembly

Transcripome sequencing across the six different tissue types (hypobranchial glands, prostate glands, capsule glands, albumen glands, mantle and foot tissues) resulted in approximately 221 million sequencing reads (Table 1). Raw sequencing reads in FASTQ format were first checked for quality using FASTQC [34] followed by removal of adapter sequences, poly-N stretches and low quality (Phred score < 20) reads using the BBDuck module of the BBMap software package (version 34_90, http://sourceforge.net/projects/bbmap), which resulted in 201 million high quality reads. Table 1 shows the number of raw and quality controlled sequencing reads for all the tissues.
High quality reads were de novo assembled into 219,437 contigs using CLC Genomics server (version 4.9, CLC Bio, Aarhus, Denmark) (Table 2). Transdecoder (version 2.0.3, http://transdecoder.github.io/) identified 76,152 contigs that contained candidate-coding regions that were used for annotation and further downstream analysis.

2.2. Transcriptome Annotation

Basic Local Alignment Search Tool (BLAST) analysis was performed using 76,152 contigs with Open Reading Frames (ORF)s against non-redundant protein database National Center for Biotechnology Information (NCBI) Nr (Mollusc specific proteins), Swiss-Prot and KEGG protein databases. Overall 28,364 contigs (~37%) had significant BLAST hits (e value 11 × 10−5). The D. orbita contig BLAST hit rate was in a similar range to other whole mollusc genome/transcriptome studies with hits ranging from 25% to 40% [35]. A total of 24,996 contigs were assigned to various cellular components, molecular function and biological process gene ontology (GO) categories, as shown in Figure 2. General cell and cell parts were the most frequent subcategories of the cellular components ontology category, while binding and cellular process was the most represented subcategory of molecular function and biological process. Analysis of KEGG pathways showed the largest number of contigs were involved in metabolic pathways (914 contigs), followed by biosynthesis of secondary metabolites (304 contigs) and microbial metabolism in diverse environment (173 contigs) (Supplementary Table S1). Previous studies have revealed a diversity of bacterial symbionts in D. orbita tissues [36,37]. However, we checked the overrepresented k-mers generated in the quality control phase of RNAseq reads and confirmed that these are mostly mitochondria RNA rather than bacterial (<0.1%).

2.3. Tryptophan Metabolism and Phenylalanine, Tyrosine, Tryptophan Biosynthetic Pathways

Specific searches in the KEGG PATHWAY database [38] were undertaken to identify genes potentially involved in the biosynthesis of Tyrian purple precursors and choline esters. Indoles like tyrindoxyl sulfate are thought to be derived from tryptophan metabolism, and we identified 28 enzyme contigs mapped to 35 reactions in the tryptophan metabolism pathway (Figure 3; note that there are multiple KEGG enzyme (EC) numbers for some enzymes, and they can occur at different positions in the pathways, generating more matches to reactions than the number of matching contigs; this applies to all other pathways below). The list of 28 mapped contigs with the KEGG orthology assignment is provided in Supplementary Table S2.
The annotated genomes of only three other molluscs are available for comparison in the KEGG PATHWAYS database, the gastropod Lottia gigantea, cephalopod Octopus bimaculoides and the bivalve Crassostrea gigas. Nearly all the genes involved in tryptophan metabolism found in D. orbita (Figure 3) were identical to those found in the other three molluscs (Figure S1). However, an important point of difference is that unlike these other molluscs, the D. orbita transcriptome contained tryptophanase (4.1.99.1; Figure 3), which converts tryptophan to indole. Consistent with the other molluscs such as L. gigantea and C. gigas, D. orbita expresses aromatic-l-amino-acid decarboxylase (4.1.1.28; Figure 3) that converts tryptophan to tryptamine. However, we did not detect a transcript for tryptophan 5-monooxygenase (1.14.16.4), which converts tryptophan to 5-hydroxy-l-tryptophan, despite the presence of aromatic-l-amino-acid decarboxylase (4.1.1.28) involved in the production of serotonin. As the other three molluscs all contain matches to tryptophan 5-monoxygenase in their genomes (Figure S1), it is possible that this gene is also present in D. orbita but was not detected in our transcriptome due to low expression. Nevertheless, it appears likely that D. orbita diverts the conversion of tryptophan away from 5-hydroxy-l-tryptophan in favour of indoles, to facilitate Tyrian purple precursor production by higher expression of the tryptophanase gene.
In the D. orbita transcriptome, we found no match in the reaction pathway for the conversion of indole to indoxyl, a precursor to indoxyl sulphate (Figure 3, 1.14.16). However, in bacteria cytochrome P450 enzymes [39,40], and/or mono- or dioxygenases [41,42], are involved in the formation of indoxyl sulfate and indigo. We found matches to two cytochrome P450 enzymes, as well as a monoxidase and several dioxygenases in the D. orbita transcriptome (Supplementary Table S2). Our recent studies have also revealed numerous Vibrio sp. that have the ability to synthesize indoles in the Tyrian purple producing tissues of D. orbita [36,37] and these may provide a supplementary source of indoles for Tyrian purple production.
In the phenylalanine, tyrosine and tryptophan biosynthetic pathway (Figure 4), there was a match in the D. orbita transcriptome to tryptophan synthase alpha chain (4.2.1.20; Figure 4), which converts tryptophan to indoles and vice versa, the initial precursors for Tyrian purple biosynthesis [3,43,44]. This tryptophan synthase was only found in the D. orbita transcriptome and not found in the annotated genomes of L. gigantea, O. bimaculoides and C. gigas. Most of the genes involved in phenylalanine, tyrosine and tryptophan biosynthesis in D. orbita were found to be different to the other three molluscs and this Neogastropoda had more matches to these aromatic amino acid biosynthesis reactions (14 genes) compared to L. gigantea (five genes), O. bimaculoides (five genes) and C. gigas (five genes) (Figure S2).

2.4. Sulfur, Cysteine and Methionine Metabolisms Pathway in Dicathais orbita

Tyrindoxyl sulphate contains a methane thiol group at the 2′ position of the indole ring (Figure 1A). Examination of the D. orbita transcripts with reference to the sulfur metabolism pathway revealed that most of the genes involved in D. orbita sulfur metabolism were found to be similar to the other three molluscs. Enzymes like dimethyl-sulfide monooxygenase (1.14.13.131; Figure 5) that produces methyl mercaptan from dimethyl disulphide were not detected in any of the molluscs, including D. orbita. Overall, the number of genes detected in the D. orbita sulfur metabolism pathway (15 genes, Figure 5) was only slightly higher than the number identified in the genome of L. gigantea (12 genes), O. bimaculoides (12 genes) and C. gigas (14 genes) (Figure S3).
The metabolism of sulfur containing amino acids provides another possible source of the methanethiol group in tyrindoxyl sulfate. Enzyme matches in the cysteine and methionine metabolism pathways (Figure 6) indicate that D. orbita has the ability to produce 3-methylthioproprionate, 3-mercaptopyruvate and thiosulfate, but no match was found to methionine-gamma-lyase (4.4.1.11; Figure 6), which converts the amino acid l-Methionine directly into methanethiol. Nevertheless, we did detect a transcript for tyrosine aminotransferase (2.6.1.5; Figure 6), which may play a role in placing the methane thiol onto the aromatic indole ring. Tyrosine aminotransferase genes were also found in the Lottia, octopus and oyster genomes. Overall, the D. orbita transcriptome had more matches to enzymes in the cysteine and methionine metabolism pathway (41 genes) compared to L. gigantea (32 genes), O. bimaculoides (30 genes) and C. gigas (33 genes) (Figure S4). This ability to metabolise sulfur from various sources is consistent with the biosynthesis of indole mercaptans in Muricidae.

2.5. Bromoperoxidase Enzymes

Tyridoxyl sulfate is a 6-brominated indole derivative (Figure 1A), and bromoperoxidase activity has been detected in the hypobranchial glands of D. orbita [45] and other Muricidae species [46]. Consequently, a search was undertaken for bromoperoxidase genes by aligning D. orbita transcripts against known bromoperoxidase genes using BLAST (e value 1 × 10−3). However, no evidence of matches to bromoperoxidase genes was found using these sequences. This is consistent with a previous study that examined the transcripts that were up-regulated in the D. orbita hypobranchial glands relative to other D. orbita tissues, which also found no matches to bromoperoxidase genes [33]. There appears to be no previous reports of bromoperoxidase coding genes or transcripts in any gastropods or other molluscs. However, recent metagenomic analyses of D. orbita hypobranchial glands revealed the presence of bacterial symbionts known to produce bromoperoxidase enzymes [36]. Our D. orbita transcriptome data therefore supports the suggestion that symbiotic bacteria play a role in the biosynthesis of Tyrian purple.

2.6. Dicathais Orbita Glycerophospholipid and Histidine Metabolism Pathway

Tyrindoxyl sulfate is stored as a choline ester of murexine, which contains a choline ester group and imidazole moiety (Figure 1B). The glycerophospholipid metabolism pathway has a role in the biosynthesis of choline esters [47] and more matches to reactions in the glycerophospholipid metabolism pathway (53 genes) were found in the D. orbita transcriptome compared to the annotated genomes of other molluscs (L. gigantea = 48, O. bimaculoides = 40 and C. gigas = 44 genes) (Figure S5). The D. orbita transcriptome included matches to choline/ethanol amine kinase (2.7.1.32; Figure 7), which produces choline from phosphocholine, as well as choline O-acetyltransferase (2.3.1.6; Figure 7) and acetylcholinesterase (3.1.1.7; Figure 7), which produce acetylcholine. Laffy et al. [33] found that acetylcholinesterase (3.1.1.7) transcripts were upregulated in the hypobranchial gland of D. orbita relative to foot tissue. The octopus O. bimaculoides and gastropod L. gigantea genomes were also found to contain matches to these enzymes for acetylcholine biosynthesis, whereas the oyster C. gigas lacks phosphocholine phosphatase. In comparison to the gastropods and the bivalve, the octopus O. bimaculoides was found to lack any matches to genes in the phosphothanolamine N-methyltransferase pathway for the production of phosphocholine from phosphoethanolamine.
Choline or acetyl choline combines with imidazole to produce the muscle relaxant murexine [48] (Figure 1B). Imidazole is a derivative of the amino acid histidine [49], and there were several matches to the histidine metabolism pathway in the D. orbita transcriptome (Figure 8). These include, diamine oxidase (1.4.3.22; Figure 8), aldehyde dehydrogenase (NAD+) (1.2.1.3; Figure 8) and monoamine oxidase (1.4.3.4; Figure 8), which convert histamine into imidazole. There was also a match to histidine ammonia-lyase (4.3.1.3; Figure 8), which converts l-histidine to urocanate, which could feasibly combine with choline ester to form murexine. All these biosynthetic enzymes are also found in the L. gigantea, O. bimaculoides and C. gigas genomes.
Previous studies of choline esters in molluscs have focused on the predatory neogastropods, and there is no record of murexine or similar derivatives being isolated from limpets or oysters. Roseghini et al. [48] found no evidence for imidazole or acryl choline esters in 27 species from eight families of herbivorous and scavenging gastropods, including three Patellidae limpets, while at least one of these compounds was found in 53 of 55 species of the predatory Muricoidae superfamily. This implies the Neogastropoda have evolved a specific murexine biosynthesis pathway and, consistent with this, the D. orbita transcriptome had more matches to enzymes in the histidine metabolism pathway (19 genes) when compared to L. gigantea (10 genes), O. bimaculoides (11 genes) and C. gigas (13 genes) (Figure S6). Specifically, the enzyme involved in imidazole biosynthesis imidazoleglycerol-phosphate dehydratase (4.2.1.19; Figure 8) was only found in the D. orbita transcriptome. Overall, it appears the neogastropod D. orbita has evolved a complex suite of metabolic capabilities that are not represented in the more primitive orthogastropod or bivalve, for which complete genome sequences are available.

3. Materials and Methods

3.1. Specimen Collection

Eighteen adult specimens of D. orbita (Table 1) were collected during low tide from the sub-tidal and intertidal rocky reefs of Flat rock, Ballina (28°84′ S and 153°60′ E), NSW, Australia. Six spawning females were collected during the breeding season August 2014 and a further six females and six males were collected after breeding season in January 2015, under the permit number F89/1171-6.0 issued by the Department of Primary Industries, NSW Government, Australia. Total RNA was extracted from the hypobranchial glands of the females collected in August 2014 and from five different tissues from female and three tissues from male D. orbita (Figure 9) collected in January 2015 (Table 1).

3.2. Transcriptome Sequencing

Three independent replicate snails were used for each tissue sample summarised in Table 1. The tissues were stabilized prior to RNA extraction in RNase free 2 ml Eppendorf tube using 600 μL of RNAlater RNA stabilization reagent (Qiagen, Chadstone, Victoria, Australia). The stabilized tissue was incubated at 4 °C overnight and stored at −80 °C, prior to extracting the total RNA. The total RNA was extracted from the RNAlater stabilized tissue using the E.Z.N.A. Mollusc RNA Kit (Omega Bio-tek, Norcross, GA, USA.) following the manufacturer’s instructions. The concentration and quality (purity and integrity) of total RNA was assessed by NanoDrop and the Agilent Bioanalyzer 2100 System (Agilent Technologies, Santa Clara, CA, USA). The total RNA extracted from three biological replicates of each tissue type (hypobranchial gland, prostate gland, albumen gland, capsule gland, mantle and foot) was pooled within the same tube for each tissue in equal masses. The pooled extracted RNA was stored at −80 °C until further used. The RNA samples were shipped to Macrogen Inc. (Seoul, Korea) for high throughput sequencing. Prior to shipping, each RNA sample was precipitated in a mixture of 2× ethanol (96%) and 0.1× sodium acetate (3 M). mRNA isolation and library construction were performed by Macrogen. The libraries were sequenced using the Illumina HiSeq 2000 platform (HCS2.2.38 version, Illumina, Seoul, Korea).

3.3. De Novo Transcriptome Assembly and Annotation

FASTQ format raw sequencing reads were checked for quality using FASTQC (version 0.10.4, http://www.bioinformatics.babraham.ac.uk/projects/fastqc) [34]. The adapter sequences, poly-N stretches and low quality reads (Phred score < 20) were removed using the BBDuck module of the BBMap software package (version 34_90, http://sourceforge.net/projects/bbmap) using command “bbduk.sh in1 = in1.fastq in2 = in2.fastq out1 = out1.fastq out2 = out2.fastq ref = adapters.fa qtrim = rl trimq = 20 ktrim = rl k = 28 hdist = 1 minlength = 20”. BLAST search of over-represented k-mers against NCBI confirmed these were mostly mitochondrial and ribosomal Dicathais obita genes, and no evidence of bacterial contamination was found. CLC Genomics Workbench, version 4.9 (CLC Bio, Aarhus, Denmark; www.clcbio.com) with the option to map reads back to contigs, automatic word size and automatic bubble size was used to de novo assemble the high quality reads. All the contigs were clustered using CD-hit-est (version v4.6.1, http://weizhongli-lab.org/cd-hit) [50], and Transdecoder (version 2.0.3, http://transdecoder.github.io/) was used to identify candidate coding regions within transcript sequences.
BLAST analysis was done against non-redundant protein database Nr (Mollusc specific proteins), Swiss-Prot and KEGG protein databases. Gene ontology analysis was performed with Interproscan (version 5.10.50, https://code.google.com/archive/p/interproscan) [51] using command “interproscan.sh -appl ProDom,PfamA,PANTHER -i longest_orfs.pep.fa -o out.txt -f TSV -goterms -iprlookup -pa” and plotted using WEGO (http://wego.genomics.org.cn/cgi-bin/wego/index.p) [52]. D. orbita transcripts were searched for bromoperoxidase gene by aligning against known bromoperoxidase genes available in the NCBI GenBank using BLAST (e value 1 × 10−3) (Supplementary Table S3).

3.4. Nucleotide Sequence Accession Number

All raw sequence data were deposited in the European nucleotide archive (ENA) with the accession numbers PRJEB12262. Assembled contigs are available from the authors upon request.

4. Conclusions

This transcriptome study of D. orbita generated over 216 million high quality reads that were de novo assembled into 219,437 contigs, of which 76,152 contigs contained candidate-coding regions that were annotated with Nr, Swiss-Prot and KEGG databases. This provides a significant new molecular resource for neogastropod molluscs, and adds to pool of genomic data for molluscs in general. Several genes that are potentially associated with Tyrian purple precursor biosynthesis in D. orbita were identified. It appears the neogastropod D. orbita has evolved a complex suite of metabolic capabilities that are not represented in the more primitive orthogastropods or bivalves, for which complete genome sequences are available.

Supplementary Materials

The following are available online at www.mdpi.com/1660-3397/14/7/135/s1, Figure S1: Tryptophan metabolism pathways for (A) Crassostrea gigas; (B) Lottia gigantean; and (C) Octopus bimaculoides showing enzyme matches in green including tryptophan 5-monoxygenase (red box), which was missing form Dicathais orbita, but no match to tryptophanase (blue box), Figure S2: Phenylalanine, tyrosine and tryptophan biosynthetic pathways for (A) Crassostrea gigas; (B) Lottia gigantean; and (C) Octopus bimaculoides with enzyme matches in green, but with no match to with tryptophan synthase highlighted in the blue box, Figure S3: Sulfur metabolism pathways for (A) Crassostrea gigas; (B) Lottia gigantea; and (C) Octopus bimaculoides, Figure S4: Cysteine and methionine metabolism pathway of (A) Crassostrea gigas; (B) Lottia gigantea; and (C) Octopus bimaculoides, Figure S5: Glycerophospholipid metabolism pathway of (A) Crassostrea gigas; (B) Lottia gigantea; and (C) Octopus bimaculoides showing enzyme matches in green with those relevant to choline ester synthesis highlighted in red (present) and blue (absent), Figure S6: Histidine metabolism pathway of (A) Crassostrea gigas; (B) Lottia gigantea; and (C) Octopus bimaculoides showing matching enzymes in green, including several enzymes that convert histidine into imidazole (red boxes) but not imidazoleglycerol-phosphate dehydratase (blue box), Table S1: Analysis of KEGG pathway showing the top 20 metabolic pathway involving the largest number of contigs in D. orbita trancriptome, Table S2: List of the 28 mapped contigs and with the KEGG orthology assignment in D. orbita trancriptome for tryptophan metabolism, Table S3: List of known bromoperoxidsae genes available in NCBI GenBank used for BLAST against D. orbita.

Acknowledgments

This work was supported by the Southern Cross University International Postgraduate Research Scholarship and School of Environment, Science and Engineering Postgraduate support for A.K.N. and a philanthropic grant to K.B.

Author Contributions

K.B. initiated the study and all authors contributed to the experimental design. A.K.N. prepared the samples for sequencing and A.B. assembled and annotated the transcriptome. All authors contributed to the data analysis and interpretation. A.B., A.K.N. and K.B. wrote the paper with editorial input from D.L.E.W.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
ENAEuropean Nucleotide Archive
NCBINational Center for Biotechnology Information
BLASTBasic Local Alignment Search Tool
ORFsOpen Reading Frames
KEGGKyoto Encyclopedia of Genes and Genomes

References

  1. Westley, C.; Benkendorff, K. Sex-specific Tyrian purple genesis: Precursor and pigment distribution in the reproductive system of the marine mollusc, Dicathais orbita. J. Chem. Ecol. 2008, 34, 44–56. [Google Scholar] [CrossRef] [PubMed]
  2. Benkendorff, K.; Rudd, D.; Nongmaithem, B.D.; Liu, L.; Young, F.; Edwards, V.; Avila, C.; Abbott, C.A. Are the traditional medical uses of Muricidae molluscs substantiated by their pharmacological properties and bioactive compounds? Mar. Drugs 2015, 13, 5237–5275. [Google Scholar] [CrossRef] [PubMed]
  3. Benkendorff, K. Natural product research in the Australian marine invertebrate Dicathais orbita. Mar. Drugs 2013, 11, 1370–1398. [Google Scholar] [CrossRef] [PubMed]
  4. Baker, J.T.; Duke, C.C. Isolation of choline and choline ester salts of tyrindoxyl sulphate from the marine molluscs Dicathais orbita and Mancinella keineri. Tetrahedron Lett. 1976, 15, 1233–1234. [Google Scholar] [CrossRef]
  5. Esmaeelian, B.; Abbott, C.A.; le Leu, R.K.; Benkendorff, K. 6-bromoisatin found in muricid mollusc extracts inhibits colon cancer cell proliferation and induces apoptosis, preventing early stage tumor formation in a colorectal cancer rodent model. Mar. Drugs 2014, 12, 17–35. [Google Scholar] [CrossRef] [PubMed]
  6. Benkendorff, K.; Bremner, J.B.; Davis, A.R. Tyrian purple precursors in the egg masses of the Australian muricid, Dicathais orbita: A possible defensive role. J. Chem. Ecol. 2000, 26, 1037–1050. [Google Scholar] [CrossRef]
  7. Esmaeelian, B.; Benkendorff, K.; Johnston, M.R.; Abbott, C.A. Purified brominated indole derivatives from Dicathais orbita induce apoptosis and cell cycle arrest in colorectal cancer cell lines. Mar. Drugs 2013, 11, 3802–3822. [Google Scholar] [CrossRef] [PubMed]
  8. Edwards, V.; Benkendorff, K.; Young, F. Marine compounds selectively induce apoptosis in female reproductive cancer cells but not in primary-derived human reproductive granulosa cells. Mar. Drugs 2012, 10, 64–83. [Google Scholar] [CrossRef] [PubMed]
  9. Westley, C.B.; McIver, C.M.; Abbott, C.A.; le Leu, R.K.; Benkendorff, K. Enhanced acute apoptotic response to azoxymethane-induced DNA damage in the rodent colonic epithelium by Tyrian purple precursors: A potential colorectal cancer chemopreventative. Cancer Biol. Ther. 2010, 9, 371–379. [Google Scholar] [CrossRef] [PubMed]
  10. Freidlander, P. Ueber den farbstoff des antiken purpura aus Murex brandaris. Chem. Ber. 1909, 42, 765–770. [Google Scholar] [CrossRef]
  11. Zhang, W.; Li, F.; Nie, L. Integrating multiple “omics” analysis for microbial biology: Application and methodologies. Microbiology 2010, 156, 287–301. [Google Scholar] [CrossRef] [PubMed]
  12. Brown, T.A. Chapter 3: Genomes. In Transcriptomes and Proteomes, 2nd ed.; John Wiley & Sons: New York, NY, USA, 2002; pp. 69–91. [Google Scholar]
  13. Adams, J. Transcriptome: Connecting the genome to gene function. Nat. Educ. 2008, 1, 195. Available online: http://www.nature.com/scitable/topicpage/transcriptome-connecting-the-genome-to-gene-function-605 (accessed on 14 May 2016). [Google Scholar]
  14. Jackson, D.J.; McDougall, C.; Green, K.; Simpson, F.; Worheide, G.; Degnan, B.M. A rapidly evolving secretome builds and patterns a sea shell. BMC Biol. 2006, 4. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Lockyer, A.E.; Spinks, J.; Kane, R.A.; Hoffmann, K.F.; Fitzpatrick, J.M.; Rollinson, D.; Noble, L.R.; Jones, C.S. Biomphalaria glabrata transcriptome: cDNA microarray profiling identifies resistant- and susceptible-specific gene expression in haemocytes from snail strains exposed to Schistosoma mansoni. BMC Genom. 2008, 9. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Knight, M.; Arican-Goktas, H.D.; Ittiprasert, W.; Odoemelam, E.C.; Miller, A.N.; Bridger, J.M. Schistosomes and snails: A molecular encounter. Front. Genet. 2014, 5, 1–7. [Google Scholar] [CrossRef] [PubMed]
  17. Moroz, L.L.; Edwards, J.R.; Puthanveettil, S.V.; Kohn, A.B.; Ha, T.; Heyland, A.; Knudsen, B.; Sahni, A.; Yu, F.; Liu, L.; et al. Neuronal Transcriptome of Aplysia: Neuronal Compartments and Circuitry. Cell 2006, 127, 1453–1467. [Google Scholar] [CrossRef] [PubMed]
  18. Sadamoto, H.; Takahashi, H.; Okada, T.; Kenmoku, H.; Toyota, M.; Asakawa, Y. De novo sequencing and transcriptome analysis of the central nervous system of mollusc Lymnaea stagnalis by deep RNA sequencing. PLoS ONE 2012, 7. [Google Scholar] [CrossRef] [PubMed]
  19. Feng, Z.P.; Zhang, Z.; van Kesteren, R.E.; Straub, V.A.; van Nierop, P.; Jin, K.; Nejatbakhsh, N.; Goldberg, J.I.; Spencer, G.E.; Yeoman, M.S.; et al. Transcriptome analysis of the central nervous system of the mollusc Lymnaea stagnalis. BMC Genom. 2009, 10. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Romero, A.; Novoa, B.; Figueras, A. Genomics, immune studies and diseases in bivalve aquaculture. J. Invertebr. Pathol. 2012, 9, 110–121. [Google Scholar]
  21. Chavez, M.J.; Valenzuela, M.V.; Nunez, A.G.; Maldonado, A.W.; Gallardo, E.C. Concholepas concholepas Ferritin H-like subunit (CcFer): Molecular characterization and single nucleotide polymorphism associated to innate immune response. Fish Shellfish Immunol. 2013, 35, 910–917. [Google Scholar] [CrossRef] [PubMed]
  22. Powell, D.; Knibb, W.; Remilton, C.; Elizur, A. De-novo transcriptome analysis of the banana shrimp (Fenneropenaeus merguiensis) and identification of genes associated with reproduction and development. Mar. Genom. 2015, 22, 71–78. [Google Scholar] [CrossRef] [PubMed]
  23. Ho, K.K.Y.; Leung, P.T.Y.; Ip, J.C.H.; Qiu, J.W.; Leung, K.M.Y. De novo transcriptomic profile in the gonadal tissues of the intertidal whelk Reishia clavigera. Mar. Pollut. Bull. 2014, 85, 499–504. [Google Scholar] [CrossRef] [PubMed]
  24. Zhang, G.; Fang, X.; Guo, X.; Li, L.; Luo, R.; Xu, F.; Yang, P.; Zhang, L.; Wang, X.; Qi, H.; et al. The oyster genome reveals stress adaptation and complexity of shell formation. Nature 2012, 490, 49–54. [Google Scholar] [CrossRef] [PubMed]
  25. Meng, J.; Zhu, Q.; Zhang, L.; Li, C.; Li, L.; She, Z.; Huang, B.; Zhang, G. Genome and transcriptome analyses provide insight into the euryhaline adaptation mechanism of Crassostrea gigas. PLoS ONE 2013, 8. [Google Scholar] [CrossRef] [PubMed]
  26. Cardenas, L.; Sanchez, R.; Gomez, D.; Fuenzalida, G.; Gallardo-Escarate, C.; Tanguy, A. Transcriptome analysis in Concholepas concholepas (Gastropoda, Muricidae): Mining and characterization of new genomic and molecular markers. Mar. Genom. 2011, 4, 197–205. [Google Scholar] [CrossRef] [PubMed]
  27. Modica, M.V.; Lombardo, F.; Franchini, P.; Oliverio, M. The venomous cocktail of the vampire snail Colubraria reticulata (Mollusca, Gastropoda). BMC Genom. 2015, 16, 441. [Google Scholar] [CrossRef] [PubMed]
  28. Terrat, Y.; Biass, D.; Dutertre, S.; Favreau, P.; Remm, M.; Stocklin, R.; Piquemal, D.; Ducancel, F. High-resolution picture of a venom gland transcriptome: Case study with the marine snail Conus consors. Toxicon 2012, 59, 34–46. [Google Scholar] [CrossRef] [PubMed]
  29. Barghi, N.; Concepcion, G.P.; Olivera, B.M.; Lluisma, A.O. Comparison of the venom peptides and their expression in closely related Conus species: Insights into adaptive post-speciation evolution of Conus exogenomes. Genome Biol. Evol. 2015, 7, 1797–1814. [Google Scholar] [CrossRef] [PubMed]
  30. Watkins, M.; Hillyard, D.R.; Olivera, B.M. Genes expressed in a turrid venom duct: Divergence and similarity to conotoxins. J. Mol. Evol. 2006, 62, 247–256. [Google Scholar] [CrossRef] [PubMed]
  31. Urbarova, I.; Karlsen, B.O.; Okkenhaug, S.; Seternes, O.M.; Johansen, S.D.; Emblem, A. Digital marine bioprospecting: Mining new neurotoxin drug candidates from the transcriptomes of cold-water sea anemones. Mar. Drugs 2012, 10, 2265–2279. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Zhang, Y.; Zhang, S.F.; Lin, L.; Wang, D.Z. Comparative transcriptome analysis of a toxin-producing dinoflagellate Alexandrium catenella and its non-toxic mutant. Mar. Drugs 2014, 12, 5698–5718. [Google Scholar] [CrossRef] [PubMed]
  33. Laffy, P.W.; Benkendorff, K.; Abbott, C.A. Suppressive subtractive hybridisation transcriptomics provides a novel insight into the functional role of the hypobranchial gland in a marine mollusc. Comp. Biochem. Physiol. D Genom. Proteom. 2013, 8, 111–122. [Google Scholar] [CrossRef] [PubMed]
  34. Andrews, S. FastQC: A quality control tool for high throughput sequence data. Available online: http://www.bioinformatics.babraham.ac.uk/projects/fastqc/ (accessed on 14 May 2016).
  35. Prentis, P.J.; Pavasovic, A. The Anadara trapezia transcriptome: A resource for molluscan physiological genomics. Mar. Genom. 2014, 18, 113–115. [Google Scholar] [CrossRef] [PubMed]
  36. Ngangbam, A.K.; Baten, A.; Waters, D.L.E.; Whalan, S.; Benkendorff, K. Characterization of bacterial communities associated with the Tyrian purple producing gland in a marine gastropod. PLoS ONE 2015, 10. [Google Scholar] [CrossRef] [PubMed]
  37. Ngangbam, A.K.; Waters, D.L.E.; Whalan, S.; Baten, A.; Benkendorff, K. Indole producing bacteria from the biosynthetic organs of Muricid mollusc could contribute to Tyrian purple production. J. Shellfish Res. 2015, 34, 443–454. [Google Scholar] [CrossRef]
  38. Kanehisa, M.; Goto, S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic. Acids Res. 2000, 28, 27–30. [Google Scholar] [CrossRef] [PubMed]
  39. Banoglu, E.; Jha, G.G.; King, R.S. Hepatic microsomal metabolism of indole to indoxyl, a precursor of indoxyl sulfate. Eur. J. Drug. Metab. Pharmacokinet. 2001, 26, 235–240. [Google Scholar] [CrossRef] [PubMed]
  40. Arora, P.K.; Sharma, A.; Bae, H. Microbial degradation of indole and its derivatives. J. Chem. 2015, 2015, 1–13. [Google Scholar] [CrossRef]
  41. O’Connor, K.E.; Hartmans, S. Indigo formation by aromatic hydrocarbon-degrading bacteria. Biotechnol. Lett. 1998, 20, 219–223. [Google Scholar] [CrossRef]
  42. Han, X.; Wang, W.; Xiao, X. Microbial biosynthesis and biotransformation of indigo and indigo-like pigments. Chin. J. Biotechnol. 2008, 24, 921–926. [Google Scholar] [CrossRef]
  43. Westley, C.; Vine, K.; Benkendorff, K. A proposed functional role for indole derivatives in reproduction and defense of the Muricidae (Neogastropoda: Mollusca). In Indirubin, the Red Shade of Indigo; Meijer, L., Guyard, N., Skaltsounis, L., Eisenbrand, G., Eds.; Life in Progress: Roscoff, France, 2006; pp. 31–44. [Google Scholar]
  44. Verhecken, A. The indole pigments of Mollusca. Ann. Soc. R. Zool. Belg. 1989, 119, 181–197. [Google Scholar]
  45. Westley, C.; Benkendorff, K. The distribution of precursors and biosynthetic enzymes required for Tyrian purple genesis in the hypobranchial gland, gonoduct, an egg masses of Dicathais orbita (Gmelin, 1791) (Neogastropoda: Muricidae). Nautilus 2009, 123, 148–153. [Google Scholar]
  46. Jannun, R.; Coe, E.L. Bromoperoxidase from the marine snail, Murex-trunculus. Comp. Biochem. Physiol. B Comp. Biochem. 1987, 88, 917–922. [Google Scholar] [CrossRef]
  47. Quastel, J.H.; Tennenbaum, M.; Wheatley, A.H. Choline ester formation in, and choline esterase activities of, tissues in vitro. Biochem. J. 1936, 30, 1668–1681. [Google Scholar] [CrossRef] [PubMed]
  48. Roseghini, M.; Severini, C.; Erspamer, G.F.; Erspamer, V. Choline esters and biogenic amines in the hypobranchial gland of 55 molluscan species of the neogastropod Muricoidea Superfamily. Toxicon 1996, 34, 33–55. [Google Scholar] [CrossRef]
  49. Sundberg, R.J.; Martin, R.B. Interactions of histidine and other imidazole derivatives with transition metal ions in chemical and biological systems. Chem. Rev. 1974, 74, 471–517. [Google Scholar] [CrossRef]
  50. Li, W.; Godzik, A. Cd-hit: A fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics 2006, 22, 1658–1659. [Google Scholar] [CrossRef] [PubMed]
  51. Zdobnov, E.M.; Apweiler, R. InterProScan—An integration platform for the signature-recognition methods in InterPro. Bioinformatics 2001, 17, 847–848. [Google Scholar] [CrossRef] [PubMed]
  52. Ye, J.; Fang, L.; Zheng, H.; Zhang, Y.; Chen, J.; Zhang, Z.; Wang, J.; Li, S.; Li, R.; Bolund, L.; et al. WEGO: A web tool for plotting GO annotations. Nucleic Acids Res. 2006, 34, W293–W297. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Tyrindoxyl sulfate (A), the ultimate Tyrian purple precursor in Dictahais orbita, is held as a salt of the choline ester murexine (B).
Figure 1. Tyrindoxyl sulfate (A), the ultimate Tyrian purple precursor in Dictahais orbita, is held as a salt of the choline ester murexine (B).
Marinedrugs 14 00135 g001
Figure 2. The proportion and number of Dicathais orbita contigs assigned to gene ontology (GO) terms from biological process, cellular component and molecular function. Biological process was the most highly represented GO category followed by cellular component and molecular function.
Figure 2. The proportion and number of Dicathais orbita contigs assigned to gene ontology (GO) terms from biological process, cellular component and molecular function. Biological process was the most highly represented GO category followed by cellular component and molecular function.
Marinedrugs 14 00135 g002
Figure 3. Tryptophan metabolism pathway with matches to Dicathais orbita contigs filled in green. The match to a tryptophanase relevant to indole biosynthesis is highlighted by the red box, whereas the tryptophan 5-monoxygenase that was not detected in our transcriptome is highlighted in a blue box.
Figure 3. Tryptophan metabolism pathway with matches to Dicathais orbita contigs filled in green. The match to a tryptophanase relevant to indole biosynthesis is highlighted by the red box, whereas the tryptophan 5-monoxygenase that was not detected in our transcriptome is highlighted in a blue box.
Marinedrugs 14 00135 g003
Figure 4. Phenylalanine, tyrosine and tryptophan biosynthetic pathways showing matches to Dicathais orbita contigs highlighted in green, with tryptophan synthase highlighted in the red box.
Figure 4. Phenylalanine, tyrosine and tryptophan biosynthetic pathways showing matches to Dicathais orbita contigs highlighted in green, with tryptophan synthase highlighted in the red box.
Marinedrugs 14 00135 g004
Figure 5. Sulfur metabolism pathway with matches to Dicathais orbita contigs highlighted in green; there was no match to dimethyl-sulfide monooxygenase in our transcriptome (blue box).
Figure 5. Sulfur metabolism pathway with matches to Dicathais orbita contigs highlighted in green; there was no match to dimethyl-sulfide monooxygenase in our transcriptome (blue box).
Marinedrugs 14 00135 g005
Figure 6. Cysteine and methionine metabolism pathways showing matches to Dicathais orbita contigs highlighted in green, including tyrosine aminotransferase (red box), but no match was found to methionine-gamma-lyase (blue box).
Figure 6. Cysteine and methionine metabolism pathways showing matches to Dicathais orbita contigs highlighted in green, including tyrosine aminotransferase (red box), but no match was found to methionine-gamma-lyase (blue box).
Marinedrugs 14 00135 g006
Figure 7. Glycerophospholipid metabolism pathway with matches to D. orbita contigs highlighted in green including choline kinase, choline O-acetyltransferase and acetylcholinesterase (red boxes) used to generate the acetyl choline moiety found in murexine.
Figure 7. Glycerophospholipid metabolism pathway with matches to D. orbita contigs highlighted in green including choline kinase, choline O-acetyltransferase and acetylcholinesterase (red boxes) used to generate the acetyl choline moiety found in murexine.
Marinedrugs 14 00135 g007
Figure 8. Histidine metabolism pathway showing matches to Dicathais orbita contigs highlighted in green, including several enzymes that convert histidine into imidazole (red boxes) and imidazoleglycerol-phosphate dehydratase (red box top pathway).
Figure 8. Histidine metabolism pathway showing matches to Dicathais orbita contigs highlighted in green, including several enzymes that convert histidine into imidazole (red boxes) and imidazoleglycerol-phosphate dehydratase (red box top pathway).
Marinedrugs 14 00135 g008
Figure 9. Dicathais orbita male (A) and female (B) tissues used for RNA extraction to generate the transcriptome.
Figure 9. Dicathais orbita male (A) and female (B) tissues used for RNA extraction to generate the transcriptome.
Marinedrugs 14 00135 g009
Table 1. Summary of the number of raw sequencing reads and the percent remaining after quality control from 14 tissue samples of Dicathais orbita.
Table 1. Summary of the number of raw sequencing reads and the percent remaining after quality control from 14 tissue samples of Dicathais orbita.
Snails (S)DescriptionRaw ReadsHigh Quality Reads
NumberPercent
S1 + S2 + S3Female hypobranchial gland 1, August, Breeding season, 201415,531,32215,100,46697.23
S4 + S5 + S6Female hypobranchial gland 2, August, Breeding season, 201415,693,38515,258,67197.23
S7 + S8 + S9Female hypobranchial gland 1, January, 201515,835,27115,425,53397.41
S10 + S11 + S12Female hypobranchial gland 2, January, 201516,457,63515,990,72497.16
S13 + S14 + S15Male hypobranchial gland 1, January, 201516,142,31715,684,92697.17
S16 + S17 + S18Male hypobranchial gland 2, January, 201517,461,00716,997,49797.35
S7 + S8 + S9Female foot 1, January, 201516,015,53515,595,46397.38
S10 + S11 + S12Female foot 2, January, 201517,057,43316,653,22291.40
S13 + S14 + S15Male foot 1, January, 201514,241,69013,885,32797.50
S16 + S17 + S18Male foot 2, January, 201515,813,36315,406,03097.42
S7 + S8 + S9Capsule gland, January, 201515,805,86715,291,49896.75
S7 + S8 + S9Albumen gland, January, 201514,442,86414,011,09997.01
S13 + S14 + S15Prostate gland, January, 201515,600,68815,113,84296.88
S10 + S11 + S12Mantle 1, January, 201516,273,55615,804,24797.12
-Total222,371,933216,218,545-
Table 2. Summary statistics of the assembled contigs using CLC Genomics de novo assembler.
Table 2. Summary statistics of the assembled contigs using CLC Genomics de novo assembler.
Contig Summary Statisticsbp (Base Pair)
Number of contigs219,437
Total assembly length117,767,308
N50608
Mean contig length537
Largest contig length12,897
Number of contigs larger than 500 bp59,144
Number of contigs larger than 1000 bp22,818

Share and Cite

MDPI and ACS Style

Baten, A.; Ngangbam, A.K.; Waters, D.L.E.; Benkendorff, K. Transcriptome of the Australian Mollusc Dicathais orbita Provides Insights into the Biosynthesis of Indoles and Choline Esters. Mar. Drugs 2016, 14, 135. https://doi.org/10.3390/md14070135

AMA Style

Baten A, Ngangbam AK, Waters DLE, Benkendorff K. Transcriptome of the Australian Mollusc Dicathais orbita Provides Insights into the Biosynthesis of Indoles and Choline Esters. Marine Drugs. 2016; 14(7):135. https://doi.org/10.3390/md14070135

Chicago/Turabian Style

Baten, Abdul, Ajit Kumar Ngangbam, Daniel L. E. Waters, and Kirsten Benkendorff. 2016. "Transcriptome of the Australian Mollusc Dicathais orbita Provides Insights into the Biosynthesis of Indoles and Choline Esters" Marine Drugs 14, no. 7: 135. https://doi.org/10.3390/md14070135

APA Style

Baten, A., Ngangbam, A. K., Waters, D. L. E., & Benkendorff, K. (2016). Transcriptome of the Australian Mollusc Dicathais orbita Provides Insights into the Biosynthesis of Indoles and Choline Esters. Marine Drugs, 14(7), 135. https://doi.org/10.3390/md14070135

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop