Next Article in Journal
Towards the Elucidation of Assimilative nasABC Operon Transcriptional Regulation in Haloferax mediterranei
Previous Article in Journal
Comparative Distribution of Repetitive Sequences in the Karyotypes of Xenopus tropicalis and Xenopus laevis (Anura, Pipidae)
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 12
Issue 5
10.3390/genes12050618
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Transcriptome Analysis Provides Insights into the Mechanism of Astaxanthin Enrichment in a Mutant of the Ridgetail White Prawn Exopalaemon carinicauda

by
Yue Jin
1,2,3,
Shihao Li
1,2,4,*,
Yang Yu
1,2,4,
Chengsong Zhang
1,2,4,
Xiaojun Zhang
1,2,4 and
Fuhua Li
1,2,4,5,*
1
Key Laboratory of Experimental Marine Biology, Institute of Oceanology, Chinese Academy of Sciences, Qingdao 266071, China
2
Laboratory for Marine Biology and Biotechnology, Qingdao National Laboratory for Marine Science and Technology, Qingdao 266000, China
3
University of Chinese Academy of Sciences, Beijing 100049, China
4
Center for Ocean Mega-Science, Chinese Academy of Sciences, Qingdao 266071, China
5
The Innovation of Seed Design, Chinese Academy of Sciences, Wuhan 430072, China
*
Authors to whom correspondence should be addressed.
Genes 2021, 12(5), 618; https://doi.org/10.3390/genes12050618
Submission received: 29 March 2021 / Revised: 17 April 2021 / Accepted: 19 April 2021 / Published: 21 April 2021
(This article belongs to the Section Animal Genetics and Genomics)
Browse Figures
Review Reports Versions Notes

Abstract

:
A mutant of the ridgetail white prawn, which exhibited rare orange-red body color with a higher level of free astaxanthin (ASTX) concentration than that in the wild-type prawn, was obtained in our lab. In order to understand the underlying mechanism for the existence of a high level of free astaxanthin, transcriptome analysis was performed to identify the differentially expressed genes (DEGs) between the mutant and wild-type prawns. A total of 78,224 unigenes were obtained, and 1863 were identified as DEGs, in which 902 unigenes showed higher expression levels, while 961 unigenes presented lower expression levels in the mutant in comparison with the wild-type prawns. Based on Gene Ontology analysis and Kyoto Encyclopedia of Genes and Genomes analysis, as well as further investigation of annotated DEGs, we found that the biological processes related to astaxanthin binding, transport, and metabolism presented significant differences between the mutant and the wild-type prawns. Some genes related to these processes, including crustacyanin, apolipoprotein D (ApoD), cathepsin, and cuticle proteins, were identified as DEGs between the two types of prawns. These data may provide important information for us to understand the molecular mechanism of the existence of a high level of free astaxanthin in the prawn.

1. Introduction

Carotenoids are tetraterpene pigments that can form a series of colors such as red, orange, and yellow throughout organisms from fungus to vertebrates. A total of 850 carotenoids have been found in nature, of which about 50 are carotenes and 800 are xanthophylls [1,2]. Carotenes are hydrocarbons, including α-carotene, β-carotene, and so on, while xanthophylls, such as lutein, zeaxanthin, and astaxanthin, contain oxygen atoms such as hydroxy [3]. Carotenoids are important nutrients for animals, acting as vitamin A providers and antioxidants [4], with essential effects on ovarian maturation and growth [5,6]. In addition, the presence of carotenoids in the feathers, skin, and muscles always endows the animals with carotenoid-based colors [3,7,8]. Carotenoids exhibit diverse forms in animals. For example, β-carotene, β-cryptoxanthin, lutein, and astaxanthin can be accumulated in birds [7,9,10]. β-carotene, zeaxanthin, 3-hydroxyechinenone, and astaxanthin are found in some insects [11,12]. Astaxanthin is enriched in salmonid flesh [13]. Pectenolone is in a rare variant of the shellfish Yesso scallop [14]. In crustaceans, astaxanthin is the main carotenoid and accumulated in the shell, the epidermis of crabs, and prawns [15,16,17].
However, carotenoids cannot be synthesized de novo in animals, which are usually accumulated from diets or modified through metabolic reactions. The accumulation of carotenoids in animals involves a series of processes, such as absorption, transport, metabolism, binding, and storage [18]. Usually, the intestine and stomach are important tissues for carotenoid absorption, and the dissolution of carotenoids is the first phase of the process of absorption [19]. Carotenoids can be emulsified into lipid droplets in the stomach, where gastric lipase hydrolyses play important roles [20]. During the process of digestion, carotenoids are incorporated with other lipids into the mixed micelles, which is necessary for their uptake by the enterocyte. Other steps such as uptake in enterocytes and incorporation for transport into chylomicrons are then followed [19,21,22]. Once carotenoids are taken up, they will be metabolized. For animals, hepatopancreas and β-carotene oxygenases are essential metabolic tissues and enzymes. Finally, these carotenoids are metabolized into other metabolites or deposited directly into the skin, eyes or muscles. So far, studies have found that some factors from the environment or the underlying genetic mechanisms may affect the accumulation of carotenoids in animals. Research has proven that a genetic basis, such as some transport proteins or enzymes, plays essential functions in these processes [18,23,24]. Scavenger receptor class B (SRB) can recognize lipoproteins carrying carotenoids, such as lutein, carotene, and zeaxanthin, and then facilitate the delivery of carotenoids into target tissues [25,26,27]. The relationship between SRB and carotenoids uptake has been demonstrated in birds, salmonids, Drosophila, and scallops [13,28,29,30]. Lipoproteins such as the high-density lipoprotein (HDL), vitellogenin, and ATP-binding cassette transporter have been demonstrated as carriers that can control the carotenoid transport [31,32,33]. Additionally, β-carotene oxygenase is involved in the cleavage of carotenoids, and its mutation results in carotenoid accumulation in chickens and scallops [34,35].
In crustaceans, several genes were found related to carotenoid binding and carotenoid-based colors. Crustacyanin, a carotenoid-binding protein that was isolated from the carapace of the lobster Homarus gammarus, was also associated with astaxanthin binding [36]. The original color of astaxanthin is orangish red, and combining astaxanthin with crustacyanin can increase the variation in color. When boiling, the denaturation of crustacyanin causes the release of astaxanthin, which turns the shell of the crustacean red [37]. In Fenneropenaeus merguiensis, sarcoplasmic calcium-binding protein, arginine kinase, and actin were considered to play a role in astaxanthin-based color formation [38]. In Halocaridina rubra, the cytochrome P450 family was thought as the possible candidate gene responsible for the production of astaxanthin [39]. In E. carinicauda, knockdown of EcBCO-like6, encoding carotenoid oxygenases, led to a variation in carotenoid levels and color changes in the hepatopancreas [40].
The wild-type ridgetail white prawn (E. carinicauda) is usually in transparent body color, with blue or red spots. Recently, we isolated a mutant with orange-red body color. The orange-red-variant prawns have a significantly higher concentration of astaxanthin compared with the wild-type prawns. Interestingly, astaxanthin in the orange-red prawns is mainly in a free form in the muscles (94%) and waste (60.4%), while it is mainly in an esterified form in the muscles (75.6%) and waste (78.4%) of the wild-type prawns [41], which exhibits a great difference with the known existing form of astaxanthin in crustaceans. However, the underlying mechanism regulating the free astaxanthin enrichment in E. carinicauda is not clear. In the present study, comparative transcriptome analysis based on Illumina high-throughput sequencing was performed on the orange-red-variant and wild-type prawns. A number of differentially expressed genes between them were identified and analyzed. The present data provided useful information to further investigate the molecular mechanisms of carotenoid enrichment and carotenoid-based body color variation in crustaceans.

2. Materials and Methods

2.1. Sample Collection

