Is Multifactorial Sex Determination in the House Fly, Musca domestica (L.), Stable Over Time?

Turnover of sex chromosomes in the Lake Tanganyika cichlid tribe Tropheini (Teleostei: Cichlidae)

Sex chromosome replacement is frequent in many vertebrate clades, including fish, frogs, and lizards. In order to understand the mechanisms responsible for sex chromosome turnover and the early stages of sex chromosome divergence, it is necessary to study lineages with recently evolved sex chromosomes. Here we examine sex chromosome evolution in a group of African cichlid fishes (tribe Tropheini) which began to diverge from one another less than 4 MYA. We have evidence for a previously unknown sex chromosome system, and preliminary indications of several additional systems not previously reported in this group. We find a high frequency of sex chromosome turnover and estimate a minimum of 14 turnovers in this tribe. We date the origin of the most common sex determining system in this tribe (XY-LG5/19) near the base of one of two major sub-clades of this tribe, about 3.4 MY ago. Finally, we observe variation in the size of one sex-determining region that suggests independent evolution of evolutionary strata in species with a shared sex-determination system. Our results illuminate the rapid rate of sex chromosome turnover in the tribe Tropheini and set the stage for further studies of the dynamics of sex chromosome evolution in this group.

Reorganization of the ancestral sex-determining regions during the evolution of trioecy in Pleodorina starrii

The coexistence of three sexual phenotypes (male, female and bisexual) in a single species, 'trioecy', is rarely found in diploid organisms such as flowering plants and invertebrates. However, trioecy in haploid organisms has only recently been reported in a green algal species, Pleodorina starrii. Here, we generated whole-genome data of the three sex phenotypes of P. starrii to reveal a reorganization of the ancestral sex-determining regions (SDRs) in the sex chromosomes: the male and bisexual phenotypes had the same "male SDR" with paralogous gene expansions of the male-determining gene MID, whereas the female phenotype had a "female SDR" with transposition of the female-specific gene FUS1 to autosomal regions. Although the male and bisexual sex phenotypes had the identical male SDR and harbored autosomal FUS1, MID and FUS1 expression during sexual reproduction differed between them. Thus, the coexistence of three sex phenotypes in P. starrii is possible.

Divergent evolution of genetic sex determination mechanisms along environmental gradients

Sex determination (SD) is a crucial developmental process, but its molecular underpinnings are very diverse, both between and within species. SD mechanisms have traditionally been categorized as either genetic (GSD) or environmental (ESD), depending on the type of cue that triggers sexual differentiation. However, mixed systems, with both genetic and environmental components, are more prevalent than previously thought. Here, we show theoretically that environmental effects on expression levels of genes within SD regulatory mechanisms can easily trigger within-species evolutionary divergence of SD mechanisms. This may lead to the stable coexistence of multiple SD mechanisms and to spatial variation in the occurrence of different SD mechanisms along environmental gradients. We applied the model to the SD system of the housefly, a global species with world-wide latitudinal clines in the frequencies of different SD systems, and found that it correctly predicted these clines if specific genes in the housefly SD system were assumed to have temperature-dependent expression levels. We conclude that environmental sensitivity of gene regulatory networks may play an important role in diversification of SD mechanisms.

X-chromosome target specificity diverged between dosage compensation mechanisms of two closely related Caenorhabditis species

An evolutionary perspective enhances our understanding of biological mechanisms. Comparison of sex determination and X-chromosome dosage compensation mechanisms between the closely related nematode species Caenorhabditis briggsae (Cbr) and Caenorhabditis elegans (Cel) revealed that the genetic regulatory hierarchy controlling both processes is conserved, but the X-chromosome target specificity and mode of binding for the specialized condensin dosage compensation complex (DCC) controlling X expression have diverged. We identified two motifs within Cbr DCC recruitment sites that are highly enriched on X: 13 bp MEX and 30 bp MEX II. Mutating either MEX or MEX II in an endogenous recruitment site with multiple copies of one or both motifs reduced binding, but only removing all motifs eliminated binding in vivo. Hence, DCC binding to Cbr recruitment sites appears additive. In contrast, DCC binding to Cel recruitment sites is synergistic: mutating even one motif in vivo eliminated binding. Although all X-chromosome motifs share the sequence CAGGG, they have otherwise diverged so that a motif from one species cannot function in the other. Functional divergence was demonstrated in vivo and in vitro. A single nucleotide position in Cbr MEX can determine whether Cel DCC binds. This rapid divergence of DCC target specificity could have been an important factor in establishing reproductive isolation between nematode species and contrasts dramatically with the conservation of target specificity for X-chromosome dosage compensation across Drosophila species and for transcription factors controlling developmental processes such as body-plan specification from fruit flies to mice.

