Evolutionary and developmental dynamics of sex-biased gene expression in common frogs with proto-Y chromosomes

The Chinese longsnout catfish genome provides novel insights into the feeding preference and corresponding metabolic strategy of carnivores

Fish show variation in feeding habits to adapt to complex environments. However, the genetic basis of feeding preference and the corresponding metabolic strategies that differentiate feeding habits remain elusive. Here, by comparing the whole genome of a typical carnivorous fish (Leiocassis longirostris Günther) with that of herbivorous fish, we identify 250 genes through both positive selection and rapid evolution, including taste receptor taste receptor type 1 member 3 (tas1r3) and trypsin We demonstrate that tas1r3 is required for carnivore preference in tas1r3-deficient zebrafish and in a diet-shifted grass carp model. We confirm that trypsin correlates with the metabolic strategies of fish with distinct feeding habits. Furthermore, marked alterations in trypsin activity and metabolic profiles are accompanied by a transition of feeding preference in tas1r3-deficient zebrafish and diet-shifted grass carp. Our results reveal a conserved adaptation between feeding preference and corresponding metabolic strategies in fish, and provide novel insights into the adaptation of feeding habits over the evolution course.

Evolution of sex-biased genes in Drosophila species with neo-sex chromosomes: Potential contribution to reducing the sexual conflict

An advantage of sex chromosomes may be the potential to reduce sexual conflict because they provide a basis for selection to operate separately on females and males. However, evaluating the relationship between sex chromosomes and sexual conflict is challenging owing to the difficulty in measuring sexual conflict and substantial divergence between species with and without sex chromosomes. We therefore examined sex-biased gene expression as a proxy for sexual conflict in three sets of Drosophila species with and without young sex chromosomes, the so-called neo-sex chromosomes. In all sets, we detected more sex-biased genes in the species with neo-sex chromosomes than in the species without neo-sex chromosomes in larvae, pupae, and adult somatic tissues but not in gonads. In particular, many unbiased genes became either female- or male-biased after linkage to the neo-sex chromosomes in larvae, despite the low sexual dimorphism. For example, genes involved in metabolism, a key determinant for the rate of development in many animals, were enriched in the genes that acquired sex-biased expression on the neo-sex chromosomes at the larval stage. These genes may be targets of sexually antagonistic selection (i.e., large size and rapid development are selected for in females but selected against in males). These results indicate that acquiring neo-sex chromosomes may have contributed to a reduction in sexual conflict, particularly at the larval stage, in Drosophila..

Gap-free X and Y chromosome assemblies of Salix arbutifolia reveal an evolutionary change from male to female heterogamety in willows, without a change in the position of the sex-determining locus

In the Vetrix clade of Salix, a genus of woody flowering plants, sex determination involves chromosome 15, but an XY system has changed to a ZW system. We studied the detailed genetic changes involved. We used genome sequencing, with chromosome conformation capture (Hi-C) and PacBio HiFi reads to assemble chromosome level gap-free X and Y of Salix arbutifolia, and distinguished the haplotypes in the 15X- and 15Y-linked regions, to study the evolutionary history of the sex-linked regions (SLRs). Our sequencing revealed heteromorphism of the X and Y haplotypes of the SLR, with the X-linked region being considerably larger than the corresponding Y region, mainly due to accumulated repetitive sequences and gene duplications. The phylogenies of single-copy orthogroups within the SLRs indicate that S. arbutifolia and Salix purpurea share an ancestral SLR within a repeat-rich region near the chromosome 15 centromere. During the change in heterogamety, the X-linked region changed to a W-linked one, while the Z was derived from the Y.

Multiple Transitions between Y Chromosome and Autosome in Tago's Brown Frog Species Complex

Sex chromosome turnover is the transition between sex chromosomes and autosomes. Although many cases have been reported in poikilothermic vertebrates, their evolutionary causes and genetic mechanisms remain unclear. In this study, we report multiple transitions between the Y chromosome and autosome in the Japanese Tago's brown frog complex. Using chromosome banding and molecular analyses (sex-linked and autosomal single nucleotide polymorphisms, SNPs, from the nuclear genome), we investigated the frogs of geographic populations ranging from northern to southern Japan of two species, Rana tagoi and Rana sakuraii (2n = 26). Particularly, the Chiba populations of East Japan and Akita populations of North Japan in R. tagoi have been, for the first time, investigated here. As a result, we identified three different sex chromosomes, namely chromosomes 3, 7, and 13, in the populations of the two species. Furthermore, we found that the transition between the Y chromosome (chromosome 7) and autosome was repeated through hybridization between two or three different populations belonging to the two species, followed by restricted chromosome introgression. These dynamic sex chromosome turnovers represent the first such findings in vertebrates and imply that speciation associated with inter- or intraspecific hybridization plays an important role in sex chromosome turnover in frogs.

