Parallel evolution in Crassostrea oysters along the latitudinal gradient is associated with variation in multiple genes involved in adipogenesis

Parallel diversification provides a proper framework for studying the role of natural selection in evolution. Yet, empirical studies from ecological 'non-model' species of invertebrates are limited at the whole genome level. Here, we presented a chromosome-scale genome assembly for Crassostrea angulata and investigated the parallel genomic evolution in oysters. Specifically, we used population genomics approaches to compare two southern-northern oyster species pairs (C. angulata-C. gigas and southern-northern C. ariakensis) along the coast of China. The estimated divergence time of C. angulata and C. gigas is earlier than that of southern and northern C. ariakensis, which aligns with the overall elevated genome-wide divergence. However, the southern-northern C. ariakensis FST profile represented more extremely divergent "islands". Combined with recent reciprocal hybridization studies, we proposed that they are currently at an early stage of speciation. These two southern-northern oyster species pairs exhibited significant repeatability in patterns of genome-wide differentiation, especially in genomic regions with extremely high and low divergence. This suggested that divergent and purifying selection has contributed to the genomic parallelism between southern and northern latitudes. Top differentiated genomic regions shared in these two oyster species pairs contained candidate genes enriched for functions in energy metabolism, especially adipogenesis, which are closely related to reproductive behaviours. These genes might be good candidates for further investigation in vivo. In conclusion, our results suggest that similar divergent selection and shared genomic features could predictably transform standing genetic variation within one species pair into differences in another.