The orange-red-variant and wild-type prawns (Figure S1) were both cultured in our laboratory. A total of 60 individuals at the intermolt stage, including 30 orange-red-variant prawns and 30 wild-type prawns, were sampled for transcriptome analysis. The prawns at the intermolt stage were selected for experiments according to Gao et al. [42]. Additionally, another 25 orange-red-variant prawns and 25 wild-type prawns at the intermolt stage were collected for quantitative real-time PCR analysis (qPCR). The average body length and body weight of the prawns were 4.38 ± 0.42 cm and 0.98 ± 0.31 g, respectively. The cephalothoraxes of these prawns were dissected and immediately frozen in liquid nitrogen. Intestine was also dissected individually and mixed with cephalothorax. Tissues were stored at −80 °C for total RNA extraction.

2.2. RNA Isolation, Sample Pooling, and cDNA Synthesis

Total RNA was extracted separately from samples by RNAiso (TaKaRa, Kyoto, Japan) according to the manufacturer’s instructions. The extracted RNA was then assessed by electrophoresis on 1% agarose gel and quantified by a NanoDrop 2000 spectrophotometer (Thermo Fisher Scientific Inc., Walthman, MA, USA). The total RNA isolated from 30 orange-red-variant prawns were randomly divided into three subgroups. Equivalent total RNA from each sample of one subgroup (containing 10 individuals) was mixed. The total RNA from the wild-type prawns was treated accordingly. Finally, three RNA sample replicates from the orange-red prawns (EcR-1, EcR-2, and EcR-3) and three replicates from the wild prawns (EcW-1, EcW-2, and EcW-3) were prepared for RNA-seq analysis.
For another 25 orange-red-variant prawns, they were randomly divided into five sample replicates with five individuals in each replicate. The 25 wild-type prawns were grouped in the same way; the total RNA from each sample replicate was extracted and assessed as mentioned above. About 1 μg RNA was then used for cDNA generation with the PrimeScript RT Reagent Kit (TaKaRa, Kyoto, Japan). According to the manufacturer’s instructions, the genomic DNA (gDNA) was firstly removed by using a gDNA Eraser, and then the first-strand cDNA was synthesized by a PrimeScript RT Enzyme with random primers.

2.3. Library Construction and Sequencing

The construction and sequencing of RNA-seq libraries were conducted according to Illumina’s RNA sample preparation protocols by Gene Denovo Biotechnology Co. (Guangzhou, China). Briefly, the mRNA from each sample was enriched using Oligo(dT) beads. The enriched mRNA was fragmented into short fragments by using fragmentation buffer. After that, these short fragments were used as templates for first-strand cDNA synthesis with random primers. The second-strand cDNA was then synthesized by DNA polymerase I, RNase H, dNTP, and buffer. The cDNA fragments were purified with the QiaQuick PCR Extraction Kit (Qiagen, Duesseldorf, Germany), end repaired, poly(A) added, and ligated to Illumina-sequencing adapters. The ligation products were size selected by agarose gel electrophoresis, PCR amplified, and sequenced using Illumina HiSeqTM 4000.

2.4. Assembly and Annotation

