Sex-Specific Selection Drives the Evolution of Alternative Splicing in Birds.

Males and females of the same species share the majority of their genomes, yet they are frequently exposed to conflicting selection pressures. Gene regulation is widely assumed to resolve these conflicting sex-specific selection pressures, and although there has been considerable focus on elucidating the role of gene expression level in sex-specific adaptation, other regulatory mechanisms have been overlooked. Alternative splicing enables different transcripts to be generated from the same gene, meaning that exons which have sex-specific beneficial effects can in theory be retained in the gene product, whereas exons with detrimental effects can be skipped. However, at present, little is known about how sex-specific selection acts on broad patterns of alternative splicing. Here, we investigate alternative splicing across males and females of multiple bird species. We identify hundreds of genes that have sex-specific patterns of splicing and establish that sex differences in splicing are correlated with phenotypic sex differences. Additionally, we find that alternatively spliced genes have evolved rapidly as a result of sex-specific selection and suggest that sex differences in splicing offer another route to sex-specific adaptation when gene expression level changes are limited by functional constraints. Overall, our results shed light on how a diverse transcriptional framework can give rise to the evolution of phenotypic sexual dimorphism.

Gene duplication to the Y chromosome in Trinidadian Guppies.

Differences in allele frequencies at autosomal genes between males and females in a population can result from two scenarios. First, unresolved sexual conflict over survival can produce allelic differentiation between the sexes. However, given the substantial mortality costs required to produce allelic differences between males and females at each generation, it remains unclear how many loci within the genome experience significant sexual conflict over survival. Alternatively, recent studies have shown that similarity between autosomal and Y sequences can create perceived allelic differences between the sexes. However, Y duplications are most likely in species with large nonrecombining regions, in part because they simply represent larger targets for duplications. We assessed the genomes of 120 wild-caught guppies, which experience extensive predation- and pathogen-induced mortality and have a relatively small ancestral Y chromosome. We identified seven autosomal genes that show allelic differences between male and female adults. Five of these genes show clear evidence of whole or partial gene duplication between the Y chromosome and the autosomes. The remaining two genes show evidence of partial homology to the Y. Overall, our findings suggest that the guppy genome experiences a very low level of unresolved sexual conflict over survival, and instead the Y chromosome, despite its small ancestral size and recent origin, may nonetheless accumulate genes with male-specific functions.

Comparison of methodological approaches to the study of young sex chromosomes: A case study in Poecilia.

Studies of sex chromosome systems at early stages of divergence are key to understanding the initial process and underlying causes of recombination suppression. However, identifying signatures of divergence in homomorphic sex chromosomes can be challenging due to high levels of sequence similarity between the X and the Y. Variations in methodological precision and underlying data can make all the difference between detecting subtle divergence patterns or missing them entirely. Recent efforts to test for X-Y sequence differentiation in the guppy have led to contradictory results. Here, we apply different analytical methodologies to the same data set to test for the accuracy of different approaches in identifying patterns of sex chromosome divergence in the guppy. Our comparative analysis reveals that the most substantial source of variation in the results of the different analyses lies in the reference genome used. Analyses using custom-made genome assemblies for the focal population or species successfully recover a signal of divergence across different methodological approaches. By contrast, using the distantly related Xiphophorus reference genome results in variable patterns, due to both sequence evolution and structural variations on the sex chromosomes between the guppy and Xiphophorus. Changes in mapping and filtering parameters can additionally introduce noise and obscure the signal. Our results illustrate how analytical differences can alter perceived results and we highlight best practices for the study of nascent sex chromosomes.

Extreme heterogeneity in sex chromosome differentiation and dosage compensation in livebearers.

