Sex-biased gene expression on the avian Z chromosome: highly expressed genes show higher male-biased expression

Single-Hit Inactivation Drove Tumor Suppressor Genes Out of the X Chromosome during Evolution

Cancer-related genes are under intense evolutionary pressure. In this study, we conjecture that X-linked tumor suppressor genes (TSG) are not protected by the Knudson's two-hit mechanism and are therefore subject to negative selection. Accordingly, nearly all mammalian species exhibited lower TSG-to-noncancer gene ratios on their X chromosomes compared with nonmammalian species. Synteny analysis revealed that mammalian X-linked TSGs were depleted shortly after the emergence of the XY sex-determination system. A phylogeny-based model unveiled a higher X chromosome-to-autosome relocation flux for human TSGs. This was verified in other mammals by assessing the concordance/discordance of chromosomal locations of mammalian TSGs and their orthologs in Xenopus tropicalis. In humans, X-linked TSGs are younger or larger in size. Consistently, pan-cancer analysis revealed more frequent nonsynonymous somatic mutations of X-linked TSGs. These findings suggest that relocation of TSGs out of the X chromosome could confer a survival advantage by facilitating evasion of single-hit inactivation.

SIGNIFICANCE

This work unveils extensive trafficking of TSGs from the X chromosome to autosomes during evolution, thus identifying X-linked TSGs as a genetic Achilles' heel in tumor suppression.

Genetic Regulation of Avian Testis Development

As in other vertebrates, avian testes are the site of spermatogenesis and androgen production. The paired testes of birds differentiate during embryogenesis, first marked by the development of pre-Sertoli cells in the gonadal primordium and their condensation into seminiferous cords. Germ cells become enclosed in these cords and enter mitotic arrest, while steroidogenic Leydig cells subsequently differentiate around the cords. This review describes our current understanding of avian testis development at the cell biology and genetic levels. Most of this knowledge has come from studies on the chicken embryo, though other species are increasingly being examined. In chicken, testis development is governed by the Z-chromosome-linked DMRT1 gene, which directly or indirectly activates the male factors, HEMGN, SOX9 and AMH. Recent single cell RNA-seq has defined cell lineage specification during chicken testis development, while comparative studies point to deep conservation of avian testis formation. Lastly, we identify areas of future research on the genetics of avian testis development.

Variable dosage compensation is associated with female consequences of an X-linked, male-beneficial mutation

Recent theory has suggested that dosage compensation mediates sexual antagonism over X-linked genes. This process relies on the assumption that dosage compensation scales phenotypic effects between the sexes, which is largely untested. We evaluated this by quantifying transcriptome variation associated with a recently arisen, male-beneficial, X-linked mutation across tissues of the field cricket Teleogryllus oceanicus, and testing the relationship between the completeness of dosage compensation and female phenotypic effects at the level of gene expression. Dosage compensation in T. oceanicus was variable across tissues but usually incomplete, such that relative expression of X-linked genes was typically greater in females. Supporting the assumption that dosage compensation scales phenotypic effects between the sexes, we found tissues with incomplete dosage compensation tended to show female-skewed effects of the X-linked allele. In gonads, where expression of X-linked genes was most strongly female-biased, ovaries-limited genes were much more likely to be X-linked than were testes-limited genes, supporting the view that incomplete dosage compensation favours feminization of the X. Our results support the expectation that sex chromosome dosage compensation scales phenotypic effects of X-linked genes between sexes, substantiating a key assumption underlying the theoretical role of dosage compensation in determining the dynamics of sexual antagonism on the X.

Evolution of dosage compensation does not depend on genomic background

Organisms have evolved various mechanisms to cope with the differences in the gene copy numbers between sexes caused by degeneration of Y and W sex chromosomes. Complete dosage compensation or at least expression balance between sexes has been reported predominantly in XX/XY systems, but rarely in ZZ/ZW systems. However, this often-reported pattern is based on comparisons of lineages where sex chromosomes evolved from nonhomologous genomic regions, potentially differing in sensitivity to differences in gene copy numbers. Here we document that two reptilian lineages (XX/XY iguanas and ZZ/ZW softshell turtles), which independently co-opted the same ancestral genomic region for the function of sex chromosomes, evolved different gene dose regulatory mechanisms. The independent co-option of the same genomic region for the role of sex chromosomes as in the iguanas and the softshell turtles offers great opportunity for testing evolutionary scenarios on sex chromosome evolution under the explicit control of the genomic background and gene identity. We show that the parallel loss of functional genes from the Y chromosome of the green anole and the W chromosome of the Florida softshell turtle led to different dosage compensation mechanisms. Our approach controlling for genetic background thus does not support that the variability in the regulation of gene dose differences is a consequence of ancestral autosomal gene content.

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.