Sex-biased gene expression across mammalian organ development and evolution

Sexually dimorphic traits are common among mammals and are specified during development through the deployment of sex-specific genetic programs. Because little is known about these programs, we investigated them using a resource of gene expression profiles in males and females throughout the development of five organs in five mammals (human, mouse, rat, rabbit, and opossum) and a bird (chicken). We found that sex-biased gene expression varied considerably across organs and species and was often cell-type specific. Sex differences increased abruptly around sexual maturity instead of increasing gradually during organ development. Finally, sex-biased gene expression evolved rapidly at the gene level, with differences between organs in the evolutionary mechanisms used, but more slowly at the cellular level, with the same cell types being sexually dimorphic across species.

Sex chromosomes as supergenes of speciation: why amphibians defy the rules?

As reflected by the two rules of speciation (Haldane's rule and the large X-/Z-effect), sex chromosomes are expected to behave like supergenes of speciation: they recombine only in one sex (XX females or ZZ males), supposedly recruit sexually antagonistic genes and evolve faster than autosomes, which can all contribute to pre-zygotic and post-zygotic isolation. While this has been mainly studied in organisms with conserved sex-determining systems and highly differentiated (heteromorphic) sex chromosomes like mammals, birds and some insects, these expectations are less clear in organismal groups where sex chromosomes repeatedly change and remain mostly homomorphic, like amphibians. In this article, we review the proposed roles of sex-linked genes in isolating nascent lineages throughout the speciation continuum and discuss their support in amphibians given current knowledge of sex chromosome evolution and speciation modes. Given their frequent recombination and lack of differentiation, we argue that amphibian sex chromosomes are not expected to become supergenes of speciation, which is reflected by the rarity of empirical studies consistent with a 'large sex chromosome effect' in frogs and toads. The diversity of sex chromosome systems in amphibians has a high potential to disentangle the evolutionary mechanisms responsible for the emergence of sex-linked speciation genes in other organisms. This article is part of the theme issue 'Genomic architecture of supergenes: causes and evolutionary consequences'.

Sex-biased gene expression and recent sex chromosome turnover

Cichlids are well known for their propensity to radiate generating arrays of morphologically and ecologically diverse species in short evolutionary time. Following this rapid evolutionary pace, cichlids show high rates of sex chromosome turnover. We here studied the evolution of sex-biased gene (SBG) expression in 14 recently diverged taxa of the Lake Tanganyika Tropheini cichlids, which show different XY sex chromosomes. Across species, sex chromosome sequence divergence predates divergence in expression between the sexes. Only one sex chromosome, the oldest, showed signs of demasculinization in gene expression and potentially contribution to the resolution of sexual conflict. SBGs in general showed high rates of turnovers and evolved mostly under drift. Sexual selection did not shape the rapid evolutionary changes of SBGs. Male-biased genes evolved faster than female-biased genes, which seem to be under more phylogenetic constraint. We found a relationship between the degree of sex bias and sequence evolution driven by sequence differences among the sexes. Consistent with other species, strong sex bias towards sex-limited expression contributes to resolving sexual conflict in cichlids. This article is part of the theme issue 'Challenging the paradigm in sex chromosome evolution: empirical and theoretical insights with a focus on vertebrates (Part II)'.

Sex-chromosome evolution in frogs: what role for sex-antagonistic genes?

Sex-antagonistic (SA) genes are widely considered to be crucial players in the evolution of sex chromosomes, being instrumental in the arrest of recombination and degeneration of Y chromosomes, as well as important drivers of sex-chromosome turnovers. To test such claims, one needs to focus on systems at the early stages of differentiation, ideally with a high turnover rate. Here, I review recent work on two families of amphibians, Ranidae (true frogs) and Hylidae (tree frogs), to show that results gathered so far from these groups provide no support for a significant role of SA genes in the evolutionary dynamics of their sex chromosomes. The findings support instead a central role for neutral processes and deleterious mutations. This article is part of the theme issue 'Challenging the paradigm in sex chromosome evolution: empirical and theoretical insights with a focus on vertebrates (Part I)'.

The Diversity and Evolution of Sex Chromosomes in Frogs

Frogs are ideal organisms for studying sex chromosome evolution because of their diversity in sex chromosome differentiation and sex-determination systems. We review 222 anuran frogs, spanning ~220 Myr of divergence, with characterized sex chromosomes, and discuss their evolution, phylogenetic distribution and transitions between homomorphic and heteromorphic states, as well as between sex-determination systems. Most (~75%) anurans have homomorphic sex chromosomes, with XY systems being three times more common than ZW systems. Most remaining anurans (~25%) have heteromorphic sex chromosomes, with XY and ZW systems almost equally represented. There are Y-autosome fusions in 11 species, and no W-/Z-/X-autosome fusions are known. The phylogeny represents at least 19 transitions between sex-determination systems and at least 16 cases of independent evolution of heteromorphic sex chromosomes from homomorphy, the likely ancestral state. Five lineages mostly have heteromorphic sex chromosomes, which might have evolved due to demographic and sexual selection attributes of those lineages. Males do not recombine over most of their genome, regardless of which is the heterogametic sex. Nevertheless, telomere-restricted recombination between ZW chromosomes has evolved at least once. More comparative genomic studies are needed to understand the evolutionary trajectories of sex chromosomes among frog lineages, especially in the ZW systems.

A bird-like genome from a frog: Mechanisms of genome size reduction in the ornate burrowing frog, Platyplectrum ornatum

The diversity of genome sizes across the tree of life is of key interest in evolutionary biology. Various correlates of variation in genome size, such as accumulation of transposable elements (TEs) or rate of DNA gain and loss, are well known, but the underlying molecular mechanisms driving or constraining genome size are poorly understood. Here, we study one of the smallest genomes among frogs characterized thus far, that of the ornate burrowing frog (Platyplectrum ornatum) from Australia, and compare it to other published frog and vertebrate genomes to examine the forces driving reduction in genome size. At ∼1.06 gigabases (Gb), the P. ornatum genome is like that of birds, revealing four major mechanisms underlying TE dynamics: reduced abundance of all major classes of TEs; increased net deletion bias in TEs; drastic reduction in intron lengths; and expansion via gene duplication of the repertoire of TE-suppressing Piwi genes, accompanied by increased expression of Piwi-interacting RNA (piRNA)-based TE-silencing pathway genes in germline cells. Transcriptomes from multiple tissues in both sexes corroborate these results and provide insight into sex-differentiation pathways in Platyplectrum Genome skimming of two closely related frog species (Lechriodus fletcheri and Limnodynastes fletcheri) confirms a reduction in TEs as a major driver of genome reduction in Platyplectrum and supports a macroevolutionary scenario of small genome size in frogs driven by convergence in life history, especially rapid tadpole development and tadpole diet. The P. ornatum genome offers a model for future comparative studies on mechanisms of genome size reduction in amphibians and vertebrates generally.

A single QTL with large effect is associated with female functional virginity in an asexual parasitoid wasp

During the transition from sexual to asexual reproduction, a suite of reproduction-related sexual traits become superfluous, and may be selected against if costly. Female functional virginity refers to asexual females resisting to mate or not fertilizing eggs after mating. These traits appear to be among the first that evolve during transitions from sexual to asexual reproduction. The genetic basis of female functional virginity remains elusive. Previously, we reported that female functional virginity segregates as expected for a single recessive locus in the asexual parasitoid wasp Asobara japonica. Here, we investigate the genetic basis of this trait by quantitative trait loci (QTL) mapping and candidate gene analyses. Consistent with the segregation of phenotypes, we found a single QTL of large effect, spanning over 4.23 Mb and comprising at least 131 protein-coding genes, of which 15 featured sex-biased expression in the related sexual species Asobara tabida. Two of the 15 sex-biased genes were previously identified to differ between related sexual and asexual population/species: CD151 antigen and nuclear pore complex protein Nup50. A third gene, hormone receptor 4, is involved in steroid hormone mediated mating behaviour. Overall, our results are consistent with a single locus, or a cluster of closely linked loci, underlying rapid evolution of female functional virginity in the transition to asexuality. Once this variant, causing rejection to mate, has swept through a population, the flanking region does not get smaller owing to lack of recombination in asexuals.

When Sex Chromosomes Recombine Only in the Heterogametic Sex: Heterochiasmy and Heterogamety in Hyla Tree Frogs

Sex chromosomes are classically predicted to stop recombining in the heterogametic sex, thereby enforcing linkage between sex-determining (SD) and sex-antagonistic (SA) genes. With the same rationale, a pre-existing sex asymmetry in recombination is expected to affect the evolution of heterogamety, for example, a low rate of male recombination might favor transitions to XY systems, by generating immediate linkage between SD and SA genes. Furthermore, the accumulation of deleterious mutations on nonrecombining Y chromosomes should favor XY-to-XY transitions (which discard the decayed Y), but disfavor XY-to-ZW transitions (which fix the decayed Y as an autosome). Like many anuran amphibians, Hyla tree frogs have been shown to display drastic heterochiasmy (males only recombine at chromosome tips) and are typically XY, which seems to fit the above expectations. Instead, here we demonstrate that two species, H. sarda and H. savignyi, share a common ZW system since at least 11 Ma. Surprisingly, the typical pattern of restricted male recombination has been maintained since then, despite female heterogamety. Hence, sex chromosomes recombine freely in ZW females, not in ZZ males. This suggests that heterochiasmy does not constrain heterogamety (and vice versa), and that the role of SA genes in the evolution of sex chromosomes might have been overemphasized.

The rise and fall of globins in the amphibia

The globin gene repertoire of gnathostome vertebrates is dictated by differential retention and loss of nine paralogous genes: androglobin, neuroglobin, globin X, cytoglobin, globin Y, myoglobin, globin E, and the α- and β-globins. We report the globin gene repertoire of three orders of modern amphibians: Anura, Caudata, and Gymnophiona. Combining phylogenetic and conserved synteny analysis, we show that myoglobin and globin E were lost only in the Batrachia clade, but retained in Gymnophiona. The major amphibian groups also retained different paralogous copies of globin X. None of the amphibian presented α^D-globin gene. Nevertheless, two clades of β-globins are present in all amphibians, indicating that the amphibian ancestor possessed two paralogous proto β-globins. We also show that orthologs of the gene coding for the monomeric hemoglobin found in the heart of Rana catesbeiana are present in Neobatrachia and Pelobatoidea species we analyzed. We suggest that these genes might perform myoglobin- and globin E-related functions. We conclude that the repertoire of globin genes in amphibians is dictated by both retention and loss of the paralogous genes cited above and the rise of a new globin gene through co-option of an α-globin, possibly facilitated by a prior event of transposition.

Hybridization and introgression between toads with different sex chromosome systems

The growing interest in the lability of sex determination in non‐model vertebrates such as amphibians and fishes has revealed high rates of sex chromosome turnovers among closely related species of the same clade. Can such lineages hybridize and admix with different sex‐determining systems, or could the changes have precipitated their speciation? We addressed these questions in incipient species of toads (Bufonidae), where we identified a heterogametic transition and characterized their hybrid zone with genome‐wide markers (RADseq). Adult and sibship data confirmed that the common toad B. bufo is female heterogametic (ZW), while its sister species the spined toad B. spinosus is male heterogametic (XY). Analysis of a fine scale transect across their parapatric ranges in southeastern France unveiled a narrow tension zone (∼10 km), with asymmetric mitochondrial and nuclear admixture over hundreds of kilometers southward and northward, respectively. The geographic extent of introgression is consistent with an expansion of B. spinosus across B. bufo’s former ranges in Mediterranean France, as also suggested by species distribution models. However, widespread cyto‐nuclear discordances (B. spinosus backrosses carrying B. bufo mtDNA) run against predictions from the dominance effects of Haldane's rule, perhaps because Y and W heterogametologs are not degenerated. Common and spined toads can thus successfully cross‐breed despite fundamental differences in their sex determination mechanisms, but remain partially separated by reproductive barriers. Whether and how the interactions of their XY and ZW genes contribute to these barriers shall provide novel insights on the debated role of labile sex chromosomes in speciation.

Differential Gene Expression between Fungal Mating Types Is Associated with Sequence Degeneration

Degenerative mutations in non-recombining regions, such as in sex chromosomes, may lead to differential expression between alleles if mutations occur stochastically in one or the other allele. Reduced allelic expression due to degeneration has indeed been suggested to occur in various sex-chromosome systems. However, whether an association occurs between specific signatures of degeneration and differential expression between alleles has not been extensively tested, and sexual antagonism can also cause differential expression on sex chromosomes. The anther-smut fungus Microbotryum lychnidis-dioicae is ideal for testing associations between specific degenerative signatures and differential expression because 1) there are multiple evolutionary strata on the mating-type chromosomes, reflecting successive recombination suppression linked to mating-type loci; 2) separate haploid cultures of opposite mating types help identify differential expression between alleles; and 3) there is no sexual antagonism as a confounding factor accounting for differential expression. We found that differentially expressed genes were enriched in the four oldest evolutionary strata compared with other genomic compartments, and that, within compartments, several signatures of sequence degeneration were greater for differentially expressed than non-differentially expressed genes. Two particular degenerative signatures were significantly associated with lower expression levels within differentially expressed allele pairs: upstream insertion of transposable elements and mutations truncating the protein length. Other degenerative mutations associated with differential expression included nonsynonymous substitutions and altered intron or GC content. The association between differential expression and allele degeneration is relevant for a broad range of taxa where mating compatibility or sex is determined by genes located in large regions where recombination is suppressed.

