Frequencies of house fly proto-Y chromosomes across populations are predicted by temperature heterogeneity within populations

Sex chromosomes often differ between closely related species and can even be polymorphic within populations. Species with multifactorial sex determination segregate for multiple different sex-determining loci within populations, making them uniquely informative of the selection pressures that drive the evolution of sex chromosomes. The house fly (Musca domestica) is a model species for studying multifactorial sex determination because male-determining genes have been identified on all six of the chromosomes, which means that any chromosome can be a "proto-Y." Natural populations of house flies also segregate for a recently derived female-determining locus, meaning house flies also have a proto-W chromosome. The different proto-Y chromosomes are distributed along latitudinal clines on multiple continents, their distributions can be explained by seasonality in temperature, and they have temperature-dependent effects on physiological and behavioral traits. It is not clear, however, how the clinal distributions interact with the effect of seasonality on the frequencies of house fly proto-Y and proto-W chromosomes across populations. To address this question, we measured the frequencies of house fly proto-Y and proto-W chromosomes across nine populations in the United States of America. We confirmed the clinal distribution along the eastern coast of North America, but it is limited to the eastern coast. In contrast, the annual mean daily temperature range predicts proto-Y chromosome frequencies across the entire continent. Our results therefore suggest that temperature heterogeneity can explain the distributions of house fly proto-Y chromosomes in a way that does not depend on the cline.