Emergence of male-biased genes on the chicken Z-chromosome: sex-chromosome contrasts between male and female heterogametic systems

Disclosing complex mutational dynamics at a Y chromosome palindrome evolving through intra- and inter-chromosomal gene conversion

The human MSY ampliconic region is mainly composed of large duplicated sequences that are organized in eight palindromes (termed P1-P8), and may undergo arm-to-arm gene conversion. Although the importance of these elements is widely recognized, their evolutionary dynamics are still nuanced. Here, we focused on the P8 palindrome, which shows a complex evolutionary history, being involved in intra- and inter-chromosomal gene conversion. To disclose its evolutionary complexity, we performed a high-depth (50×) targeted next-generation sequencing of this element in 157 subjects belonging to the most divergent lineages of the Y chromosome tree. We found a total of 72 polymorphic paralogous sequence variants that have been exploited to identify 41 Y-Y gene conversion events that occurred during recent human history. Through our analysis, we were able to categorize P8 arms into three portions, whose molecular diversity was modelled by different evolutionary forces. Notably, the outer region of the palindrome is not involved in any gene conversion event and evolves exclusively through the action of mutational pressure. The inner region is affected by Y-Y gene conversion occurring at a rate of 1.52 × 10-5 conversions/base/year, with no bias towards the retention of the ancestral state of the sequence. In this portion, GC-biased gene conversion is counterbalanced by a mutational bias towards AT bases. Finally, the middle region of the arms, in addition to intra-chromosomal gene conversion, is involved in X-to-Y gene conversion (at a rate of 6.013 × 10-8 conversions/base/year) thus being a major force in the evolution of the VCY/VCX gene family.

Sex-biased gene expression and recent sex chromosome turnover

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

Sex chromosome evolution among amniotes: is the origin of sex chromosomes non-random?

Sex chromosomes are a great example of a convergent evolution at the genomic level, having evolved dozens of times just within amniotes. An intriguing question is whether this repeated evolution was random, or whether some ancestral syntenic blocks have significantly higher chance to be co-opted for the role of sex chromosomes owing to their gene content related to gonad development. Here, we summarize current knowledge on the evolutionary history of sex determination and sex chromosomes in amniotes and evaluate the hypothesis of non-random emergence of sex chromosomes. The current data on the origin of sex chromosomes in amniotes suggest that their evolution is indeed non-random. However, this non-random pattern is not very strong, and many syntenic blocks representing putatively independently evolved sex chromosomes are unique. Still, repeatedly co-opted chromosomes are an excellent model system, as independent co-option of the same genomic region for the role of sex chromosome offers a great opportunity for testing evolutionary scenarios on the sex chromosome evolution under the explicit control for the genomic background and gene identity. Future studies should use these systems more to explore the convergent/divergent evolution of sex chromosomes. This article is part of the theme issue 'Challenging the paradigm in sex chromosome evolution: empirical and theoretical insights with a focus on vertebrates (Part II)'.

Transcriptomic and Epigenetic Preservation of Genetic Sex Identity in Estrogen-feminized Male Chicken Embryonic Gonads

Whereas in ovo exposure of genetically male (ZZ) chicken embryos to exogenous estrogens temporarily feminizes gonads at the time of hatching, the morphologically ovarian ZZ-gonads (FemZZs for feminized ZZ gonads) are masculinized back to testes within 1 year. To identify the feminization-resistant "memory" of genetic male sex, FemZZs showing varying degrees of feminization were subjected to transcriptomic, DNA methylome, and immunofluorescence analyses. Protein-coding genes were classified based on their relative mRNA expression across normal ZZ-testes, genetically female (ZW) ovaries, and FemZZs. We identified a group of 25 genes that were strongly expressed in both ZZ-testes and FemZZs but dramatically suppressed in ZW-ovaries. Interestingly, 84% (21/25) of these feminization-resistant testicular marker genes, including the DMRT1 master masculinizing gene, were located in chromosome Z. Expression of representative marker genes of germline cells (eg, DAZL or DDX4/VASA) was stronger in FemZZs than normal ZZ-testes or ZW-ovaries. We also identified 231 repetitive sequences (RSs) that were strongly expressed in both ZZ-testes and FemZZs, but these RSs were not enriched in chromosome Z. Although 94% (165/176) of RSs exclusively expressed in ZW-ovaries were located in chromosome W, no feminization-inducible RS was detected in FemZZs. DNA methylome analysis distinguished FemZZs from normal ZZ- and ZW-gonads. Immunofluorescence analysis of FemZZ gonads revealed expression of DMRT1 protein in medullary SOX9+ somatic cells and apparent germline cell populations in both medulla and cortex. Taken together, our study provides evidence that both somatic and germline cell populations in morphologically feminized FemZZs maintain significant transcriptomic and epigenetic memories of genetic sex.