Once recombination is halted between the X and Y chromosomes, sex chromosomes begin to differentiate and transition to heteromorphism. While there is a remarkable variation across clades in the degree of sex chromosome divergence, far less is known about the variation in sex chromosome differentiation within clades. Here, we combined whole-genome and transcriptome sequencing data to characterize the structure and conservation of sex chromosome systems across Poeciliidae, the livebearing clade that includes guppies. We found that the Poecilia reticulata XY system is much older than previously thought, being shared not only with its sister species, Poecilia wingei, but also with Poecilia picta, which diverged roughly 20 million years ago. Despite the shared ancestry, we uncovered an extreme heterogeneity across these species in the proportion of the sex chromosome with suppressed recombination, and the degree of Y chromosome decay. The sex chromosomes in P. reticulata and P. wingei are largely homomorphic, with recombination in the former persisting over a substantial fraction. However, the sex chromosomes in P. picta are completely nonrecombining and strikingly heteromorphic. Remarkably, the profound degradation of the ancestral Y chromosome in P. picta is counterbalanced by the evolution of functional chromosome-wide dosage compensation in this species, which has not been previously observed in teleost fish. Our results offer important insight into the initial stages of sex chromosome evolution and dosage compensation.

Multi-Copy Gene Family Evolution on the Avian W Chromosome.

The sex chromosomes often follow unusual evolutionary trajectories. In particular, the sex-limited chromosomes frequently exhibit a small but unusual gene content in numerous species, where many genes have undergone massive gene amplification. The reasons for this remain elusive with a number of recent studies implicating meiotic drive, sperm competition, genetic drift, and gene conversion in the expansion of gene families. However, our understanding is primarily based on Y chromosome studies as few studies have systematically tested for copy number variation on W chromosomes. Here, we conduct a comprehensive investigation into the abundance, variability, and evolution of ampliconic genes on the avian W. First, we quantified gene copy number and variability across the duck W chromosome. We find a limited number of gene families as well as conservation in W-linked gene copy number across duck breeds, indicating that gene amplification may not be such a general feature of sex chromosome evolution as Y studies would initially suggest. Next, we investigated the evolution of HINTW, a prominent ampliconic gene family hypothesized to play a role in female reproduction and oogenesis. In particular, we investigated the factors driving the expansion of HINTW using contrasts between modern chicken and duck breeds selected for different female-specific selection regimes and their wild ancestors. Although we find the potential for selection related to fecundity in explaining small-scale gene amplification of HINTW in the chicken, purifying selection seems to be the dominant mode of evolution in the duck. Together, this challenges the assumption that HINTW is key for female fecundity across the avian phylogeny.

ngsLD: evaluating linkage disequilibrium using genotype likelihoods.

MOTIVATION: Linkage disequilibrium (LD) measures the correlation between genetic loci and is highly informative for association mapping and population genetics. As many studies rely on called genotypes for estimating LD, their results can be affected by data uncertainty, especially when employing a low read depth sequencing strategy. Furthermore, there is a manifest lack of tools for the analysis of large-scale, low-depth and short-read sequencing data from non-model organisms with limited sample sizes. RESULTS: ngsLD addresses these issues by estimating LD directly from genotype likelihoods in a fast, reliable and user-friendly implementation. This method makes use of the full information available from sequencing data and provides accurate estimates of linkage disequilibrium patterns compared with approaches based on genotype calling. We conducted a case study to investigate how LD decays over physical distance in two avian species. AVAILABILITY AND IMPLEMENTATION: The methods presented in this work were implemented in C/C and are freely available for non-commercial use from https://github.com/fgvieira/ngsLD. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Signature of sexual conflict is actually conflict resolved.

There has been substantial interest of late in using population genetic methods to study sexual conflict, where an allele increases the fitness of one sex at some cost to the other (Mank, 2017). Population genomic scans for sexual conflict offer an important advance given the difficulties of identifying antagonistic alleles from more traditional methods, and could greatly increase our understanding of the extent and loci of sexual conflict. This is particularly true for studies in natural populations, for which obtaining accurate fitness measurements for each sex can be challenging. In this issue of Molecular Ecology, Bissegger, Laurentino, Roesti, and Berner (2019) present a cautionary tale about how to interpret these population genomic data.

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.

Phenotypic sexual dimorphism is associated with genomic signatures of resolved sexual conflict.

