Skip to main content

Genome-wide DNA methylation profile changes associated with shell colouration in the Yesso scallop (Patinopecten yessoensis) as measured by whole-genome bisulfite sequencing

Abstract

Background

Mollusca, a phylum of highly rich species, possess vivid shell colours, but the underlying molecular mechanism remains to be elucidated. DNA methylation, one of the most common epigenetic modifications in eukaryotes, is believed to play a vital role in various biological processes. However, analysis of the effects of DNA methylation on shell colouration has rarely been performed in molluscs, limiting the current knowledge of the molecular mechanism of shell colour formation.

Results

In the present study, to reveal the role of epigenetic regulation in shell colouration, WGBS, the “gold standard” of DNA methylation analysis, was first performed on the mantle tissues of Yesso scallops (Patinopecten yessoensis) with different shell colours (brown and white), and DNA methylomes at single-base resolution were generated. About 3% of cytosines were methylated in the genome of the Yesso scallop. A slight increase in mCG percentage and methylation level was found in brown scallops. Sequence preference of nearby methylated cytosines differed between high and low methylation level sites and between the brown- and white-shelled scallops. DNA methylation levels varied among the different genomic regions; all the detected regions in the brown group exhibited higher methylation levels than the white group. A total of 41,175 DMRs (differentially methylated regions) were detected between brown and white scallops. GO functions and pathways associated with the biosynthesis of melanin and porphyrins were significantly enriched for DMRs, among which several key shell colour-related genes were identified. Further, different correlations between mRNA expression levels and DNA methylation status were found in these genes, suggesting that DNA methylation regulates shell colouration in the Yesso scallop.

Conclusions

This study provides genome-wide DNA methylation landscapes of Yesso scallops with different shell colours, offering new insights into the epigenetic regulatory mechanism underlying shell colour.

Peer Review reports

Background

Colour is a common and important trait in many organisms, and is often associated with camouflage, aposematism, mate attraction, thermoregulation and immunity [13]. Mollusca is the second-largest phylum of animals on the planet, comprising highly rich species. Molluscs are known for their colourful shells, which have attracted collectors and scientists for hundreds of years. Compared with plants and vertebrates, there is limited knowledge of the molecular mechanism underlying shell colouration in molluscs due to its complexity. Generally, the shell colour is related to the presence of shell pigments, including melanin, carotenoids, porphyrins and bile pigments, which have been identified in some molluscs [4]. Though environmental conditions and food intake can sometimes affect shell pigmentation [5, 6], many breeding studies have demonstrated that the shell colour trait is under genetic control [710], and some shell colour-related genetic markers have been detected [1114]. Benefiting from next-generation sequencing, extensive transcriptome analysis has been performed in molluscan species, and some pathways (such as melanogenesis, the heme biosynthesis pathway and the porphyrin and chlorophyll metabolism pathway) associated with shell colour have been identified [1521]. However, the underlying regulation mechanism remains to be elucidated.

Epigenetic modification is an important gene regulatory system affecting gene expression without altering nucleotide sequences [2224]. DNA methylation is one of the most common epigenetic modifications in eukaryotes and is believed to play a vital role in various biological processes, including transcription regulation, embryogenesis, cellular differentiation, genomic imprinting, X-chromosome inactivation, transposon silencing and so on [2531]. In plants and some vertebrates, DNA methylation is also involved in the regulation of colouration, such as flower colours [32], eggshell colours [33] and animal skin and coat colours [34]; this suggests that it is also likely involved in the regulation of shell colour. However, to date, no published studies have addressed this hypothesis.

Though limited, an increasing number of DNA methylation analyses have been performed in molluscs due to the important regulation role of DNA methylation; however, most have focused on oysters. For example, global DNA methylation analysis by whole-genome bisulfite sequencing (WGBS) in the Pacific oyster Crassostrea gigas indicated that DNA methylation plays a complex role in regulating genome activity and is involved in the regulation of development processes and stress and environmental responses [3538]. In the pearl oyster Pinctada fucata martensii, DNA methylomes of biomineralization-related tissues were analysed by methylated DNA immunoprecipitation sequencing (MeDIP-Seq) to explore the regulatory mechanism underlying shell formation [39]. In the Yesso scallop Patinopecten yessoensis, dynamic DNA methylation levels were detected during gametogenesis and early development [40]. Further, different DNA methylation levels have been detected in Yesso scallops with different adductor muscle colours using the MethylRAD technique [41]. However, there are few published analyses of the effects of DNA methylation on shell colouration in molluscs, limiting the current knowledge of the molecular mechanism of shell colour formation.

The Yesso scallop (P. yessoensis) is a large and old (dating back to ~ 350 Ma) bivalve animal that can be found living along the coastlines of northern Japan, the Russian Far East and northern Korea. It is one of the most important aquaculture shellfish in Asian countries and is consumed worldwide. The shell colour of the Yesso scallop represents a polymorphism whereby most individuals possess brown left valves while a small percentage of individuals have white left valves. Unlike some other molluscs with complex shell colours, the shell colour pattern of the Yesso scallop is unambiguous and relatively simple, making it a good material for shell colour study. Due to the important evolutionary and economic position of the Yesso scallop, much research has been conducted on this species, including genome and transcriptome sequencing [42], epigenetics analyses [14, 40] and genetic breeding for shell colour [41]; these studies provide valuable genetic resources for research on the biological problems of the Yesso scallop. WGBS, which is considered the “gold standard” of DNA methylation analysis, is the most effective method to study genomic methylation patterns; it can generate methylomes at single-base resolution without a preference [43]. In the present study, WGBS analysis was first performed on the mantle tissues of Yesso scallops with two different shell colours (brown and white) in order to generate DNA methylomes at single-base resolution for this species. The general characteristics of DNA methylation in the two shell colour Yesso scallops were analysed and DNA methylation changes between the two kinds of scallops were examined. The aim was to explore the impact of DNA methylation on shell colouration and provide new genome-wide evidence for understanding the epigenetic regulatory mechanism underlying the shell colour trait.