Strong variation in frequencies of male and female determiners between neighboring housefly populations

Sex-determination mechanisms evolve rapidly and vary between species. Occasionally, polymorphic systems are found, like in the housefly. Studying the dynamics and stability of such systems can provide a better understanding of the evolution of sex-determination systems. In the housefly, dominant male-determining loci (M) can lie not only on the Y chromosome (M^Y ), but also on autosomes (M^A ) or the X chromosome (M^X ). M enforces male development by inhibiting the female-determining gene transformer (tra). A mutant tra allele, tra^D , is insensitive to M and is a dominant female determiner. M^Y prevails at high latitudes and polymorphic M loci together with tra^D at low latitudes. To get more insight into the stability and frequencies of these sex determiners with mutually exclusive dominance, we investigated 5 regional Spanish populations. We found strong variation among populations. Two populations with hemizygous M^III were found, 2 contained homozygous M^X with additional hemizygous M^I and M^II in 1 population. One population contained homozygous and hemizygous M^X with additionally hemizygous M^II . All females in populations with homozygous M, had tra^D , whereas no tra^D was found in populations without homozygous M. Our results indicate locally stable systems may either harbor a single hemizygous M and no tra^D , corresponding to a male heterogametic system, or homozygous and/or multiple M and heterozygous tra^D , reminiscent of a female heterogametic system. They support that M loci can accumulate in the presence of a dominant female determiner. Limited migration between populations may contribute to the stability of these systems.

Thermal tolerance and preference are both consistent with the clinal distribution of house fly proto-Y chromosomes

Selection pressures can vary within localized areas and across massive geographical scales. Temperature is one of the best studied ecologically variable abiotic factors that can affect selection pressures across multiple spatial scales. Organisms rely on physiological (thermal tolerance) and behavioral (thermal preference) mechanisms to thermoregulate in response to environmental temperature. In addition, spatial heterogeneity in temperatures can select for local adaptation in thermal tolerance, thermal preference, or both. However, the concordance between thermal tolerance and preference across genotypes and sexes within species and across populations is greatly understudied. The house fly, Musca domestica, is a well-suited system to examine how genotype and environment interact to affect thermal tolerance and preference. Across multiple continents, house fly males from higher latitudes tend to carry the male-determining gene on the Y chromosome, whereas those from lower latitudes usually have the male determiner on the third chromosome. We tested whether these two male-determining chromosomes differentially affect thermal tolerance and preference as predicted by their geographical distributions. We identify effects of genotype and developmental temperature on male thermal tolerance and preference that are concordant with the natural distributions of the chromosomes, suggesting that temperature variation across the species range contributes to the maintenance of the polymorphism. In contrast, female thermal preference is bimodal and largely independent of congener male genotypes. These sexually dimorphic thermal preferences suggest that temperature-dependent mating dynamics within populations could further affect the distribution of the two chromosomes. Together, the differences in thermal tolerance and preference across sexes and male genotypes suggest that different selection pressures may affect the frequencies of the male-determining chromosomes across different spatial scales.

Evolutionary transition to XY sex chromosomes associated with Y-linked duplication of a male hormone gene in a terrestrial isopod

