Sequence diversity patterns suggesting balancing selection in partially sex‐linked genes of the plant Silene latifolia are not generated by demographic history or gene flow

DNA sequence diversity in genes in the partially sex‐linked pseudoautosomal region (PAR) of the sex chromosomes of the plant Silene latifolia is higher than expected from within‐species diversity of other genes. This could be the footprint of sexually antagonistic (SA) alleles that are maintained by balancing selection in a PAR gene (or genes) and affect polymorphism in linked genome regions. SA selection is predicted to occur during sex chromosome evolution, but it is important to test whether the unexpectedly high sequence polymorphism could be explained without it, purely by the combined effects of partial linkage with the sex‐determining region and the population's demographic history, including possible introgression from Silene dioica. To test this, we applied approximate Bayesian computation‐based model choice to autosomal sequence diversity data, to find the most plausible scenario for the recent history of S. latifolia and then to estimate the posterior density of the most relevant parameters. We then used these densities to simulate variation to be expected at PAR genes. We conclude that an excess of variants at high frequencies at PAR genes should arise in S. latifolia populations only for genes with strong associations with fully sex‐linked genes, which requires closer linkage with the fully sex‐linked region than that estimated for the PAR genes where apparent deviations from neutrality were observed. These results support the need to invoke selection to explain the S. latifolia PAR gene diversity, and encourage further work to test the possibility of balancing selection due to sexual antagonism.     

The potential for sexually antagonistic polymorphism in different genome regions

Sex differences in the fitness effects of alleles at a single locus (intralocus sexual antagonism, or SA) have several evolutionary consequences. Among the consequences of SA, polymorphisms at genes partially linked to the sex-determining region of the sex chromosome pair potentially drive the evolution of suppressed recombination between the sex chromosomes. Understanding the conditions under which SA polymorphism can exist at such pseudo-autosomal (or PAR) loci should increase understanding of the evolution of recombination between sex chromosome pairs, and can help predict when we may expect potentially empirically detectable allele frequency differences between the sexes. Models so far published have concluded that PAR genes can maintain SA polymorphisms over a wider range of selection coefficients than autosomal ones, but have used restrictive assumptions. We expand the modeling of SA alleles at a single locus with the full range of degrees of linkage to the male-specific region, to include strong or weak selection and the possibility of different dominance coefficients in the two sexes. We confirm the previous major conclusion that SA polymorphisms are generally maintained in a larger region of parameter space if the locus is in the PAR than if it is autosomal.     

Male-biased mutation, sex linkage, and the rate of adaptive evolution

An interaction between sex-linked inheritance and sex-biased mutation rates may affect the rate of adaptive evolution. Males have much higher mutation rates than females in several vertebrate and plant taxa. When evolutionary rates are limited by the supply of favorable new mutations, then genes will evolve faster when located on sex chromosomes that spend more time in males. For mutations with additive effects, Y-linked genes evolve fastest, followed by Z-linked genes, autosomal genes, X-linked genes, and finally W-linked and cytoplasmic genes. This ordering can change when mutations show dominance. The predicted differences in substitution rates may be detectable at the molecular level. Male-biased mutation could cause adaptive changes to accumulate more readily on certain kinds of chromosomes and favor animals with Z-W sex determination to have rapidly evolving male sexual displays.     

Testing for the Footprint of Sexually Antagonistic Polymorphisms in the Pseudoautosomal Region of a Plant Sex Chromosome Pair

The existence of sexually antagonistic (SA) polymorphism is widely considered the most likely explanation for the evolution of suppressed recombination of sex chromosome pairs. This explanation is largely untested empirically, and no such polymorphisms have been identified, other than in fish, where no evidence directly implicates these genes in events causing loss of recombination. We tested for the presence of loci with SA polymorphism in the plant Silene latifolia, which is dioecious (with separate male and female individuals) and has a pair of highly heteromorphic sex chromosomes, with XY males. Suppressed recombination between much of the Y and X sex chromosomes evolved in several steps, and the results in Bergero et al. (2013) show that it is still ongoing in the recombining or pseudoautosomal, regions (PARs) of these chromosomes. We used molecular evolutionary approaches to test for the footprints of SA polymorphisms, based on sequence diversity levels in S. latifolia PAR genes identified by genetic mapping. Nucleotide diversity is high for at least four of six PAR genes identified, and our data suggest the existence of polymorphisms maintained by balancing selection in this genome region, since molecular evolutionary (HKA) tests exclude an elevated mutation rate, and other tests also suggest balancing selection. The presence of sexually antagonistic alleles at a locus or loci in the PAR is suggested by the very different X and Y chromosome allele frequencies for at least one PAR gene.     

