Evolutionary strata on the chicken Z chromosome: implications for sex chromosome evolution

The human X chromosome exhibits four "evolutionary strata," interpreted to represent distinct steps in the process whereby recombination became arrested between the proto X and proto Y. To test if this is a general feature of sex chromosome evolution, we studied the Z-W sex chromosomes of birds, which have female rather than male heterogamety and evolved from a different autosome pair than the mammalian X and Y. Here we analyze all five known gametologous Z-W gene pairs to investigate the "strata" hypothesis in birds. Comparisons of the rates of synonymous substitution and intronic divergence between Z and W gametologs reveal the presence of at least two evolutionary strata spread over the p and q arms of the chicken Z chromosome. A phylogenetic analysis of intronic sequence data from different avian lineages indicates that Z-W recombination ceased in the oldest stratum (on Zq; CHD1Z, HINTZ, and SPINZ) 102-170 million years ago (MYA), before the split of the Neoaves and Eoaves. However, recombination continued in the second stratum (on Zp; UBAP2Z and ATP5A1Z) until after the divergence of extant avian orders, with Z and W diverging 58-85 MYA. Our data suggest that progressive and stepwise cessation of recombination is a general feature behind sex chromosome evolution.     

Molecular evolution of the avian CHD1 genes on the Z and W sex chromosomes

Genes shared between the nonrecombining parts of the two types of sex chromosomes offer a potential means to study the molecular evolution of the same gene exposed to different genomic environments. We have analyzed the molecular evolution of the coding sequence of the first pair of genes found to be shared by the avian Z (present in both sexes) and W (female-specific) sex chromosomes, CHD1Z and CHD1W. We show here that these two genes evolve independently but are highly conserved at nucleotide as well as amino acid levels, thus not indicating a female-specific role of the CHD1W gene. From comparisons of sequence data from three avian lineages, the frequency of nonsynonymous substitutions (K(a)) was found to be higher for CHD1W (1.55 per 100 sites) than for CHD1Z (0.81), while the opposite was found for synonymous substitutions (K(s), 13.5 vs. 22.7). We argue that the lower effective population size and the absence of recombination on the W chromosome will generally imply that nonsynonymous substitutions accumulate faster on this chromosome than on the Z chromosome. The same should be true for the Y chromosome relative to the X chromosome in XY systems. Our data are compatible with a male-biased mutation rate, manifested by the faster rate of neutral evolution (synonymous substitutions) on the Z chromosome than on the female-specific W chromosome.     

Male-Biased Mutation Rates Revealed from Z and W Chromosome-Linked ATP Synthase α-Subunit (ATP5A1) Sequences in Birds

Whether the mutation rate differs between sexes has been a matter of discussion for years. Molecular analyses of mammals have indicated that males mutate more often than females, as manifested by the faster rate of neutral sequence evolution on the Y chromosome than on the X chromosome. However, these observations can as well be interpreted as specific reduction of the X chromosome mutation rate, which would be adaptive because of reducing the number of slightly deleterious recessive mutations exposed in hemizygote males. Recently, data from birds have suggested that vertebrate mutation rates may indeed be male-biased. In birds, females are the heterogametic sex (ZW), and analyses of the Z-linked CHD1Z gene have shown that it evolves faster than its W-linked and thus female-specific homologue, CHD1W. We have now studied the second avian gene known to exist in a copy on the nonrecombining regions of both the Z and the W chromosome, viz., the ATP synthase α-subunit (ATP5A1). In independent comparisons of three pairs of bird species from divergent lineages, intron sequences of the Z-linked copy (ATP5A1Z) were consistently found to evolve faster than the W-linked copy (ATP5A1W). From these data we calculated male-to-female mutation rate ratios (α) of 1.8, 2.3, and 5.0 in Galliform, Anseriform, and Ciconiiform lineages, respectively. Therefore, this study provides independent support for a male-biased mutation rate in birds. Received: 15 July 1999 / Accepted: 5 January 2000     

Retroposon insertions and the chronology of avian sex chromosome evolution

