Skip to main content

A de novo transcriptome of the noble scallop, Chlamys nobilis, focusing on mining transcripts for carotenoid-based coloration

Abstract

Background

The noble scallop Chlamys nobilis Reeve displays polymorphism in shell and muscle colors. Previous research showed that the orange scallops with orange shell and muscle had a significantly higher carotenoid content than the brown ones with brown shell and white muscle. There is currently a need to identify candidate genes associated with carotenoid-based coloration.

Results

In the present study, 454 GS-FLX sequencing of noble scallop transcriptome yielded 1,181,060 clean sequence reads, which were assembled into 49,717 isotigs, leaving 110,158 reads as the singletons. Of the 159,875 unique sequences, 11.84% isotigs and 9.35% singletons were annotated. Moreover, 3,844 SSRs and over 120,000 high confidence variants (SNPs and INDELs) were identified. Especially, one class B scavenge receptor termed SRB-like-3 was discovered to express only in orange scallops and absent in brown ones, suggesting a significant association with high carotenoid content. Down-regulation of SRB-like-3 mRNA by RNA interference remarkably decreased blood carotenoid, providing compelling evidence that SRB-like-3 is an ideal candidate gene controlling carotenoid deposition and determining orange coloration.

Conclusion

Transcriptome analysis of noble scallop reveals a novel scavenger receptor significantly associated with orange scallop rich in carotenoid content. Our findings pave the way for further functional elucidation of this gene and molecular basis of carotenoid deposition in orange scallop.

Background

Carotenoids are bright yellow and red pigments that are responsible for some coloration found in animals [1]. Carotenoids also play important physiological roles such as acting as antioxidants in the immune system [2,3]. Unlike other pigments types such as melanins, carotenoids cannot be synthesized by animals and must be acquired through diet [4]. There are a number of factors (such as food source, seasonal change) that potentially limit the ability of animals to deposit carotenoids in their body tissues [5,6]. Although carotenoid traits have often been shown to be condition-dependent, carotenoid coloration and accumulation is also dependent on underlying genetic mechanisms. Animals preferentially deposit certain carotenoids over others, and are able to enzymatically convert and cleave dietary carotenoids into other derived forms [7], implying strongly the involvement of genes encoding appropriate carotenoid-binding and transport proteins or enzymes participating carotenoid metabolism.

The SRB (scavenger receptor class B) is first identified as playing a role in the uptake of lutein [8], carotene [9], zeaxanthin and xanthophylls [10], and lycopene [11]. A SRB homologue, ninaD, is essential for cellular uptake of carotenoids in Drosophila and a mutation in this gene results in carotenoid-free and thus a vitamin A deficient phenotype [12]. Two recently cloned genes, Cameo2 and SCRB15 of CD36 (Cluster Determinant 36), which are homologous to SRB, have been shown to be involved in the selective transport of lutein and β-carotene, respectively, into the silk gland of Bombyx mori [13,14]. StAR (steroidogenic acute regulatory)/MLN64 (metastatic lymph node 64) are members of the StAR domain family that are involved in the intracellular transport of cholesterol for the initiation of steroidogenesis [15]. StAR isolated in the macula of primate retina could selectively bind lutein with high affinity [16]. The B. mori carotenoid-binding protein (CBP) is an orthologue of vertebrate MLN64, and is involved in the transport of lutein [17]. The BCMO (β,β-carotene-15,15’-monooxygenase)/BCDO(β,β-carotene-9,9’-oxygenase) is involved in the enzymatic cleavage of carotenoids [18]. Loss-of-function mutation in BCMO results in hypercarotenemia [19]. Carotenoids are cleaved to form colorless apo-carotenoid derivatives in chickens with white skin, while yellow-skinned chickens presumably have one or more cis-acting regulatory mutations in BCDO, resulting in a yellow coloration in the skin because of deposition of uncleaved carotenoids [20]. Other genes involved in the transport and binding of carotenoids are Niemann Pick C1-like 1 (NPC1L1) [21], ATP-binding cassette sub-family G member 5 (ABCG5) [22], Glutathione S-transferase Pi1 (GST) [23] and crustacyanin [24]. Lastly, intestinal transcription factor (ISX) [25] and retinoic acid receptor (RAR)/retinoid X receptor (RXR) [26] are important transcription factors that regulate the expression of genes (such as BCMO and SRB) involved in carotenoid deposition.

The noble scallop Chlamys nobilis Reeve, an important aquaculture bivalve in China, displays conspicuous polymorphism in shell color (such as orange, orange-purple, brown, etc.) and difference in muscle color (such as orange and white). The orange scallops have carotenoid-based orange mantle and adductor muscle due to high presence of carotenoids. Our previous work showed that the orange scallops with orange shell and muscle had a significantly higher carotenoid content than the brown ones with brown shell and white muscle [27]. By establishing different scallop lines, both shell color and muscle color have been confirmed to be control by at least two loci, with one locus showing dominance epistasis to the other [28,29]. Therefore, the carotenoid-based orange coloration in muscle is likely due to differential expression of one or a few genes at the site of carotenoid deposition.