Sex‐antagonistic genes,XY recombination and feminized Y chromosomes

The canonical model of sex‐chromosome evolution predicts that sex‐antagonistic (SA) genes play an instrumental role in the arrest of XY recombination and ensuing Y chromosome degeneration. Although this model might account for the highly differentiated sex chromosomes of birds and mammals, it does not fit the situation of many lineages of fish, amphibians or nonavian reptiles, where sex chromosomes are maintained homomorphic through occasional XY recombination and/or high turnover rates. Such situations call for alternative explanatory frameworks. A crucial issue at stake is the effect of XY recombination on the dynamics of SA genes and deleterious mutations. Using individual‐based simulations, we show that a complete arrest of XY recombination actually benefits females, not males. Male fitness is maximized at different XY recombination rates depending on SA selection, but never at zero XY recombination. This should consistently favour some level of XY recombination, which in turn generates a recombination load at sex‐linked SA genes. Hill–Robertson interferences with deleterious mutations also impede the differentiation of sex‐linked SA genes, to the point that males may actually fix feminized phenotypes when SA selection and XY recombination are low. We argue that sex chromosomes might not be a good localization for SA genes, and sex conflicts seem better solved through the differential expression of autosomal genes.     

Conflict theory of genomic imprinting in mammals

In mammals, some embryonic genes are expressed differently depending on whether they are inherited from the sperm or egg, a phenomenon known as genomic imprinting. The information on the parental origin is transmitted by an epigenetic mark. Both the molecular mechanisms and evolutionary processes of genomic imprinting have been studied extensively. Here, I illustrate the simplest evolutionary dynamics of imprinting evolution based on the “conflict theory,” by considering the evolution of a gene encoding an embryonic growth factor controlling the maternal resource supply. It demonstrates that (a) the autosomal genes controlling placenta development to modify maternal resource acquisition may evolve a strong asymmetry of gene expression, provided the mother has some chance of accepting multiple males. (b) The genomic imprinting may not evolve if there is a small fraction of recessive deleterious mutations on the gene. (c) The growth-enhancing genes should evolve to paternally expressed, while the growth-suppressing genes should evolve to maternally expressed. (d) The X-linked genes also evolve genomic imprinting, but the main evolutionary force is the sex difference in the optimal embryonic size. I discuss other aberrations that can be explained by the modified versions of the basic model.     

SEX LINKAGE,SEX‐SPECIFIC SELECTION,AND THE ROLE OF RECOMBINATION IN THE EVOLUTION OF SEXUALLY DIMORPHIC GENE EXPRESSION

Sex‐biased genes—genes that are differentially expressed within males and females—are nonrandomly distributed across animal genomes, with sex chromosomes and autosomes often carrying markedly different concentrations of male‐ and female‐biased genes. These linkage patterns are often gene‐ and lineage‐dependent, differing between functional genetic categories and between species. Although sex‐specific selection is often hypothesized to shape the evolution of sex‐linked and autosomal gene content, population genetics theory has yet to account for many of the gene‐ and lineage‐specific idiosyncrasies emerging from the empirical literature. With the goal of improving the connection between evolutionary theory and a rapidly growing body of genome‐wide empirical studies, we extend previous population genetics theory of sex‐specific selection by developing and analyzing a biologically informed model that incorporates sex linkage, pleiotropy, recombination, and epistasis, factors that are likely to vary between genes and between species. Our results demonstrate that sex‐specific selection and sex‐specific recombination rates can generate, and are compatible with, the gene‐ and species‐specific linkage patterns reported in the genomics literature. The theory suggests that sexual selection may strongly influence the architectures of animal genomes, as well as the chromosomal distribution of fixed substitutions underlying sexually dimorphic traits.     

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.     

RECENT GENE‐CAPTURE ON THE UV SEX CHROMOSOMES OF THE MOSS CERATODON PURPUREUS

