Abstract
Understanding the combined effects of environmental heterogeneity and evolutionary processes on marine populations is a primary goal of seascape genomic approaches. Here, we utilized genomic approaches to identify local adaptation signatures in Littoraria flava, a widely distributed marine gastropod in the tropical West Atlantic population. We also performed molecular evolution analyses to investigate potential selective signals across the genome. After obtaining 6,298 and 16,137 single nucleotide polymorphisms derived from genotyping-by-sequencing and RNA sequencing, respectively, 69 from genotyping-by-sequencing (85 specimens) and four from RNA sequencing (40 specimens) candidate single nucleotide polymorphisms were selected and further evaluated. The correlation analyses support different evolutionary pressures over transcribed and non-transcribed regions. Thus, single nucleotide polymorphisms within transcribed regions could account for the genotypic and possibly phenotypic divergences in periwinkles. Our molecular evolution tests based on synonymous and non-synonymous ratio (kN/kS) showed that genotype divergences containing putative adaptive single nucleotide polymorphisms arose mainly from synonymous and/or UTR substitutions rather than polymorphic proteins. The distribution of genotypes across different localities seems to be influenced by marine currents, pH, and temperature variations, suggesting that these factors may impact the species dispersion. The combination of RNA sequencing and genotyping-by-sequencing derived datasets provides a deeper understanding of the molecular mechanisms underlying selective forces responses on distinct genomic regions and could guide further investigations on seascape genomics for non-model species.