Results

Genome-wide DNA methylation landscape in the Yesso scallop

To investigate the overall methylation patterns of the Yesso scallop and detect those sites associated with shell pigmentation, WGBS was performed on the mantle tissues of two different shell colour populations (brown and white) of Yesso scallops. A total of 147.54 Gb and 150.25 Gb raw bases were generated for the two groups, respectively, with an average of 49.18 Gb and 50.08 Gb, respectively, for each replicate, representing > 30 × the reference genome [44]. Detailed description of the sequencing data is provided in Table 1. The raw data have been deposited into the NCBI SRA database with the accession number PRJNA695315.

Table 1 Data summary of BS-Seq for the two shell colour Yesso scallops; data represents three replicates

The distributions of the DNA methylation sites in the whole genome of the Yesso scallop were separately analysed in the brown and white individuals (Fig. 1a, Table S1). Overall, approximately 3.07 and 2.94% of all genomic C sites were methylated in the brown and white groups, respectively. Methylation was found to exist in three sequence contexts: CG, CHG and CHH (where H represents A, C or T). The genome of the brown group exhibited 24.35% (mCG), 0.12% (mCHG) and 0.09% (mCHH) methylation on the total sequenced CG, CHG and CHH sites, respectively. Accordingly, the white group exhibited 22.71, 0.12 and 0.09% methylation in the CG, CHG and CHH contexts, respectively, with a slight decrease in CG methylation percentage compared with the brown group. The proportions of mCG, mCHG and mCHH on the total mC sites were analysed (Fig. 1b, Table S1), which reflects the distributions of mCs in the three sequence contexts. Similar rates were found in the brown and white groups, with methylcytosine most often found at CG sites (~ 97%) and much less frequently at CHG (~ 1%) and CHH (~ 2%) sites.

Fig. 1
figure 1

DNA methylation landscape of Yesso scallops with different shell colours. A Density plot of 5-methylvytosine in mCG, mCHG and mCHH sequence contexts. The first circle (outermost) represents the 30 longest scaffolds in the genome of the Yesso scallop; circles 2–5 represent the ratios of each sequence context (i.e., mC/covered C, mCG/covered CG, mCHG/covered CHG, mCHH/covered CHH); circles 6–9 represent the average methylation levels of each sequence context. B The proportion of mC in the three sequence contexts

On the whole, the methylation levels of the two shell colour groups shared similar trends, with high levels in CG sites and low levels in CHG and CHH sites (Fig. 2, Table S1). The average methylation levels in the genome of the brown group were 2.34, 17.96, 0.18 and 0.16% for all detected C, CG, CHG and CHH sites, and 71.33, 72.24, 40.61 and 39.98% for the mC, mCG, mCHG and mCHH sites, respectively. Moreover, the average methylation levels in the genome of the white group were 2.17, 16.39, 0.15 and 0.12% for all detected C, CG, CHG and CHH sites, and 69.94, 70.86, 39.49 and 39.12% for the mC, mCG, mCHG and mCHH sites, respectively. Overall, all types of sites in the white group exhibited slightly decreased methylation levels compared with those of the brown group.

Fig. 2
figure 2

The distribution of DNA methylation levels of all mC sites and different sequence contexts in Yesso scallops with two different shell colours

Sequence preferences analysis for DNA methylation

To explore the sequence preference of methylated sites, Logo Plots [45] were utilized to analyse the sequence information of each methylated site and its nearby regions (a total of 9-mer sequences containing a methylated site). For the CG context, sites with a methylation level > 75% were defined as high methylation level sites and the others were defined as low methylation level sites. For the CHG and CHH contexts, sites with a methylation level > 25% were defined as high methylation level sites and the others were defined as low methylation level sites. The sequence preference around methylated C sites varied among the three types of sequence contexts, and obvious differences were found between the high and low methylation level sites; differences also existed between the two shell colour groups (Supplementary Figure S1). In the CG context, the bases around the high methylation level sites were distributed evenly, while for the low methylation level sites, mC frequently occurred at the ACGA sequence for the brown group and ACGT sequence for the white group (Figure S1A). In the CHG context, mC in the high methylation level sites often appeared at the ACAGA sequence for the brown group and TCAGA sequence for the white group, while mC in the low methylation level sites appeared at the ACTGA and ACTGT sequences for the brown and white groups, respectively (Figure S1B). In the CHH context, mC frequently occurred at the CAT sequence regardless of the high or low methylation level status of the region (Figure S1C); however, differences in the composition of the nearby bases in the CHH context were also observed between the brown and white groups, e.g., cytosines preferred to be methylated at the CATA sequence for all the methylated sites in the brown group, while preference was for the CATT sequence in the white group.

DNA methylation pattern in different genomic regions

To study the DNA methylation profiles of different genomic regions of the Yesso scallop and explore the differences between the brown and white groups, the DNA methylation levels of different genomic regions were compared for each sequence context of the two kinds of scallops. As shown in Fig. 3, DNA methylation levels varied among the different genomic regions; the CG context showed the highest methylation level among the three sequence contexts in each genomic region, while the CHH context showed the lowest. The methylation level of the CG context in both of the shell colour groups was highest in the CDS, followed by the exon, and was lowest in the intergenic regions (Fig. 3a). Similar patterns were found in the other two contexts, but the differences between the different genomic regions were not as obvious as those in the CG context (Fig. 3b and c). The methylation levels within transcribed regions, including the 5’UTR, exon, intron and 3’UTR, and 2 kb upstream (i.e., promoter region) and downstream regions were also evaluated (Fig. 3d, f). In the CG context, the methylation level in the promoter region was the lowest and gradually declined toward the 5’UTR, then rose high in the 5’UTR. 5’UTR and 3’UTR had almost the same methylation levels, with methylation levels slightly higher than the intron regions. The exon regions exhibited the highest methylation level. The methylation level dropped sharply at the start of the 2-kb downstream regions and then rose back (Fig. 3d). A similar methylation pattern was also found in the CHG and CHH contexts, but the average methylation levels within the transcribed regions were comparable (Fig. 3e and f). Moreover, similar methylation patterns were observed in the two shell colour scallops, but all the detected regions in the brown group had higher methylation levels than those in the white group. The differences in the methylation levels of cytosines between the two kinds of scallops may be correlated with shell colouration.