Sex chromosomes evolve from ordinary autosomes through the expansion and subsequent degeneration of a region of suppressed recombination that is inherited through one sex. Here we investigate the relative timing of these processes in the UV sex chromosomes of the moss Ceratodon purpureus using molecular population genetic analyses of eight newly discovered sex‐linked loci. In this system, recombination is suppressed on both the female‐transmitted (U) sex chromosome and the male‐transmitted (V) chromosome. Genes on both chromosomes therefore should show the deleterious effects of suppressed recombination and sex‐limited transmission, while purifying selection should maintain homologs of genes essential for both sexes on both sex chromosomes. Based on analyses of eight sex‐linked loci, we show that the nonrecombining portions of the U and V chromosomes expanded in at least two events (～0.6–1.3 MYA and ～2.8–3.5 MYA), after the divergence of C. purpureus from its dioecious sister species, Trichodon cylindricus and Cheilothela chloropus. Both U‐ and V‐linked copies showed reduced nucleotide diversity and limited population structure, compared to autosomal loci, suggesting that the sex chromosomes experienced more recent selective sweeps that the autosomes. Collectively these results highlight the dynamic nature of gene composition and molecular evolution on nonrecombining portions of the U and V sex chromosomes.     

EVOLUTIONARY FEEDBACKS BETWEEN REPRODUCTIVE MODE AND MUTATION RATE EXACERBATE THE PARADOX OF SEX

Evolutionary theory suggests that low mutation rates should favor the persistence of asexuals. Additionally, given the observation that most nonneutral mutations are deleterious, asexuality may strengthen selection for reduced mutation rates. This reciprocal relationship raises the possibility of a positive feedback loop between sex and mutation rate. We explored the consequences of this evolutionary feedback with an individual‐based model in which a sexual population is continually challenged by the introduction of asexual clones. We found that asexuals were more likely to spread in a population when mutation rates were able to evolve relative to a model in which mutation rates were held constant. In fact, under evolving mutation rates, asexuals were able to spread to fixation even when sexuals faced no cost of sex whatsoever. The added success of asexuals was the result of their ability to evolve lower mutation rates and thereby slow the process of mutation accumulation that otherwise limited their spread. Given the existence of ample mutation rate variation in natural populations, our findings show that the evolutionary feedback between sex and mutation rate may intensify the “paradox of sex,” supporting the argument that deleterious mutation accumulation alone is likely insufficient to overcome the reproductive advantage of asexual competitors in the short term.     

The horse pseudoautosomal region (PAR): characterization and comparison with the human, chimp and mouse PARs

The pseudoautosomal region (PAR) is a genomic segment on mammalian sex chromosomes where sequence homology mimics that seen between autosomal homologues. The region is essential for pairing and proper segregation of sex chromosomes during male meiosis. As yet, only human/chimp and mouse PARs have been characterized. The two groups of species differ dramatically in gene content and size of the PAR and therefore do not provide clues about the likely evolution and constitution of PAR among mammals. Here we characterize the equine PAR by i) isolating and arranging 71 BACs containing 129 markers (110 STS and 19 genes) into two contigs spanning the region, ii) precisely localizing the pseudoautosomal boundary (PAB), and iii) describing part of the contiguous X- and Y-specific regions. We also report the discovery of an approximately 200 kb region in the middle of the PAR that is present in the male-specific region of the Y (MSY) as well. Such duplication is a novel observation in mammals. Further, comparison of the equine PAR with the human counterpart shows that despite containing orthologs from an additional 1 Mb region beyond the human PAR1, the equine PAR is around 0.9 Mb smaller than the size of the human PAR. We theorize that the PAR varies in size and gene content across evolutionarily closely as well as distantly related mammals. Although striking differences like those observed between human and mouse may be rare, variations similar to those seen between horse and human may be prevalent among mammals.     

Sex‐biased transcriptome divergence along a latitudinal gradient