Evolutionary Dynamics of Sex Chromosomes of Paleognathous Birds

Standard models of sex chromosome evolution propose that recombination suppression leads to the degeneration of the heterogametic chromosome, as is seen for the Y chromosome in mammals and the W chromosome in most birds. Unlike other birds, paleognaths (ratites and tinamous) possess large nondegenerate regions on their sex chromosomes (PARs or pseudoautosomal regions). It remains unclear why these large PARs are retained over >100 Myr, and how this retention impacts the evolution of sex chromosomes within this system. To address this puzzle, we analyzed Z chromosome evolution and gene expression across 12 paleognaths, several of whose genomes have recently been sequenced. We confirm at the genomic level that most paleognaths retain large PARs. As in other birds, we find that all paleognaths have incomplete dosage compensation on the regions of the Z chromosome homologous to degenerated portions of the W (differentiated regions), but we find no evidence for enrichments of male-biased genes in PARs. We find limited evidence for increased evolutionary rates (faster-Z) either across the chromosome or in differentiated regions for most paleognaths with large PARs, but do recover signals of faster-Z evolution in tinamou species with mostly degenerated W chromosomes, similar to the pattern seen in neognaths. Unexpectedly, in some species, PAR-linked genes evolve faster on average than genes on autosomes, suggested by diverse genomic features to be due to reduced efficacy of selection in paleognath PARs. Our analysis shows that paleognath Z chromosomes are atypical at the genomic level, but the evolutionary forces maintaining largely homomorphic sex chromosomes in these species remain elusive.

Insights into Avian Incomplete Dosage Compensation: Sex-Biased Gene Expression Coevolves with Sex Chromosome Degeneration in the Common Whitethroat

Non-recombining sex chromosomes (Y and W) accumulate deleterious mutations and degenerate. This poses a problem for the heterogametic sex (XY males; ZW females) because a single functional gene copy often implies less gene expression and a potential imbalance of crucial expression networks. Mammals counteract this by dosage compensation, resulting in equal sex chromosome expression in males and females, whereas birds show incomplete dosage compensation with significantly lower expression in females (ZW). Here, we study the evolution of Z and W sequence divergence and sex-specific gene expression in the common whitethroat (Sylvia communis), a species within the Sylvioidea clade where a neo-sex chromosome has been formed by a fusion between an autosome and the ancestral sex chromosome. In line with data from other birds, females had lower expression than males at the majority of sex-linked genes. Results from the neo-sex chromosome region showed that W gametologs have diverged functionally to a higher extent than their Z counterparts, and that the female-to-male expression ratio correlated negatively with the degree of functional divergence of these gametologs. We find it most likely that sex-linked genes are being suppressed in females as a response to W chromosome degradation, rather than that these genes experience relaxed selection, and thus diverge more, by having low female expression. Overall, our data of this unique avian neo-sex chromosome system suggest that incomplete dosage compensation evolves, at least partly, through gradual accumulation of deleterious mutations at the W chromosome and declining female gene expression.

Evolution of Sex Chromosome Dosage Compensation in Animals: A Beautiful Theory, Undermined by Facts and Bedeviled by Details

Many animals with genetic sex determination harbor heteromorphic sex chromosomes, where the heterogametic sex has half the gene dose of the homogametic sex. This imbalance, if reflected in the abundance of transcripts or proteins, has the potential to deleteriously disrupt interactions between X-linked and autosomal loci in the heterogametic sex. Classical theory predicts that molecular mechanisms will evolve to provide dosage compensation that recovers expression levels comparable to ancestral expression prior to sex chromosome divergence. Such dosage compensating mechanisms may also, secondarily, result in balanced sex-linked gene expression between males and females. However, numerous recent studies addressing sex chromosome dosage compensation (SCDC) in a diversity of animals have yielded a surprising array of patterns concerning dosage compensation in the heterogametic sex, as well as dosage balance between sexes. These results substantially contradict longstanding theory, catalyzing both novel perspectives and new approaches in dosage compensation research. In this review, we summarize the theory, analytical approaches, and recent results concerning evolutionary patterns of SCDC in animals. We also discuss methodological challenges and discrepancies encountered in this research, which often underlie conflicting results. Finally, we discuss what outstanding questions and opportunities exist for future research on SCDC.

Gene expression in the brain of a migratory songbird during breeding and migration

BACKGROUND

We still have limited knowledge about the underlying genetic mechanisms that enable migrating species of birds to navigate the globe. Here we make an attempt to get insight into the genetic architecture controlling this complex innate behaviour. We contrast the gene expression profiles of two closely related songbird subspecies with divergent migratory phenotypes. In addition to comparing differences in migratory strategy we include a temporal component and contrast patterns between breeding adults and autumn migrating juvenile birds of both subspecies. The two willow warbler subspecies, Phylloscopus trochilus trochilus and P. t. acredula, are remarkably similar both in phenotype and genotype and have a narrow contact zone in central Scandinavia. Here we used a microarray gene chip representing 23,136 expressed sequence tags (ESTs) from the zebra finch Taeniopygia guttata to identify mRNA level differences in willow warbler brain tissue in relation to subspecies and season.

