Unraveling the genotype by environment interaction in a thermosensitive fish with a polygenic sex determination system

In most animals, sex determination occurs at conception, when sex chromosomes are segregated following Mendelian laws. However, in multiple reptiles and fishes, this genetic sex can be overridden by external factors after fertilization or birth. In some species, the genetic sex may also be governed by multiple genes, further limiting our understanding of sex determination in such species. We used the European sea bass (Dicentrarchus labrax) as a model and combined genomic (using a single nucleotide polymorphism chip) and transcriptomic (RNA-Sequencing) approaches to thoroughly depict this polygenic sex determination system and its interaction with temperature. We estimated genetic sex tendency (eGST), defined as the estimated genetic liability to become a given sex under a liability threshold model for sex determination, which accurately predicts the future phenotypic sex. We found evidence that energetic pathways, concerning the regulation of lipids and glucose, are involved in sex determination and could explain why females tend to exhibit higher energy levels and improved growth compared to males. Besides, early exposure to high-temperature up-regulated sox3, followed by sox9a in individuals with intermediate eGST, but not in individuals showing highly female-biased eGST, providing the most parsimonious explanation for temperature-induced masculinization. This gonadal state was maintained likely by DNA methylation and the up-regulation of several genes involved in histone modifications, including jmjd1c. Overall, we describe a sex determination system resulting from continuous genetic and environmental influences in an animal. Our results provide significant progress in our understanding of the mechanisms underlying temperature-induced masculinization in fish.

Polygenic sex determination produces modular sex polymorphism in an African cichlid fish

For many vertebrates, a single genetic locus initiates a cascade of developmental sex differences in the gonad and throughout the organism, resulting in adults with two, phenotypically distinct sexes. Species with polygenic sex determination (PSD) have multiple interacting sex determination alleles segregating within a single species, allowing for more than two genotypic sexes, and scenarios where sex genotype at a given locus can be decoupled from gonadal sex. Here we investigate the effects of PSD on secondary sexual characteristics in the cichlid fish Metriaclima mbenjii, where one female (W) and one male (Y) sex determination allele interact to produce siblings with four possible sex classes: ZZXX females, ZWXX females, ZWXY females, and ZZXY males. We find that PSD in M. mbenjii produces an interplay of sex-linkage and sex-limitation resulting in modular variation in morphological and behavioral traits. Further, the evolution or introgression of a novel sex determiner creates additional axes of phenotypic variation for varied traits, including genital morphology, craniofacial morphology, gastrointestinal morphology, and home tank behaviors. In contrast to single-locus sex determination, which broadly results in sexual dimorphism, polygenic sex determination can induce higher-order sexual polymorphism. The modularity of secondary sexual characteristics produced by PSD provides novel context for understanding the evolutionary causes and consequences of maintenance, gain, or loss of sex determination alleles in populations.

The maintenance of polygenic sex determination depends on the dominance of fitness effects which are predictive of the role of sexual antagonism

In species with polygenic sex determination, multiple male- and female-determining loci on different proto-sex chromosomes segregate as polymorphisms within populations. The extent to which these polymorphisms are at stable equilibria is not yet resolved. Previous work demonstrated that polygenic sex determination is most likely to be maintained as a stable polymorphism when the proto-sex chromosomes have opposite (sexually antagonistic) fitness effects in males and females. However, these models usually consider polygenic sex determination systems with only two proto-sex chromosomes, or they do not broadly consider the dominance of the alleles under selection. To address these shortcomings, I used forward population genetic simulations to identify selection pressures that can maintain polygenic sex determination under different dominance scenarios in a system with more than two proto-sex chromosomes (modeled after the house fly). I found that overdominant fitness effects of male-determining proto-Y chromosomes are more likely to maintain polygenic sex determination than dominant, recessive, or additive fitness effects. The overdominant fitness effects that maintain polygenic sex determination tend to have proto-Y chromosomes with sexually antagonistic effects (male-beneficial and female-detrimental). In contrast, dominant fitness effects that maintain polygenic sex determination tend to have sexually antagonistic multi-chromosomal genotypes, but the individual proto-sex chromosomes do not have sexually antagonistic effects. These results demonstrate that sexual antagonism can be an emergent property of the multi-chromosome genotype without individual sexually antagonistic chromosomes. My results further illustrate how the dominance of fitness effects has consequences for both the likelihood that polygenic sex determination will be maintained as well as the role sexually antagonistic selection is expected to play in maintaining the polymorphism.

Temperature-dependent effects of house fly proto-Y chromosomes on gene expression act independently of the sex determination pathway

Sex determination, the developmental process by which sexually dimorphic phenotypes are established, evolves fast. Species with polygenic sex determination, in which master regulatory genes are found on multiple different proto-sex chromosomes, are informative models to study the evolution of sex determination. House flies are such a model system, with male determining loci possible on all six chromosomes and a female-determiner on one of the chromosomes as well. The distributions of the two most common male-determining proto-Y chromosomes across natural populations suggests that temperature variation is an important selection pressure responsible for maintaining polygenic sex determination in this species. To test that hypothesis, we used RNA-seq to identify temperature-dependent effects of the proto-Y chromosomes on gene expression. We find no evidence for ecologically meaningful temperature-dependent expression of sex determining genes between male genotypes, but we identified hundreds of other genes whose expression depends on the interaction between proto-Y chromosome genotype and temperature. Notably, genes with genotype-by-temperature interactions on expression are not enriched on the proto-sex chromosomes. Moreover, there is no evidence that temperature-dependent expression is driven by chromosome-wide expression divergence between the proto-Y and proto-X alleles. Therefore, if temperature-dependent gene expression is responsible for differences in phenotypes and fitness of proto-Y genotypes across house fly populations, these effects are driven by a small number of temperature-dependent alleles on the proto-Y chromosomes.

The maintenance of polygenic sex determination depends on the dominance of fitness effects which are predictive of the role of sexual antagonism

In species with polygenic sex determination, multiple male- and/or female-determining loci on different proto-sex chromosomes segregate as polymorphisms within populations. The extent to which these polymorphisms are stable equilibria is not yet resolved. Previous work demonstrated that polygenic sex determination is most likely to be maintained as a stable polymorphism when the proto-sex chromosomes have opposite (sexually antagonistic) fitness effects in males and females. However, these models usually consider polygenic sex determination systems with only two proto-sex chromosomes, or they do not broadly consider the dominance of the variants under selection. To address these shortcomings, I used forward population genetic simulations to identify selection pressures that can maintain polygenic sex determination under different dominance scenarios in a system with more than two proto-sex chromosomes (modeled after the house fly). I found that overdominant fitness effects of male-determining proto-Y chromosomes in males are more likely to maintain polygenic sex determination than dominant, recessive, or additive fitness effects. I also found that additive fitness effects that maintain polygenic sex determination have the strongest signatures of sexually antagonistic selection, but there is also some evidence for sexually antagonism when fitness effects of proto-Y chromosomes are dominant or recessive. More generally, these results suggest that the expected effect of sexually antagonistic selection on the maintenance of genetic variation in natural populations will depend on whether the alleles are sex-linked and the dominance of their fitness effects.

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.