Animals
This study was performed on “Bohai Red” scallop, which is a new bay scallop strain selected from the hybrids between the bay scallop and Peruvian scallop. It is widely cultured in northern China(C et al., 2017). The scallops were conditioned to breed and mature in the scallop hatchery of Yantai Spring-Sea AquaSeed Co., Ltd in Laizhou, Shandong, and were cultured in Yangma Island, Mouping, Shandong.
Transcriptome sequencing and analysis
In total, six scallops at the same stage of development were dissected(three black mantle and three white mantle), and their mantle was collected, with using normal-colouredcolored (white) mantle tissuesscallop samples as controls, the transcriptome of black(treated group) and white(control group) mantle tissues was compared. Total RNA was extracted from the tissue(TIANGEN, Beijing, China) and quantified using Nanodrop 2000. The RNA integrity was confirmed using 1% agarose gel electrophoresis. Library construction and high-throughput sequencing were performed using Mouse Genome Informatics (MGI).
Raw data obtained from the sequencing were filtered by removing adapters, unknown nucleotides (> 10%), and low-quality reads (Q-value ≤ 20) and bases (> 50%) using SOAPnuke(Chen et al., 2018). Raw reads were mapped to the scallop genome (BioSample accessions: SAMN08322131(Liu et al., 2020)) using HISAT2(Kim et al., 2015) and Bowtie2(Langmead, 2010). The results obtained from Bowtie2 were statically analyzed using RSEM(Dewey and Bo, 2011), and the gene expression was determined using fragments per kilobase per million bases(Trapnell et al., 2010).
Differentially expressed genes (DEGs) were identified using DESeq2(Love et al., 2014) with |log2fold change (FC)| > 1 and false discovery rate (FDR) < 0.05. Further, the DEGs were subjected to Gene Ontology (GO)(Ashburner et al., 2000) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses(Kanehisa et al., 2004, Kanehisa, 2002, Chen et al., 2011).
Metabolomic sequencing and analysis
To explore the differentially expressed metabolites, six scallops at the same stage of development were dissected(three black mantle and three white mantle) to collect their mantle which will divided into two groups(black and white). The mantle sample put into a centrifuge tube. To this, 400 µL methanol and some steel beads were added and vortexed for 3 min. If the sample was not completely broken, more steel beads were added and vortexed for 3 min. Further, the sample was sonicated for 10 min on an ice-water bath, followed by incubation at − 20°C for 30 min. Further, the sample was centrifuged at 12,000 rpm and 4°C for 10 min. The supernatants were harvested for subsequent analysis.
Data transformation raw files were converted into mzML format using Proteowizard. Peak extraction, alignment, and retention time correction were performed using the XCMS program. Subsequently, the identified metabolite information was obtained by searching the laboratory’s self-built database, the public database (Metlin), and metDNA. Statistical analysis was performed using R statistical program. Based on OPLS-DA analysis, differentially expressed metabolites were initially screened. Further screening was performed based on p-value < 0.05 and FC ≥ 2 and ≤ 0.5. The differentially expressed metabolites were functionally annotated using KEGG pathway enrichment analysis(Kanehisa et al., 2004).
Estimation of melanin
The black component and its amount in the mantle of scallops were assessed using Shellfish Melanin ELISA Test Kit (Shanghai Yuanjie Biotechnology Centre, Shanghai, China). The mantle obtained from scallops was divided into two groups: black mantle and white mantle, each comprising three replicates. The mantle tissue was homogenized with a high speed tissue grinder(OSE-Y30,TIANGEN,Beijing) in an appropriate volume of physiological saline solution and centrifuged at 3000 rpm for 10 min. The supernatant was collected for subsequent experiments.
ELISA was performed using standard melanin samples at 0, 3, 6, 12, 24, and 48 mg/mL, along with samples and blank wells(0 mg/mL standard melanin samples) as controls, using a 96-well plate. In each standard and sample well, 100 µL of horseradish peroxidase (HRP)-labeled detection antibody was added, and the wells were sealed with sealing membranes and incubated at 37°C for 60 min. After removing the liquid from the wells, the plates were tapped on absorbent paper. To each well, wash buffer was added, incubated for 1 min, and removed by flicking. The plates were tapped dry on an absorbent paper. This washing process was repeated five times. Subsequently, in each well, 50 µL of substrates A and B were added and incubated at 37°C for 15 min in dark. Finally, 50 µL of stop solution was added to each well, and the OD values were determined at 450 nm within 15 min. A standard calibration curve was plotted with the concentrations of the standard samples as the abscissa and their corresponding OD values as the ordinate. Linear regression analysis was performed, and the concentrations of the respective samples were calculated based on the equation derived from the calibration curve.