Sex determination, gonadal sex differentiation, and plasticity in vertebrate species

A diverse array of sex determination (SD) mechanisms, encompassing environmental to genetic, have been found to exist among vertebrates, covering a spectrum from fixed SD mechanisms (mammals) to functional sex change in fishes (sequential hermaphroditic fishes). A major landmark in vertebrate SD was the discovery of the SRY gene in 1990. Since that time, many attempts to clone an SRY ortholog from nonmammalian vertebrates remained unsuccessful, until 2002, when DMY/dmrt1by was discovered as the SD gene of a small fish, medaka. Surprisingly, however, DMY/dmrt1by was found in only 2 species among more than 20 species of medaka, suggesting a large diversity of SD genes among vertebrates. Considerable progress has been made over the last 3 decades, such that it is now possible to formulate reasonable paradigms of how SD and gonadal sex differentiation may work in some model vertebrate species. This review outlines our current understanding of vertebrate SD and gonadal sex differentiation, with a focus on the molecular and cellular mechanisms involved. An impressive number of genes and factors have been discovered that play important roles in testicular and ovarian differentiation. An antagonism between the male and female pathway genes exists in gonads during both sex differentiation and, surprisingly, even as adults, suggesting that, in addition to sex-changing fishes, gonochoristic vertebrates including mice maintain some degree of gonadal sexual plasticity into adulthood. Importantly, a review of various SD mechanisms among vertebrates suggests that this is the ideal biological event that can make us understand the evolutionary conundrums underlying speciation and species diversity.

The Female-Specific W Chromosomes of Birds Have Conserved Gene Contents but Are Not Feminized

Sex chromosomes are unique genomic regions with sex-specific or sex-biased inherent patterns and are expected to be more frequently subject to sex-specific selection. Substantial knowledge on the evolutionary patterns of sex-linked genes have been gained from the studies on the male heterogametic systems (XY male, XX female), but the understanding of the role of sex-specific selection in the evolution of female-heterogametic sex chromosomes (ZW female, ZZ male) is limited. Here we collect the W-linked genes of 27 birds, covering the three major avian clades: Neoaves (songbirds), Galloanserae (chicken), and Palaeognathae (ratites and tinamous). We find that the avian W chromosomes exhibit very conserved gene content despite their independent evolution of recombination suppression. The retained W-linked genes have higher dosage-sensitive and higher expression level than the lost genes, suggesting the role of purifying selection in their retention. Moreover, they are not enriched in ancestrally female-biased genes, and have not acquired new ovary-biased expression patterns after becoming W-linked. They are broadly expressed across female tissues, and the expression profile of the W-linked genes in females is not deviated from that of the homologous Z-linked genes. Together, our new analyses suggest that female-specific positive selection on the avian W chromosomes is limited, and the gene content of the W chromosomes is mainly shaped by purifying selection.

Evolutionary Patterns of Sex-Biased Genes in Three Species of Haplodiploid Insects

Females and males often differ obviously in morphology and behavior, and the differences between sexes are the result of natural selection and/or sexual selection. To a great extent, the differences between the two sexes are the result of differential gene expression. In haplodiploid insects, this phenomenon is obvious, since males develop from unfertilized zygotes and females develop from fertilized zygotes. Whiteflies of the Bemisia tabaci species complex are typical haplodiploid insects, and some species of this complex are important pests of many crops worldwide. Here, we report the transcriptome profiles of males and females in three species of this whitefly complex. Between-species comparisons revealed that non-sex-biased genes display higher variation than male-biased or female-biased genes. Sex-biased genes evolve at a slow rate in protein coding sequences and gene expression and have a pattern of evolution that differs from those of social haplodiploid insects and diploid animals. Genes with high evolutionary rates are more related to non-sex-biased traits-such as nutrition, immune system, and detoxification-than to sex-biased traits, indicating that the evolution of protein coding sequences and gene expression has been mainly driven by non-sex-biased traits.

Recessive Z-linked lethals and the retention of haplotype diversity in a captive butterfly population