Fig. 3
figure 3

DNA methylation levels in different genomic regions of Yesso scallops with different shell colours. A-C Distributions of DNA methylation levels on different genome elements separately for the CG, CHG and CHH contexts. D-F Distributions of DNA methylation levels among gene features separately for the CG, CHG and CHH contexts

DMR analysis between the brown and white Yesso scallops

Differentially methylated regions (DMRs) were examined between the brown and white groups to assess the role of methylation in colouration. A total of 41,175 DMRs were identified (Supplementary Table S2). To analyse the function of these DMRs, GO (Gene Ontology term) and KEGG (Kyoto Encyclopaedia of Genes and Genomes) enrichment analyses were separately performed on the DMR-associated genes. A total of 1132 GO terms were significantly enriched (FDR <  0.05, gene number > 3), covering a variety of functions (Supplementary Table S3). Many of the enriched GO terms were associated with histone modification, such as histone phosphorylation, acetylation, deacetylation, methylation, demethylation and so on; chromosome and DNA structure changes, such as chromatin remodelling, binding and silencing, DNA helicase activity, duplex unwinding and so on; and transcription regulation, such as different transcription factors activities (Table 2). More importantly, several pigmentation-related GO terms were enriched, i.e., pigmentation, melanocyte differentiation, melanosome organization, melanosome transport, endosome to melanosome transport, porphyrin-containing compound biosynthetic process and heme biosynthetic process, which are closely related to melanin and porphyrins formation. KEGG enrichment analysis showed that 39 pathways were significantly enriched at a threshold of FDR <  0.05 and gene number > 3; these pathways were related to RNA biosynthesis, degradation and transport, fatty acid biosynthesis, degradation and metabolism, glycan biosynthesis, carbon metabolism and so on (Fig. 4). It should be noted that the pathway of ‘porphyrin and chlorophyll metabolism’ (KO00860), which was closely associated with shell colour formation, was also significantly enriched with a threshold of p <  0.05 (Supplementary Table S4).

Table 2 GO enrichment analysis of DMRs
Fig. 4
figure 4