Sex chromosomes are highly variable in some taxonomic groups, but the evolutionary mechanisms underlying this diversity are not well understood. In terrestrial isopod crustaceans, evolutionary turnovers in sex chromosomes are frequent, possibly caused by Wolbachia, a vertically-transmitted endosymbiont causing male-to-female sex reversal. Here, we use surgical manipulations and genetic crosses, plus genome sequencing, to examine sex chromosomes in the terrestrial isopod Trachelipus rathkei. Although an earlier cytogenetics study suggested a ZZ/ZW sex chromosome system in this species, we surprisingly find multiple lines of evidence that in our study population, sex is determined by an XX/XY system. Consistent with a recent evolutionary origin for this XX/XY system, the putative male-specific region of the genome is small. The genome shows evidence of Y-linked duplications of the gene encoding the androgenic gland hormone, a major component of male sexual differentiation in isopods. Our analyses also uncover sequences horizontally acquired from past Wolbachia infections, consistent with the hypothesis that Wolbachia may have interfered with the evolution of sex determination in T. rathkei. Overall, these results provide evidence for the co-occurrence of multiple sex chromosome systems within T. rathkei, further highlighting the relevance of terrestrial isopods as models for the study of sex chromosome evolution.

Temperature-dependent effects of house fly proto-Y chromosomes on gene expression could be responsible for fitness differences that maintain polygenic sex determination

Sex determination, the developmental process by which sexually dimorphic phenotypes are established, evolves fast. Evolutionary turnover in a sex determination pathway may occur via selection on alleles that are genetically linked to a new master sex determining locus on a newly formed proto-sex chromosome. 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 and sex chromosomes. 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 two most common male-determining proto-Y chromosomes form latitudinal clines on multiple continents, suggesting that temperature variation is an important selection pressure responsible for maintaining polygenic sex determination in this species. Temperature-dependent fitness effects could be manifested through temperature-dependent gene expression differences across proto-Y chromosome genotypes. These gene expression differences may be the result of cis regulatory variants that affect the expression of genes on the proto-sex chromosomes, or trans effects of the proto-Y chromosomes on genes elswhere in the genome. We used RNA-seq to identify genes whose expression depends on proto-Y chromosome genotype and temperature in adult male house flies. We found no evidence for ecologically meaningful temperature-dependent expression differences of sex determining genes between male genotypes, but we were probably not sampling an appropriate developmental time-point to identify such effects. In contrast, we identified many other genes whose expression depends on the interaction between proto-Y chromosome genotype and temperature, including genes that encode proteins involved in reproduction, metabolism, lifespan, stress response, and immunity. Notably, genes with genotype-by-temperature interactions on expression were not enriched on the proto-sex chromosomes. Moreover, there was no evidence that temperature-dependent expression is driven by chromosome-wide cis-regulatory 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 that may have trans effects on the expression of genes on other chromosomes.

Gene-Level, but Not Chromosome-Wide, Divergence between a Very Young House Fly Proto-Y Chromosome and Its Homologous Proto-X Chromosome

X and Y chromosomes are usually derived from a pair of homologous autosomes, which then diverge from each other over time. Although Y-specific features have been characterized in sex chromosomes of various ages, the earliest stages of Y chromosome evolution remain elusive. In particular, we do not know whether early stages of Y chromosome evolution consist of changes to individual genes or happen via chromosome-scale divergence from the X. To address this question, we quantified divergence between young proto-X and proto-Y chromosomes in the house fly, Musca domestica. We compared proto-sex chromosome sequence and gene expression between genotypic (XY) and sex-reversed (XX) males. We find evidence for sequence divergence between genes on the proto-X and proto-Y, including five genes with mitochondrial functions. There is also an excess of genes with divergent expression between the proto-X and proto-Y, but the number of genes is small. This suggests that individual proto-Y genes, but not the entire proto-Y chromosome, have diverged from the proto-X. We identified one gene, encoding an axonemal dynein assembly factor (which functions in sperm motility), that has higher expression in XY males than XX males because of a disproportionate contribution of the proto-Y allele to gene expression. The upregulation of the proto-Y allele may be favored in males because of this gene's function in spermatogenesis. The evolutionary divergence between proto-X and proto-Y copies of this gene, as well as the mitochondrial genes, is consistent with selection in males affecting the evolution of individual genes during early Y chromosome evolution.

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 (PSD), 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 PSD 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 PSD 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 PSD 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 PSD than dominant, recessive, or additive fitness effects. The overdominant fitness effects that maintain PSD tend to have proto-Y chromosomes with sexually antagonistic effects (male-beneficial and female-detrimental). In contrast, dominant fitness effects that maintain PSD 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 PSD will be maintained as well as the role sexually antagonistic selection is expected to play in maintaining the polymorphism.