Sex chromosomes are predicted to harbour elevated levels of sexually antagonistic variation due to asymmetries in the heritability of recessive traits in the homogametic versus heterogametic sex. This evolutionary dynamic may manifest as high recessive load specifically affecting the homogametic sex, and the retention of haplotype diversity in small populations. We tested the hypothesis that the Z chromosome in the butterfly Bicyclus anynana carries a high inbred load for male fertility and viability. Homozygosity of Z chromosome blocks was produced by daughter-father backcrosses, and inferred from marker loci positioned via a linkage map. Male sterility was, in general, unrelated to homozygosity in any region of the Z, but there was an extreme deficit of homozygous males within a 2 cM interval in all families. In contrast, no corresponding skew in Z genotype was detected in their (hemizygous) sisters. The same pattern was observed in historically inbred lines, indicating a high frequency of recessive lethals in the ancestral population. Allele-frequency changes between 1993 and 2006 (70 generations at N_e ~ 160) show that, despite the loss of many haplotypes, diversity was retained significantly above the neutral expectation. Effective overdominance in the lethal region can account for this effect locally but not in other parts of the chromosome, that are also associated with persistent linkage disequilibrium. These unexpected patterns suggest the operation of other factors, such as epistatic selection, recombination suppression, assortative mating and meiotic drive. Our results highlight the role of balancing selection in maintaining the inbred load and linked genetic diversity.

Expression Profile of Chicken Sex Chromosome Gene BTF3 is Linked to Gonadal Phenotype

In birds, the female is heterogametic (ZW) and the male homogametic (ZZ). The small W chromosome comprises only 28 protein coding genes (homologues to Z chromosome counterparts) and a number of repeat regions. Here, we report our analysis of one of these genes, BTF3 (basic transcription factor 3), which exhibits differential expression during gonadogenesis. We measured RNA levels of both Z and W homologues and BTF3 protein levels in male and female gonads during development of the chicken embryo. In addition, BTF3 RNA and protein levels were compared in female gonads (ovary) and in female gonads following treatment to induce sex reversal (testis). Combined BTF3 RNA levels were higher in female gonads than male gonads, while BTF3-Z was expressed at similar levels in males and females. Surprisingly, BTF3 protein levels were higher in male gonads than female gonads at embryonic day 6 (E6), suggesting translational rather than transcriptional regulation. BTF3 protein was expressed in both somatic and germ cells and was restricted to the medulla of the developing ovary in females and the sex cords of the developing testis in males. In addition, in gonadal sex-reversed females, RNA and protein levels of BTF3 were similar to those normally found in male gonads, suggesting that BTF3 expression correlated with the gonadal phenotype.

Gonadal and Endocrine Analysis of a Gynandromorphic Chicken

Birds have a ZZ male and ZW female sex chromosome system. The relative roles of genetics and hormones in regulating avian sexual development have been revealed by studies on gynandromorphs. Gynandromorphs are rare bilateral sex chimeras, male on one side of the body and female on the other. We examined a naturally occurring gynandromorphic chicken that was externally male on the right side of the body and female on the left. The bird was diploid but with a mix of ZZ and ZW cells that correlated with the asymmetric sexual phenotype. The male side was 96% ZZ, and the female side was 77% ZZ and 23% ZW. The gonads of this bird at sexual maturity were largely testicular. The right gonad was a testis, with SOX9+ Sertoli cells, DMRT1+ germ cells, and active spermatogenesis. The left gonad was primarily testicular, but with some peripheral aromatase-expressing follicles. The bird had low levels of serum estradiol and high levels of testosterone, as expected for a male. Despite the low percentage of ZW cells on that side, the left side had female sex-linked feathering, smaller muscle mass, smaller leg and spur, and smaller wattle than the male side. This indicates that these sexually dimorphic structures must be at least partly independent of sex steroid effects. Even a small percentage of ZW cells appears sufficient to support female sexual differentiation. Given the lack of chromosome-wide dosage compensation in birds, various sexually dimorphic features may arise due to Z-gene dosage differences between the sexes.

Sex chromosomes and speciation in birds and other ZW systems