In recent years, transcriptome analysis has been widely recognized as a very useful tool to identify candidate genes underlying molecular mechanisms. In the present study, we first sequenced and assembled the transcriptome of noble scallop C. nobilis using a GS-FLX 454 platform. Second, we quantified the expression of genes that are homologous to known carotenoid candidate genes in the adductor muscle, which actively deposits carotenoids in scallop, and investigated whether differential expression was associated with carotenoid content variation in orange scallops versus brown scallops. Our goals were to 1) generate a transcriptome database useful for functional genetic studies of C. nobilis; and 2) identify candidate transcripts involved in carotenoid-based coloration or carotenoid deposition.

Methods

Ethics statement

The scallops used in this study were taken from Nan’ao Marine Biology Station of Shantou University, located at Nan’ao island of Shantou, Guangdong, China. No specific permits were required for the described field studies, as the sampling locations were not privately owned or protected in any way. These field studies also did not include endangered or protected species. The animals were processed according to “the Regulations for the Administration of Affairs Concerning Experimental Animals” established by the Guangdong Provincial Department of Science and Technology on the Use and Care of Animals.

Sample collection and preparation

In noble scallop, orange color was dominant to brown color. Color segregation occurred when crossing two orange scallops [28,29]. Both orange and brown scallops used in the present study were from a line of F2 generation produced by continuous crossing orange parental scallops (Figure 1A). A total of 20 orange scallops (rich in carotenoids) and 20 brown scallops (lack of carotenoids) at 14-month old were randomly chosen. Average shell size in length, height and width for the orange was 67.28 ± 4.22 mm, 72.92 ± 3.73 mm and 24.06 ± 1.66 mm, and for the brown was 66.80 ± 3.09 mm, 71.69 ± 3.73 mm and 23.88 ± 1.15 mm, respectively. Tissues including the gonad, mantle, gill and adductor muscle (Figure 1C) were sampled and homogenized with the QiaShredder (Qiagen, Germany) for total RNA extraction using the Qiagen RNeasy (Qiagen, Germany) kit. mRNA was then purified by using the Qiagen Oligotex mRNA purification kit. Equal amounts of mRNA from four tissues were pooled for either an “orange” or a “brown” scallop sample. From the two pooled samples, about 600 ng mRNA was used for cDNA generation with the SMART cDNA synthesis kit (Clontech Laboratories, USA). Quality control in each extraction step was investigated using gel electrophoresis and nanodrop spectrophotometry (Peqlab, Germany). Both orange and brown scallop cDNAs were further checked with a Bioanalyzer 2100 (Agilent Technologies, USA). As a result, two cDNA libraries (one for orange scallop and the other for brown scallop) with an average length of 400 bp were generated according to the manufacturers’ protocol and sequenced on a 454 Genome Sequencer system (Roche Life Sciences, USA) with FLX and Titanium chemistry.

Figure 1
figure1

Line of the C. nobilis was used in transcriptome sequencing. (A) Both parents are orange scallop. (B) Offspring with orange and brown coloration segregation obtained by crossing two orange scallops. (C) Four kinds of tissues: gonad (a), mantle (b), gill (c), and adductor muscle (d). Male scallop has a white or lighter orange color gonad, while female scallop has a heavier orange color gonad.

Sequence assembly and functional annotation

All sequence reads taken directly from the 454 GS-FLX sequencer were run through the sff file program (Newbler v2.6, Roche) to remove sequencing adapters A and B. Barcodes were removed by Seqclean (Lastest86_64) program and poor sequence data were further cleaned by Lucy v1.20 program (–m 50 –e 0.03 0.03 –w 30 0.03 10 0.1 –b 4 0.03). Sequences with homopolymers of a single nucleotide occupying 60% of the read and those less than 50 nucleotides in length were discarded. Trimmed sequences from orange or brown scallop were mixed and then assembled de novo using the default parameters of Newbler v2.6 (Roche). All C. nobilis EST (expressed sequence tags) sequences were submitted to NCBI Sequence Read Archive under Accession No. SRX253988. ESTs that did not form isotigs (singletons) and isotigs resulting from the assembly of multiple sequences were referred to as unique sequences. These unique sequences were translated into six reading frames and used as a query to search the public databases including Non-redundant protein database (Nr) and Swiss-Prot database (Swiss-Prot). All unique sequences were sequentially compared using BlastX (cut-off E-value of 1e-5) with the sequences in two public protein databases (Nr and Swiss-Prot). Once a sequence had a blast hit in one of the databases, a description was built from the description of that hit. Additionally, Gene Ontology (GO) terms were deduced from the blast results using Blast2GO, and sorted into the immediate subcategories for ‘molecular function’, ‘cellular component’ and ‘biological process’.

Identification of EST-SSR motifs and EST-SNPs