Genetic evaluation and characterization of behavioral resistance to imidacloprid in the house fly

Insecticide resistance in pest populations is an increasing problem in both urban and rural settings due to over-application of insecticides and lack of rotation among insecticidal chemical classes. The house fly (Musca domestica L.) is a cosmopolitan pest fly species implicated in the transmission of numerous pathogens. The evolution of insecticide resistance long has been documented in house flies, with resistance reported to all major insecticide classes. House fly resistance to imidacloprid, the most widely used neonicotinoid insecticide available for fly control, has evolved in field populations through both physiological and behavioral mechanisms. Previous studies have characterized and mapped the genetic changes that confer physiological resistance to imidacloprid, but no study have examined the genetics involved in behavioral resistance to imidacloprid to date. In the current study, several approaches were utilized to characterize the genetics and inheritance of behavioral resistance to imidacloprid in the house fly. These include behavioral observation analyses, preference assays, and the use of genetic techniques for the identification of house fly chromosome(s) carrying factors. Behavioral resistance was mapped to autosomes 1 and 4. Inheritance of resistance was shown to be neither fully dominant nor recessive. Factors on autosomes 1 and 4 independently conferred contact-dependent avoidance of imidacloprid and a feeding preference for sugar alone or for sugar with dinotefuran, another neonicotinoid insecticide, over imidacloprid. This study serves as the first linkage analysis of a behavioral trait in the house fly, and provides new avenues for research regarding inherited behavior in the house fly and other animals.

Network architecture and sex chromosome turnovers: Do epistatic interactions shape patterns of sex chromosome replacement?

Recent studies have revealed an astonishing diversity of sex chromosomes in many vertebrate lineages, prompting questions about the mechanisms of sex chromosome turnover. While there is considerable population genetic theory about the evolutionary forces promoting sex chromosome replacement, this theory has not yet been integrated with our understanding of the molecular and developmental genetics of sex determination. Here, we review recent data to examine four questions about how the structure of gene networks influences the evolution of sex determination. We argue that patterns of epistasis, arising from the structure of genetic networks, may play an important role in regulating the rates and patterns of sex chromosome replacement.

How to identify sex chromosomes and their turnover

Although sex is a fundamental component of eukaryotic reproduction, the genetic systems that control sex determination are highly variable. In many organisms the presence of sex chromosomes is associated with female or male development. Although certain groups possess stable and conserved sex chromosomes, others exhibit rapid sex chromosome evolution, including transitions between male and female heterogamety, and turnover in the chromosome pair recruited to determine sex. These turnover events have important consequences for multiple facets of evolution, as sex chromosomes are predicted to play a central role in adaptation, sexual dimorphism, and speciation. However, our understanding of the processes driving the formation and turnover of sex chromosome systems is limited, in part because we lack a complete understanding of interspecific variation in the mechanisms by which sex is determined. New bioinformatic methods are making it possible to identify and characterize sex chromosomes in a diverse array of non-model species, rapidly filling in the numerous gaps in our knowledge of sex chromosome systems across the tree of life. In turn, this growing data set is facilitating and fueling efforts to address many of the unanswered questions in sex chromosome evolution. Here, we synthesize the available bioinformatic approaches to produce a guide for characterizing sex chromosome system and identity simultaneously across clades of organisms. Furthermore, we survey our current understanding of the processes driving sex chromosome turnover, and highlight important avenues for future research.

Hemimetabolous insects elucidate the origin of sexual development via alternative splicing

Insects are the only known animals in which sexual differentiation is controlled by sex-specific splicing. The doublesex transcription factor produces distinct male and female isoforms, which are both essential for sex-specific development. dsx splicing depends on transformer, which is also alternatively spliced such that functional Tra is only present in females. This pathway has evolved from an ancestral mechanism where dsx was independent of tra and expressed and required only in males. To reconstruct this transition, we examined three basal, hemimetabolous insect orders: Hemiptera, Phthiraptera, and Blattodea. We show that tra and dsx have distinct functions in these insects, reflecting different stages in the changeover from a transcription-based to a splicing-based mode of sexual differentiation. We propose that the canonical insect tra-dsx pathway evolved via merger between expanding dsx function (from males to both sexes) and narrowing tra function (from a general splicing factor to dedicated regulator of dsx).

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.