Intralocus sexual conflict, where an allele benefits one sex at the expense of the other, has an important role in shaping genetic diversity of populations through balancing selection. However, the potential for mating systems to exert balancing selection through sexual conflict on the genome remains unclear. Furthermore, the nature and potential for resolution of sexual conflict across the genome has been hotly debated. To address this, we analysed de novo transcriptomes from six avian species, chosen to reflect the full range of sexual dimorphism and mating systems. Our analyses combine expression and population genomic statistics across reproductive and somatic tissue, with measures of sperm competition and promiscuity. Our results reveal that balancing selection is weakest in the gonad, consistent with the resolution of sexual conflict and evolutionary theory that phenotypic sex differences are associated with lower levels of ongoing conflict. We also demonstrate a clear link between variation in sexual conflict and levels of genetic variation across phylogenetic space in a comparative framework. Our observations suggest that this conflict is short-lived, and is resolved via the decoupling of male and female gene expression patterns, with important implications for the role of sexual selection in adaptive potential and role of dimorphism in facilitating sex-specific fitness optima.

Recent Sex Chromosome Divergence despite Ancient Dioecy in the Willow Salix viminalis.

Sex chromosomes can evolve when recombination is halted between a pair of chromosomes, and this can lead to degeneration of the sex-limited chromosome. In the early stages of differentiation sex chromosomes are homomorphic, and even though homomorphic sex chromosomes are very common throughout animals and plants, we know little about the evolutionary forces shaping these types of sex chromosomes. We used DNA- and RNA-Seq data from females and males to explore the sex chromosomes in the female heterogametic willow, Salix viminalis, a species with ancient dioecy but with homomorphic sex chromosomes. We detected no major sex differences in read coverage in the sex determination (SD) region, indicating that the W region has not significantly degenerated. However, single nucleotide polymorphism densities in the SD region are higher in females compared with males, indicating very recent recombination suppression, followed by the accumulation of sex-specific single nucleotide polymorphisms. Interestingly, we identified two female-specific scaffolds that likely represent W-chromosome-specific sequence. We show that genes located in the SD region display a mild excess of male-biased expression in sex-specific tissue, and we use allele-specific gene expression analysis to show that this is the result of masculinization of expression on the Z chromosome rather than degeneration of female-expression on the W chromosome. Together, our results demonstrate that insertion of small DNA fragments and accumulation of sex-biased gene expression can occur before the detectable decay of the sex-limited chromosome.

Slow evolution of sex-biased genes in the reproductive tissue of the dioecious plant Salix viminalis.

The relative rate of evolution for sex-biased genes has often been used as a measure of the strength of sex-specific selection. In contrast to studies in a wide variety of animals, far less is known about the molecular evolution of sex-biased genes in plants, particularly in dioecious angiosperms. Here, we investigate the gene expression patterns and evolution of sex-biased genes in the dioecious plant Salix viminalis. We observe lower rates of sequence evolution for male-biased genes expressed in the reproductive tissue compared to unbiased and female-biased genes. These results could be partially explained by the lower codon usage bias for male-biased genes leading to elevated rates of synonymous substitutions compared to unbiased genes. However, the stronger haploid selection in the reproductive tissue of plants, together with pollen competition, would also lead to higher levels of purifying selection acting to remove deleterious variation. Future work should focus on the differential evolution of haploid- and diploid-specific genes to understand the selective dynamics acting on these loci.

Sexual selection drives evolution and rapid turnover of male gene expression.

The profound and pervasive differences in gene expression observed between males and females, and the unique evolutionary properties of these genes in many species, have led to the widespread assumption that they are the product of sexual selection and sexual conflict. However, we still lack a clear understanding of the connection between sexual selection and transcriptional dimorphism, often termed sex-biased gene expression. Moreover, the relative contribution of sexual selection vs. drift in shaping broad patterns of expression, divergence, and polymorphism remains unknown. To assess the role of sexual selection in shaping these patterns, we assembled transcriptomes from an avian clade representing the full range of sexual dimorphism and sexual selection. We use these species to test the links between sexual selection and sex-biased gene expression evolution in a comparative framework. Through ancestral reconstruction of sex bias, we demonstrate a rapid turnover of sex bias across this clade driven by sexual selection and show it to be primarily the result of expression changes in males. We use phylogenetically controlled comparative methods to demonstrate that phenotypic measures of sexual selection predict the proportion of male-biased but not female-biased gene expression. Although male-biased genes show elevated rates of coding sequence evolution, consistent with previous reports in a range of taxa, there is no association between sexual selection and rates of coding sequence evolution, suggesting that expression changes may be more important than coding sequence in sexual selection. Taken together, our results highlight the power of sexual selection to act on gene expression differences and shape genome evolution.