The vast majority of extant birds possess highly differentiated Z and W sex chromosomes. Nucleotide sequence data from gametologs (homologs on opposite sex chromosomes) suggest that this divergence occurred throughout early bird evolution via stepwise cessation of recombination between identical sex chromosomal regions. Here, we investigated avian sex chromosome differentiation from a novel perspective, using retroposon insertions and random insertions/deletions for the reconstruction of gametologous gene trees. Our data confirm that the CHD1Z/CHD1W genes differentiated in the ancestor of the neognaths, whereas the NIPBLZ/NIPBLW genes diverged in the neoavian ancestor and independently within Galloanserae. The divergence of the ATP5A1Z/ATP5A1W genes in galloanserans occurred independently in the chicken, the screamer, and the ancestor of duck-related birds. In Neoaves, this gene pair differentiated in each of the six sampled representatives, respectively. Additionally, three of our investigated loci can be utilized as universal, easy-to-use independent tools for molecular sexing of Neoaves or Neognathae.     

The ZW sex chromosomes of <Emphasis Type="Italic">Gekko hokouensis</Emphasis> (Gekkonidae,Squamata) represent highly conserved homology with those of avian species

Populations of the gecko lizard Gekko hokouensis (Gekkonidae, Squamata) on Okinawajima Island and a few other islands of the Ryukyu Archipelago, Japan, have the morphologically differentiated sex chromosomes, the acrocentric Z chromosome and the subtelocentric W chromosome, although the continental representative of this species reportedly shows no sex chromosome heteromorphism. To investigate the origin of sex chromosomes and the process of sex chromosomal differentiation in this species, we molecularly cloned the homologues of six chicken Z-linked genes and mapped them to the metaphase chromosomes of the Okinawajima sample. They were all localized to the Z and W chromosomes in the order ACO1/IREBP–RPS6–DMRT1–CHD1–GHR–ATP5A1, indicating that the origin of ZW chromosomes in G. hokouensis is the same as that in the class Aves, but is different from that in the suborder Ophidia. These results suggest that in reptiles the origin of sex chromosomes varies even within such a small clade as the order Squamata, employing a variety of genetic sex determination. ACO1/IREBP, RPS6, and DMRT1 were located on the Z long arm and the W short arm in the same order, suggesting that multiple rearrangements have occurred in this region of the W chromosome, where genetic differentiation between the Z and W chromosomes has been probably caused by the cessation of meiotic recombination.     

Evolution of the Spindlin Gene in Birds: Independent Cessation of the Recombination of Sex Chromosomes at the Spindlin Locus in Neognathous Birds and Tinamous, a Palaeognathous Avian Family

Tinamous (Aves, Palaeognathae, Tinamiformes) are primitive birds, generally considered to be the sister group to the ratites. Tinamous possess a W sex-chromosome, intermediate in heterochromatization between the largely euchromatic W chromosome of the ratites and the highly condensed W chromosome of the neognathous birds. Of the four genes which are known to have diverged copies on the neognathous avian W and Z chromosome (ATP5A1, CHD1, PKC and SPIN) only the spindlin gene has W- and Z-chromosomal forms in the tinamiformes. This paper describes experiments which show that the sequences of these forms are more similar to each other and to the homologous undifferentiated spindlin gene sequences in the ratite genome than to the W or Z forms of the spindlin gene in other, neognathous species. This suggests that cessation of recombination at the spindlin locus of the ancestral W and Z chromosomes of the paleognathous tinamiformes and the neognathous avian species were independent events.     

Fast Accumulation of Nonsynonymous Mutations on the Female-Specific W Chromosome in Birds