Theory and empirical patterns suggest a disproportionate role for sex chromosomes in evolution and speciation. Focusing on ZW sex determination (females ZW, males ZZ; the system in birds, many snakes, and lepidopterans), I review how evolutionary dynamics are expected to differ between the Z, W and the autosomes, discuss how these differences may lead to a greater role of the sex chromosomes in speciation and use data from birds to compare relative evolutionary rates of sex chromosomes and autosomes. Neutral mutations, partially or completely recessive beneficial mutations, and deleterious mutations under many conditions are expected to accumulate faster on the Z than on autosomes. Sexually antagonistic polymorphisms are expected to arise on the Z, raising the possibility of the spread of preference alleles. The faster accumulation of many types of mutations and the potential for complex evolutionary dynamics of sexually antagonistic traits and preferences contribute to a role for the Z chromosome in speciation. A quantitative comparison among a wide variety of bird species shows that the Z tends to have less within-population diversity and greater between-species differentiation than the autosomes, likely due to both adaptive evolution and a greater rate of fixation of deleterious alleles. The W chromosome also shows strong potential to be involved in speciation, in part because of its co-inheritance with the mitochondrial genome. While theory and empirical evidence suggest a disproportionate role for sex chromosomes in speciation, the importance of sex chromosomes is moderated by their small size compared to the whole genome.

The Z chromosome is enriched for sperm proteins in two divergent species of Lepidoptera

Genes that promote sexual conflict, such as those with a sex-limited fitness benefit, are expected to accumulate differentially on sex chromosomes relative to autosomes. Few tests of this hypothesis exist for male homogametic (ZZ) taxa, however, and most use RNA expression data to identify such genes. Here, we employ a different identification method by using proteomic analysis of sperm cells to identify genes with a sex-limited benefit. We tested for a bias in genomic location of sperm protein genes in two species of Lepidoptera. An excess of sperm protein genes was identified on the Z chromosomes of both the Carolina sphinx moth (Manduca sexta) and the monarch butterfly (Danaus plexippus). Taking into consideration a Z-autosome fusion in monarchs, we discover that the ancestrally sex-linked portion of the genome is the source of this enrichment, while the newly sex-linked portion still appears similar to autosomes in relative abundance of sperm protein genes. Together, these results point to an enrichment of male-beneficial genes on the Z chromosome and demonstrate the usefulness of proteomic datasets in sexual conflict research.

A sex-chromosome inversion causes strong overdominance for sperm traits that affect siring success

Male reproductive success depends on the competitive ability of sperm to fertilize the ova, which should lead to strong selection on sperm characteristics. This raises the question of how heritable variation in sperm traits is maintained. Here we show that in zebra finches (Taeniopygia guttata) nearly half of the variance in sperm morphology is explained by an inversion on the Z chromosome with a 40% allele frequency in the wild. The sperm of males that are heterozygous for the inversion had the longest midpieces and the highest velocity. Furthermore, such males achieved the highest fertility and the highest siring success, both within-pair and extra-pair. Males homozygous for the derived allele show detrimental sperm characteristics and the lowest siring success. Our results suggest heterozygote advantage as the mechanism that maintains the inversion polymorphism and hence variance in sperm design and in fitness.

RNA sequencing reveals a complete but an unconventional type of dosage compensation in the domestic silkworm Bombyx mori

Sex chromosomal dose difference between sexes is often normalized by a gene regulatory mechanism called dosage compensation (DC). Studies indicate that DC mechanisms are generally effective in XY rather than ZW systems. However, DC studies in lepidopterans (ZW system) gave bewildering results. In Manduca sexta, DC was complete and in Plodia interpunctella, it was incomplete. In Heliconius species, dosage was found to be partly incomplete. In domesticated silkmoth Bombyx mori, DC studies have yielded contradictory results thus far, showing incomplete DC based on microarray data and a possible existence of DC based on recent reanalysis of same data. In this study, analysis of B. mori sexed embryos (78, 96 and 120 h) and larval heads using RNA sequencing suggest an onset of DC at 120 h. The average Z-linked expression is substantially less than autosomes, and the male-biased Z-linked expression observed at initial stages (78 and 96 h) gets almost compensated at 120 h embryonic stage and perfectly compensated in heads. Based on these findings, we suggest a complete but an unconventional type of DC, which may be achieved by reduced Z-linked expression in males (ZZ). To our knowledge, this is the first next-generation sequencing report showing DC in B. mori, clarifying the previous contradictions.

Evolution of sex-biased gene expression in a dioecious plant