Sperm competition shapes gene expression and sequence evolution in the ocellated wrasse.

Gene expression differences between males and females often underlie sexually dimorphic phenotypes, and the expression levels of genes that are differentially expressed between the sexes are thought to respond to sexual selection. Most studies on the transcriptomic response to sexual selection treat sexual selection as a single force, but postmating sexual selection in particular is expected to specifically target gonadal tissue. The three male morphs of the ocellated wrasse (Symphodus ocellatus) make it possible to test the role of postmating sexual selection in shaping the gonadal transcriptome. Nesting males hold territories and have the highest reproductive success, yet we detected feminization of their gonadal gene expression compared to satellite males. Satellite males are less brightly coloured and experience more intense sperm competition than nesting males. In line with postmating sexual selection affecting gonadal gene expression, we detected a more masculinized expression profile in satellites. Sneakers are the lowest quality males and showed both de-masculinization and de-feminization of gene expression. We also detected higher rates of gene sequence evolution of male-biased genes compared to unbiased genes, which could at least in part be explained by positive selection. Together, these results reveal the potential for postmating sexual selection to drive higher rates of gene sequence evolution and shape the gonadal transcriptome profile.

Positive Selection Underlies Faster-Z Evolution of Gene Expression in Birds.

The elevated rate of evolution for genes on sex chromosomes compared with autosomes (Fast-X or Fast-Z evolution) can result either from positive selection in the heterogametic sex or from nonadaptive consequences of reduced relative effective population size. Recent work in birds suggests that Fast-Z of coding sequence is primarily due to relaxed purifying selection resulting from reduced relative effective population size. However, gene sequence and gene expression are often subject to distinct evolutionary pressures; therefore, we tested for Fast-Z in gene expression using next-generation RNA-sequencing data from multiple avian species. Similar to studies of Fast-Z in coding sequence, we recover clear signatures of Fast-Z in gene expression; however, in contrast to coding sequence, our data indicate that Fast-Z in expression is due to positive selection acting primarily in females. In the soma, where gene expression is highly correlated between the sexes, we detected Fast-Z in both sexes, although at a higher rate in females, suggesting that many positively selected expression changes in females are also expressed in males. In the gonad, where intersexual correlations in expression are much lower, we detected Fast-Z for female gene expression, but crucially, not males. This suggests that a large amount of expression variation is sex-specific in its effects within the gonad. Taken together, our results indicate that Fast-Z evolution of gene expression is the product of positive selection acting on recessive beneficial alleles in the heterogametic sex. More broadly, our analysis suggests that the adaptive potential of Z chromosome gene expression may be much greater than that of gene sequence, results which have important implications for the role of sex chromosomes in speciation and sexual selection.

Masculinization of gene expression is associated with exaggeration of male sexual dimorphism.

Gene expression differences between the sexes account for the majority of sexually dimorphic phenotypes, and the study of sex-biased gene expression is important for understanding the genetic basis of complex sexual dimorphisms. However, it has been difficult to test the nature of this relationship due to the fact that sexual dimorphism has traditionally been conceptualized as a dichotomy between males and females, rather than an axis with individuals distributed at intermediate points. The wild turkey (Meleagris gallopavo) exhibits just this sort of continuum, with dominant and subordinate males forming a gradient in male secondary sexual characteristics. This makes it possible for the first time to test the correlation between sex-biased gene expression and sexually dimorphic phenotypes, a relationship crucial to molecular studies of sexual selection and sexual conflict. Here, we show that subordinate male transcriptomes show striking multiple concordances with their relative phenotypic sexual dimorphism. Subordinate males were clearly male rather than intersex, and when compared to dominant males, their transcriptomes were simultaneously demasculinized for male-biased genes and feminized for female-biased genes across the majority of the transcriptome. These results provide the first evidence linking sexually dimorphic transcription and sexually dimorphic phenotypes. More importantly, they indicate that evolutionary changes in sexual dimorphism can be achieved by varying the magnitude of sex-bias in expression across a large proportion of the coding content of a genome.