All EST sequences were searched for SSR motifs using the MISA (MIcroSAtellite identification) program (http://pgrc.ipk-gatersleben.de/misa/). Default settings were employed to detect perfect di-, tri-, tetra-, penta-, and hexa-nucleotide motifs (including compound motifs). To be assigned, di-nucleotide SSRs (Simple Sequence Repeats) required a minimum of 6 repeats, and all other SSR types needed a minimum of 5 repeats. Two neighboring SSRs with the maximum interruption no more than 100 nucleotides were considered as a compound SSR.

Multiple nucleotide sequence alignments of isotigs identified among the EST libraries were undertaken to identify putative SNPs. Since few reference sequences were available, SNPs were identified as superimposed nucleotide peaks where 2 or more reads contained polymorphisms at the variant allele. SNPs were identified using default parameters in gsMapper v2.3 (Roche) to align isotigs from two color datasets. In addition, only an overall transition vs transversion (Ts/Tv) ratio was calculated across the dataset. Perl script modules linked to the primer modeling software Primer3 were used to design PCR primers flanking for each unique SNP region identified.

Data mining of transcripts with putative function involved in carotenoid Deposition

From public databases, we compiled a dataset of the 15 known gene involved in carotenoid deposition were collected (Table 1). The amino-acid sequences of the known carotenoid deposition genes, covering carotenoid absorption, transport and cleavage, were used to search (tBlastn) for homologues in 454-derived sequences. Those sequences with scores more than or equal to 100 and E values less than or equal to 1e-10 were clustered to develop unigenes, and all of the unigenes were considered as candidate transcripts. The resulting unigenes were in turn used to search the GenBank databases by BlastX to confirm their putative carotenoid-related functions.

Table 1 Carotenoid-related candidate gene

mRNA expression of selected candidate transcripts in orange and brown scallop

Expression of selected transcripts was investigated in adductor muscle from 6 orange scallop or 6 brown scallops at 14-month old, and two technical replicates were performed for each scallop. All scallops used in this experiment were from a F2 generation as described above, and cultured in the same cage. Total RNA was extracted and quality and quantity determined using a nanodrop spectrophotometer. 1 μg mRNA was used to synthesize cDNA by PrimeScript RT reagent kit with gDNA Eraser (TaKaRa). Quantitative real-time RT-PCR was conducted in a LightCycler®480 System using the SYBR Premix Ex Taq II qRT-PCR Kit (TaKaRa). Each assay was performed with β-actin mRNA as the internal control. The real-time PCR program was 95°C for 30s, followed by 40 cycles of 95°C for 5 s, and 60°C for 30s according to the instructions of the manufacturer. Dissociation analysis of amplification products was performed at the end of each PCR reaction to confirm that only one PCR product was amplified and detected. The comparative CT method (2-ΔΔCT method) was used to analyze the expression level of each candidate genes. All data were given in terms of relative mRNA expressed as means ± SE. The data were subjected to analysis of one-way ANOVA, and p-values smaller than 0.05 were considered statistically significant.

Detecting presence of SRB (scavenger receptor class B)-3-like and measurement of total carotenoid content in scallops

Four scallop lines derived from orange parents, which have color segregation of orange and brown, were chosen to performed this experiment. In total, 80 scallops (40 orange and 40 brown), derived from 4 lines produced by crossing two orange scallops in the Spring of 2012, were used to detect the presence of SRB-3-like in the blood and determine total carotenoid content in the adductor muscle. Presence of SRB-3-like in the blood was detected using primers S3F1: CGATTTTGGAACGGTAACAGTAACTTGGA and S3R1: ATGGATTGACTGATGTGAGATGT. PCR amplification product was confirmed by sequencing. Total carotenoid content in the adductor muscle was determined using the method of Zheng et al. [27].

dsRNA synthesis

SRB-like-3 gene was amplified through PCR with noble scallop cDNA as template and 1Fi and 1Ri as primers (Table 2). The PCR products were separated, purified, ligated with vector pMD-18 T (Takara), and transformed into DH5α E. coli cell. The plasmid was extracted using MiniBEST Plasmid Purification Kit Ver.4.0 (Takara) according to the manufacturer’s protocol.

Table 2 Primer sequences used for dsRNA synthesis

For dsRNA synthesis, SRB-like-3 was amplified by PCR with the primers 2Fi and 2Ri (containing T7 promoter) using the recombinant plasmid pMD-18 T-SRB as the template and (Table 2). Similarly, for dsRNA synthesis of EGFP gene [30], plasmid pEGFP-N1 was used as the template for PCR using EGFPF and EGFPR as the primers. Quantity and quality of the DNA fragments were assessed by nanodrop spectrophotometry and electrophoresis in 1.0% agarose gel. dsRNA was synthesized in vitro using MEGAscript RNAi Kit (Life Technology) following the manufacturer’s protocol. After being incubated at 75°C for 5 min, dsRNA was cooled to room temperature, digested with DNase and RNase, and purified.

RNAi (RNA interference) assay

Forty orange scallops were used, and each of them was injected with 40 μg dsRNA of SRB-like-3 or EGFP gene (as a control) into the adductor muscle. Scallops were labeled and placed in a cage. The blank group was injected with Rnase-free water. Five individuals were sampled at 3, 6, 12, and 24 h for each group. Adductor muscle muscle, blood and intestine were subjected to total RNA extraction. Real-Time PCR was performed as described above with 2 technological replicates for each sample.

Effect of dsRNA on carotenoid deposition in the blood and adductor muscle

Orange scallops were randomly chosen, and 20 of them were injected with dsRNA of SRB-like-3 or EGFP gene, 5 of them were injected with RNA-free water as the blank group. 24 h later, they were injected again. 12 h after the second injection, 5 scallops from Rnase-free water group, 10 scallops from dsSRB-like-3 group, and 10 scallops from dsEGFP group were sampled. 1 ml blood from each scallop was freeze-dried and added with 0.5 ml acetone to extract caroteoid for about 2-4 h at darkness. Caroteoids from adductor muscle were extracted according to method by Zheng et al. [27]. The samples were always under N2 until measurement of absorption at 480 nm to determine their carotenoid content.

Results and discussion

Roche 454 GS-FLX sequencing and isotigs assembly

A total of 1,416,522 raw reads were obtained from the 454 GS-FLX sequencer, including 1,181,060 clean reads (averaging 308 bp in length) after adaptor trimming, size-selection and quality control (Table 3). The size distribution of raw reads and clean reads are shown in Figure 2A and B.

Table 3 Summary statistics for EST and de novo assembly
Figure 2
figure2

Overview of the C. nobilis transcriptome sequencing and assembly. (A) Size distribution of 454 raw reads. (B) Size distribution of 454 reads after removal of adaptor and short sequences. (C) Log-log plot showing the dependence of isotigs lengths on the number of reads assembled into each isotigs. (D) Size distribution of isotigs.

Sequences that passed basic quality standards were clustered and assembled de novo (Newbler v2.60; Roche). Overall, approximately 91% (1,070,902) reads were assembled into 49,717 isotigs, and the others (110,158) remained as singletons (Table 1). Sequencing coverage of isotigs is shown in Figure 2C with an average 7-fold coverage. The size distribution of isotigs is shown in Figure 2D, which ranges from 50 to 7,102 bp with an average of 580 bp. The percentage of reads assembled de novo is similar to that found in other studies [31-33]. The large numbers of unique sequences (singletons and isotigs) in this study are likely due to the extensive diversity in the initial RNA samples as mentioned above. Different organs and sexes, and sequence variants in individuals are known to produce extensive alternatively spliced transcripts, resulting in misalignments and incorrect assembly between reads arising from the same genomic region [34].

Annotation of the transcriptome

All isotigs and singletons were subjected to CD-hit program (version 4.5.6) to remove redundant sequences, leaving 111,670 unique sequences (46,284 isotigs, 65,386 singletons). The annotation for unique sequences from C. nobiliswas based on sequence similarity searches against public databases. These databases included NCBI Nr and SwissProt. About 21.19% of all sequences (11.84% isotigs and 9.35% singletons) identified a homologue mostly in both or at least one of the databases at e-value 1e-5 (Additional file 1: Table S1). Because the significance of sequence similarity depends in part on the length of the query sequence, short unique sequences frequently cannot be matched to known genes [31,35]. The proportion of sequences with matches in public databases was greater for the longer assembled sequences. Namely, 24.12% matches were recorded for sequences ≥ 300 bp, 45.29% for those longer than 1 kb, but only 5.79% for those short sequences (<300 bp) (Table 4). The percentage of sequences with annotation information in this study was considerably low (approximately 17.43%). The poor annotation efficiency may be due to the insufficient sequences in public databases for phylogenetically close species to date [31,36].

Table 4 Summary of annotation of the C. nobilis transcriptome

Gene ontology assignments

Gene Ontology (GO) [37] could provide a dynamic, controlled vocabulary and hierarchical relationships for the representation of information on molecular function, cellular component and biological process, allowing a coherent annotation of gene products. Of annotated unique sequences in Nr or SwissProt database, 10,051 unique sequences were assigned to one or more GO terms (Additional file 2: Table S2). A total of 4,031 GO terms were obtained, with 21,182 unique sequences for biological processes, 24,348 unique sequences for cellular components and 15,122 unique sequences for molecular function (Figure 3). Of those sequences for biological processes, the major ones were cellular process (27.88%) and metabolic process (26.51%). For cellular component, the most represented categories were cell (28.55%) and cell part (28.55%). Regarding molecular functions, binding (44.95%) was the most represented sequences according to GO terms, followed by catalytic activity (33.79%). Similar results were also found in other species such as the Yesso scallop (Patinopecten yessoensis) [31] and the freshwater prawn (Macrobrachium rosenbergii) [33]. These GO annotations can provide a comprehensive information on transcript functions of C. nobilis.

Figure 3
figure3

Functional annotation of assembled sequences based on gene ontology (GO) categorization. GO analysis was performed at the level 2 for three main categories (cellular component, molecular function and biological process).

Putative molecular markers

Selection based on molecular markers is well known to be vital in shellfish aquaculture industry [38]. 454 sequencing may provide an excellent source for mining and development of these markers for C. nobilis, because few genetic markers are currently available. A total of 3,259 unique sequences were found to possess 3,479 SSRs, of which 276 (7.93%) unique sequences annotated in Nr and/or Swissport database were considered as priority candidates for maker development (Table 5, Additional file 3: Table S3). The most frequent repeat motifs were di-nucleotide repeats (68.12%), followed by tri-nucleotide repeats (25.29%) and tetra-nucleotide repeats (5.81%). Penta-nucleotides repeats and Hexa-nucleotides repeats only accounted for 0.78%. Among the di-nucleotide repeat classes, TA was the most frequent dimer motif (47.93%). With regard to tri-nucleotide repeats, TGG (6.48%) was the most common motif, followed by GAT (6.02%) and TGA (5.00%).

Table 5 Summary of simple sequence repeat (SSR) nucleotide classes among different nucleotide types found in C. nobilis sequences

SNPs in C. nobilis EST isotigs were identified using the ssahaSNP program (http://www.sanger.ac.uk). Of the 71,719 SNPs detected, 43,433 were putative transitions (Ts) and 28,286 were putative transversions (Tv), giving a mean Ts: Tv ratio of 1.54: 1.00 across the transcriptome (Figure 4, Additional file 4: Table S4). The Ts: Tv ratio can help to identify genes affected by selection. A total of 96,320 INDELs across the transcriptome were detected. However, much caution must be paid to those INDELs because of technical problems associated with 454 pyrosequencing [39].

Figure 4
figure4

Classification of single nucleotide polymorphisms (SNPs) identified from 454 sequences. The overall frequency of these SNP types in C. nobilis transcriptome is one per 278 bp.

The overall frequency of all types of SNPs in the transcriptome, including INDELs, was one per 278 bp. Of the predicted SNPs, including INDELs, 122,927 (73.15%) were identified from isotigs covered by ten or more reads, suggesting the majority of SNPs identified in this study were covered at sufficient sequencing depth and more likely represent ‘true’ SNPs [40]. Among the SNPs, 53,831 (32.03%) were identified from isotigs with annotation information.

Twenty five of these predicted SNPs were randomly selected for validation using PCR and Sanger sequencing, and 17 of these tests (68%) were successful (Additional file 5: Table S5). The result here confirmed that the majority of computationally predicted SNPs from the 454 transcriptome sequences would benefit us in our future genetic markers development.

Identification of carotenoid-based coloration transcripts from 454 sequences

Full length protein sequences of the 15 known genes responsible for carotenoid absorption, binding or carotenoids cleavage were used to perform tBlastn searches against 454-derived sequences. A total of 44 isotigs and 20 singletons from 454-derived sequences were identified with scores more than or equal to 100 and E values less than or equal to 1e-10, and 48 non-redundant sequences were developed and used to identify their putative functions by BlastX searches against the GenBank databases (Additional file 6: Table S6). After BlastX searches against Nr database, 26 transcripts were considered as tentative C. nobilis carotenoid-related transcripts (Additional file 6: Table S6).

Quantitative real-time reverse transcription PCR was performed for each of the 26 candidate transcripts to determine their relative levels of expression in the adductor muscle taken from orange scallop and brown scallop. However, we failed to find any transcripts with significant difference in expression between the orange and white adductor muscle (Figure 5). We speculated sequence variations might exist in these transcripts.

Figure 5
figure5

Comparison of the expression level of 26 selected tentative carotenoid deposition transcripts in orange and brown scallop adductor muscle. 6 scallops were used in the experiment and each expression analysis was also performed in two independent experiments. Significant difference was performed by one-way ANOVA test (P < 0.05).

Members of three gene families, SRB/CD36, StAR/MLN64, and BCMO/BCDO have been implicated in uptake and deposition of carotenoids in animal tissues, providing plausible candidates for carotenoids accumulation. To find out whether sequence mutations existed in transcripts that belonged to these three families, we cloned 4 SRB-like genes, 2 STAR-like genes and 2 BCMO-like genes based on transcriptome data by RACE PCR (Additional file 7: Table S7), and screened possible mutation sites of CDS. Several missense mutations were found, which, however, showed no correlation with carotenoids accumulation (data not shown).

Nonetheless, a scavenge receptor gene termed SRB-like-3 was identified, which showed a high sequence similarity with SRB-like-2 (Figure 6). PCR was performed using primers (S3F1 and S3R1 shown in Figure 6) for SRB-like-3. Interestingly, SRB-like-3 was only detected in orange scallop and absent in brown scallop, suggesting that SRB-like-3 is significantly associated with high carotenoid content and possibly an important candidate gene in carotenoid deposition (Table 6). Tissue expression profile (Figure 7) showed that SRB-like-3 was highly expressed in the gonad (having highest carotenoid content) and intestine (where carotenoids were absorbed), indicating that SRB-like-3 mainly functions in these two tissues.

Figure 6
figure6

Alignment of SRB-like-2 and SRB-like-3. Primers S3F1 and S3R1 give special PCR amplification of SRB-like-3.

Table 6 Progeny testing of the four lines in scallop C. nobilis
Figure 7
figure7

Tissue expression profile of SRB-like-3 in orange scallop. Different letter means significant difference by one-way ANOVA test (P < 0.05).

mRNA expression of SRB-lile-3 after RNAi

To find out whether SRB-like-3 was involved in carotenoid deposition, dsRNA of SRB-like-3 (dsSRB-like-3) was synthesized and injected into orange scallop adductor muscle. Results of Real-Time PCR showed that the mRNA level of SRB-like-3 in tested tissue was down-regulated by dsSRB-like-3 (Figure 8). The expression level of SRB-like-3 mRNA was significantly suppressed by 67%, 48% and 45% in the adductor muscle, blood and intestine, respectively, when compared to that of the dsEGFP injected group 24 h after injection.

Figure 8
figure8

mRNA expression of SRB-like-3 after injection of dsRNA. Letter indicates comparison of the same tissue. Different letter means significant difference by one-way ANOVA test (P < 0.05).

Carotenoid content in the blood and adductor muscle after RNAi

Total caronoid content in the blood was measured after RNAi. Color density of the dsSRB-like-3 group is lighter than that of dsEGFP group (Figure 9), implying that the dsSRB-like-3 group had a relative low carotenoid content. Carotenoid content in the blood was measured, showing that the dsSRB-like-3 group indeed had a significantly lower carotenoid content when compared to that of the dsEGFP group or blank group (Table 7), while there was no statistical difference between the dsEGFP group and the RNase-free water group. Carotenoid content in the adductor muscle showed no remarkable difference among the dsSRB-like-3, dsEGFP and Rnase-free water group (Table 8). Anyhow, our result providing compelling evidence that SRB-like-3 was a candidate gene that was at least responsible for blood carotenoid content, implying that SRB-like-3 might take part in absorption of carotenoid to blood.

Figure 9
figure9

Carotenoids extration from the blood (dsEGFP: dsRNA of EGFP; dsSRB-like-3: dsRNA of SRB-like-3).

Table 7 Carotenoid content (CC) of blood
Table 8 Carotenoid content (CC) of adductor

Conclusion

Here we documented a large-scale, multi-organ transcriptome for the noble scallop C. nobilis, which has the unique characterization of carotenoid accumulation but few molecular knowledge has been available. Our findings provide a nearly complete description of the expressed genes, which is a substantial contribution to the existing sequence resources for this species. Application of these resources will greatly enhance future genetic and genomic studies on scallop and other mollusks. The description of the expressed genes and their functions was illustrated according to annotation and GO assignment. 3,844 SSRs and over 120,000 high confidence variants (SNPs and INDELs) were identified that can be useful for mapping and QTLs in this scallop and related species. The most important point is that a scavenge receptor termed SRB-like-3 is only expressed in orange scallop but absent in brown scallop, significantly associated with high carotenoid content, suggesting SRB-like-3 is possibly a candidate gene responsible for carotenoid deposition in orange scallop. Results from RNAi study of this gene provides convincing evidence that SRB-like-3 is involved in carotenoid deposition in blood.

References

  1. 1.

    Matsuno T. Aquatic animal carotenoids. Fisheries Sci. 2001;67(5):771–83.

    Article  CAS  Google Scholar 

  2. 2.

    Chew BP, Park JS. Carotenoid action on the immune response. The J Nutr. 2004;134(1):257S–61.

    CAS  Google Scholar 

  3. 3.

    Hughes DA. Dietary carotenoids and human immune function. Nutrition. 2001;17(10):823–7.

    Article  CAS  PubMed  Google Scholar 

  4. 4.

    Walsh N, Dale J, McGraw KJ, Pointer MA, Mundy NI. Candidate genes for carotenoid coloration in vertebrates and their expression profiles in the carotenoid-containing plumage and bill of a wild bird. Proc Biol Sci. 2012;279(1726):58–66.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  5. 5.

    Chimsung N, Lall S, Tantikitti C, Verlhac-Trichet V, Milley J. Effects of dietary cholesterol on astaxanthin transport in plasma of Atlantic salmon (Salmo salar). Comp Biochem Phys B. 2013;165(1):73–81.

    Article  CAS  Google Scholar 

  6. 6.

    Olsen RE, Kiessling A, Milley JE, Ross NW, Lall SP. Effect of lipid source and bile salts in diet of Atlantic salmon, Salmo salar L., on astaxanthin blood levels. Aquaculture. 2005;250(3):804–12.

    Article  CAS  Google Scholar 

  7. 7.

    Maoka T. Recent progress in structural studies of carotenoids in animals and plants. Arch Biochem Biophys. 2009;483(2):191–5.

    Article  CAS  PubMed  Google Scholar 

  8. 8.

    Reboul E, Abou L, Mikail C, Ghiringhelli O, Andre M, Portugal H, et al. Lutein transport by caco-2 TC-7 cells occurs partly by a facilitated process involving the scavenger receptor class B type I (SR-BI). J Biol Chem. 2005;387(Pt 2):455–61.

    CAS  Google Scholar 

  9. 9.

    van Bennekum A, Werder M, Thuahnai ST, Han CH, Duong P, Williams DL, et al. Class B scavenger receptor-mediated intestinal absorption of dietary beta-carotene and cholesterol. Biochemistry-us. 2005;44(11):4517–25.

    Article  Google Scholar 

  10. 10.

    During A, Doraiswamy S, Harrison EH. Xanthophylls are preferentially taken up compared with beta-carotene by retinal cells via a SRBI-dependent mechanism. J Lipid Res. 2008;49(8):1715–24.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  11. 11.

    Moussa M, Landrier JF, Reboul E, Ghiringhelli O, Comera C, Collet X, et al. Lycopene absorption in human intestinal cells and in mice involves scavenger receptor class B type I but not Niemann-Pick C1-like 1. J Nutr. 2008;138(8):1432–6.

    CAS  PubMed  Google Scholar 

  12. 12.

    Kiefer C, Sumser E, Wernet MF, von Lintig J. A class B scavenger receptor mediates the cellular uptake of carotenoids in Drosophila. Proc Natl Acad Sci U S A. 2002;99(16):10581–6.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  13. 13.

    Sakudoh T, Iizuka T, Narukawa J, Sezutsu H, Kobayashi I, Kuwazaki S, et al. A CD36-related transmembrane protein Is coordinated with an intracellular lipid-binding protein in selective carotenoid transport for cocoon coloration. J Biol Chem. 2010;285:7739–51.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  14. 14.

    Sakudoh T, Kuwazaki S, Iizuka T, Narukawa J, Yamamoto K, Uchino K, et al. CD36 homolog divergence is responsible for the selectivity of carotenoid species migration to the silk gland of the silkworm Bombyx mori. J Lipid Res. 2013;54(2):482–95.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  15. 15.

    Watari H, Arakane F, Moog-Lutz C, Kallen CB, Tomasetto C, Gerton GL, et al. MLN64 contains a domain with homology to the steroidogenic acute regulatory protein (StAR) that stimulates steroidogenesis. Proc Natl Acad Sci U S A. 1997;94(16):8462–7.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  16. 16.

    Li B, Vachali P, Frederick JM, Bernstein PS. Identification of StARD3 as a lutein-binding protein in the macula of the primate retina. Biochemistry-us. 2011;50(13):2541–9.

    Article  CAS  Google Scholar 

  17. 17.

    Sakudoh T, Tsuchida K, Kataoka H. BmStart1, a novel carotenoid-binding protein isoform from Bombyx mori, is orthologous to MLN64, a mammalian cholesterol transporter. Biochem Bioph Res Co. 2005;336:1125–35.

    Article  CAS  Google Scholar 

  18. 18.

    Kiefer C, Hessel S, Lampert JM, Vogt K, Lederer MO, Breithaupt DE, et al. Identification and characterization of a mammalian enzyme catalyzing the asymmetric oxidative cleavage of provitamin A. J Biol Chem. 2001;276:14110–6.

    CAS  PubMed  Google Scholar 

  19. 19.

    Lindqvist A, Sharvill J, Sharvill DE, Andersson S. Loss-of-function mutation in carotenoid 15,15'-monooxygenase identified in a patient with hypercarotenemia and hypovitaminosis A. J Nutr. 2007;137(11):2346–50.

    CAS  PubMed  Google Scholar 

  20. 20.

    Eriksson J, Larson G, Gunnarsson U, Bed’hom B, Tixier-Boichard M, Stromstedt L, et al. Identification of the Yellow skin gene reveals a hybrid origin of the domestic chicken. PLoS Genetics. 2008;4(2):e1000010.

    Article  PubMed Central  PubMed  Google Scholar 

  21. 21.

    During A, Dawson HD, Harrison EH. Carotenoid transport is decreased and expression of the lipid transporters SR-BI, NPC1L1, and ABCA1 is downregulated in Caco-2 cells treated with ezetimibe. J Nutr. 2005;135(10):2305–12.

    CAS  PubMed  Google Scholar 

  22. 22.

    Herron KL, McGrane MM, Waters D, Lofgren IE, Clark RM, Ordovas JM, et al. The ABCG5 polymorphism contributes to individual responses to dietary cholesterol and carotenoids in eggs. J Nutr. 2006;136(5):1161–5.

    CAS  PubMed  Google Scholar 

  23. 23.

    Bhosale P, Larson AJ, Frederick JM, Southwick K, Thulin CD, Bernstein PS. Identification and characterization of a Pi isoform of glutathione S-transferase (GSTP1) as a zeaxanthin-binding protein in the macula of the human eye. J Biol Chem. 2004;279(47):49447–54.

    Article  CAS  PubMed  Google Scholar 

  24. 24.

    Wade NM, Tollenaere A, Hall MR, Degnan BM. Evolution of a novel carotenoid-binding protein responsible for crustacean shell color. Mol Biol Evol. 2009;26(8):1851–64.

    Article  CAS  PubMed  Google Scholar 

  25. 25.

    Lobo GP, Hessel S, Eichinger A, Noy N, Moise AR, Wyss A, et al. ISX is a retinoic acid-sensitive gatekeeper that controls intestinal beta, beta-carotene absorption and vitamin A production. Faseb J. 2010;24(6):1656–66.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  26. 26.

    Xiao JH, Durand B, Chambon P, Voorhees JJ. Endogenous retinoic acid receptor (RAR)-retinoid X receptor (RXR) heterodimers are the major functional forms regulating retinoid-responsive elements in adult human keratinocytes. Binding of ligands to RAR only is sufficient for RAR-RXR heterodimers to confer ligand-dependent activation of hRAR beta 2/RARE (DR5). J Biol Chem. 1995;270(7):3001–11.

    Article  CAS  PubMed  Google Scholar 

  27. 27.

    Zheng H, Liu H, Zhang T, Wang S, Sun Z, Liu W, et al. Total carotenoid differences in scallop tissues of Chlamys nobilis (Bivalve: Pectinidae) with regard to gender and shell colour. Food Chem. 2010;122(4):1164–7.

    Article  CAS  Google Scholar 

  28. 28.

    Zheng H, Zhang T, Sun Z, Liu W, Liu H. Inheritance of shell colours in the noble scallop Chlamys nobilis (Bivalve: Pectinidae). Aquac Res. 2012;44(8):1229–35.

    Article  Google Scholar 

  29. 29.

    Liu H, Zheng H, Sun Z, Zhang Q, Wang D. Effect of artificial selection on total carotenoid content in the noble scallop Chlamys nobilis. 18th International Pectinid Workshop 2011. p. 46–47 (Abstract).

  30. 30.

    You Y, Huan P, Liu B. RNAi assay in primary cells: a new method for gene function analysis in marine bivalve. Mol Biol Rep. 2012;39(8):8209–16.

    Article  CAS  PubMed  Google Scholar 

  31. 31.

    Hou R, Bao Z, Wang S, Su H, Li Y, Du H, et al. Transcriptome sequencing and de novo analysis for yesso scallop (Patinopecten yessoensis) using 454 GS FLX. PLoS One. 2011;6(6):e21560.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  32. 32.

    Du H, Bao Z, Hou R, Wang S, Su H, Yan J, et al. Transcriptome sequencing and characterization for the sea cucumber Apostichopus japonicus (Selenka, 1867). PLoS One. 2012;7(3):e33311.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  33. 33.

    Jung H, Lyons RE, Dinh H, Hurwood DA, McWilliam S, Mather PB. Transcriptomics of a giant freshwater prawn (Macrobrachium rosenbergii): de novo assembly, annotation and marker discovery. PLoS One. 2011;6(12):e27938.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  34. 34.

    Kim E, Magen A, Ast G. Different levels of alternative splicing among eukaryotes. Nucleic Acids Res. 2007;35(1):125–31.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  35. 35.

    Novaes E, Drost DR, Farmerie WG, Pappas Jr GJ, Grattapaglia D, Sederoff RR, et al. High-throughput gene and SNP discovery in Eucalyptus grandis, an uncharacterized genome. BMC Genomics. 2008;9:312.

    Article  PubMed Central  PubMed  Google Scholar 

  36. 36.

    Gilbert JA, Clark MS, Thorne MAS, Toullec J-Y, Meng Y, Guan LL, et al. Antarctic krill 454 pyrosequencing reveals chaperone and stress transcriptome. PLoS One. 2011;6(1):e15919.

    Article  Google Scholar 

  37. 37.

    Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, et al. Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000;25(1):25–9.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  38. 38.

    Zhan A, Hu J, Hu X, Hui M, Wang M, Peng W, et al. Construction of microsatellite based linkage maps and identification of size related quantitative trait loci for Zhikong scallop (Chlamys farreri). Anim Genet. 2009;40(6):821–31.

    Article  CAS  PubMed  Google Scholar 

  39. 39.

    Hale MC, McCormick CR, Jackson JR, Dewoody JA. Next-generation pyrosequencing of gonad transcriptomes in the polyploid lake sturgeon (Acipenser fulvescens): the relative merits of normalization and rarefaction in gene discovery. BMC Genomics. 2009;10:203.

    Article  PubMed Central  PubMed  Google Scholar 

  40. 40.

    Zakas C, Schult N, McHugh D, Jones KL, Wares JP. Transcriptome analysis and SNP development can resolve population differentiation of Streblospio benedicti, a developmentally dimorphic marine annelid. PLoS One. 2012;7(2):e31613.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

Download references

Acknowledgments

We are grateful to Dr. Chiju Wei for his careful revision and many constructive comments. Funding for this research was provided by National Natural Science Foundation of China (41076107, 31372528), Ministry of Education of P.R. China (20114402110001), China Modern Agro-industry Technology Research System (CARS-48), National Basic Research Program of China (973 Program, No. 2010CB126402), Knowledge Innovation Program of Deep Sea Science and Engineering, Chinese Academy of Sciences (SIDSSE-QN-201407), Department of Education (2050205-95), Department of Science & Technology (2013B020503061) and Oceanic and Fisheries Administrator (B201300B06) of Guangdong Province, China.

Author information

Affiliations

Authors

Corresponding authors

Correspondence to Huaiping Zheng or Guofan Zhang.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

Conceived and designed the experiments: HL, HZheng. Performed the experiments: HL, HZhang. Analyzed the data: HL, LD, FM, YW. Contributed reagents/materials/analysis tools: HL, HZhang, ZG, GZ. Wrote the paper: HL, HZheng, WL, SW, SL. All authors read and approved the final manuscript.

Additional files

Additional file 1: Table S1.

Sequences with significant BlastX matches against Swiss-Prot and Nr database.

Additional file 2: Table S2.

Gene Ontology of C. nobilis unique sequences.

Additional file 3: Table S3.

Details of SSR motifs in C. nobilis unique sequences.

Additional file 4: Table S4.

Details of SNPs and Indels in C. nobilis isotigs.

Additional file 5: Table S5.

Validation of predicted SNPs by PCR and re-sequencing.

Additional file 6: Table S6.

Candidate transcripts.

Additional file 7: Table S7.

Information of 8 candidate carotenoid deposition genes.

Rights and permissions

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Liu, H., Zheng, H., Zhang, H. et al. A de novo transcriptome of the noble scallop, Chlamys nobilis, focusing on mining transcripts for carotenoid-based coloration. BMC Genomics 16, 44 (2015). https://doi.org/10.1186/s12864-015-1241-x

Download citation

Keywords

  • Chlamys nobilis
  • Transcriptome sequencing
  • Carotenoid coloration
  • Candidate genes