RESULTS

Out of the 22,109 EST probe sets that remained after filtering poorly binding probes, we found 11,898 (51.8 %) probe sets that could be reliably and uniquely matched to a total of 6,758 orthologous zebra finch genes. The two subspecies showed very similar levels of gene expression with less than 0.1 % of the probe sets being significantly differentially expressed. In contrast, 3,045 (13.8 %) probe sets were found to be differently regulated between samples collected from breeding adults and autumn migrating juvenile birds. The genes found to be differentially expressed between seasons appeared to be enriched for functional roles in neuronal firing and neuronal synapse formation.

CONCLUSIONS

Our results show that only few genes are differentially expressed between the subspecies. This suggests that the different migration strategies of the subspecies might be governed by few genes, or that the expression patterns of those genes are time-structured or tissue-specific in ways, which our approach fails to uncover. Our findings will be useful in the planning of new experiments designed to unravel the genes involved in the migratory program of birds.

Characterization of the genome and transcriptome of the blue tit Cyanistes caeruleus: polymorphisms, sex-biased expression and selection signals

Decoding genomic sequences and determining their variation within populations has potential to reveal adaptive processes and unravel the genetic basis of ecologically relevant trait variation within a species. The blue tit Cyanistes caeruleus--a long-time ecological model species--has been used to investigate fitness consequences of variation in mating and reproductive behaviour. However, very little is known about the underlying genetic changes due to natural and sexual selection in the genome of this songbird. As a step to bridge this gap, we assembled the first draft genome of a single blue tit, mapped the transcriptome of five females and five males to this reference, identified genomewide variants and performed sex-differential expression analysis in the gonads, brain and other tissues. In the gonads, we found a high number of sex-biased genes, and of those, a similar proportion were sex-limited (genes only expressed in one sex) in males and females. However, in the brain, the proportion of female-limited genes within the female-biased gene category (82%) was substantially higher than the proportion of male-limited genes within the male-biased category (6%). This suggests a predominant on-off switching mechanism for the female-limited genes. In addition, most male-biased genes were located on the Z-chromosome, indicating incomplete dosage compensation for the male-biased genes. We called more than 500,000 SNPs from the RNA-seq data. Heterozygote detection in the single reference individual was highly congruent between DNA-seq and RNA-seq calling. Using information from these polymorphisms, we identified potential selection signals in the genome. We list candidate genes which can be used for further sequencing and detailed selection studies, including genes potentially related to meiotic drive evolution. A public genome browser of the blue tit with the described information is available at http://public-genomes-ngs.molgen.mpg.de.

Sex differences and rapid estrogen signaling: A look at songbird audition

The actions of estrogens have been associated with brain differentiation and sexual dimorphism in a wide range of vertebrates. Here we consider the actions of brain-derived 'neuroestrogens' in the forebrain and the accompanying differences and similarities observed between males and females in a variety of species. We summarize recent evidence showing that baseline and fluctuating levels of neuroestrogens within the auditory forebrain of male and female zebra finches are largely similar, and that neuroestrogens enhance auditory representations in both sexes. With a comparative perspective we review evidence that non-genomic mechanisms of neuroestrogen actions are sexually differentiated, and we propose a working model for nonclassical estrogen signaling via the MAPK intracellular signaling cascade in the songbird auditory forebrain that is informed by the way sex differences may be compensated. This view may lead to a more comprehensive understanding of how sex influences estradiol-dependent modulation of sensorimotor representations.

Complete dosage compensation and sex-biased gene expression in the moth Manduca sexta

Sex chromosome dosage compensation balances homogametic sex chromosome expression with autosomal expression in the heterogametic sex, leading to sex chromosome expression parity between the sexes. If compensation is incomplete, this can lead to expression imbalance and sex-biased gene expression. Recent work has uncovered an intriguing and variable pattern of dosage compensation across species that includes a lack of complete dosage compensation in ZW species compared with XY species. This has led to the hypothesis that ZW species do not require complete compensation or that complete compensation would negatively affect their fitness. To date, only one study, a study of the moth Bombyx mori, has discovered evidence for complete dosage compensation in a ZW species. We examined another moth species, Manduca sexta, using high-throughput sequencing to survey gene expression in the head tissue of males and females. We found dosage compensation to be complete in M. sexta with average expression between the Z chromosome in males and females being equal. When genes expressed at very low levels are removed by filtering, we found that average autosome expression was highly similar to average Z expression, suggesting that the majority of genes in M. sexta are completely dosage compensated. Further, this compensation was accompanied by sex-specific gene expression associated with important sexually dimorphic traits. We suggest that complete dosage compensation in ZW species might be more common than previously appreciated and linked to additional selective processes, such as sexual selection. More ZW and lepidopteran species should now be examined in a phylogenetic framework, to understand the evolution of dosage compensation.