Following cessation of recombination during sex chromosome evolution, the nonrecombining sex chromosome is affected by a number of degenerative forces, possibly resulting in the fixation of deleterious mutations. This might take place because of weak selection against recessive or partly recessive deleterious mutations due to permanent heterozygosity of nonrecombining chromosomes. Furthermore, population genetic processes, such as selective sweeps, background selection, and Muller’s ratchet, result in a reduction in N_e, which increase the likelihood of fixation of deleterious mutations. Theory thus predicts that nonrecombining genes should show increased levels of nonsynonymous (d_N) to synonymous substitutions (d_S). We tested this in an avian system by estimating the ratio between d_N and d_S in six gametologous gene pairs located on the Z chromosome and the nonrecombining, female-specific W chromosome. In comparisons, we found a significantly higher d_N/d_S ratio for the W-linked than the Z-linked copy in three of the investigated genes. In a concatenated alignment of all six genes, the d_N/d_S ratio was six times higher for W-linked than Z-linked genes. By using human and mouse as outgroup in maximum likelihood analyses, W-linked genes were found to evolve differently compared with their Z-linked gametologues and outgroup sequences. This seems not to be a consequence of functional diversification because d_N/d_S ratios between gametologous gene copies were consistently low. We conclude that deleterious mutations are accumulating at a high rate on the avian W chromosome, probably as a result of the lack of recombination in this female-specific chromosome. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users. [Reviewing Editor: Dr. Deborah Charlesworth]     

Is the rate of insertion and deletion mutation male biased?: Molecular evolutionary analysis of avian and primate sex chromosome sequences.

The rate of mutation for nucleotide substitution is generally higher among males than among females, likely owing to the larger number of DNA replications in spermatogenesis than in oogenesis. For insertion and deletion (indel) mutations, data from a few human genetic disease loci indicate that the two sexes may mutate at similar rates, possibly because such mutations arise in connection with meiotic crossing over. To address origin- and sex-specific rates of indel mutation we have conducted the first large-scale molecular evolutionary analysis of indels in noncoding DNA sequences from sex chromosomes. The rates are similar on the X and Y chromosomes of primates but about twice as high on the avian Z chromosome as on the W chromosome. The fact that indels are not uncommon on the nonrecombining Y and W chromosomes excludes meiotic crossing over as the main cause of indel mutation. On the other hand, the similar rates on X and Y indicate that the number of DNA replications (higher for Y than for X) is also not the main factor. Our observations are therefore consistent with a role of both DNA replication and recombination in the generation of short insertion and deletion mutations. A significant excess of deletion compared to insertion events is observed on the avian W chromosome, consistent with gradual DNA loss on a nonrecombining chromosome.     

The evolution of sex chromosomes in organisms with separate haploid sexes

The evolution of dimorphic sex chromosomes is driven largely by the evolution of reduced recombination and the subsequent accumulation of deleterious mutations. Although these processes are increasingly well understood in diploid organisms, the evolution of dimorphic sex chromosomes in haploid organisms (U/V) has been virtually unstudied theoretically. We analyze a model to investigate the evolution of linkage between fitness loci and the sex‐determining region in U/V species. In a second step, we test how prone nonrecombining regions are to degeneration due to accumulation of deleterious mutations. Our modeling predicts that the decay of recombination on the sex chromosomes and the addition of strata via fusions will be just as much a part of the evolution of haploid sex chromosomes as in diploid sex chromosome systems. Reduced recombination is broadly favored, as long as there is some fitness difference between haploid males and females. The degeneration of the sex‐determining region due to the accumulation of deleterious mutations is expected to be slower in haploid organisms because of the absence of masking. Nevertheless, balancing selection often drives greater differentiation between the U/V sex chromosomes than in X/Y and Z/W systems. We summarize empirical evidence for haploid sex chromosome evolution and discuss our predictions in light of these findings.     

The origin and differentiation of the heteromorphic sex chromosomes Z, W, X, and Y in the frog Rana rugosa, inferred from the sequences of a sex-linked gene, ADP/ATP translocase

Sex chromosomes of the Japanese frog Rana rugosa are heteromorphic in the male (XX/XY) or in the female (ZZ/ZW) in two geographic forms, whereas they are still homomorphic in both sexes in two other forms (Hiroshima and Isehara types). To make clear the origin and differentiation mechanisms of the heteromorphic sex chromosomes, we isolated a sex-linked gene, ADP/ATP translocase, and constructed a phylogenetic tree of the genes derived from the sex chromosomes. The tree shows that the Hiroshima gene diverges first, and the rest form two clusters: one includes the Y and Z genes and the other includes the X, W, and Isehara genes. The Hiroshima gene shares more sequence similarity with the Y and Z genes than with the X, W, and Isehara genes. This suggests that the Y and Z sex chromosomes originate from the Hiroshima type, whereas the X and W chromosomes originate from the Isehara-type sex chromosome. Thus, we infer that hybridization between two ancestral forms, with the Hiroshima-type sex chromosome in one and the Isehara-type sex chromosome in the other, was the primary event causing differentiation of the heteromorphic sex chromosomes.      