Sex‐dependent gene expression is likely an important genomic mechanism that allows sex‐specific adaptation to environmental changes. Among Drosophila species, sex‐biased genes display remarkably consistent evolutionary patterns; male‐biased genes evolve faster than unbiased genes in both coding sequence and expression level, suggesting sex differences in selection through time. However, comparatively little is known of the evolutionary process shaping sex‐biased expression within species. Latitudinal clines offer an opportunity to examine how changes in key ecological parameters also influence sex‐specific selection and the evolution of sex‐biased gene expression. We assayed male and female gene expression in Drosophila serrata along a latitudinal gradient in eastern Australia spanning most of its endemic distribution. Analysis of 11 631 genes across eight populations revealed strong sex differences in the frequency, mode and strength of divergence. Divergence was far stronger in males than females and while latitudinal clines were evident in both sexes, male divergence was often population specific, suggesting responses to localized selection pressures that do not covary predictably with latitude. While divergence was enriched for male‐biased genes, there was no overrepresentation of X‐linked genes in males. By contrast, X‐linked divergence was elevated in females, especially for female‐biased genes. Many genes that diverged in D. serrata have homologs also showing latitudinal divergence in Drosophila simulans and Drosophila melanogaster on other continents, likely indicating parallel adaptation in these distantly related species. Our results suggest that sex differences in selection play an important role in shaping the evolution of gene expression over macro‐ and micro‐ecological spatial scales.     

Sex ratio polymorphism: The impact of mutation and drift on evolution

This paper addresses the question, which sex ratio will evolve in a population that is subject to mutation and drift. The problem is analyzed using a simulation model as well as analytical methods. A detailed simulation model for the evolution of a population's allele distribution shows that for the sex ratio game a wide spectrum of different population states may evolve from on the one hand a monomorphic state with one predominant allele and with all other alleles suppressed by the forces of selection, to on the other hand a polymorphism determined by recurrent mutations. Which of these states will evolve depends on the population size, the mating system and the rate of mutations. For the sex ratio game the evolutionary stable strategy (ESS), as defined by evolutionary game theory, can only predict the population sex ratio but not the underlying stable population state. A comparison of different approaches to the problem shows that false predictions of the stable population states might result from two simplifying assumptions that are fairly common in evolutionary biology: a) it is assumed that mutations are rare events and there is never more than one mutant gene present in a population at any one time; b) a deterministic relationship is assumed between the fitness assigned to an individual's strategy and the individual's contribution to the gene pool of future generations.     

EVOLUTIONARILY STABLE SEX RATIOS AND MUTATION LOAD

Frequency‐dependent selection should drive dioecious populations toward a 1:1 sex ratio, but biased sex ratios are widespread, especially among plants with sex chromosomes. Here, we develop population genetic models to investigate the relationships between evolutionarily stable sex ratios, haploid selection, and deleterious mutation load. We confirm that when haploid selection acts only on the relative fitness of X‐ and Y‐bearing pollen and the sex ratio is controlled by the maternal genotype, seed sex ratios evolve toward 1:1. When we also consider haploid selection acting on deleterious mutations, however, we find that biased sex ratios can be stably maintained, reflecting a balance between the advantages of purging deleterious mutations via haploid selection, and the disadvantages of haploid selection on the sex ratio. Our results provide a plausible evolutionary explanation for biased sex ratios in dioecious plants, given the extensive gene expression that occurs across plant genomes at the haploid stage.     

The Human Pseudoautosomal Region (PAR): Origin, Function and Future

The pseudoautosomal regions (PAR1 and PAR2) of the human X and Y chromosomes pair and recombine during meiosis. Thus genes in this region are not inherited in a strictly sex-linked fashion. PAR1 is located at the terminal region of the short arms and PAR2 at the tips of the long arms of these chromosomes. To date, 24 genes have been assigned to the PAR1 region. Half of these have a known function. In contrast, so far only 4 genes have been discovered in the PAR2 region. Deletion of the PAR1 region results in failure of pairing and male sterility. The gene SHOX (short stature homeobox-containing) resides in PAR1. SHOX haploinsufficiency contributes to certain features in Turner syndrome as well as the characteristics of Leri-Weill dyschondrosteosis. Only two of the human PAR1 genes have mouse homologues. These do not, however, reside in the mouse PAR1 region but are autosomal. The PAR regions seem to be relics of differential additions, losses, rearrangements and degradation of the X and Y chromosome in different mammalian lineages. Marsupials have three homologues of human PAR1 genes in their autosomes, although, in contrast to mouse, do not have a PAR region at all. The disappearance of PAR from other species seems likely and this region will only be rescued by the addition of genes to both X and Y, as has occurred already in lemmings. The present review summarizes the current understanding of the evolution of PAR and provides up-to-date information about individual genes residing in this region.     

