Genome-wide association study reveals genetic variations associated with ocean acidification resilience in Yesso scallop Patinopecten yessoensis.

Ocean acidification (OA), which refers to a gradual decrease in seawater pH due to the absorption of atmospheric carbon dioxide, profoundly affects the growth, development and survival of bivalves. Relatively limited studies have assessed the resilience of bivalve to OA. In the present study, Patinopecten yessoensis, an economically and ecologically significant species, were exposed to low pH (pH = 7.5) for 4 weeks. Forty-seven scallops that died in the first week were considered pH-sensitive population, and 20 that were alive at the end of the experiment were considered pH-tolerant population. A genome-wide association study was conducted to identify the genomic loci associated the resilience of P. yessoensis to OA. Twenty-one single nucleotide polymorphisms were significantly associated with resilience, which were distributed in 11 linkage groups. Within the linkage disequilibrium block region (± 300 kb) surrounding the 21 SNPs, 193 candidate genes were successfully identified. Particularly, five associated SNPs were directly located on five genes, including SP24, CFDH, 5HTR3, HSDL1 and ZFP346. The GO enrichment and KEGG pathway analyses showed that the molecular response of P. yessoensis to OA mainly involved neural signal transmission, energy metabolism and redox reaction. Candidate genes were expressed during larval development and in adult tissues. Furthermore, the expression of 30 candidate genes changed significantly under low pH stress in the mantle. Our results reveal certain SNPs and candidate genes that could elucidate the different responses of P. yessoensis to OA. The genetic variations indicated molecular resilience in P. yessoensis populations, which may enable adaptation to future acidification stress.