Digest: Climate plays marginal role for homomorphic sex chromosome differentiation in common frogs†

In systems with early stage sex-chromosome evolution, climate gradients can largely explain changes in the sex-determining systems (i.e., genetic or environmental factors). However, in the common frog Rana temporaria, Phillips et al. found that phylogeography, rather than elevation (used as a proxy for climate), was associated with homomorphic sex-chromosome differentiation levels.

No evidence that Y-chromosome differentiation affects male fitness in a Swiss population of common frogs

The canonical model of sex-chromosome evolution assigns a key role to sexually antagonistic (SA) genes on the arrest of recombination and ensuing degeneration of Y chromosomes. This assumption cannot be tested in organisms with highly differentiated sex chromosomes, such as mammals or birds, owing to the lack of polymorphism. Fixation of SA alleles, furthermore, might be the consequence rather than the cause of recombination arrest. Here we focus on a population of common frogs (Rana temporaria) where XY males with genetically differentiated Y chromosomes (nonrecombinant Y haplotypes) coexist with both XY° males with proto-Y chromosomes (only differentiated from X chromosomes in the immediate vicinity of the candidate sex-determining locus Dmrt1) and XX males with undifferentiated sex chromosomes (genetically identical to XX females). Our study finds no effect of sex-chromosome differentiation on male phenotype, mating success or fathering success. Our conclusions rejoin genomic studies that found no differences in gene expression between XY, XY° and XX males. Sexual dimorphism in common frogs might result more from the differential expression of autosomal genes than from sex-linked SA genes. Among-male variance in sex-chromosome differentiation seems better explained by a polymorphism in the penetrance of alleles at the sex locus, resulting in variable levels of sex reversal (and thus of X-Y recombination in XY females), independent of sex-linked SA genes.

