Parallel genetic changes underlie integrated craniofacial traits in an adaptive radiation of trophic specialist pupfishes

Many factors such as divergence time, shared standing genetic variation, frequency of introgression, and mutation rates can influence the likelihood of whether populations adapt to similar environments via parallel or non-parallel genetic changes. However, the frequency of parallel vs non-parallel genetic changes resulting in parallel phenotypic evolution is still unknown. In this study, we used a QTL mapping approach to investigate the genetic basis of highly divergent craniofacial traits between scale- and snail-eating trophic specialist species across similar hypersaline lake environments in an adaptive radiation of pupfishes endemic to San Salvador Island, Bahamas. We raised F2 intercrosses of scale- and snail-eaters from two different lake populations of sympatric specialists, estimated linkage maps, scanned for significant QTL for 30 skeletal and craniofacial traits, and compared the location of QTL between lakes to quantify parallel and non-parallel genetic changes. We found strong support for parallel genetic changes in both lakes for five traits in which we detected a significant QTL in at least one lake. However, many of these shared QTL affected different, but highly correlated craniofacial traits in each lake, suggesting that pleiotropy and trait integration should not be neglected when estimating rates of parallel evolution. We further observed a 23-52% increase in adaptive introgression within shared QTL, suggesting that introgression may be important for parallel evolution. Overall, our results suggest that the same genomic regions contribute to parallel integrated craniofacial phenotypes across lakes. We also highlight the need for more expansive searches for shared QTL when testing for parallel evolution.