Functional and population genomic divergence within and between two species of killifish adapted to different osmotic niches

Adaptation to salinity affects species distributions, promotes speciation, and guides many evolutionary patterns in fishes. To uncover the basis of a complex trait like osmoregulation, genome-level analyses are sensible. We combine population genomic scans with genome expression profiling to discover candidate genes and pathways associated with divergence between osmotic environments. We compared transcriptome sequence divergence between multiple freshwater and saltwater populations of the rainwater killifish, Lucania parva. We also compared sequence divergence between L. parva and its sister species, Lucania goodei, a freshwater specialist. We found highly differentiated single nucleotide polymorphisms (SNPs) between freshwater and saltwater L. parva populations in cell junction and ion transport genes, including V-type H+ ATPase. Between species, we found divergence in reproduction and osmotic stress genes. Genes that were differentially expressed between species during osmotic acclim...