BackgroundShell color polymorphisms of Mollusca have contributed to development of evolutionary biology and population genetics, while the genetic bases and molecular mechanisms underlying shell pigmentation are poorly understood. The Pacific oyster (Crassostrea gigas) is one of the most important farmed oysters worldwide. Through successive family selection, four shell color variants (white, golden, black and partially pigmented) of C. gigas have been developed. To elucidate the genetic mechanisms of shell coloration in C. gigas and facilitate the selection of elite oyster lines with desired coloration patterns, differentially expressed genes (DEGs) were identified among the four shell color variants by RNA-seq.ResultsDigital gene expression generated over fifteen million reads per sample, producing expression data for 28,027 genes. A total number of 2,645 DEGs were identified from pair-wise comparisons, of which 432, 91, 43 and 39 genes specially were up-regulated in white, black, golden and partially pigmented shell of C. gigas, respectively. Three genes of Abca1, Abca3 and Abcb1 which belong to the ATP-binding cassette (ABC) transporters super-families were significantly associated with white shell formation. A tyrosinase transcript (CGI_10008737) represented consistent up-regulated pattern with golden coloration. We proposed that white shell variant of C. gigas could employ “endocytosis” to down-regulate notch level and to prevent shell pigmentation.ConclusionThis study discovered some potential shell coloration genes and related molecular mechanisms by the RNA-seq, which would provide foundational information to further study on shell coloration and assist in selective breeding in C. gigas.