Tsetse fly evolution, genetics and the trypanosomiases - A review

This reviews work published since 2007. Relative efforts devoted to the agents of African trypanosomiasis and their tsetse fly vectors are given by the numbers of PubMed accessions. In the last 10 years PubMed citations number 3457 for Trypanosoma brucei and 769 for Glossina. The development of simple sequence repeats and single nucleotide polymorphisms afford much higher resolution of Glossina and Trypanosoma population structures than heretofore. Even greater resolution is offered by partial and whole genome sequencing. Reproduction in T. brucei sensu lato is principally clonal although genetic recombination in tsetse salivary glands has been demonstrated in T. b. brucei and T. b. rhodesiense but not in T. b. gambiense. In the past decade most genetic attention was given to the chief human African trypanosomiasis vectors in subgenus Nemorhina e.g., Glossina f. fuscipes, G. p. palpalis, and G. p. gambiense. The chief interest in Nemorhina population genetics seemed to be finding vector populations sufficiently isolated to enable efficient and long-lasting suppression. To this end estimates were made of gene flow, derived from F_ST and its analogues, and Ne, the size of a hypothetical population equivalent to that under study. Genetic drift was greater, gene flow and Ne typically lesser in savannah inhabiting tsetse (subgenus Glossina) than in riverine forms (Nemorhina). Population stabilities were examined by sequential sampling and genotypic analysis of nuclear and mitochondrial genomes in both groups and found to be stable. Gene frequencies estimated in sequential samplings differed by drift and allowed estimates of effective population numbers that were greater for Nemorhina spp than Glossina spp. Prospects are examined of genetic methods of vector control. The tsetse long generation time (c. 50 d) is a major contraindication to any suggested genetic method of tsetse population manipulation. Ecological and modelling research convincingly show that conventional methods of targeted insecticide applications and traps/targets can achieve cost-effective reduction in tsetse densities.

Drift-Induced Selection Between Male and Female Heterogamety

Evolutionary transitions between male and female heterogamety are common in both vertebrates and invertebrates. Theoretical studies of these transitions have found that, when all genotypes are equally fit, continuous paths of intermediate equilibria link the two sex chromosome systems. This observation has led to a belief that neutral evolution along these paths can drive transitions, and that arbitrarily small fitness differences among sex chromosome genotypes can determine the system to which evolution leads. Here, we study stochastic evolutionary dynamics along these equilibrium paths. We find non-neutrality, both in transitions retaining the ancestral pair of sex chromosomes, and in those creating a new pair. In fact, substitution rates are biased in favor of dominant sex determining chromosomes, which fix with higher probabilities than mutations of no effect. Using diffusion approximations, we show that this non-neutrality is a result of "drift-induced selection" operating at every point along the equilibrium paths: stochastic jumps off the paths return with, on average, a directional bias in favor of the dominant segregating sex chromosome. Our results offer a novel explanation for the observed preponderance of dominant sex determining genes, and hint that drift-induced selection may be a common force in standard population genetic systems.

The Reversible Sex of Gonochoristic Fish: Insights and Consequences

Fish sex reversal is a means to understand sex determination and differentiation, but it is also used to control sex in aquaculture. This review discusses sex reversal in gonochoristic fish, with the coexistence of genetic and environmental influences. The different periods of fish sensitivity to sex reversal treatments are presented with the mechanisms implicated. The old players of sex differentiation are revisited with transcriptome data and loss of function studies following hormone- or temperature-induced sex reversal. We also discuss whether cortisol is the universal mediator of sex reversal in fish due to its implication in ovarian meiosis and 11KT increase. The large plasticity in fish for sex reversal is also evident in the brain, with a reversibility existing even in adulthood. Studies on epigenetics are presented, since it links the environment, gene expression, and sex reversal, notably the association of DNA methylation in sex reversal. Manipulations with exogenous factors reverse the primary sex in many fish species under controlled conditions, but several questions arise on whether this can occur under wild conditions and what is the ecological significance. Cases of sex reversal in wild fish populations are shown and their fitness and future perspectives are discussed.