Separate sexes and sex-biased gene expression have repeatedly evolved in animals and plants, but the underlying changes in gene expression remain unknown. Here, we studied a pair of plant species, one in which separate sexes and sex chromosomes evolved recently and one which maintained hermaphrodite flowers resembling the ancestral state, to reconstruct expression changes associated with the evolution of dioecy. We found that sex-biased gene expression has evolved in autosomal and sex-linked genes in the dioecious species. Most expression changes relative to hermaphrodite flowers occurred in females rather than males, with higher and lower expression in females leading to female-biased and male-biased expression, respectively. Expression changes were more common in genes located on the sex chromosomes than the autosomes and led to feminization of the X chromosome and masculinization of the Y chromosome. Our results support a scenario in which sex-biased gene expression evolved during the evolution of dioecy to resolve intralocus sexual conflicts over the allocation of resources.

Quantitative Mass Spectrometry Reveals Partial Translational Regulation for Dosage Compensation in Chicken

There is increasing evidence that dosage compensation is not a ubiquitous feature following sex chromosome evolution, especially not in organisms where females are the heterogametic sex, like in birds. Even when it occurs, compensation can be incomplete and limited to dosage-sensitive genes. However, previous work has mainly studied transcriptional regulation of sex-linked genes, which may not reflect expression at the protein level. Here, we used liquid chromatography-tandem mass spectrometry to detect and quantify expressed levels of more than 2,400 proteins in ten different tissues of male and female chicken embryos. For comparison, transcriptome sequencing was performed in the same individuals, five of each sex. The proteomic analysis revealed that dosage compensation was incomplete, with a mean male-to-female (M:F) expression ratio of Z-linked genes of 1.32 across tissues, similar to that at the RNA level (1.29). The mean Z chromosome-to-autosome expression ratio was close to 1 in males and lower than 1 in females, consistent with partly reduced Z chromosome expression in females. Although our results exclude a general mechanism for chromosome-wide dosage compensation at translation, 30% of all proteins encoded from Z-linked genes showed a significant change in the M:F ratio compared with the corresponding ratio at the RNA level. This resulted in a pattern where some genes showed balanced expression between sexes and some close to 2-fold higher expression in males. This suggests that proteomic analyses will be necessary to reveal a more complete picture of gene regulation and sex chromosome evolution.

Transcriptome assemblies for studying sex-biased gene expression in the guppy, Poecilia reticulata

BACKGROUND

Sexually dimorphic phenotypes are generally associated with differential gene expression between the sexes. The study of molecular evolution and genomic location of these differentially expressed, or sex-biased, genes is important for understanding inter-sexual divergence under sex-specific selection pressures. Teleost fish provide a unique opportunity to examine this divergence in the presence of variable sex-determination mechanisms of recent origin. The guppy, Poecilia reticulata, displays sexual dimorphism in size, ornaments, and behavior, traits shaped by natural and sexual selection in the wild.

RESULTS

To gain insight into molecular mechanisms underlying the guppy's sexual dimorphism, we assembled a reference transcriptome combining genome-independent as well as genome-guided assemblies and analyzed sex-biased gene expression between different tissues of adult male and female guppies. We found tissue-associated sex-biased expression of genes related to pigmentation, signal transduction, and spermatogenesis in males; and growth, cell-division, extra-cellular matrix organization, nutrient transport, and folliculogenesis in females. While most sex-biased genes were randomly distributed across linkage groups, we observed accumulation of ovary-biased genes on the sex linkage group, LG12. Both testis-biased and ovary-biased genes showed a significantly higher rate of non-synonymous to synonymous substitutions (dN/dS) compared to unbiased genes. However, in somatic tissues only female-biased genes, including those co-expressed in multiple tissues, showed elevated ratios of non-synonymous substitutions.

CONCLUSIONS

Our work identifies a set of annotated gene products that are candidate factors affecting sexual dimorphism in guppies. The differential genomic distribution of gonad-biased genes provides evidence for sex-specific selection pressures acting on the nascent sex chromosomes of the guppy. The elevated rates of evolution of testis-biased and female-biased genes indicate differing evolution under distinct selection pressures on the reproductive versus non-reproductive tissues.

Transcriptome sequencing reveals the character of incomplete dosage compensation across multiple tissues in flycatchers