Variation in promiscuity and sexual selection drives avian rate of Faster-Z evolution.

Higher rates of coding sequence evolution have been observed on the Z chromosome relative to the autosomes across a wide range of species. However, despite a considerable body of theory, we lack empirical evidence explaining variation in the strength of the Faster-Z Effect. To assess the magnitude and drivers of Faster-Z Evolution, we assembled six de novo transcriptomes, spanning 90 million years of avian evolution. Our analysis combines expression, sequence and polymorphism data with measures of sperm competition and promiscuity. In doing so, we present the first empirical evidence demonstrating the positive relationship between Faster-Z Effect and measures of promiscuity, and therefore variance in male mating success. Our results from multiple lines of evidence indicate that selection is less effective on the Z chromosome, particularly in promiscuous species, and that Faster-Z Evolution in birds is due primarily to genetic drift. Our results reveal the power of mating system and sexual selection in shaping broad patterns in genome evolution.

W chromosome expression responds to female-specific selection.

The W chromosome is predicted to be subject to strong female-specific selection stemming from its female-limited inheritance and therefore should play an important role in female fitness traits. However, the overall importance of directional selection in shaping the W chromosome is unknown because of the powerful degradative forces that act to decay the nonrecombining sections of the genome. Here we greatly expand the number of known W-linked genes and assess the expression of the W chromosome after >100 generations of different female-specific selection regimens in different breeds of chicken and in the wild ancestor, the Red Jungle Fowl. Our results indicate that female-specific selection has a significant effect on W chromosome gene-expression patterns, with a strong convergent pattern of up-regulation associated with increased female-specific selection. Many of the transcriptional changes in the female-selected breeds are the product of positive selection, suggesting that selection is an important force in shaping the evolution of gene expression on the W chromosome, a finding consistent with both the importance of the W chromosome in female fertility and the haploid nature of the W. Taken together, these data provide evidence for the importance of the sex-limited chromosome in a female heterogametic species and show that sex-specific selection can act to preserve sex-limited chromosomes from degrading forces.

The evolution of gene expression and the transcriptome-phenotype relationship.

Changes in gene expression underlie the adaptive evolution in many complex phenotypes, and the recent increase in the availability of multi-species comparative transcriptome data has made it possible to scan whole transcriptomes for loci that have experienced adaptive changes in expression. However, despite the increase in data availability, current models of gene expression evolution often do not account for the complexities and inherent noise associated with transcriptome data. Additionally, in contrast to current models of gene sequence evolution, models of transcriptome evolution often lack the sophistication to effectively determine whether transcriptional differences between species or within a clade are the result of neutral or adaptive processes. In this review, we discuss the tools, methods and models that define our current understanding of the relationship between gene expression and complex phenotype evolution. Our goal is to summarize what we know about the evolution of global gene expression patterns underlying complex traits, as well to identify some of the questions that remain to be answered.

Recent progress in understanding the genomic architecture of sexual conflict.

Genomic conflict between the sexes over shared traits is widely assumed to be resolved through the evolution of sex-biased expression and the subsequent emergence of sexually dimorphic phenotypes. However, while there is support for a broad relationship between genome-wide patterns of expression level and sexual conflict, recent studies suggest that sex differences in the nature and strength of interactions between loci are instead key to conflict resolution. Furthermore, the advent of new technologies for measuring and perturbing expression means we now have much more power to detect genomic signatures of sexual conflict. Here, we review our current understanding of the genomic architecture of sexual conflict in the light of these new studies and highlight the potential for novel approaches to address outstanding knowledge gaps.