A comparative view on sex determination in medaka

In fish, an amazing variety of sex determination mechanisms are known, ranging from hermaphroditism to gonochorism and from environmental to genetic sex determination. This makes fish especially suited for studying sex determination from the evolutionary point of view. In several fish groups, different sex determination mechanisms are found in closely related species, and evolution of this process is still ongoing in recent organisms. The medaka (Oryzias latipes) has an XY-XX genetic sex determination system. The Y-chromosome in this species is at an early stage of evolution. The molecular differences between X and Y are only very subtle and the Y-specific segment is very small. The sex-determining region has accumulated duplicated sequences from elsewhere in the genome, leading to recombinational isolation. The region contains a candidate for the male sex-determining gene named dmrt1bY. This gene arose through duplication of an autosomal chromosome fragment of linkage group 9. While all other genes degenerated, dmrt1bY is the only functional gene in the Y-specific region. The duplication leading to dmrt1bY occurred recently during evolution of the genus Oryzias. This suggests that different genes might be the master sex-determining gene in other fish.     

Analysis of phage resistance in Staphylococcus aureus SA003 reveals different binding mechanisms for the closely related Twort-like phages ɸSA012 and ɸSA039

We have previously generated strains of Staphylococcus aureus SA003 resistant to its specific phage ɸSA012 through a long-term coevolution experiment. However, the DNA mutations responsible for the phenotypic change of phage resistance are unknown. Whole-genome analysis revealed eight genes that acquired mutations: six point mutations (five missense mutations and one nonsense mutation) and two deletions. Complementation of the phage-resistant strains by the wild-type alleles showed that five genes were linked to phage adsorption of ɸSA012, and two mutated host genes were linked to the inhibition of post-adsorption. Unlike ɸSA012, infection by ɸSA039, a close relative of ɸSA012, onto early coevolved phage-resistant SA003 (SA003R2) was impaired drastically. Here, we identified that ɸSA012 and ɸSA039 adsorb to the cell surface S. aureus SA003 through a different mechanism. ɸSA012 requires the backbone of wall teichoic acids (WTA), while ɸSA039 requires both backbone and the β-GlcNAc residue. In silico analysis of the ɸSA039 genome revealed that several proteins in the tail and baseplate region were different from ɸSA012. The difference in tail and baseplate proteins might be the factor for specificity difference between ɸSA012 and ɸSA039.

     

Nuclear–mitochondrial epistasis: a gene's eye view of genomic conflict

We use population genetic models to investigate the cooperative and conflicting synergistic fitness effects between genes from the nucleus and the mitochondrion. By varying fitness parameters, we examine the scope for conflict relative to cooperation among genomes and the utility of the “gene's eye view” analytical approach, which is based on the marginal average fitness of specific alleles. Because sexual conflict can maintain polymorphism of mitochondrial haplotypes, we can explore two types of evolutionary conflict (genomic and sexual) with one epistatic model. We find that the nuclear genetic architecture (autosomal, X‐linked, or Z‐linked) and the mating system change the regions of parameter space corresponding to the evolution by sexual and genomic conflict. For all models, regardless of conflict or cooperation, we find that population mean fitness increases monotonically as evolution proceeds. Moreover, we find that the process of gene frequency change with positive, synergistic fitnesses is self‐accelerating, as the success of an allele in one genome or in one sex increases the frequency of the interacting allele upon which its success depends. This results in runaway evolutionary dynamics caused by the positive intergenomic associations generated by selection. An inbreeding mating system tends to further accelerate these runaway dynamics because it maintains favorable host–symbiont or male–female gene combinations. In contrast, where conflict predominates, the success of an allele in one genome or in one sex diminishes the frequency of the corresponding allele in the other, resulting in considerably slower evolutionary dynamics. The rate of change of mean fitness is also much faster with positive, synergistic fitnesses and much slower where conflict is predominant. Consequently, selection rapidly fixes cooperative gene combinations, while leaving behind a slowing evolving residue of conflicting gene combinations at mutation–selection balance. We discuss how an emphasis on marginal fitness averages may obscure the interdependence of allelic fitness across genomes, making the evolutionary trajectories appear independent of one another when they are not.