Sex chromosome divergence, which follows the cessation of recombination and degeneration of the sex-limited chromosome, can cause a reduction in expression level for sex-linked genes in the heterozygous sex, unless some mechanisms of dosage compensation develops to counter the reduction in gene dose. Because large-scale perturbations in expression levels arising from changes in gene dose might have strong deleterious effects, the evolutionary response should be strong. However, in birds and in at least some other female heterogametic organisms, wholesale sex chromosome dosage compensation does not seem to occur. Using RNA-seq of multiple tissues and individuals, we investigated male and female expression levels of Z-linked and autosomal genes in the collared flycatcher, a bird for which a draft genome sequence recently has been reported. We found that male expression of Z-linked genes was on average 50% higher than female expression, although there was considerable variation in the male-to-female ratio among genes. The ratio for individual genes was well correlated among tissues and there was also a correlation in the extent of compensation between flycatcher and chicken orthologs. The relative excess of male expression was positively correlated with expression breadth, expression level, and number of interacting proteins (protein connectivity), and negatively correlated with variance in expression. These observations lead to a model of compensation occurring on a gene-by-gene basis, supported by an absence of clustering of genes on the Z chromosome with respect to the extent of compensation. Equal mean expression level of autosomal and Z-linked genes in males, and 50% higher expression of autosomal than Z-linked genes in females, is compatible with that partial compensation is achieved by hypertranscription from females' single Z chromosome. A comparison with male-to-female expression ratios in orthologous Z-linked genes of ostriches, where Z-W recombination still occurs, suggests that male-biased expression of Z-linked genes is a derived trait after avian sex chromosome divergence.

Masculinization of the x chromosome in the pea aphid

Evolutionary theory predicts that sexually antagonistic mutations accumulate differentially on the X chromosome and autosomes in species with an XY sex-determination system, with effects (masculinization or feminization of the X) depending on the dominance of mutations. Organisms with alternative modes of inheritance of sex chromosomes offer interesting opportunities for studying sexual conflicts and their resolution, because expectations for the preferred genomic location of sexually antagonistic alleles may differ from standard systems. Aphids display an XX/X0 system and combine an unusual inheritance of the X chromosome with the alternation of sexual and asexual reproduction. In this study, we first investigated theoretically the accumulation of sexually antagonistic mutations on the aphid X chromosome. Our results show that i) the X is always more favourable to the spread of male-beneficial alleles than autosomes, and should thus be enriched in sexually antagonistic alleles beneficial for males, ii) sexually antagonistic mutations beneficial for asexual females accumulate preferentially on autosomes, iii) in contrast to predictions for standard systems, these qualitative results are not affected by the dominance of mutations. Under the assumption that sex-biased gene expression evolves to solve conflicts raised by the spread of sexually antagonistic alleles, one expects that male-biased genes should be enriched on the X while asexual female-biased genes should be enriched on autosomes. Using gene expression data (RNA-Seq) in males, sexual females and asexual females of the pea aphid, we confirm these theoretical predictions. Although other mechanisms than the resolution of sexual antagonism may lead to sex-biased gene expression, we argue that they could hardly explain the observed difference between X and autosomes. On top of reporting a strong masculinization of the aphid X chromosome, our study highlights the relevance of organisms displaying an alternative mode of sex chromosome inheritance to understanding the forces shaping chromosome evolution.

Lack of dosage compensation accompanies the arrested stage of sex chromosome evolution in ostriches

Sex chromosome evolution is usually seen as a process that, once initiated, will inevitably progress toward an advanced stage of degeneration of the nonrecombining chromosome. However, despite evidence that avian sex chromosome evolution was initiated >100 Ma, ratite birds have been trapped in an arrested stage of sex chromosome divergence. We performed RNA sequencing of several tissues from male and female ostriches and assembled the transcriptome de novo. A total of 315 Z-linked genes fell into two categories: those that have equal expression level in the two sexes (for which Z–W recombination still occurs) and those that have a 2-fold excess of male expression (for which Z–W recombination has ceased). We suggest that failure to evolve dosage compensation has constrained sex chromosome divergence in this basal avian lineage. Our results indicate that dosage compensation is a prerequisite for, not only a consequence of, sex chromosome evolution.

Network analysis of functional genomics data: application to avian sex-biased gene expression