The chicken (Gallus gallus) Z chromosome contains at least three nonlinear evolutionary strata

Birds have female heterogamety with Z and W sex chromosomes. These evolved from different autosomal precursor chromosomes than the mammalian X and Y. However, previous work has suggested that the pattern and process of sex chromosome evolution show many similarities across distantly related organisms. Here we show that stepwise restriction of recombination between the protosex chromosomes of birds has resulted in regions of the chicken Z chromosome showing discrete levels of divergence from W homologs (gametologs). The 12 genes analyzed fall into three levels of estimated divergence values, with the most recent divergence (d_S = 0.18–0.21) displayed by 6 genes in a region on the Z chromosome corresponding to the interval 1–11 Mb of the assembled genome sequence. Another 4 genes show intermediate divergence (d_S = 0.27–0.38) and are located in the interval 16–53 Mb. Two genes (at positions 42 and 50 Mb) with higher d_S values are located proximal to the most distal of the 4 genes with intermediate divergence, suggesting an inversion event. The distribution of genes and their divergence indicate at least three evolutionary strata, with estimated times for cessation of recombination between Z and W of 132–150 (stratum 1), 71–99 (stratum 2), and 47–57 (stratum 3) million years ago. An inversion event, or some other form of intrachromosomal rearrangement, subsequent to the formation of strata 1 and 2 has scrambled the gene order to give rise to the nonlinear arrangement of evolutionary strata currently seen on the chicken Z chromosome. These observations suggest that the progressive restriction of recombination is an integral feature of sex chromosome evolution and occurs also in systems of female heterogamety.     

A trans-species missense SNP in Amhr2 is associated with sex determination in the tiger pufferfish, Takifugu rubripes (fugu)

Heterogametic sex chromosomes have evolved independently in various lineages of vertebrates. Such sex chromosome pairs often contain nonrecombining regions, with one of the chromosomes harboring a master sex-determining (SD) gene. It is hypothesized that these sex chromosomes evolved from a pair of autosomes that diverged after acquiring the SD gene. By linkage and association mapping of the SD locus in fugu (Takifugu rubripes), we show that a SNP (C/G) in the anti-Müllerian hormone receptor type II (Amhr2) gene is the only polymorphism associated with phenotypic sex. This SNP changes an amino acid (His/Asp384) in the kinase domain. While females are homozygous (His/His384), males are heterozygous. Sex in fugu is most likely determined by a combination of the two alleles of Amhr2. Consistent with this model, the medaka hotei mutant carrying a substitution in the kinase domain of Amhr2 causes a female phenotype. The association of the Amhr2 SNP with phenotypic sex is conserved in two other species of Takifugu but not in Tetraodon. The fugu SD locus shows no sign of recombination suppression between X and Y chromosomes. Thus, fugu sex chromosomes represent an unusual example of proto-sex chromosomes. Such undifferentiated X-Y chromosomes may be more common in vertebrates than previously thought.     

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

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

Genetic mapping in a natural population of collared flycatchers (Ficedula albicollis): conserved synteny but gene order rearrangements on the avian Z chromosome

Data from completely sequenced genomes are likely to open the way for novel studies of the genetics of nonmodel organisms, in particular when it comes to the identification and analysis of genes responsible for traits that are under selection in natural populations. Here we use the draft sequence of the chicken genome as a starting point for linkage mapping in a wild bird species, the collared flycatcher - one of the most well-studied avian species in ecological and evolutionary research. A pedigree of 365 flycatchers was established and genotyped for single nucleotide polymorphisms in 23 genes selected from (and spread over most of) the chicken Z chromosome. All genes were also found to be located on the Z chromosome in the collared flycatcher, confirming conserved synteny at the level of gene content across distantly related avian lineages. This high degree of conservation mimics the situation seen for the mammalian X chromosome and may thus be a general feature in sex chromosome evolution, irrespective of whether there is male or female heterogamety. Alternatively, such unprecedented chromosomal conservation may be characteristic of most chromosomes in avian genome evolution. However, several internal rearrangements were observed, meaning that the transfer of map information from chicken to nonmodel bird species cannot always assume conserved gene orders. Interestingly, the rate of recombination on the Z chromosome of collared flycatchers was only approximately 50% that of chicken, challenging the widely held view that birds generally have high recombination rates.     