KEGG pathway enrichment analysis of DMRs (ranked by FDR; FDR < 0.05; KEGG pathways got from: https://www.kegg.jp/kegg/kegg1.html)

Genes involved in the biosynthesis pathway of melanin and porphyrins were identified among DMRs, and are listed in Table 3. MITF (microphthalmia-associated transcription factor) is the primary regulator of melanin synthesis. There were three DMRs detected in the gene of MITF, and they all exhibited significantly higher methylation levels in the white group than the brown group. Gene expression analysis for the MITF gene showed a higher expression level in the brown group than that in the white group, though this difference was not statistically significant (p ≥ 0.05) (Fig. 5). The genes ALAS (5-aminolevulinate synthase), ALAD (delta-aminolevulinic acid dehydratase), UROS (uroporphyrinogen-III synthase), UROD (uroporphyrinogen decarboxylase), CPOX (coproporphyrinogen-III oxidase), PPOX (protoporphyrinogen oxidase) and FECH (ferrochelatase), all from the heme biosynthesis pathway, are involved in the large pathway porphyrin and chlorophyll metabolism. Seven of eight genes in the heme biosynthesis pathway were differentially methylated between the brown and white groups. The methylation levels of all the DMRs in the genes ALAS, UROS and PPOX were higher in the brown group than in the white group; while the methylation levels of DMRs in the ALAD and CPOX genes were lower in the brown group than the white group; and different methylation sites in the UROD and FECH genes showed different methylation level changes between the two shell colour groups. Gene expression analysis for the above seven genes showed that aside from CPOX, which had a significantly higher expression level in the white group (p < 0.05), the other six genes had significantly higher expression levels in the brown group than the white (p < 0.05) (Fig. 5).

Table 3 DMRs involved in genes associated with shell colouration
Fig. 5
figure 5

Relative expression levels of shell colour-related genes in the mantle tissues of Yesso scallops with different shell colours. Vertical bars represent the mean ± S.E. (N = 3). Asterisks indicate significant differences (p < 0.05)

Discussion

DNA methylation is one of the most important epigenetic modifications in eukaryotes, but little is known about its role in shell colour formation. In the present study, genome-wide DNA methylation analysis was first performed on the mantle tissues of Yesso scallops with different shell colours using WGBS. The aim of this study was to explore the potential epigenetic regulatory mechanism of the shell colour trait.

DNA methylation is evolutionarily conserved, but markedly different methylation patterns exist among kingdoms [46]. Generally, vertebrates exhibit high DNA methylation levels while the genomes of invertebrates are relatively less methylated [47]. For example, approximately 4–6% of genomic cytosines are methylated in humans [48, 49], and a similar ratio was also found in a marine organism sea squirt from Urochordata (4.07%) [46]; however, quite a small percentage of cytosines (~ 0.1–0.2%) are methylated in the whole genomes of some insects [50, 51, 52]. In the limited available studies, different methylation levels have been found in marine invertebrates; for example, ~ 1.4% of genomic cytosines are methylated in the sea anemone [47], while ~ 3.5% were detected in the sea cucumber [53]. Studies of molluscs have primarily focused on oysters. In the present study, about 3% of genomic cytosines were found to be methylated in the Yesso scallop, slightly higher than that in the Pacific oyster (~ 1.95%) [38] and the Eastern oyster (~ 2.7%) [54]. Relatively low mapping rate of BS-seq sequences was found in the present study (~ 45.5–50.2%), which was probably due to the relatively high polymorphism of molluscan species or sequence differences between the studied individuals and the reference genome. Similar results have also been found in Pacific oyster, such as ~ 48.4–49.3% mapping rate for the BS-seq sequences of the wild Pacific oyster by Wang et al. [38], and ~ 53% in the study by Olson and Roberts [37].

Few studies investigating the function of DNA methylation in shell colouration of molluscs have been carried out. In the current study, whole-genome methylation patterns and levels of the Yesso scallop were first analysed and compared between brown and white shell colour individuals. Overall, the DNA methylation landscapes of the two groups were similar, but slight increases in the mCG percentage and methylation level were detected in the brown group. The sequence preference of methylated cytosines is not always consistent among different species [48, 55]. In the present study, the nearby sequences of the mCG, mCHG and mCHH contexts varied across different methylation level sites (high and low), and differences were also observed between the two shell colour groups. DNA methylation levels varied among different genomic regions, with the highest in CDS; this is similar to the analysis in the Pacific oyster by Wang et al. [38]. In the transcribed regions, exons showed higher methylation levels than introns, consistent with the result observed in the sea cucumber [53], but different to most other species including plants [55, 56], mammals [48, 49] and also some molluscs, such as the Pacific oyster [36], eastern oyster [54] and pearl oyster [39], which showed higher methylation levels in introns. Further, all detected regions in the brown individuals exhibited higher methylation levels than those in the white group. All of the above DNA methylation differences between the two shell colour groups may reflect a mechanism that regulates the process of shell colour formation.

To further assess the function of DNA methylation in shell colouration, DMRs were analysed between the brown and white Yesso scallops. The results suggested that DNA methylation alterations could regulate the formation or transportation processes of melanin and porphyrins, thereby affecting shell colouration. Melanin is one of the most common pigments in animals, and has been identified in many molluscan species, including the Yesso scallop [4, 57]. In the present study, many melanin-related GO functions were significantly enriched in the DMRs between the brown and white Yesso scallops, which indicates that DNA methylation plays an important role in regulating the biosynthesis of shell melanin. MITF is a key regulation gene in the melanogenesis pathway [17, 1758, 59]. The study by Sun et al. detected a significantly higher expression level of MITF in brown Yesso scallops compared to white individuals [17]. Though not statistically significant, a higher expression level of MITF was found in the brown scallops in the present study, and a negative correlation was observed between the mRNA expression and DNA methylation levels of the MITF gene. Similar findings have also been reported in a study of duck (Anas platyrhynchos) plumage colouration, where the DNA methylation level was negatively correlated with the mRNA expression of MITF-M [60]. Thus, the results of the present study suggest that the expression of MITF could be negatively regulated by DNA methylation, thereby affecting the biosynthesis of shell melanin.

Porphyrin is another pigment in molluscan shells, and has also been found in some gastropods and bivalves [4, 61, 62]. In vertebrates, porphyrins are synthesized as by-products via an evolutionarily ancient heme biosynthesis pathway, which includes eight genes controlling eight enzymatic reactions [63, 64]. The pathway has also been identified in some molluscan species, such as marine snails Clanculus margaritarius and C. pharaonius, and has been shown to be related to shell colouration [18]. Transcriptome analyses in the Manila clam (Ruditapes philippinarum) [20], hard clam (Mercenaria mercenaria) [21] and Yesso scallop [15, 19] have identified that the ‘porphyrin and chlorophyll metabolism’ (KO00860) pathway, which is a large pathway containing the heme biosynthesis pathway, is involved in shell colour formation. In the present study, GO functions, including porphyrin and heme biosynthetic processes, and the porphyrin and chlorophyll metabolism pathway were also significantly enriched in the DMRs between the two shell colour scallops, indicating that the biosynthesis process of porphyrins may be regulated by DNA methylation. Seven genes involved in this process were found to be differentially methylated, but complex methylation level changes were observed. Specifically, there were positive correlations between mRNA expression levels and DNA methylation status for the ALAS, UROS, CPOX and PPOX genes, a negative correlation was detected for the ALAD gene, while different correlations were observed for different DNA methylation sites in the UROD and FECH genes. These results probably indicate that different DNA methylation modulation mechanisms exist for different genes involved in shell colouration. Further comprehensive studies are needed to verify this assertion.

Conclusions

In conclusion, WGBS analysis was first performed in Yesso scallops with different shell colours, and DNA methylomes at single-base resolution were generated. Characteristics of the DNA methylation patterns and levels of the Yesso scallop were analysed and compared between the two shell colour groups. DMR analysis between the two kinds of scallops indicated that DNA methylation regulates genes and pathways involved in melanin and porphyrin formation, thereby affecting shell colouration. Moreover, different correlations between mRNA expression levels and DNA methylation status were found in different shell colour-related genes. This study provides new genome-wide evidence for understanding the epigenetic regulatory mechanism underlying the shell colour trait.

Materials and methods

Scallop collection

Two-year-old Yesso scallops with brown shell colors were collected from the Dalian Zhangzidao Sea area (Liaoning, China). The scallops were acclimated in the laboratory for 1 week and then sampled following the procesures in [65]. Briefly, filtered and aerated seawater was maintained at approximately 8 °C, which is within the optimum temperature range for their growth. The mantle tissues of the left valves were sampled for DNA extraction. The tissues were immediately frozen in liquid nitrogen and stored in an ultra-low temperature freezer (− 80 °C). To avoid the effect of sex on DNA methylation, scallops with the same sex (male) were chosen in the present study. No specific permits were needed for the described field studies as the Yesso scallops utilized in this study were cultured marine animals provided for market sale, and they are not part of an endangered or protected species. All experiments were conducted following the regulations of the local and central government.

DNA extraction, library construction and sequencing

Genomic DNA of the mantle tissues was extracted using a standard phenol-chloroform method. The DNA quality and purity were tested by 1% agarose gel electrophoresis and the NanoDrop2000 spectrophotometer (Thermo Scientific, Wilmington, DE, USA). A total of 5.2 μg genomic DNA spiked with 26 ng lambda DNA was fragmented to 200 ~ 300 bp by sonication with Covaris S220 Focused-ultra sonicator (Covaris, Woburn, MA, USA). Lambda DNA was used as an unmethylated control for calculating the bisulphite conversion rate. End-repair, 3′-end dA addition and cytosine-methylated adaptor ligation were subsequently performed. Then, the DNA fragments were treated twice with bisulphite by the EN DNA methylation Gold™ Kit (Zymo Research, Irvine, CA, USA). Through bisulphite treatment, unmethylated C changed to U while methylated C kept unchanged. The libraries were constructed using TruSeq DNA Methylation Kit (Illumina, San Diego, CA, USA) according to the manufacturer’s instructions. Finally, all the libraries were subjected to 150 bp paired-end sequencing on the Illumina Hiseq X TEN genomic sequencing platform. Three biological replicates separately for each kind of shell color scallops were utilized for WGBS libraries construction and sequencing.

Sequencing data analysis

Raw sequencing data was processed using fastp software [66] to remove reads containing adapters or ambiguous “N” nucleotides and low-quality reads. The remaining high-quality clean reads were used for the subsequent analyses. Then, the clean reads were aligned to the Yesso scallop genome [44] using Bismark software [67] with the default parameters. Conversions of C to T and G to A were separately made in the reference genome and sequence reads before alignment, and the bisulfite-converted genome was indexed by Bowtie2 [68]. The base distribution nearby the methylated cytosine was analyzed by Seqlogo [45], including the frequency of different methylated cytosine sequence contexts, the Logo Plots for 9-mer sequence containing the methylated cytosine, and sequence features of different sequence contexts.

Methylation sites were identified using MethylKit software [69] with the settings that at least 10 reads covering a site (Covered C) and at last one read supported methylated C. Binomial test was used to determine a real methylation site with the FDR ≤ 0.05. Differentially methylated regions (DMRs) were subsequently identified by MethylKit with an unadjusted p value threshold of 0.05 and an absolute cutoff of 10% between the two shell color groups. GO and KEGG pathway (https://www.kegg.jp/kegg/kegg1.html) enrichment analysis of DMR related genes were performed respectively using R software based on the hypergeometric distribution algorithm with the following formulas:

$$ p=1-\sum \limits_{f-0}^{m=1}\frac{\left(\begin{array}{l}M\\ {}i\end{array}\right)\left(\begin{array}{l}N-M\\ {}n-i\end{array}\right)}{\left(\begin{array}{l}N\\ {}n\end{array}\right)} $$
$$ Enrichment=\raisebox{1ex}{$\frac{m}{n}$}\!\left/ \!\raisebox{-1ex}{$\frac{M}{N}$}\right. $$

where N is the number of all genes with GO or KEGG annotations, n is the number of genes related to DMR involved in N, M is the number of genes with a specific GO or KEGG annotation, and m is the number of genes related to DMR involved in M. GO functions or pathways with the FDR lower than 0.05 were considered to be significantly enriched.

Quantitative real-time PCR analysis

qRT-PCR analysis was carried out mainly according to the method mentioned in [65]. Total RNA was isolated from each sample using an RNAprep pure tissue kit (Tiangen). The quantity and quality of total RNA were determined by using the Micro-Spectrophotometer NV3000 (Vastech) and agarose gel electrophoresis. First strand cDNA was synthesized using PrimeScript RT reagent Kit (TaKaRa) following the manufacturer’s protocol. All of the cDNA products were diluted to 200 ng/ul for use as the templates in quantitative real-time PCR. qRT-PCR was conducted using the FastStart Essential DNA Green Master kit (Roche) on a Roche Light Cycler 96 System (Roche) The running program was as follows: 95 °C for 10 min, followed by 40 cycles of 95 °C for 10 s and 60 °C for 30 s, then 95 °C for 10 s and 65 °C for 1 min, finally 97 °C for 1 s and 37 °C for 30 s. The β-actin gene was selected as a reference gene. All the primers used for qRT-PCR were designed by Primer Premier 5.0 software, and the primer sequences and products length were listed in Supplementary Table S5. The specificity of the primers was assessed by alignment with the Yesso scallop genome by BLASTN with an e-value of 1e-10. Melting curve analysis was also performed to verify that each primer set amplified a single product. Three technical replicates were performed for each reaction and three individuals for each type of shell color strain were used as biological replicates.

The data from real-time PCR was analysed by the 2-ΔΔCT method [70] to evaluate the relative gene expression levels of the mantle tissues in brown and white Yesso scallops. Statistical analysis of the data was performed with SPSS (version 19.0) software using the independent T-Test. Differences were considered significant at p < 0.05.

Availability of data and materials

The raw data of WGBS have been deposited into the NCBI SRA database with the accession number of PRJNA695315. The datasets generated during this study are included in the article and its supplementary information files.

Abbreviations

WGBS:

Whole-genome bisulfite sequencing

DMR:

Differentially methylated region

MITF :

Microphthalmia-associated transcription factor

ALAS :

5-aminolevulinate synthase

ALAD :

Delta-aminolevulinic acid dehydratase

UROS :

Uroporphyrinogen-III synthase

UROD :

Uroporphyrinogen decarboxylase

CPOX :

Coproporphyrinogen-III oxidase

PPOX :

Protoporphyrinogen oxidase

FECH :

Ferrochelatase

GO:

Gene Ontology term

KEGG:

Kyoto Encyclopedia of Genes and Genomes

References

  1. Caro T. The adaptive significance of coloration in mammals. BioScience. 2005;55(2):125–36.

    Article  Google Scholar 

  2. Hoekstra H. Genetics, development and evolution of adaptive pigmentation in vertebrates. Heredity. 2006;97(3):222–34.

    Article  CAS  PubMed  Google Scholar 

  3. Hubbard J, Uy J, Hauber M, Hoekstra H, Safran R. Vertebrate pigmentation: from underlying genes to adaptive function. Trends Genet. 2010;26(5):231–9.

    Article  CAS  PubMed  Google Scholar 

  4. Williams S. Molluscan shell colour. Biol Rev. 2016;92(2):1039–58.

    Article  PubMed  Google Scholar 

  5. Lindberg D, Pearse J. Experimental manipulation of Shell color and morphology of the limpets Lottia-Asmi (Middendorff) and Lottia-Digitalis (Rathke) (Mollusca, Patellogastropoda). J Exp Mar Biol Ecol. 1990;140(3):173–85.

    Article  Google Scholar 

  6. Zheng H, Liu H, Zhang T, Wang S, Sun Z, 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 

  7. Kobayashi T, Kawahara I, Hasakura O, Kijima A. Genetic control of bluish shell color variation in the Pacific abalone Haliotis discus hannai. J Shellfish Res. 2004;23(4):1153–6.

    Google Scholar 

  8. Liu X, Wu F, Zhao H, Zhang G, Guo X. A novel shell color variant of the Pacific abalone Haliotis discus hannai Ino subject to genetic control and dietary influence. J Shellfish Res. 2009;28(2):419–24.

    Article  CAS  Google Scholar 

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

  10. Richards P, Liu M, Lowe N, Davey J, Blaxter M, et al. RAD-Seq derived markers flank the shell colour and banding loci of the Cepaea nemoralis supergene. Mol Ecol. 2013;22(11):3077–89.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Qin Y, Liu X, Zhang H, Zhang G, Guo X. Identification and mapping of amplified fragment length polymorphism markers linked to shell color in bay scallop, Argopecten irradians irradians (Lamarck, 1819). Mar Biotechnol. 2007;9(1):66–73.

    Article  CAS  Google Scholar 

  12. Zou K, Zhang D, Guo H, Zhu C, Li M, et al. A preliminary study for identification of candidate AFLP markers under artificial selection for shell color in pearl oyster Pinctada fucata. Gene. 2014;542(1):8–15.

    Article  CAS  PubMed  Google Scholar 

  13. Zhong X, Li Q, Guo X, Yu H, Kong L. QTL mapping for glycogen content and shell pigmentation in the Pacific oyster Crassostrea gigas using microsatellites and SNPs. Aquac Int. 2014;22(6):1877–89.

    Article  CAS  Google Scholar 

  14. Zhao L, Li Y, Li Y, Yu J, Liao H, et al. A genome-wide association study identifies the genomic region associated with shell color in yesso scallop, Patinopecten yessoensis. Mar Biotechnol. 2017;19(3):301–9.

    Article  CAS  Google Scholar 

  15. Ding J, Zhao L, Chang Y, Zhao W, Du Z, et al. Transcriptome sequencing and characterization of Japanese scallop Patinopecten yessoensis from different shell color lines. PLoS One. 2015;10(2):e0116406.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  16. Sun X, Yang A, Wu B, Zhou L, Liu Z. Characterization of the mantle transcriptome of Yesso scallop (Patinopecten yessoensis): identification of genes potentially involved in biomineralization and pigmentation. PLoS One. 2015;10(4):e0122967.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  17. Sun X, Liu Z, Zhou L, Wu B, Dong Y, et al. Integration of next generation sequencing and EPR analysis to uncover molecular mechanism underlying shell color variation in scallops. PLoS One. 2016;11(8):e0161876.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  18. Williams S, Lockyer A, Dyal P, Nakano T, Churchill C, et al. Colorful seashells: identification of haem pathway genes associated with the synthesis of porphyrin shell color in marine snails. Ecol Evol. 2017;7(23):10379–97.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Mao J, Zhang W, Wang X, Song J, Yin D, et al. Histological and expression differences among different mantle regions of the yesso scallop (Patinopecten yessoensis) provide insights into the molecular mechanisms of biomineralization and pigmentation. Mar Biotechnol. 2019;21(5):683–96.

    Article  CAS  Google Scholar 

  20. Nie H, Jiang K, Jiang L, Huo Z, Ding J, et al. Transcriptome analysis reveals the pigmentation related genes in four different shell color strains of the Manila clam Ruditapes philippinarum. Genomics. 2020;112(2):2011–20.

    Article  CAS  PubMed  Google Scholar 

  21. Hu Z, Song H, Zhou C, Yu Z, Yang M, et al. De novo assembly transcriptome analysis reveals the preliminary molecular mechanism of pigmentation in juveniles of the hard clam Mercenaria mercenaria. Genomics. 2020;112(5):3636–47.

    Article  CAS  PubMed  Google Scholar 

  22. Wolffe A, Matzke M. Epigenetics: regulation through repression. Science. 1999;286:481–6.

    Article  CAS  PubMed  Google Scholar 

  23. Jablonka E, Lamb M. The changing concept of epigenetics. Ann N Y Acad Sci. 2002;981(1):82–96.

    Article  PubMed  Google Scholar 

  24. Turner B. Epigenetic responses to environmental change and their evolutionary implications. Philos Trans R Soc Lond B Biol Sci. 2009;364(1534):3403–18.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Li E, Bestor T, Jaenisch R. Targeted mutation of the DNA methyltransferase gene results in embryonic lethality. Cell. 1992;69(6):915–26.

    Article  CAS  PubMed  Google Scholar 

  26. Hsieh C. Dependence of transcriptional repression on CpG methylation density. Mol Cell Biol. 1994;14(8):5487–94.

    CAS  PubMed  PubMed Central  Google Scholar 

  27. Okano M, Bell D, Haber D, Li E. DNA methyltransferases Dnmt3a and Dnmt3b are essential for de novo methylation and mammalian development. Cell. 1999;99(3):247–57.

    Article  CAS  PubMed  Google Scholar 

  28. Maloisel L, Rossignol J. Suppression of crossing-over by DNA methylation in Ascobolus. Genes Dev. 1998;12(9):1381–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Zilberman D, Gehring M, Tran R, Ballinger T, Henikoff S. Genome-wide analysis of Arabidopsis thaliana DNA methylation uncovers an interdependence between methylation and transcription. Nat Genet. 2007;39(1):61–9.

    Article  CAS  PubMed  Google Scholar 

  30. Bell A, Felsenfeld G. Methylation of a CTCF-dependent boundary controls imprinted expression of the Igf2 gene. Nature. 2000;405(6785):482–5.

    Article  CAS  PubMed  Google Scholar 

  31. Csankovszki G, Nagy A, Jaenisch R. Synergism of Xist RNA, DNA methylation, and histone Hypoacetylation in maintaining X chromosome inactivation. J Cell Biol. 2001;153(4):773–84.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Jiang L, Zhang M, Ma K. Whole-genome DNA methylation associated with differentially expressed genes regulated anthocyanin biosynthesis within flower color chimera of ornamental tree Prunus mume. Forests. 2020;11(1):90.

    Article  Google Scholar 

  33. Li Z, Ren T, Li W, Zhou Y, Han R, et al. Association between the methylation statuses at CpG sites in the promoter region of the SLCO1B3, RNA expression and color change in blue eggshells in Lushi chickens. Front Genet. 2019;10:161.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Jung K, Kim S, Yoon J. Differential methylation of melanin-related epigenetic genes during brindle cattle growth. J Anim Health Prod. 2020;8(2):80–8.

    Article  Google Scholar 

  35. Gavery M, Roberts S. DNA methylation patterns provide insight into epigenetic regulation in the Pacific oyster (Crassostrea gigas). BMC Genomics. 2010;11(1):1–9.

    Article  CAS  Google Scholar 

  36. Gavery M, Roberts S. Predominant intragenic methylation is associated with gene expression characteristics in a bivalve mollusc. PeerJ. 2013;1:e215.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  37. Olson C, Roberts S. Genome-wide profiling of DNA methylation and gene expression in Crassostrea gigas male gametes. Front Physiol. 2014;5:224.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Wang X, Li Q, Lian J, Li L, Jin L, et al. Genome-wide and single-base resolution DNA methylomes of the Pacific oyster Crassostrea gigas provide insight into the evolution of invertebrate CpG methylation. BMC Genomics. 2014;15(1):1–12.

    Article  CAS  Google Scholar 

  39. Zhang J, Luo S, Gu Z, Deng Y, Jiao Y. Genome-wide DNA methylation analysis of mantle edge and mantle central from pearl oyster Pinctada fucata martensii. Mar Biotechnol. 2020;22(3):380–90.

    Article  CAS  Google Scholar 

  40. Li Y, Zhang L, Li Y, Li W, Guo Z, et al. Dynamics of DNA methylation and DNMT expression during gametogenesis and early development of scallop Patinopecten yessoensis. Mar Biotechnol. 2019;21(2):196–205.

    Article  CAS  Google Scholar 

  41. Wang S, Lv J, Zhang L, Dou J, Sun Y, et al. MethylRAD: a simple and scalable method for genome-wide DNA methylation profiling using methylation-dependent restriction enzymes. Open Biol. 2015;5(11):150130.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

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

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Harris R, Wang T, Coarfa C, Nagarajan R, Hong C, et al. Comparison of sequencing-based methods to profile DNA methylation and identification of monoallelic epigenetic modifications. Nat Biotechnol. 2010;28(10):1097–105.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Wang S, Zhang J, Jiao W, Li J, Xun X, et al. Scallop genome provides insights into evolution of bilaterian karyotype and development. Nat Ecol Evol. 2017;1(5):0120.

    Article  Google Scholar 

  45. Wagih O. Ggseqlogo: a versatile R package for drawing sequence logos. Bioinformatics. 2017;33(22):3645–7.

    Article  CAS  PubMed  Google Scholar 

  46. Feng S, Cokus S, Zhang X, Chen P, Bostick M, et al. Conservation and divergence of methylation patterning in plants and animals. Proc Natl Acad Sci. 2010;107(19):8689–94.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Zemach A, McDaniel I, Silva P, Zilberman D. Genome-wide evolutionary analysis of eukaryotic DNA methylation. Science. 2010;328(5980):916–9.

    Article  CAS  PubMed  Google Scholar 

  48. Lister R, Pelizzola M, Dowen R, Hawkins R, Hon G, et al. Human DNA methylomes at base resolution show widespread epigenomic differences. Nature. 2009;462(7271):315–22.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Li Y, Zhu J, Tian G, Li N, Li Q, et al. The DNA methylome of human peripheral blood mononuclear cells. PLoS Biol. 2010;8(11):e1000533.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  50. Xiang H, Zhu J, Chen Q, Dai F, Li X, et al. Single base-resolution methylome of the silkworm reveals a sparse epigenomic map. Nat Biotechnol. 2010;28(5):516.

    Article  CAS  PubMed  Google Scholar 

  51. Lyko F, Foret S, Kucharski R, Wolf S, Falckenhayn C, et al. The honey bee epigenomes: differential methylation of brain DNA in queens and workers. PLoS Biol. 2010;8(11):e1000506.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  52. Bonasio R, Li Q, Lian J, Mutti N, Jin L, et al. Genome-wide and caste-specific DNA methylomes of the ants Camponotus floridanus and Harpegnathos saltator. Curr Biol. 2012;22(19):1755–64.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Yang Y, Zheng Y, Sun L, Chen M. Genome-wide DNA methylation signatures of sea cucumber Apostichopus japonicus during environmental induced aestivation. Genes. 2020;11(9):1020.

    Article  CAS  PubMed Central  Google Scholar 

  54. Venkataraman Y, Downey-Wall A, Ries J, Westfield I, White S, et al. General DNA methylation patterns and environmentally-induced differential methylation in the eastern oyster (Crassostrea virginica). Front Mar Sci. 2020;7:225.

    Article  Google Scholar 

  55. Lu X, Wang X, Chen X, Shu N, Wang J, et al. Single-base resolution methylomes of upland cotton (Gossypium hirsutum L.) reveal epigenome modifications in response to drought stress. BMC Genomics. 2017;18(1):1–14.

    Article  Google Scholar 

  56. Xu J, Zhou S, Gong X, Song Y, van Nocker S, et al. Single-base methylome analysis reveals dynamic epigenomic differences associated with water deficit in apple. Plant Biotechnol J. 2018;16(2):672–87.

    Article  CAS  PubMed  Google Scholar 

  57. Sun X, Wu B, Zhou L, Liu Z, Dong Y, et al. Isolation and characterization of melanin pigment from yesso scallop Patinopecten yessoensis. J Ocean Univ China. 2017;16(2):279–84.

    Article  CAS  Google Scholar 

  58. Zhang S, Wang H, Yu J, Jiang F, Yue X, et al. Identification of a gene encoding microphthalmia-associated transcription factor and its association with shell color in the clam Meretrix petechialis. Comp Biochem Physiol B: Biochem Mol Biol. 2018;225:75–83.

    Article  CAS  Google Scholar 

  59. Mao J, Zhang X, Zhang W, Tian Y, Wang X, et al. Genome-wide identification, characterization and expression analysis of the MITF gene in Yesso scallops (Patinopecten yessoensis) with different shell colors. Gene. 2019;688:155–62.

    Article  CAS  PubMed  Google Scholar 

  60. Lin R, Lin W, Zhou S, et al. Integrated analysis of mRNA expression, CpG island methylation, and polymorphisms in the MITF gene in ducks (Anas platyrhynchos). Biomed Res Int. 2019;2019:8512467.

    Article  PubMed  PubMed Central  Google Scholar 

  61. Williams S, Ito S, Wakamatsu K, Goral T, Edwards N, et al. Identification of shell colour pigments in marine snails Clanculus pharaonius and C. margaritarius (Trochoidea; Gastropoda). PLoS One. 2016;11(7):e0156664.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  62. Verdes A, Cho W, Hossain M, Brennan P, Hanley D, et al. Nature’s palette: characterization of shared pigments in colorful avian and mollusk shells. PLoS One. 2015;10(12):e0143545.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  63. Ajioka R, Phillips J, Kushner J. Biosynthesis of heme in mammals. Biochimica et Biophysica Acta. 2006;1763(7):723–36.

    Article  CAS  PubMed  Google Scholar 

  64. Heinemann I, Jahn M, Jahn D. The biochemistry of heme biosynthesis. Arch Biochem Biophys. 2008;474(2):238–51.

    Article  CAS  PubMed  Google Scholar 

  65. Xing Q, Yu Q, Dou H, Wang J, Li R, et al. Genome-wide identification, characterization and expression analyses of two TNFRs in Yesso scallop (Patinopecten yessoensis) provide insight into the disparity of responses to bacterial infections and heat stress in bivalves. Fish Shellfish Immunol. 2016;52:44–56.

    Article  CAS  PubMed  Google Scholar 

  66. Chen S, Zhou Y, Chen Y, Gu J. Fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics. 2018;34(17):i884–90.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  67. Krueger F, Andrews S. Bismark: a flexible aligner and methylation caller for bisulfite-Seq applications. Bioinformatics. 2011;27(11):1571–2.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  68. Langmead B, Salzberg S. Fast gapped-read alignment with bowtie 2. Nat Methods. 2012;9(4):357.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  69. Akalin A, Kormaksson M, Li S, Garrett-Bakelman F, Figueroa M, et al. methylKit: a comprehensive R package for the analysis of genome-wide DNA methylation profiles. Genome Biol. 2012;13(10):1–9.

    Article  Google Scholar 

  70. Livak K, Schmittgen T. Analysis of relative gene expression data using real-time quantitative PCR and the 2− ΔΔCT method. Methods. 2001;25(4):402–8.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This project was supported by Dalian High-level Talent Innovation Support Program (2019RQ143), Scientific research fund of Liaoning Provincial Department of Education (LJKZ0705), Dalian Science and Technology Innovation Grant (2018J11CY029), the university innovation team and innovation talent support plan of Liaoning Province (LT2019003) and the National Natural Science Foundation of China (31,702,342 and 42,076,101).

Author information

Authors and Affiliations

Authors

Contributions

JM and YC conceived and designed the experiments; CY, HS, YW and MG prepared the samples; CY and HS performed the experiments; JM, CY, XW and YT analyzed the data; JM and CY wrote the paper; JD, ZH and YC reviewed the paper; YC, JM, HZ and JD acquired the funding. The authors read and approved the final manuscript.

Corresponding authors

Correspondence to Junxia Mao or Yaqing Chang.

Ethics declarations

Ethics approval and consent to participate

All experiments were conducted following the regulations of the local and central government.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Figure S1.

Sequence preferences for DNA methylation in different sequence contexts in two shell colour Yesso scallops. A. Sequence preference analysis in CG context. For this type, sites with methylation level >75% were defined as high methylation level sites and others were low methylation level sites. B. Sequence preference analysis in CHG context. C. Sequence preferences in CHH context. For CHG and CHH context, sites with methylation level >25% were defined as high methylation level sites and others were low methylation level sites. Table S1. Statistical results of methylated cytosines in different sequence contexts. Table S5. Primers for quantitative real-time PCR used in this study.

Additional file 2: Table S2.

DMRs detected between the brown and white Yesso scallops. Table S3. All GO terms significantly enriched for differentially methylated regions. Table S4. All KEGG pathways significantly enriched for differentially methylated regions (p < 0.05).

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Yuan, C., Mao, J., Sun, H. et al. Genome-wide DNA methylation profile changes associated with shell colouration in the Yesso scallop (Patinopecten yessoensis) as measured by whole-genome bisulfite sequencing. BMC Genomics 22, 740 (2021). https://doi.org/10.1186/s12864-021-08055-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-021-08055-6

Keywords