Gene expression analysis is often used to investigate the molecular and functional underpinnings of a phenotype. However, differential expression of individual genes is limited in that it does not consider how the genes interact with each other in networks. To address this shortcoming we propose a number of network-based analyses that give additional functional insights into the studied process. These were applied to a dataset of sex-specific gene expression in the chicken gonad and brain at different developmental stages. We first constructed a global chicken interaction network. Combining the network with the expression data showed that most sex-biased genes tend to have lower network connectivity, that is, act within local network environments, although some interesting exceptions were found. Genes of the same sex bias were generally more strongly connected with each other than expected. We further studied the fates of duplicated sex-biased genes and found that there is a significant trend to keep the same pattern of sex bias after duplication. We also identified sex-biased modules in the network, which reveal pathways or complexes involved in sex-specific processes. Altogether, this work integrates evolutionary genomics with systems biology in a novel way, offering new insights into the modular nature of sex-biased genes.

Trade-off between selection for dosage compensation and masculinization on the avian Z chromosome

Following the suppression of recombination, gene expression levels decline on the sex-limited chromosome, and this can lead to selection for dosage compensation in the heterogametic sex to rebalance average expression from the X or Z chromosome with average autosomal expression. At the same time, due to their unequal pattern of inheritance in males and females, the sex chromosomes are subject to unbalanced sex-specific selection, which contributes to a nonrandom distribution of sex-biased genes compared to the remainder of the genome. These two forces act against each other, and the relative importance of each is currently unclear. The Gallus gallus Z chromosome provides a useful opportunity to study the importance and trade-offs between sex-specific selection and dosage compensation in shaping the evolution of the genome as it shows incomplete dosage compensation and is also present twice as often in males than females, and therefore predicted to be enriched for male-biased genes. Here, we refine our understanding of the evolution of the avian Z chromosome, and show that multiple strata formed across the chromosome over ∼130 million years. We then use this evolutionary history to examine the relative strength of selection for sex chromosome dosage compensation vs. the cumulative effects of masculinizing selection on gene expression. We find that male-biased expression increases over time, indicating that selection for dosage compensation is relatively less important than masculinizing selection in shaping Z chromosome gene expression.

Segmental dataset and whole body expression data do not support the hypothesis that non-random movement is an intrinsic property of Drosophila retrogenes

Background

Several studies in Drosophila have shown excessive movement of retrogenes from the X chromosome to autosomes, and that these genes are frequently expressed in the testis. This phenomenon has led to several hypotheses invoking natural selection as the process driving male-biased genes to the autosomes. Metta and Schlötterer (BMC Evol Biol 2010, 10:114) analyzed a set of retrogenes where the parental gene has been subsequently lost. They assumed that this class of retrogenes replaced the ancestral functions of the parental gene, and reported that these retrogenes, although mostly originating from movement out of the X chromosome, showed female-biased or unbiased expression. These observations led the authors to suggest that selective forces (such as meiotic sex chromosome inactivation and sexual antagonism) were not responsible for the observed pattern of retrogene movement out of the X chromosome.

Results

We reanalyzed the dataset published by Metta and Schlötterer and found several issues that led us to a different conclusion. In particular, Metta and Schlötterer used a dataset combined with expression data in which significant sex-biased expression is not detectable. First, the authors used a segmental dataset where the genes selected for analysis were less testis-biased in expression than those that were excluded from the study. Second, sex-biased expression was defined by comparing male and female whole-body data and not the expression of these genes in gonadal tissues. This approach significantly reduces the probability of detecting sex-biased expressed genes, which explains why the vast majority of the genes analyzed (parental and retrogenes) were equally expressed in both males and females. Third, the female-biased expression observed by Metta and Schlötterer is mostly found for parental genes located on the X chromosome, which is known to be enriched with genes with female-biased expression. Fourth, using additional gonad expression data, we found that autosomal genes analyzed by Metta and Schlötterer are less up regulated in ovaries and have higher chance to be expressed in meiotic cells of spermatogenesis when compared to X-linked genes.

Conclusions

The criteria used to select retrogenes and the sex-biased expression data based on whole adult flies generated a segmental dataset of female-biased and unbiased expressed genes that was unable to detect the higher propensity of autosomal retrogenes to be expressed in males. Thus, there is no support for the authors’ view that the movement of new retrogenes, which originated from X-linked parental genes, was not driven by selection. Therefore, selection-based genetic models remain the most parsimonious explanations for the observed chromosomal distribution of retrogenes.

Evidence of a paucity of genes that interact with the mitochondrion on the X in mammals