The raw reads were filtered by removing adaptors. Low-quality reads containing more than 10% of unknown nucleotides or more than 50% of low-quality bases (Q ≤ 10) were removed. The high-quality reads were mapped to a ribosome RNA (rRNA) sequence to remove residual rRNA reads. The clean reads were then assembled into unigenes using the Trinity program, which combines three components: Inchworm, Chrysalis, and Butterfly [43]. For annotation analysis, unigenes were aligned by BLASTx (E < 10−5) to databases, including the National Center for Biotechnology Information non-redundant protein (Nr) database (http://www.ncbi.nlm.nih.gov, accessed on 19 November 2018), the Swiss-Prot protein database (http://www.expasy.ch/sprot, accessed on 19 November 2018), the Kyoto Encyclopedia of Genes and Genomes (KEGG) database (http://www.genome.jp/kegg, accessed on 19 November 2018), and the Cluster of Orthologous Groups (COG/KOG) database (http://www.ncbi.nlm.nih.gov/COG, accessed on 19 November 2018). Protein functional annotations could then be obtained according to the best alignment results.

2.5. Analysis of Differentially Expressed Unigenes

Gene abundances were calculated and normalized to reads per kb per million (RPKM) reads [44]. To identify differentially expressed genes (DEGs) across samples or groups, the edgeR package was used. Genes with a fold change ≥2 and a false discovery rate (FDR) <0.05 in the two groups were deemed DEGs. DEGs were then subjected to enrichment analysis of GO functions and KEGG pathways. The calculated p-value went through FDR correction, taking an FDR ≤0.05 as a threshold. Only GO terms and KEGG pathways meeting this condition could be defined as significantly enriched terms and pathways in DEGs.

2.6. Multiple Sequence Alignment and Phylogenetic Analysis

Differentially expressed lipoproteins encoding genes were used for further phylogenetic analysis. Complete open reading frame (ORF) regions and amino acid sequences were deduced using ORFfinder (https://www.ncbi.nlm.nih.gov/orffinder/, accessed on 8 April 2021). Conserved protein domains were predicted with SMART databases (http://smart.embl-heidelberg.de/, accessed on 8 April 2021) [45]. Amino acid sequences of lipoproteins from different species for phylogenetic analysis were obtained from the GenBank™/EBI or UniProt databases [46]. Multiple sequence alignment was performed using ClustalW (https://myhits.sib.swiss/cgi-bin/clustalw, accessed on 8 April 2021), and phylogenetic analysis was performed using MEGA7.0 software [47]. The phylogenetic tree was constructed by the neighbor-joining (NJ) distance algorithm [48]. A bootstrapping test was adopted with 1000 replications.
Mutation analysis was performed for the differentially expressed crustacyanin genes. Clean reads from the wild-type prawns and the orange-red-variant prawns were used to assemble their transcriptome data separately. The sequences of crustacyanin from these two transcriptomes were then extracted. Deduced amino acid sequences were obtained using ORFfinder, and multiple sequence alignment was performed using ClustalW.

2.7. Validation by Quantitative Real-Time PCR

In order to verify the accuracy of RNA-seq data, nine unigenes were selected to detect their expression levels in orange-red prawns and wild prawns by quantitative real-time PCR (qPCR) analysis using the Eppendorf Mastercycler ep realplex (Eppendorf, Hamburg, Germany). The 18S rRNA gene was used as an internal reference. Information of the primers is listed in Table 1. The qPCR program was set as follows: one cycle of 94 °C for 2 min, followed by 40 cycles of: 15 s at 95 °C, 15 s at annealing temperature for each pair of primers, and 20 s at 72 °C, and a final 20 min melting curve step was added to test the specificity of these products. All samples were performed in quadruplicates. Relative gene expression levels were calculated following the comparative Ct method with the analysis formula 2-ΔΔCt [49,50,51]. The qPCR results were compared with the transcriptome data to detect the expression correlation of the selected genes.

2.8. Statistical Analyses

Statistical significance between orange-red-variant prawns and wild-type prawns was analyzed by using SPSS (version 20). All data obtained from the studies were presented as means ± SE, and the results were analyzed by one-way analysis of variance (ANOVA). Significant differences at p < 0.05 were labeled with asterisks.

3. Results

3.1. Basic Information of the Transcriptome

In total, 76558224, 75407822, 78269646, 84068372, 72076762, and 70924248 raw reads were obtained for each sample by using Illumina HiSeqTM 4000. After filtration, 73335584 (95.79%), 72788448 (96.53%), 75141044 (96.00%), 80832778 (96.15%), 69700338 (96.7%), and 68597552 (96.72%) high-quality clean reads were obtained, respectively (Table 2). The filtered high-quality clean reads were then mapped to the reference transcriptome, and 87.67–88.34% filtered reads were mapped. Finally, these clean reads were assembled and then clustered into 78,224 unigenes with an average length of 824 bp and half of the total assembly (N50) length of 1565 bp. A total of 75,947 existed in the EcR group, and 76,247 existed in the EcW group. In total, 19,461 unigenes were annotated in the database, including 19,285 (24.7%) unigenes with an annotation in the Nr database, 13,808 (17.7%) unigenes in the Swiss-Prot database, 5934 (7.59%) unigenes in the GO database, 12,732 (16.3%) unigenes in the KOG database, and 9150 (11.7%) unigenes in the KEGG database.

3.2. Identification and Verification of Differentially Expressed Genes (DEGs)

A total of 1863 unigenes were identified as DEGs between the orange-red prawns (EcR) and the wild-type prawns (EcW), including 902 highly expressed DEGs and 961 low-expressed DEGs in EcR compared with the wild-type prawns. Nine unigenes (Table 1) were selected for quantitative real-time PCR analysis. Among them, Unigene0059050 encoded an esterase, which could hydrolyze esterified compounds. Unigene0036561 was annotated as an apolipoprotein, which might be related to the process of transport. Unigene0054180 was annotated as crustacyanin. Unigene0037069, Unigene0038920, and Unigene0043447 were the genes coding cuticle proteins. These genes might be related to transport and binding functions. The transcriptome data showed that most selected genes showed significantly lower expression levels in the orange-red prawns than in the wild-type prawns (Figure 1A). The qPCR results showed highly consistent expression profiles of these genes compared to the transcriptome data (Figure 1B), which suggested that the transcriptome data were reliable.

3.3. GO Term and KEGG Pathway Enrichment Analysis of DEGs

In order to understand which biological processes and signaling pathways might participate in the formation of orange-red color in E. carinicauda, GO and KEGG analyses were performed on the DEGs. All DEGs were sorted into 29 functional groups (Figure 2). For the group EcW and EcR, the significant GO molecular function terms were mainly concentrated into categories including “catalytic activity”, “binding”, “transporter activity”, and “structural molecule activity”. Of the biological-process-related genes, most were involved in a “cellular process”, a “single-organism process”, and a “metabolic process”. “Membrane” is the most significantly enriched GO term of cellular components. Results also revealed that the top five terms where DEGs had a much higher expression level in the group EcR included “catalytic activity”, “cellular process”, “binding”, “single-organism process”, and “metabolic process”. Of the low-expressed DEGs, they were mainly enriched in three terms, such as “catalytic activity”, “metabolic process”, and “binding”. In general, “catalytic activity”, “metabolic process”, and “binding” were the terms with the most differentially expressed unigenes.
Only 29 identified DEGs were enriched in pathways by KEGG analysis. The top 20 pathways were shown in Figure 3A. However, most of the enriched pathways only contained one unigene with a Q-value greater than 0.1. The “Phagosome” and “Lysosome” were the most significantly enriched pathways (Q < 0.05).

3.4. Functional Classification of DEGs

In order to further identify genes related to astaxanthin enrichment in orange-red prawns, DEGs with an annotation were categorized according to their reported functions. Among the 1863 DEGs, 279 unigenes were annotated in the Nr database, including 139 highly expressed DEGs (Table S1) and 140 low-expressed DEGs (Table S2) in the EcR prawns compared with the wild-type prawns.
Notably, according to the transcriptome analysis, a total of 31 unigenes encoding cuticular proteins or related proteins, which were with similar domains or functional annotations, exhibited higher expression levels in EcW than in EcR (Table 3). These genes were classified into three groups based on their domains. Among them, the deduced amino acid sequences of seven unigenes, including Unigene0017103, Unigene0032773, Unigene0035142, Unigene0037069, Unigene0038920, Unigene0043447 and Unigene0073146 contained a cuticle_1 domain. Unigene0003529, Unigene0019577, Unigene0019752, Unigene0019787, Unigene0022558, Unigene0029471, Unigene0039421, and Unigene00056168 encoded proteins containing a chitin_bind_4 domain, while Unigene0061823 and Unigene0001656 encoded proteins containing a CBM 14 domain.
Binding is an important process for carotenoid metabolism. GO analysis also revealed “binding” was the main term with much more differentially expressed unigenes. Except for cuticular proteins, many other structural genes in the “binding” term, such as tubulin, fibrillin, and fibronectin, were significantly higher in the orange-red prawns compared with those in the wild-type prawns. Additionally, as shown in Table 4, carotenoid-binding protein crustacyanin and hemocyanin genes also have different expression levels between the group EcR and EcW. The result of the sequence alignment revealed that crustacyanin in the wild-type prawns and the orange-red-variant prawns shared 94.85% sequence identity, and both shared a relative similarity with crustacyanin subunit A3 protein in the lobster (Figure 4). Several amino acid residues (residues 6, 87, 92, 104–106, 108, 118–119, 124–125, 129, 131, 134, 144, 191, and 193) showed differences between the EcR and EcW groups. Interestingly, significant changes in the region (PYEVIETDYETYSC) were found, which was the lipocalin consensus region (Figure 4, green colored).
Genes related to transport and metabolism were also analyzed. DEGs including vitellogenin, ApoD, and StAR-related lipid transfer protein were related to the transport process (Table 4). Genes annotated as vitellogenin were all with higher expression levels in the orange-red-variant prawns. Phylogenetic analysis showed that they belonged to apolipocrustacein (apoCr) rather than vitellogenin in crustaceans (Figure 5). By contrast, the three apolipoprotein genes showed different patterns, some of which were with higher expression levels, while the others were with lower expressions in the variant. The result of sequence alignment revealed that Unigene0036561 and Unigene0044620 shared a 65.63% sequence similarity, while they shared 40.61% and 37.31% sequence similarities with Unigene0028810, respectively (Figure S2). Phylogenetic analysis also demonstrated that Unigene0036561 and Unigene0044620 were clustered together, which were all with lower expression levels in the orange-red-variant prawns. In addition, most pathways in metabolism were downregulated in orange-red prawns (Figure 3B). Among them, “Lysosome” was the most significantly enriched pathway, and DEGs such as cathepsin B were included in this pathway (Table 5).

4. Discussion

Carotenoid absorption, transport, and metabolism are several crucial steps of carotenoid-dependent coloration. In the present study, the expression profiles of many genes related to carotenoid absorption, transport, metabolism, binding, and deposition were compared between orange-red prawns and wild-type prawns. Lipid digestion and micelle formation are important processes of carotenoid absorption before its uptake in enterocytes [21,52]. After the dissolution of carotenoids, fat-soluble carotenoids are incorporated with other lipids into mixed micelles, which are composed of phospholipids, cholesterol, lipid digestion products, and bile salts [53]. It seems that only free forms of carotenoids are absorbed by the intestinal mucosa, suggesting that esterified forms are first hydrolyzed, and esterification may affect the absorption efficiency [19,54]. The enzymes involved in the hydrolysis of carotenoid esters are probably lipase and esterase, as they have been shown to hydrolyze carotenoid esters [55]. However, genes encoding these kinds of enzymes were hardly found in the DEGs of the present transcriptome, except for Unigene0059050, which encoded a carboxylesterase (esterase FE4) and exhibited higher expression levels in EcW than in EcR. In addition to carotenoid absorption, previous studies have suggested that the incorporation of carotenoids into the carotenoid-binding protein complex and subsequent deposition require a step for the hydrolysis of carotenoid esters [16]. Therefore, the transcriptome data showed that there was no obvious difference in the carotenoid absorption process between EcW and EcR, and the high expression of one carboxylesterase encoding gene in EcW might be related to other carotenoid metabolic processes.
The lipophilicity of carotenoids determines that their transport to various tissues via circulatory and lymphatic systems is intimately linked to lipoproteins. Once loaded with carotenoids, apolipoproteins, components of lipoproteins, can help absorbed carotenoids target various tissues [56,57,58]. Apolipoprotein D (ApoD), a member of the lipocalin family, was once considered a candidate gene associated with carotenoid transport [59]. In Macrobrachium rosenbergii, knockdown of the MLC gene, a member of the lipocalin family gene, resulted in a body color shift from blue to orange-red [60]. Additionally, astaxanthin can be transported by HDL and vitellogenin in the chum salmon [33,61]. StARD3 was also reported to be involved in carotenoid transport [62]. The differential expression of these genes in the orange-red prawns and the wild-type prawns indicated that the transport process might be important for carotenoid enrichment in the orange-red prawns.
Metabolism and utilization in target tissues will greatly contribute to the in vivo accumulation of astaxanthin. Phagocytosis is a process whereby cells engulf large particles, which involves engulfment and internalization [63]. Once a target is bound and the necessary signals are activated, these nascent phagosomes initiate the cellular degradation processes, leading to the eventual destruction of the phagocytosed particle. Usually, there are interactions between phagocytes, lysosomes, and endosomes [63,64]. Lysosome, which contains more than 60 enzymes, is the main organelle for substance degradation and metabolism [65]. In the orange-red prawns, the “lysosome” pathway was mainly enriched with related genes at lower expression levels compared with those in the wild-type prawns. Compared with these two pathways in the wild-type prawns, differences between “phagocytosis” and “lysosome” processes might affect the utilization level and accumulation of astaxanthin in the orange-red prawns, which would also contribute to the orange-red body color of the orange-red prawns.
Another process, the binding of carotenoids to a specific protein, is usually related to body coloration and an essential step for carotenoid deposition. In the present study, plenty of genes related to carotenoid binding were identified as DEGs in the transcriptome data. In invertebrates, crustacyanin is proven to bind and stabilize astaxanthin by forming carotene–protein complexes [66]. Astaxanthin is in the free state within the protein complex [67,68]. Once the connection between them is broken, astaxanthin will dissociate from crustacyanin and reveal its natural color. A lower expression level of crustacyanin might lead to a limitation in the binding and subsequent utilization of free astaxanthin in the orange-red prawns. Furthermore, we found that several amino acid residues in the lipocalin consensus region (PYEVIETDYETYSC) were different in a crustacyanin between the wild-type and orange-red-variant prawns. In the crustacyanin homolog A3 protein, several β strands, which are the components of β sheets, are included in this region, which is essential for carotenoid binding [68]. This indicated that a change of these residues might affect the structure and function of the crustacyanin lipocalin domain, leading to a different astaxanthin-binding ability in the orange-red-variant prawns. However, it still needs further investigation to reveal whether they contribute to the binding function with astaxanthin. Additionally, a group of genes encoding cuticular proteins or related proteins showed lower expression levels in the orange-red prawns compared with those in the wild-type prawns. Cuticle proteins are a kind of structural protein involved in the formation of the crustacean exoskeleton. These proteins bind to chitin and play functions in cuticle calcification through different types of domains, including cuticle_1, chitin_bind_4, and CBM 14 [69]. Cuticle proteins also protect insects from adverse environmental stresses, such as UV, insecticides, and infection by pathogens [70,71]. Additionally, in the silkworm, deficiency of one cuticle protein, BmorCPH24, could disrupt Wnt1 function and lead to abnormal body coloration [72]. A lower expression of cuticular proteins in the orange-red prawns might affect their body color by influencing other pigment-related genes, pathways, or the layer of pigment cells in the cuticle.

5. Conclusions

In conclusion, the present study identified some genes that differentially expressed between the orange-red prawns and the wild-type prawns through comparative transcriptome analysis. GO term enrichment analysis and KEGG pathway enrichment analysis, as well as further investigation of annotated DEGs, revealed the main differences in the biological processes related to astaxanthin binding and transport and metabolism between two types of prawns. The present data provided insights into understanding the molecular mechanism of the body coloration of the orange-red prawns.

Supplementary Materials

The following are available online at: https://www.mdpi.com/article/10.3390/genes12050618/s1, Figure S1: Figure of the wild-type prawn (up) and the orange-red-variant prawn (down), Figure S2: Alignment of apolipoprotein D amino acid sequences, Table S1: Information of DEGs with higher expression levels in EcR prawns, Table S2: Information of DEGs with lower expression levels in EcR prawns, Table S3: Information of lipoproteins used in phylogenetic analysis.

Author Contributions

Contributed conception and design of the study, S.L. and F.L.; conducted the investigation, performed all experiments, and wrote the manuscript. Y.J.; collected data and performed the statistical analysis for the work, Y.Y., C.Z., and X.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the General Program of the National Natural Science Foundation of China (31872552), the National Key R&D Program of China (2018YFD0900103), and the Major Project of Agriculture New Variety of Jiangsu Province (PZCZ201747).

Institutional Review Board Statement

All animals were handled and treated according to the alpaca guidelines approved by the Animal Ethics Committee [2020(37)] at the Institute of Oceanology, Chinese Academy of Sciences (Qingdao, Shandong, China).

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Britton, G.; Liaaen-Jensen, S.; Pfander, H. Carotenoids Handbook; Birkhäuser: Basel, Switzerland, 2004. [Google Scholar]
  2. Maoka, T. Recent progress in structural studies of carotenoids in animals and plants. Arch. Biochem. Biophys. 2009, 483, 191–195. [Google Scholar] [CrossRef]
  3. Maoka, T. Carotenoids as natural functional pigments. J. Nat. Med. 2020, 74, 1–16. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Rao, A.V.; Rao, L.G. Carotenoids and human health. Pharm. Res. 2007, 55, 207–216. [Google Scholar] [CrossRef] [PubMed]
  5. Dall, W.S.; Smith, D.M.; Moore, L.E. Carotenoids in the tiger prawn Penaeus esculentus during ovarian maturation. Mar. Biol. 1995, 123, 435–441. [Google Scholar] [CrossRef]
  6. Liñán-Cabello, M.A.; Paniagua-Michel, J. Induction factors derived from carotenoids and vitamin A during the ovarian maturation of Litopenaeus vannamei. Aquac. Int. 2004, 12, 583–592. [Google Scholar] [CrossRef]
  7. Stradi, R.; Pini, E.; Celentano, G. Carotenoids in bird plumage: The complement of red pigments in the plumage of wild and captive bullfinch (Pyrrhula pyrrhula). Comp. Biochem. Physiol. B Biochem. Mol. Biol. 2001, 128, 529–535. [Google Scholar] [CrossRef]
  8. Lehnert, S.J.; Christensen, K.A.; Vandersteen, W.E.; Sakhrani, D.; Pitcher, T.E.; Heath, J.W.; Koop, B.F.; Heath, D.D.; Devlin, R.H. Carotenoid pigmentation in salmon: Variation in expression at BCO2-l locus controls a key fitness trait affecting red coloration. Proc. Biol. Sci. 2019, 286, 20191588. [Google Scholar] [CrossRef] [Green Version]
  9. Lopes, R.J.; Johnson, J.D.; Toomey, M.B.; Ferreira, M.S.; Araujo, P.M.; Melo-Ferreira, J.; Andersson, L.; Hill, G.E.; Corbo, J.C.; Carneiro, M. Genetic Basis for Red Coloration in Birds. Curr. Biol. 2016, 26, 1427–1434. [Google Scholar] [CrossRef] [Green Version]
  10. Mundy, N.I.; Stapley, J.; Bennison, C.; Tucker, R.; Twyman, H.; Kim, K.W.; Burke, T.; Birkhead, T.R.; Andersson, S.; Slate, J. Red Carotenoid Coloration in the Zebra Finch Is Controlled by a Cytochrome P450 Gene Cluster. Curr. Biol. 2016, 26, 1435–1440. [Google Scholar] [CrossRef] [Green Version]
  11. Zhu, L.; Zhang, Y.Q. Identification and analysis of the pigment composition and sources in the colored cocoon of the silkworm, Bombyx mori, by HPLC-DAD. J. Insect Sci. 2014, 14, 31. [Google Scholar] [CrossRef] [Green Version]
  12. Wybouw, N.; Kurlovs, A.H.; Greenhalgh, R.; Bryon, A.; Kosterlitz, O.; Manabe, Y.; Osakabe, M.; Vontas, J.; Clark, R.M.; Van Leeuwen, T. Convergent evolution of cytochrome P450s underlies independent origins of keto-carotenoid pigmentation in animals. Proc. Biol. Sci. 2019, 286, 20191039. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Sundvold, H.; Helgeland, H.; Baranski, M.; Omholt, S.W.; Vage, D.I. Characterisation of a novel paralog of scavenger receptor class B member I (SCARB1) in Atlantic salmon (Salmo salar). BMC Genet. 2011, 12, 52. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Li, N.; Hu, J.J.; Wang, S.; Cheng, J.; Hu, X.L.; Lu, Z.Y.; Lin, Z.J.; Zhu, W.M.; Bao, Z.M. Isolation and identification of the main carotenoid pigment from the rare orange muscle of the Yesso scallop. Food Chem. 2010, 118, 616–619. [Google Scholar] [CrossRef]
  15. Coral-Hinostroza, G.N.; Bjerkeng, B. Astaxanthin from the red crab langostilla (Pleuroncodes planipes): Optical R/S isomers and fatty acid moieties of astaxanthin esters. Comp. Biochem. Physiol. B Biochem. Mol. Biol. 2002, 133, 437–444. [Google Scholar] [CrossRef]
  16. Wade, N.; Goulter, K.C.; Wilson, K.J.; Hall, M.R.; Degnan, B.M. Esterified astaxanthin levels in lobster epithelia correlate with shell colour intensity: Potential role in crustacean shell colour formation. Comp. Biochem. Physiol. B Biochem. Mol. Biol. 2005, 141, 307–313. [Google Scholar] [CrossRef] [PubMed]
  17. Nm, S.; Narayan, B.; Mahendrakar, N. Carotenoids in different body components of Indian shrimps. J. Sci. Food Agric. 2005, 85, 167–172. [Google Scholar]
  18. Toews, D.P.L.; Hofmeister, N.R.; Taylor, S.A. The evolution and genetics of carotenoid processing in animals. Trends Genet. 2017, 33, 171–182. [Google Scholar] [CrossRef]
  19. Reboul, E. Absorption of vitamin A and carotenoids by the enterocyte: Focus on transport proteins. Nutrients 2013, 5, 3563–3581. [Google Scholar] [CrossRef] [Green Version]
  20. Carriere, F.; Barrowman, J.A.; Verger, R.; Laugier, R. Secretion and contribution to lipolysis of gastric and pancreatic lipases during a test meal in humans. Gastroenterology 1993, 105, 876–888. [Google Scholar] [CrossRef]
  21. van Het Hof, K.H.; West, C.E.; Weststrate, J.A.; Hautvast, J.G. Dietary factors that affect the bioavailability of carotenoids. J. Nutr. 2000, 130, 503–506. [Google Scholar] [CrossRef]
  22. Parker, R.S. Absorption, metabolism, and transport of carotenoids. FASEB J. 1996, 10, 542–551. [Google Scholar] [CrossRef] [PubMed]
  23. Pointer, M.A.; Prager, M.; Andersson, S.; Mundy, N.I. A novel method for screening a vertebrate transcriptome for genes involved in carotenoid binding and metabolism. Mol. Ecol. Resour. 2012, 12, 149–159. [Google Scholar] [CrossRef] [PubMed]
  24. Walsh, N.; Dale, J.; McGraw, K.J.; Pointer, M.A.; Mundy, N.I. 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, 58–66. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Reboul, E.; Abou, L.; Mikail, C.; Ghiringhelli, O.; Andre, M.; Portugal, H.; Jourdheuil-Rahmani, D.; Amiot, M.J.; Lairon, D.; Borel, P. Lutein transport by Caco-2 TC-7 cells occurs partly by a facilitated process involving the scavenger receptor class B type I (SR-BI). Biochem. J. 2005, 387, 455–461. [Google Scholar] [CrossRef] [PubMed]
  26. During, A.; Doraiswamy, S.; Harrison, E.H. Xanthophylls are preferentially taken up compared with β-carotene by retinal cells via a SRBI-dependent mechanism. J. Lipid Res. 2008, 49, 1715–1724. [Google Scholar] [CrossRef] [Green Version]
  27. van Bennekum, A.; Werder, M.; Thuahnai, S.T.; Han, C.H.; Duong, P.; Williams, D.L.; Wettstein, P.; Schulthess, G.; Phillips, M.C.; Hauser, H. Class B scavenger receptor-mediated intestinal absorption of dietary β-carotene and cholesterol. Biochemistry 2005, 44, 4517–4525. [Google Scholar] [CrossRef]
  28. Toomey, M.B.; Lopes, R.J.; Araujo, P.M.; Johnson, J.D.; Gazda, M.A.; Afonso, S.; Mota, P.G.; Koch, R.E.; Hill, G.E.; Corbo, J.C.; et al. High-density lipoprotein receptor SCARB1 is required for carotenoid coloration in birds. Proc. Natl. Acad. Sci. USA 2017, 114, 5219–5224. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Kiefer, C.; Sumser, E.; Wernet, M.F.; Von Lintig, J. A class B scavenger receptor mediates the cellular uptake of carotenoids in Drosophila. Proc. Natl. Acad. Sci. USA 2002, 99, 10581–10586. [Google Scholar] [CrossRef] [Green Version]
  30. Liu, H.; Zheng, H.; Zhang, H.; Deng, L.; Liu, W.; Wang, S.; Meng, F.; Wang, Y.; Guo, Z.; Li, S.; et al. A de novo transcriptome of the noble scallop, Chlamys nobilis, focusing on mining transcripts for carotenoid-based coloration. BMC Genom. 2015, 16, 44. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Tsuchida, K.; Arai, M.; Tanaka, Y.; Ishihara, R.; Ryan, R.O.; Maekawa, H. Lipid transfer particle catalyzes transfer of carotenoids between lipophorins of Bombyx mori. Insect Biochem. Mol. Biol. 1998, 28, 927–934. [Google Scholar] [CrossRef]
  32. During, A.; Dawson, H.D.; Harrison, E.H. 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, 2305–2312. [Google Scholar] [CrossRef]
  33. Ando, S.; Hatano, M. Bilirubin-binding protein in the serum of spawning-migrating chum salmon, Oncorhynchus keta: Its identity with carotenoid-carrying lipoprotein. Fish. Physiol. Biochem. 1988, 5, 69–78. [Google Scholar] [CrossRef]
  34. Eriksson, J.; Larson, G.; Gunnarsson, U.; Bed’hom, B.; Tixier-Boichard, M.; Stromstedt, L.; Wright, D.; Jungerius, A.; Vereijken, A.; Randi, E.; et al. Identification of the yellow skin gene reveals a hybrid origin of the domestic chicken. PLoS Genet. 2008, 4, e1000010. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Li, X.; Wang, S.; Xun, X.; Zhang, M.; Wang, S.; Li, H.; Zhao, L.; Fu, Q.; Wang, H.; Li, T.; et al. A carotenoid oxygenase is responsible for muscle coloration in scallop. Biochim. Biophys. Acta Mol. Cell Biol. Lipids 2019, 1864, 966–975. [Google Scholar] [CrossRef]
  36. Weesie, R.J.; Askin, D.; Jansen, F.J.; de Groot, H.J.; Lugtenburg, J.; Britton, G. Protein-chromophore interactions in α-crustacyanin, the major blue carotenoprotein from the carapace of the lobster, Homarus gammarus. A study by 13C magic angle spinning NMR. FEBS Lett. 1995, 362, 34–38. [Google Scholar] [CrossRef] [Green Version]
  37. Zagalsky, P.F. β-Crustacyanin, the blue-purple carotenoprotein of lobster carapace: Consideration of the bathochromic shift of the protein-bound astaxanthin. Acta Cryst. D Biol. Cryst. 2003, 59, 1529–1531. [Google Scholar] [CrossRef] [Green Version]
  38. Ertl, N.G.; Elizur, A.; Brooks, P.; Kuballa, A.V.; Anderson, T.A.; Knibb, W.R. Molecular characterisation of colour formation in the prawn Fenneropenaeus merguiensis. PLoS ONE 2013, 8, e56920. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Weaver, R.J.; Gonzalez, B.K.; Santos, S.R.; Havird, J.C. Red Coloration in an Anchialine Shrimp: Carotenoids, Genetic Variation, and Candidate Genes. Biol. Bull. 2020, 238, 119–130. [Google Scholar] [CrossRef] [PubMed]
  40. Jin, Y.; Yu, Y.; Zhang, C.; Li, S.; Zhang, X.; Li, F. Characterization and function analysis of the β-carotene oxygenase-like genes in carotenoids metabolism of the ridgetail white prawn Exopalaemon carinicauda. Front. Physiol. 2020, 11, 745. [Google Scholar] [CrossRef] [PubMed]
  41. Zhang, C.; Su, F.; Li, S.; Yu, Y.; Xiang, J.; Liu, J.; Li, F. Isolation and identification of the main carotenoid pigment from a new variety of the ridgetail white prawn Exopalaemon carinicauda. Food Chem. 2018, 269, 450–454. [Google Scholar] [CrossRef]
  42. Gao, Y.; Zhang, X.; Wei, J.; Sun, X.; Yuan, J.; Li, F.; Xiang, J. Whole Transcriptome Analysis Provides Insights into Molecular Mechanisms for Molting in Litopenaeus vannamei. PLoS ONE 2015, 10, e0144350. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Grabherr, M.G.; Haas, B.J.; Yassour, M.; Levin, J.Z.; Thompson, D.A.; Amit, I.; Adiconis, X.; Fan, L.; Raychowdhury, R.; Zeng, Q.; et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat. Biotechnol 2011, 29, 644–652. [Google Scholar] [CrossRef] [Green Version]
  44. Mortazavi, A.; Williams, B.A.; McCue, K.; Schaeffer, L.; Wold, B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat. Methods 2008, 5, 621–628. [Google Scholar] [CrossRef] [PubMed]
  45. Letunic, I.; Doerks, T.; Bork, P. SMART: Recent updates, new developments and status in 2015. Nucleic Acids Res. 2015, 43, D257–D260. [Google Scholar] [CrossRef]
  46. Avarre, J.C.; Lubzens, E.; Babin, P.J. Apolipocrustacein, formerly vitellogenin, is the major egg yolk precursor protein in decapod crustaceans and is homologous to insect apolipophorin II/I and vertebrate apolipoprotein B. BMC Evol. Biol. 2007, 7, 3. [Google Scholar] [CrossRef] [PubMed]
  47. Kumar, S.; Stecher, G.; Tamura, K. MEGA7: Molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol. Biol. Evol. 2016, 33, 1870–1874. [Google Scholar] [CrossRef] [Green Version]
  48. Zhang, W.; Sun, Z. Random local neighbor joining: A new method for reconstructing phylogenetic trees. Mol. Phylogenet. Evol. 2008, 47, 117–128. [Google Scholar] [CrossRef]
  49. Bustin, S.A.; Benes, V.; Garson, J.A.; Hellemans, J.; Huggett, J.; Kubista, M.; Mueller, R.; Nolan, T.; Pfaffl, M.W.; Shipley, G.L.; et al. The MIQE guidelines: Minimum information for publication of quantitative real-time PCR experiments. Clin. Chem. 2009, 55, 611–622. [Google Scholar] [CrossRef] [Green Version]
  50. Johnson, G.; Nour, A.A.; Nolan, T.; Huggett, J.; Bustin, S. Minimum information necessary for quantitative real-time PCR experiments. Methods Mol. Biol. 2014, 1160, 5–17. [Google Scholar]
  51. Livak, K.; Schmittgen, T. Analysis of Relative Gene Expression Data Using Real-Time Quantitative PCR and the 2−ΔΔCt Method. Methods 2000, 25, 402–408. [Google Scholar] [CrossRef]
  52. Chimsung, N.; Lall, S.P.; Tantikitti, C.; Verlhac-Trichet, V.; Milley, J.E. Effects of dietary cholesterol on astaxanthin transport in plasma of Atlantic salmon (Salmo salar). Comp. Biochem. Physiol. B Biochem. Mol. Biol. 2013, 165, 73–81. [Google Scholar] [CrossRef]
  53. Deming, D.M.; Erdman, J.W. Mammalian carotenoid absorption and metabolism. Pure Appl. Chem. 1999, 71, 2213–2223. [Google Scholar] [CrossRef] [Green Version]
  54. Borel, P. Factors affecting intestinal absorption of highly lipophilic food microconstituents (fat-soluble vitamins, carotenoids and phytosterols). Clin. Chem. Lab. Med. 2003, 41, 979–994. [Google Scholar] [CrossRef]
  55. Breithaupt, D.E.; Bamedi, A.; Wirt, U. Carotenol fatty acid esters: Easy substrates for digestive enzymes? Comp. Biochem. Physiol. B Biochem. Mol. Biol. 2002, 132, 721–728. [Google Scholar] [CrossRef]
  56. Bhosale, P.; Bernstein, P.S. Vertebrate and invertebrate carotenoid-binding proteins. Arch. Biochem. Biophys. 2007, 458, 121–127. [Google Scholar] [CrossRef] [Green Version]
  57. Parker, R.S. Carotenoids in human blood and tissues. J. Nutr. 1989, 119, 101–104. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Hata, T.R.; Scholz, T.A.; Ermakov, I.V.; McClane, R.W.; Khachik, F.; Gellermann, W.; Pershing, L.K. Non-invasive raman spectroscopic detection of carotenoids in human skin. J. Investig. Derm. 2000, 115, 441–448. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Rassart, E.; Bedirian, A.; Do Carmo, S.; Guinard, O.; Sirois, J.; Terrisse, L.; Milne, R. Apolipoprotein D. Biochim. Biophys. Acta 2000, 1482, 185–198. [Google Scholar] [CrossRef]
  60. Yang, F.; Wang, M.R.; Ma, Y.G.; Ma, W.M.; Yang, W.J. Prawn lipocalin: Characterization of a color shift induced by gene knockdown and ligand binding assay. J. Exp. Zool. A Ecol. Genet. Physiol. 2011, 315, 562–571. [Google Scholar] [CrossRef] [PubMed]
  61. Ando, S.; Takeyama, T.; Hatano, M. Transport Associated with Serum Vitellogenin of Carotenoid in Chum Salmon (Oncorhynchus keta). Agric. Biol. Chem. 1986, 50, 557–563. [Google Scholar] [CrossRef] [Green Version]
  62. Li, B.; Vachali, P.; Frederick, J.M.; Bernstein, P.S. Identification of StARD3 as a lutein-binding protein in the macula of the primate retina. Biochemistry 2011, 50, 2541–2549. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Shandala, T.; Lim, C.; Sorvina, A.; Brooks, D.A. A Drosophila model to image phagosome maturation. Cells 2013, 2, 188–201. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Trimble, W.S.; Coppolino, M.G. Phagosome Maturation. In Molecular Mechanisms of Phagocytosis; Rosales, C., Ed.; Springer: Boston, MA, USA, 2005; pp. 133–150. [Google Scholar]
  65. Xu, H.; Ren, D. Lysosomal physiology. Annu. Rev. Physiol. 2015, 77, 57–80. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Cheesman, D.; Lee, W.; Zagalsky, P. Carotenoproteins in Invertebrates. Biol. Rev. 2008, 42, 131–160. [Google Scholar] [CrossRef]
  67. Chayen, N.E.; Cianci, M.; Grossmann, J.G.; Habash, J.; Helliwell, J.R.; Nneji, G.A.; Raftery, J.; Rizkallah, P.J.; Zagalsky, P.F. Unravelling the structural chemistry of the colouration mechanism in lobster shell (vol D59, pg 2072, 2003). Acta Cryst. D 2005, 61, 1312. [Google Scholar] [CrossRef] [Green Version]
  68. Cianci, M.; Rizkallah, P.J.; Olczak, A.; Raftery, J.; Chayen, N.E.; Zagalsky, P.F.; Helliwell, J.R. The molecular basis of the coloration mechanism in lobster shell: β-crustacyanin at 3.2-angstrom resolution. Proc. Natl. Acad. Sci. USA 2002, 99, 9795–9800. [Google Scholar] [CrossRef] [Green Version]
  69. Kuballa, A.V.; Merritt, D.J.; Elizur, A. Gene expression profiling of cuticular proteins across the moult cycle of the crab Portunus pelagicus. BMC Biol. 2007, 5, 45. [Google Scholar] [CrossRef] [Green Version]
  70. Xiong, G.; Tong, X.; Yan, Z.; Hu, H.; Duan, X.; Li, C.; Han, M.; Lu, C.; Dai, F. Cuticular protein defective Bamboo mutant of Bombyx mori is sensitive to environmental stresses. Pestic. Biochem. Physiol. 2018, 148, 111–115. [Google Scholar] [CrossRef]
  71. Asano, T.; Taoka, M.; Shinkawa, T.; Yamauchi, Y.; Isobe, T.; Sato, D. Identification of a cuticle protein with unique repeated motifs in the silkworm, Bombyx mori. Insect Biochem. Mol. Biol. 2013, 43, 344–351. [Google Scholar] [CrossRef]
  72. Xiong, G.; Tong, X.; Gai, T.; Li, C.; Qiao, L.; Monteiro, A.; Hu, H.; Han, M.; Ding, X.; Wu, S.; et al. Body Shape and Coloration of Silkworm Larvae Are Influenced by a Novel Cuticular Protein. Genetics 2017, 207, 1053–1066. [Google Scholar] [CrossRef] [Green Version]
Genes 12 00618 g001 550
Figure 1. Expression profiles of nine selected unigenes between EcR and EcW prawns. The x-axis represents the names of selected unigenes. Columns represent the RPKM value from the transcriptome result (A) and the means of relative expression levels from the qPCR result (B). The significant difference in the expression levels between the EcR group and the EcW group was labeled with asterisks (*) at p < 0.05.
Figure 1. Expression profiles of nine selected unigenes between EcR and EcW prawns. The x-axis represents the names of selected unigenes. Columns represent the RPKM value from the transcriptome result (A) and the means of relative expression levels from the qPCR result (B). The significant difference in the expression levels between the EcR group and the EcW group was labeled with asterisks (*) at p < 0.05.
Genes 12 00618 g001
Genes 12 00618 g002 550
Figure 2. GO annotation of DEGs. DEGs with a GO annotation were divided into three major categories: biological process, cellular component, and molecular function. Red columns represented the DEGs with higher expression levels in the EcR group, while green columns represented the DEGs with lower expression levels in the EcR group.
Figure 2. GO annotation of DEGs. DEGs with a GO annotation were divided into three major categories: biological process, cellular component, and molecular function. Red columns represented the DEGs with higher expression levels in the EcR group, while green columns represented the DEGs with lower expression levels in the EcR group.
Genes 12 00618 g002
Genes 12 00618 g003 550
Figure 3. Top 20 enriched pathways (A) and pathways related to metabolism (B). The circle size and filled portions represented the number of differentially expressed genes in each pathway. The statistical significance was colored according to Q-values.
Figure 3. Top 20 enriched pathways (A) and pathways related to metabolism (B). The circle size and filled portions represented the number of differentially expressed genes in each pathway. The statistical significance was colored according to Q-values.
Genes 12 00618 g003
Genes 12 00618 g004 550
Figure 4. Alignment of crustacyanin amino acid sequences. Unigene0054180-EcR and Unigene0054180-EcW represented the crustacyanin in the orange-red-variant prawns and the wild-type prawns. A1 and A3 represented the crustacyanin subunit proteins in the lobster. Amino acid sequences of lipocalin consensus regions were colored, and different amino acid residues between the orange-red-variant prawns and the wild-type prawns were highlighted in the red box. Asterisks (*) indicated positions which had a single, fully conserved residue, while one black dot (.) and two black dots (:) indicated that the residue sites were less conserved but their physicochemical properties were similar or very similar.
Figure 4. Alignment of crustacyanin amino acid sequences. Unigene0054180-EcR and Unigene0054180-EcW represented the crustacyanin in the orange-red-variant prawns and the wild-type prawns. A1 and A3 represented the crustacyanin subunit proteins in the lobster. Amino acid sequences of lipocalin consensus regions were colored, and different amino acid residues between the orange-red-variant prawns and the wild-type prawns were highlighted in the red box. Asterisks (*) indicated positions which had a single, fully conserved residue, while one black dot (.) and two black dots (:) indicated that the residue sites were less conserved but their physicochemical properties were similar or very similar.
Genes 12 00618 g004
Genes 12 00618 g005 550
Figure 5. Evolutionary relationships of lipoproteins from E. carinicauda and other species. All DEGs with higher expression levels in the orange-red-variant prawns were marked with red font in the tree, while DEGs with lower expression levels were marked with green font. apoCr, apolipocrustacein; MTP, microsomal triglyceride transfer protein; apoB, vertebrate apolipoprotein B; Vtg, vitellogenin; apoLp-II/I, insect apolipophorin II/I precursor. Accession numbers for these sequences were indicated in Table S3.
Figure 5. Evolutionary relationships of lipoproteins from E. carinicauda and other species. All DEGs with higher expression levels in the orange-red-variant prawns were marked with red font in the tree, while DEGs with lower expression levels were marked with green font. apoCr, apolipocrustacein; MTP, microsomal triglyceride transfer protein; apoB, vertebrate apolipoprotein B; Vtg, vitellogenin; apoLp-II/I, insect apolipophorin II/I precursor. Accession numbers for these sequences were indicated in Table S3.
Genes 12 00618 g005
Table
Table 1. Information of primers used for real-time PCR.
Table 1. Information of primers used for real-time PCR.
Gene IDDescriptionSequence (5′ to 3′)Tm (°C)
Unigene0059050Esterase FE4F: GTTGGACCAAATTTCGGCTCTT
R: TGGAGAGAGCGCATGAAAGTTA		59
Unigene0054180Crustacyanin subunit AF: TCCTGAGGGAACGCACGAGA
R: TGGAGAAGACGAAGGCGAACT		58
Unigene0036561Apolipoprotein DF: CTGGGCACAGATTACGAGAACT
R: GGTGGGTACGATGGAAGAGG 		64
Unigene0037069Cuticular protein 34F: CAGAAGCCACGAGCCAGAAG
R: GGGTCCTACGAGGATGATGTTAT		57
Unigene0038920Cuticle proproteinF: CGTAACTGGCACCTGGATGA
R: AGAGGATAGCGTGGGCTGGA		59
Unigene0043447Cuticular protein 34F: TCGTCAGAGCCGATGGAAAC
R: TGGGAGTGCCATCATCAGTG		58
Unigene0030347NoneF: CTATGCGAATGAATAAGATGAGGAG
R: CCAGCGTACCAAGTAATACTGAAA		57
Unigene0033842Retrovirus-related Pol polyproteinF: CGATTCACTGTCCCCACTACTC
R: GTCTATTTCCTTGATGCTCTTACCA		56
Unigene0047842Lathosterol oxidaseF: AGTGTAAGAATCCACCAATGCC
R: GCCAAGTTTCCAATGAGACAGC		58
18S rRNA18S rRNAF: TATACGCTAGTGGAGCTGGAA
R: GGGGAGGTAGTGACGAAAAAT		55
Note: F—forward primer; R—reverse primer.
Table
Table 2. Summary of sequencing and assembly of the transcriptome.
Table 2. Summary of sequencing and assembly of the transcriptome.
SampleRaw ReadsClean ReadsMapping RatioUnigene NumberRatio
EcW-1765582247333558464071854 (87.71%)73,57694.06%
EcW-2754078227278844863707819 (87.67%)74,17494.82%
EcW-3782696467514104466104983 (88.14%)74,00694.61%
EcR-1840683728083277871088041 (88.17%)73,69794.21%
EcR-2720767626970033860751470 (87.70%)73,49793.96%
EcR-3709242486859755260474700 (88.34%)72,32092.45%
EcW230235692221265076193884656 (87.84%)76,24797.47%
EcR227069382219130668192314211 (88.07%)75,94797.09%
All457305074440395744387591019 (88.01%)77,07998.54%
Note: EcW: the wild prawn group, consisting of three replicates (EcW-1, EcW-2, and EcW-3). EcR: the orange-red prawn group, consisting of three replicates (EcR-1, EcR-2, and EcR-3). Mapping Ratio = Mapped Reads number/the rRNA-removed high-quality clean reads. Ratio = gene number expressed in each sample/all reference gene numbers (78,224).
Table
Table 3. Information of cuticular proteins in DEGs.
Table 3. Information of cuticular proteins in DEGs.
Genelog2 Ratio (EcR/EcW)DescriptionSpecies
Unigene0039421−2.98Calcification-associated peptide−2Procambarus clarkii
Unigene0037661−2.56Calcification-associated peptide−2P. clarkii
Unigene0019577−6.61Calcified cuticle protein CP14.1 Callinectes sapidus
Unigene0010114−3.86Calcified cuticle protein CP19.0 isoform AC. sapidus
Unigene0072453−6.58Calcified cuticle protein CP19.0 isoform BC. sapidus
Unigene0042544−5.85Calcified cuticle protein CP19.0 isoform BC. sapidus
Unigene0038920−3.89Cuticle proprotein, partial Palaemon varians
Unigene0032773−5.13Cuticle protein 1Cherax quadricarinatus
Unigene0056168−4.36Cuticle protein 18.6, isoform B Lepeophtheirus salmonis
Unigene0016948−4.23Cuticle protein BD2 Portunus pelagicus
Unigene0035595−2.25Cuticle protein, partial Daphnia magna
Unigene0003529−1.74Cuticular proteinTenebrio molitor
Unigene0073146−9.77Cuticular protein 34 Eriocheir sinensis
Unigene0035142−5.96Cuticular protein 34E. sinensis
Unigene0037069−5.30Cuticular protein 34E. sinensis
Unigene0017103−3.76Cuticular protein 34E. sinensis
Unigene0043447−2.56Cuticular protein 34E. sinensis
Unigene0061823−1.05Cuticular protein analogous to peritrophins 1-F, partialD. magna
Unigene0019752−1.54Cuticular protein RR-2 motif 78 precursorBombyx mori
Unigene0029471−2.05DD5Marsupenaeus japonicus
Unigene0019787−1.97DD5M. japonicus
Unigene0029470−1.69DD5M. japonicus
Unigene0044640−2.05DD9A, partialM. japonicus
Unigene0001656−1.55Obstructor BLocusta migratoria
Unigene0050714−1.99Post-molt protein 1C. quadricarinatus
Unigene0022558−1.80Strongly chitin-binding protein-1P. clarkii
Unigene0046800−1.19Pupal cuticle proteinCephus cinctus
Unigene0008712−1.12Pupal cuticle proteinC. cinctus
Unigene0023218−1.42Cuticle protein 6 Blaberus craniifer
Unigene0064526−1.29Cuticle protein 6 B. craniifer
Unigene0025406−3.35Cuticle protein CP1876Cancer pagurus
Table
Table 4. Transport- and binding-related DEGs.
Table 4. Transport- and binding-related DEGs.
Functional CategoriesGeneIDlog2 Ratio (EcR/EcW)DescriptionSpecies
TransportUnigene00655501.57Vitellogenin 2Pandalopsis japonica
Unigene00637523.65VitellogeninPandalus hypsinotus
Unigene00167791.04VitellogeninE. carinicauda
Unigene00025421.52StAR-related lipid transfer protein 3Penaeus monodon
Unigene00288103.54Apolipoprotein DZootermopsis nevadensis
Unigene0036561−2.46Apolipoprotein DD. magna
Unigene0044620−1.26Apolipoprotein DD. magna
BindingUnigene00299182.56TubulinTetrahymena thermophila
Unigene00734212.13Tubulin/FtsZ familyT. thermophila
Unigene00681541.79Tubulin α-1 chainTrichinella murrelli
Unigene00137641.07Tubulin α-1 chain-likeHalyomorpha halys
Unigene00260121.92Fibronectin-likeCiona intestinalis
Unigene00220241.53Fibrillin-1Exaiptasia pallida
Unigene00365971.06HemocyaninE. carinicauda
Unigene00082141.15HemocyaninE. carinicauda
Unigene0027904−1.80HemocyaninMacrobrachium nipponense
Unigene0054180−1.85Crustacyanin subunit AF. merguiensis
Table
Table 5. Differentially expressed genes in Lysosome pathway.
Table 5. Differentially expressed genes in Lysosome pathway.
GeneIDlog2 Ratio (EcR/EcW)DescriptionSpecies
Unigene00334041.02Cathepsin BPandalus borealis
Unigene0049424−1.25CrustapainPandalus borealis
Unigene0018100−2.52Cathepsin L2Litopenaeus vannamei
Unigene00245152.31Cathepsin LD. magna
Unigene0066106−3.44Heparan-α-glucosaminide N-acetyltransferaseD. magna
Unigene0022872−1.11N-acetyl-β-d-glucosaminidaseM. nipponense
Unigene0017001−2.85Ecdysteroid regulated-like proteinLitopenaeus vannamei
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Jin, Y.; Li, S.; Yu, Y.; Zhang, C.; Zhang, X.; Li, F. Transcriptome Analysis Provides Insights into the Mechanism of Astaxanthin Enrichment in a Mutant of the Ridgetail White Prawn Exopalaemon carinicauda. Genes 2021, 12, 618. https://doi.org/10.3390/genes12050618

AMA Style

Jin Y, Li S, Yu Y, Zhang C, Zhang X, Li F. Transcriptome Analysis Provides Insights into the Mechanism of Astaxanthin Enrichment in a Mutant of the Ridgetail White Prawn Exopalaemon carinicauda. Genes. 2021; 12(5):618. https://doi.org/10.3390/genes12050618

Chicago/Turabian Style

Jin, Yue, Shihao Li, Yang Yu, Chengsong Zhang, Xiaojun Zhang, and Fuhua Li. 2021. "Transcriptome Analysis Provides Insights into the Mechanism of Astaxanthin Enrichment in a Mutant of the Ridgetail White Prawn Exopalaemon carinicauda" Genes 12, no. 5: 618. https://doi.org/10.3390/genes12050618

APA Style

Jin, Y., Li, S., Yu, Y., Zhang, C., Zhang, X., & Li, F. (2021). Transcriptome Analysis Provides Insights into the Mechanism of Astaxanthin Enrichment in a Mutant of the Ridgetail White Prawn Exopalaemon carinicauda. Genes, 12(5), 618. https://doi.org/10.3390/genes12050618

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