Phylogeography, more than elevation, accounts for sex chromosome differentiation in Swiss populations of the common frog (Rana temporaria)

Sex chromosomes in vertebrates range from highly heteromorphic (as in most birds and mammals) to strictly homomorphic (as in many fishes, amphibians, and nonavian reptiles). Reasons for these contrasted evolutionary trajectories remain unclear, but species such as common frogs with polymorphism in the extent of sex chromosome differentiation may potentially deliver important clues. By investigating 92 common frog populations from a wide range of elevations throughout Switzerland, we show that sex chromosome differentiation strongly correlates with alleles at the candidate sex-determining gene Dmrt1. Y-specific Dmrt1 haplotypes cluster into two main haplogroups, Y_A and Y_B , with a phylogeographic signal that parallels mtDNA haplotypes: Y_A populations, with mostly well-differentiated sex chromosomes, occur primarily south of the main alpine ridge that bisects Switzerland, whereas Y_B populations, with mostly undifferentiated (proto-)sex chromosomes, occur north of this ridge. Elevation has only a marginal effect, opposing previous suggestions of a major role for climate on sex chromosome differentiation. The Y-haplotype effect might result from differences in the penetrance of alleles at the sex-determining locus (such that sex reversal and ensuing X-Y recombination are more frequent in Y_B populations), and/or fixation of an inversion on Y_A (as supported by the empirical observation that Y_A haplotypes might not recombine in XY_A females).

