2.3 Fixation index (FST) calculation
For each dataset, we used ANGSD to generate unbiased allele frequency spectra of our populations (AFS) and from that, calculated the weightedF ST for each gene (Korneliussen et al., 2014). While assigning SNP loci to genes, we extended the annotated gene boundaries one-thousand bases at the start and end of the gene, to capture the upstream regulatory region and the downstream untranslated RNA region that might be missing in the gene annotation.