Minimal effects of proto-Y chromosomes on house fly gene expression in spite of evidence that selection maintains stable polygenic sex determination

Sex determination is the developmental process by which organismal sex is established. Sex determination evolves fast, often due to changes in the master regulators at the top of the pathway. In addition, some species are polymorphic for multiple different master regulators within natural populations. Understanding the forces that maintain this polygenic sex determination can be informative of the factors that drive the evolution of sex determination. The house fly, Musca domestica, is a well-suited model to those ends because natural populations harbor male-determining loci on each of the six chromosomes and a bi-allelic female-determiner. Multiple lines of evidence suggest that natural selection maintains polygenic sex determination in house fly. However, previous work found that there are very few sequence differences between proto-Y chromosomes and their homologous X chromosomes. This suggests that there is not much genetic variation upon which natural selection could act to maintain polygenic sex determination in house fly. To address this paradox, we performed RNA-seq experiments that examine the effects of the two most common proto-Y chromosomes on gene expression. We find that the proto-Y chromosomes do indeed have a relatively minor effect on gene expression, as expected based on the minimal X-Y sequence differences. Despite these minimal gene expression differences, we identify some patterns that are consistent with sex-specific selection acting on phenotypic effects of proto-Y chromosomes. Our results suggest that, if natural selection maintains polygenic sex determination in house fly, the phenotypic differences under selection are minor and possibly depend on ecological contexts that were not tested in our experimental design.