Multi-Tissue Transcriptomes Yield Information on High-Altitude Adaptation and Sex-Determination in Scutiger cf. sikimmensis

The Himalayas are one of earth’s hotspots of biodiversity. Among its many cryptic and undiscovered organisms, including vertebrates, this complex high-mountain ecosystem is expected to harbour many species with adaptations to life in high altitudes. However, modern evolutionary genomic studies in Himalayan vertebrates are still at the beginning. Moreover, in organisms, like most amphibians with relatively high DNA content, whole genome sequencing remains bioinformatically challenging and no complete nuclear genomes are available for Himalayan amphibians. Here, we present the first well-annotated multi-tissue transcriptome of a Greater Himalayan species, the lazy toad Scutiger cf. sikimmensis (Anura: Megophryidae). Applying Illumina NextSeq 500 RNAseq to six tissues, we obtained 41.32 Gb of sequences, assembled to ~111,000 unigenes, translating into 54362 known genes as annotated in seven functional databases. We tested 19 genes, known to play roles in anuran and reptile adaptation to high elevations, and potentially detected diversifying selection for two (TGS1, SENP5) in Scutiger. Of a list of 37 genes, we also identify 27 candidate genes for sex determination or sexual development, all of which providing the first such data for this non-model megophryid species. These transcriptomes will serve as a valuable resource for further studies on amphibian evolution in the Greater Himalaya as a biodiversity hotspot.