Comparison of the X and Y chromosome organization in Silene latifolia

Here we compare gene orders on the Silene latifolia sex chromosomes. On the basis of the deletion mapping results (11 markers and 23 independent Y chromosome deletion lines used), we conclude that a part of the Y chromosome (covering a region corresponding to at least 23.9 cM on the X chromosome) has been inverted. The gradient in silent-site divergence suggests that this inversion took place after the recombination arrest in this region. Because recombination arrest events followed by Y chromosome rearrangements also have been found in the human Y chromosome, this process seems to be a general evolutionary pathway.     

Phylogeny and sex chromosome evolution of Palaeognathae

Many paleognaths (ratites and tinamous) have a pair of homomorphic ZW sex chromosomes in contrast to the highly differentiated sex chromosomes of most other birds. To understand the evolutionary causes for the different tempos of sex chromosome evolution, we produced female genomes of 12 paleognathous species and reconstructed the phylogeny and the evolutionary history of paleognathous sex chromosomes. We uncovered that Palaeognathae sex chromosomes had undergone stepwise recombination suppression and formed a pattern of “evolutionary strata”. Nine of the 15 studied species' sex chromosomes have maintained homologous recombination in their long pseudoautosomal regions extending more than half of the entire chromosome length. We found that in the older strata, the W chromosome suffered more serious functional gene loss. Their homologous Z-linked regions, compared with other genomic regions, have produced an excess of species-specific autosomal duplicated genes that evolved female-specific expression, in contrast to their broadly expressed progenitors. We speculate such “defeminization” of Z chromosome with underrepresentation of female-biased genes and slow divergence of sex chromosomes of paleognaths might be related to their distinctive mode of sexual selection targeting females rather than males, which evolved in their common ancestors.     

Steps in the evolution of heteromorphic sex chromosomes

We review some recently published results on sex chromosomes in a diversity of species. We focus on several fish and some plants whose sex chromosomes appear to be 'young', as only parts of the chromosome are nonrecombining, while the rest is pseudoautosomal. However, the age of these systems is not yet very clear. Even without knowing what proportions of their genes are genetically degenerate, these cases are of great interest, as they may offer opportunities to study in detail how sex chromosomes evolve. In particular, we review evidence that recombination suppression occurs progressively in evolutionarily independent cases, suggesting that selection drives loss of recombination over increasingly large regions. We discuss how selection during the period when a chromosome is adapting to its role as a Y chromosome might drive such a process.     

Investigation of pseudoautosomal and bordering regions in avian Z and W chromosomes with the use of large insert genomic BAC clones

To study pseudoautosomal and bordering regions in the avian Z and W chromosomes, we used seven BAC clones from genomic libraries as DNA probes of fragments of different gametologs of the ATP5A1 gene located close to the proximal border of the pseudoautosomal region (PAR) of sex chromosomes of domestic chicken and Japanese quail. Localization of BAC clones TAM31-b100C09, TAM31-b99N01, TAM31-b27P16, and TAM31-b95L18 in the short arm of Z chromosomes of domestic chicken and Japanese quail (region Zp23-p22) and localization of the BAC clones CHORI-261-CH46G16, CHORI-261-CH33F10, and CHORI-261-CH64F22 on W chromosomes of these species and in the short arm of Z chromosomes (region Zp23-p22) were determined by fluorescence in situ hybridization with the use of W-specific probes. The difference in the localization of the BAC clones on the Z and W chromosomes is probably explained by divergence of the nucleotide sequences of different sex chromosomes located beyond the pseudoautosomal region.