Reference genes for quantitative gene expression studies in multiple avian species

Quantitative real-time PCR (qPCR) rapidly and reliably quantifies gene expression levels across different experimental conditions. Selection of suitable reference genes is essential for meaningful normalization and thus correct interpretation of data. In recent years, an increasing number of avian species other than the chicken has been investigated molecularly, highlighting the need for an experimentally validated pan-avian primer set for reference genes. Here we report testing a set for 14 candidate reference genes (18S, ABL, GAPDH, GUSB, HMBS, HPRT, PGK1, RPL13, RPL19, RPS7, SDHA, TFRC, VIM, YWHAZ) on different tissues of the mallard (Anas platyrhynchos), domestic chicken (Gallus gallus domesticus), common crane (Grus grus), white-tailed eagle (Haliaeetus albicilla), domestic turkey (Meleagris gallopavo f. domestica), cockatiel (Nymphicus hollandicus), Humboldt penguin (Sphenicus humboldti), ostrich (Struthio camelus) and zebra finch (Taeniopygia guttata), spanning a broad range of the phylogenetic tree of birds. Primer pairs for six to 11 genes were successfully established for each of the nine species. As a proof of principle, we analyzed expression levels of 10 candidate reference genes as well as FOXP2 and the immediate early genes, EGR1 and CFOS, known to be rapidly induced by singing in the avian basal ganglia. We extracted RNA from microbiopsies of the striatal song nucleus Area X of adult male zebra finches after they had sang or remained silent. Using three different statistical algorithms, we identified five genes (18S, PGK1, RPS7, TFRC, YWHAZ) that were stably expressed within each group and also between the singing and silent conditions, establishing them as suitable reference genes. In conclusion, the newly developed pan-avian primer set allows accurate normalization and quantification of gene expression levels in multiple avian species.

X chromosome regulation of autosomal gene expression in bovine blastocysts

Although X chromosome inactivation in female mammals evolved to balance the expression of X chromosome and autosomal genes in the two sexes, female embryos pass through developmental stages in which both X chromosomes are active in somatic cells. Bovine blastocysts show higher expression of many X genes in XX than XY embryos, suggesting that X inactivation is not complete. Here, we reanalyzed bovine blastocyst microarray expression data from a network perspective with a focus on interactions between X chromosome and autosomal genes. Whereas male-to-female ratios of expression of autosomal genes were distributed around a mean of 1, X chromosome genes were clearly shifted towards higher expression in females. We generated gene coexpression networks and identified a major module of genes with correlated gene expression that includes female-biased X genes and sexually dimorphic autosomal genes for which the sexual dimorphism is likely driven by the X genes. In this module, expression of X chromosome genes correlates with autosome genes, more than the expression of autosomal genes with each other. Our study identifies correlated patterns of autosomal and X-linked genes that are likely influenced by the sexual imbalance of X gene expression when X inactivation is inefficient.

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.

RNA sequencing reveals sexually dimorphic gene expression before gonadal differentiation in chicken and allows comprehensive annotation of the W-chromosome

Background

Birds have a ZZ male: ZW female sex chromosome system and while the Z-linked DMRT1 gene is necessary for testis development, the exact mechanism of sex determination in birds remains unsolved. This is partly due to the poor annotation of the W chromosome, which is speculated to carry a female determinant. Few genes have been mapped to the W and little is known of their expression.

Results

We used RNA-seq to produce a comprehensive profile of gene expression in chicken blastoderms and embryonic gonads prior to sexual differentiation. We found robust sexually dimorphic gene expression in both tissues pre-dating gonadogenesis, including sex-linked and autosomal genes. This supports the hypothesis that sexual differentiation at the molecular level is at least partly cell autonomous in birds. Different sets of genes were sexually dimorphic in the two tissues, indicating that molecular sexual differentiation is tissue specific. Further analyses allowed the assembly of full-length transcripts for 26 W chromosome genes, providing a view of the W transcriptome in embryonic tissues. This is the first extensive analysis of W-linked genes and their expression profiles in early avian embryos.

Conclusion

Sexual differentiation at the molecular level is established in chicken early in embryogenesis, before gonadal sex differentiation. We find that the W chromosome is more transcriptionally active than previously thought, expand the number of known genes to 26 and present complete coding sequences for these W genes. This includes two novel W-linked sequences and three small RNAs reassigned to the W from the Un_Random chromosome.