Sexual Dimorphism in the Age of Genomics: How, When, Where

In mammals, sex chromosomes start to program autosomal gene expression and epigenetic patterns very soon after fertilization. Yet whether the resulting sex differences are perpetuated throughout development and how they connect to the sex-specific expression patterns in adult tissues is not known. There is a dearth of information on the timing and continuity of sex biases during development. It is also unclear whether sex-specific selection operates during embryogenesis. On the other hand, there is mounting evidence that all adult tissues exhibit sex-specific expression patterns, some of which are independent of hormonal influence and due to intrinsic regulatory effects of the sex chromosome constitution. There are many diseases with origins during embryogenesis that also exhibit sex biases. Epigenetics has provided us with viable mechanisms to explain how the genome stores the memory of developmental events. We propose that some of these marks can be traced back to the sex chromosomes, which interact with the autosomes and establish sex-specific epigenetic features soon after fertilization. Sex-biased epigenetic marks that linger after reprograming may reveal themselves at the transcriptional level at later developmental stages and possibly, throughout the lifespan. Detailed molecular information on the ontogeny of sex biases would also elucidate the sex-specific selective pressures operating on embryos and how compensatory mechanisms evolved to resolve sexual conflict.

Minimal Effects of Proto-Y Chromosomes on House Fly Gene Expression in Spite of Evidence that Selection Maintains Stable Polygenic Sex Determination

Sex determination, the developmental process by which organismal sex is established, evolves fast, often due to changes in the master regulators at the top of the pathway. Additionally, in species with polygenic sex determination, multiple different master regulators segregate as polymorphisms. Understanding the forces that maintain 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 biallelic female determiner. To investigate how natural selection maintains polygenic sex determination in the house fly, we assayed the phenotypic effects of proto-Y chromosomes by performing mRNA-sequencing experiments to measure gene expression in house fly males carrying different proto-Y chromosomes. We find that the proto-Y chromosomes have similar effects as a nonsex-determining autosome. In addition, we created sex-reversed males without any proto-Y chromosomes and they had nearly identical gene expression profiles as genotypic males. Therefore, the proto-Y chromosomes have a minor effect on male gene expression, consistent with previously described minimal X-Y sequence differences. Despite these minimal differences, we find evidence for a disproportionate effect of one proto-Y chromosome on male-biased expression, which could be partially responsible for fitness differences between males with different proto-Y chromosome genotypes. Therefore our results suggest that, if natural selection maintains polygenic sex determination in house fly via gene expression differences, the phenotypes under selection likely depend on a small number of genetic targets.