Mitochondria are essential organelles whose replication, development, and physiology are dependent upon coordinated gene interactions with both the mitochondrial and the nuclear genomes. The evolution of coadapted (CA) nuclear-mitochondrial gene combinations would be facilitated if such nuclear genes were located on the X-chromosome instead of on the autosomes because of the increased probability of cotransmission. Here, we test the prediction of the CA hypothesis by investigating the chromosomal distribution of nuclear genes that interact with mitochondria. Using the online genome database BIOMART, we compared the density of genes that have a mitochondrion cellular component annotation across chromosomes in 16 vertebrates. We find a strong and highly significant genomic pattern against the CA hypothesis: nuclear genes interacting with the mitochondrion are significantly underrepresented on the X-chromosome in mammals but not in birds. We interpret our findings in terms of sexual conflict as a mechanism that may generate the observed pattern. Our finding extends single-gene theory for the evolution of sexually antagonistic genes to nuclear-mitochondrial gene combinations.

Expression reduction in mammalian X chromosome evolution refutes Ohno's hypothesis of dosage compensation

Susumu Ohno proposed in 1967 that, during the origin of mammalian sex chromosomes from a pair of autosomes, per-allele expression levels of X-linked genes were doubled to compensate for the degeneration of their Y homologs. This conjecture forms the foundation of the current evolutionary model of sex chromosome dosage compensation, but has been tested in mammals only indirectly via a comparison of expression levels between X-linked and autosomal genes in the same genome. The test results have been controversial, because examinations of different gene sets led to different conclusions that either support or refute Ohno's hypothesis. Here we resolve this uncertainty by directly comparing mammalian X-linked genes with their one-to-one orthologs in species that diverged before the origin of the mammalian sex chromosomes. Analyses of RNA sequencing data and proteomic data provide unambiguous evidence for expression halving (i.e., no change in per-allele expression level) of X-linked genes during evolution, with the exception of only ∼5% of genes that encode members of large protein complexes. We conclude that Ohno's hypothesis is rejected for the vast majority of genes, reopening the search for the evolutionary force driving the origin of chromosome-wide X inactivation in female mammals.

Retrogenes moved out of the z chromosome in the silkworm

Previous studies on organisms with well-differentiated X and Y chromosomes, such as Drosophila and mammals, consistently detected an excess of genes moving out of the X chromosome and gaining testis-biased expression. Several selective evolutionary mechanisms were shown to be associated with this nonrandom gene traffic, which contributed to the evolution of the X chromosome and autosomes. If selection drives gene traffic, such traffic should also exist in species with Z and W chromosomes, where the females are the heterogametic sex. However, no previous studies on gene traffic in species with female heterogamety have found any nonrandom chromosomal gene movement. Here, we report an excess of retrogenes moving out of the Z chromosome in an organism with the ZW sex determination system, Bombyx mori. In addition, we showed that those "out of Z" retrogenes tended to have ovary-biased expression, which is consistent with the pattern of non-retrogene traffic recently reported in birds and symmetrical to the retrogene movement in mammals and fruit flies out of the X chromosome evolving testis functions. These properties of gene traffic in the ZW system suggest a general role for the heterogamety of sex chromosomes in determining the chromosomal locations and the evolution of sex-biased genes.

Disentangling the relationship between sex-biased gene expression and X-linkage

X chromosomes are preferentially transmitted through females, which may favor the accumulation of X-linked alleles/genes with female-beneficial effects. Numerous studies have shown that genes with sex-biased expression are under- or over-represented on the X chromosomes of a wide variety of organisms. The patterns, however, vary between different animal species, and the causes of these differences are unresolved. Additionally, genes with sex-biased expression tend to be narrowly expressed in a limited number of tissues, and narrowly expressed genes are also non-randomly X-linked in a taxon-specific manner. It is therefore unclear whether the unique gene content of the X chromosome is the result of selection on genes with sex-biased expression, narrowly expressed genes, or some combination of the two. To address this problem, we measured sex-biased expression in multiple Drosophila species and at different developmental time points. These data were combined with available expression measurements from Drosophila melanogaster and mouse to reconcile the inconsistencies in X-chromosome content among taxa. Our results suggest that most of the differences between Drosophila and mammals are confounded by disparate data collection/analysis approaches as well as the correlation between sex bias and expression breadth. Both the Drosophila and mouse X chromosomes harbor an excess of genes with female-biased expression after controlling for the confounding factors, suggesting that the asymmetrical transmission of the X chromosome favors the accumulation of female-beneficial mutations in X-linked genes. However, some taxon-specific patterns remain, and we provide evidence that these are in part a consequence of constraints imposed by the dosage compensation mechanism in Drosophila.