Candidate genes for shell colour polymorphism in Cepaea nemoralis

The characteristic ground colour and banding patterns on shells of the land snail Cepaea nemoralis form a classic study system for genetics and adaptation. We use RNAseq analysis to identify candidate genes underlying this polymorphism. We sequenced cDNA from the body and the mantle (the shell-producing tissue) of four individuals of two phenotypes and produced a de novo transcriptome of 147,397 contigs. Differential expression analysis identified a set of 1,961 transcripts that were upregulated in mantle tissue. Sequence variant analysis resulted in a set of 2,592 transcripts with single nucleotide polymorphisms (SNPs) that differed consistently between the phenotypes. Combining these results yielded a set of 197 candidate transcripts, of which 38 were annotated. Four of these transcripts are involved in production of the shell’s nacreous layer. Comparison with morph-associated RAD-tags from a published study yielded seven transcripts that were annotated as metallothionein, a protein that is thought to inhibit the production of melanin in melanocytes. These results thus provide an excellent starting point for the elucidation of the genetic regulation of the Cepaea nemoralis shell colour polymorphism.