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     

Observation of a ZZW female in a natural population: implications for avian sex determination

Avian sex determination is chromosomal; however, the underlying mechanisms are not yet understood. There is no conclusive evidence for either of two proposed mechanisms: a dominant genetic switch or a dosage mechanism. No dominant sex-determining gene on the female-specific W chromosome has been found. Birds lack inactivation of one of the Z chromosomes in males, but seem to compensate for a double dose of Z-linked genes by other mechanisms. Recent studies showing female-specific expression of two genes may support an active role of the W chromosome. To resolve the question of avian sex determination the investigation of birds with a 2A: ZZW or 2A: ZO genotype would be decisive. Here, we report the case of an apparent 2A: ZZW great reed warbler (Acrocephalus arundinaceus) female breeding in a natural population, which was detected using Z-linked microsatellites. Our data strongly suggest a role of W-linked genes in avian sex determination.     

Gene conversion drives the evolution of HINTW, an ampliconic gene on the female-specific avian W chromosome

The HINTW gene on the female-specific W chromosome of chicken and other birds is amplified and present in numerous copies. Moreover, as HINTW is distinctly different from its homolog on the Z chromosome (HINTZ), is a candidate gene in avian sex determination, and evolves rapidly under positive selection, it shows several common features to ampliconic and testis-specific genes on the mammalian Y chromosome. A phylogenetic analysis within galliform birds (chicken, turkey, quail, and pheasant) shows that individual HINTW copies within each species are more similar to each other than to gene copies of related species. Such convergent evolution is most easily explained by recurrent events of gene conversion, the rate of which we estimated at 10(-6)-10(-5) per site and generation. A significantly higher GC content of HINTW than of other W-linked genes is consistent with biased gene conversion increasing the fixation probability of mutations involving G and C nucleotides. Furthermore, and as a likely consequence, the neutral substitution rate is almost twice as high in HINTW as in other W-linked genes. The region on W encompassing the HINTW gene cluster is not covered in the initial assembly of the chicken genome, but analysis of raw sequence reads indicates that gene copy number is significantly higher than a previous estimate of 40. While sexual selection is one of several factors that potentially affect the evolution of ampliconic, male-specific genes on the mammalian Y chromosome, data from HINTW provide evidence that gene amplification followed by gene conversion can evolve in female-specific chromosomes in the absence of sexual selection. The presence of multiple and highly similar copies of HINTW may be related to protein function, but, more generally, amplification and conversion offers a means to the avoidance of accumulation of deleterious mutations in nonrecombining chromosomes.     

Adaptive molecular evolution of HINTW,a female-specific gene in birds

It is well established that many genes on the male-specific Y chromosome of organisms such as mammals are involved in male reproduction and may evolve rapidly because of positive selection on male reproductive traits. In contrast, very little is known about the function and evolution of W-linked genes restricted to the female genome of organisms with female heterogamety. For birds (males ZZ, females ZW), only one W-linked gene (HINTW) is sufficiently different from its Z-linked homolog to indicate a female-specific function. Here, we report that HINTW shows evidence of adaptive molecular evolution, implying strong positive selection for new functional properties in female birds. Moreover, because HINTW is expressed in the gonads of female birds just before sexual differentiation and is thus a candidate for sex determination, it suggests adaptive evolution related to female development. This provides the first example of Darwinian evolution of a gene restricted to the female genome of any organism. Given that HINTW exists in multiple copies on W, similar to some testis-specific genes amplified on mammalian Y, avian HINTW may thus potentially represent a female parallel to the organization and evolution of Y chromosome genes involved in male reproduction and development.     

Deleterious Mutations Accumulate Faster in Allopolyploid Than Diploid Cotton (Gossypium) and Unequally between Subgenomes

Whole-genome duplication (polyploidization) is among the most dramatic mutational processes in nature, so understanding how natural selection differs in polyploids relative to diploids is an important goal. Population genetics theory predicts that recessive deleterious mutations accumulate faster in allopolyploids than diploids due to the masking effect of redundant gene copies, but this prediction is hitherto unconfirmed. Here, we use the cotton genus (Gossypium), which contains seven allopolyploids derived from a single polyploidization event 1–2 Million years ago, to investigate deleterious mutation accumulation. We use two methods of identifying deleterious mutations at the nucleotide and amino acid level, along with whole-genome resequencing of 43 individuals spanning six allopolyploid species and their two diploid progenitors, to demonstrate that deleterious mutations accumulate faster in allopolyploids than in their diploid progenitors. We find that, unlike what would be expected under models of demographic changes alone, strongly deleterious mutations show the biggest difference between ploidy levels, and this effect diminishes for moderately and mildly deleterious mutations. We further show that the proportion of nonsynonymous mutations that are deleterious differs between the two coresident subgenomes in the allopolyploids, suggesting that homoeologous masking acts unequally between subgenomes. Our results provide a genome-wide perspective on classic notions of the significance of gene duplication that likely are broadly applicable to allopolyploids, with implications for our understanding of the evolutionary fate of deleterious mutations. Finally, we note that some measures of selection (e.g., dN/dS, π_N/π_S) may be biased when species of different ploidy levels are compared.     

Male-Driven Evolution Among Eoaves? A Test of the Replicative Division Hypothesis in a Heterogametic Female (ZW) System

Because avian females are heterogametic, the reverse of mammals, avian sex chromosomes undergo significantly different patterns and numbers of DNA replications than do those in mammals. This makes the W (female-specific) and the Z chromosomes an excellent model system for the study of the replicative division hypothesis, which purports that DNA substitution rate is determined by the number of germline replications. The sex-specific chromosome in birds (the W) is predicted to change at the slowest rate of all avian chromosomes because it undergoes the fewest rounds of replication per unit of evolutionary time. Using published data on gametogenesis from a variety of sources, we estimated the ratio of male-to-female germline replications (c) in galliforms and anseriforms to be approximately 4.4. The value of c should predict the value of the ratio of male-to-female mutation rates (α_m) if the replicative division hypothesis is true. Homologous DNA sequences including an intron and parts of two exons of the CHD gene were obtained from the W and the Z chromosomes in ostrich, sage grouse, canvasback duck, tundra swan, and snow goose. The exons show significantly different nucleotide composition from the introns, and the W-linked exons show evidence of relaxed constraint. The Z-linked intron is diverging ≈ 3.1 times faster than the W-linked intron. From this, α_m was calculated to be approximately 4.1, with a confidence interval of 3.1 to 5.1. The data support the idea that the number of replicative divisions is a major determinant of substitution rate in the Eoavian genome. Received: 19 January 1999 / Accepted: 8 June 1999     

SEX‐CHROMOSOME TURNOVERS INDUCED BY DELETERIOUS MUTATION LOAD

In sharp contrast with mammals and birds, many cold‐blooded vertebrates present homomorphic sex chromosomes. Empirical evidence supports a role for frequent turnovers, which replace nonrecombining sex chromosomes before they have time to decay. Three main mechanisms have been proposed for such turnovers, relying either on neutral processes, sex‐ratio selection, or intrinsic benefits of the new sex‐determining genes (due, e.g., to linkage with sexually antagonistic mutations). Here, we suggest an additional mechanism, arising from the load of deleterious mutations that accumulate on nonrecombining sex chromosomes. In the absence of dosage compensation, this load should progressively lower survival rate in the heterogametic sex. Turnovers should occur when this cost outweighs the benefits gained from any sexually antagonistic genes carried by the nonrecombining sex chromosome. We use individual‐based simulations of a Muller's ratchet process to test this prediction, and investigate how the relevant parameters (effective population size, strength and dominance of deleterious mutations, size of nonrecombining segment, and strength of sexually antagonistic selection) are expected to affect the rate of turnovers.     

The evolutionary landscape of the Mycobacterium tuberculosis genome

Mycobacterium tuberculosis is one of the most deadly human pathogens. The major mechanism for the adaptations of M. tuberculosis is nucleotide substitution. Previous studies have relied on the nonsynonymous-to-synonymous substitution rate (d_N/d_S) ratio as a measurement of selective constraint based on the assumed selective neutrality of synonymous substitutions. However, this assumption has been shown to be untrue in many cases. In this study, we used the substitution rate in intergenic regions (d_i) of the M. tuberculosis genome as the neutral reference, and conducted a genome-wide profiling for d_i, d_S, and the rate of insertions/deletions (indel rate) as compared with the genome of M. canettii using a 50 kb sliding window. We demonstrate significant variations in all of the three evolutionary measurements across the M. tuberculosis genome, even for regions in close vicinity. Furthermore, we identified a total of 233 genes with their d_S deviating significantly from d_i within the same window. Interestingly, d_S also varies significantly in some of the windows, indicating drastic changes in mutation rate and/or selection pressure within relatively short distances in the M. tuberculosis genome. Importantly, our results indicate that selection on synonymous substitutions is common in the M. tuberculosis genome. Therefore, the d_N/d_S ratio test must be applied carefully for measuring selection pressure on M. tuberculosis genes.     

Life history and the male mutation bias

Abstract If DNA replication is a major cause of mutation, then those life-history characters, which are expected to affect the number of male germline cell divisions, should also affect the male to female mutation bias (a_m). We tested this hypothesis by comparing several clades of bird species, which show variation both in suitable life-history characters (generation time as measured by age at first breeding and sexual selection as measured by frequency of extrapair paternity) and in a_m, which was estimated by comparing Z-linked and W-linked substitution rates in gametologous introns. a_m differences between clades were found to positively covary with both generation time and sexual selection, as expected if DNA replication causes mutation. The effects of extrapair paternity frequency on a_m suggests that increased levels of sexual selection cause higher mutation rates, which offers an interesting solution to the paradox of the loss of genetic variance associated with strong directional sexual selection. We also used relative rate tests to examine whether the observed differences in a_m between clades were due to differences in W-linked or Z-linked substitution rates. In one case, a significant difference in a_m between two clades was shown to be due to W-linked rates and not Z-linked rates, a result that suggests that mutation rates are not determined by replication alone.     

Tempo of Degeneration Across Independently Evolved Nonrecombining Regions

Recombination is beneficial over the long term, allowing more effective selection. Despite long-term advantages of recombination, local recombination suppression can evolve and lead to genomic degeneration, in particular on sex chromosomes. Here, we investigated the tempo of degeneration in nonrecombining regions, that is, the function curve for the accumulation of deleterious mutations over time, leveraging on 22 independent events of recombination suppression identified on mating-type chromosomes of anther-smut fungi, including newly identified ones. Using previously available and newly generated high-quality genome assemblies of alternative mating types of 13 Microbotryum species, we estimated degeneration levels in terms of accumulation of nonoptimal codons and nonsynonymous substitutions in nonrecombining regions. We found a reduced frequency of optimal codons in the nonrecombining regions compared with autosomes, that was not due to less frequent GC-biased gene conversion or lower ancestral expression levels compared with recombining regions. The frequency of optimal codons rapidly decreased following recombination suppression and reached an asymptote after ca. 3 Ma. The strength of purifying selection remained virtually constant at d_N/d_S = 0.55, that is, at an intermediate level between purifying selection and neutral evolution. Accordingly, nonsynonymous differences between mating-type chromosomes increased linearly with stratum age, at a rate of 0.015 per My. We thus develop a method for disentangling effects of reduced selection efficacy from GC-biased gene conversion in the evolution of codon usage and we quantify the tempo of degeneration in nonrecombining regions, which is important for our knowledge on genomic evolution and on the maintenance of regions without recombination.     

Hens,cocks and avian sex determination: A quest for genes on Z or W?

The sex of an individual is generally determined genetically by genes on one of the two sex chromosomes. In mammals, for instance, the presence of the male-specific Y chromosome confers maleness, whereas in Drosophila melanogaster and Caenorhabditis elegans it is the number of X chromosomes that matters. For birds (males ZZ, females ZW), however, the situation remains unclear. The recent discovery that the Z-linked DMRT1 gene, which is conserved across phyla as a gene involved in sexual differentiation, is expressed early in male development suggests that it might be the number of Z chromosomes that regulate sex in birds. On the other hand, the recent identification of the first protein unique to female birds, encoded by the W-linked PKCIW gene, and the observation that it is expressed early in female gonads, suggests that the W chromosome plays a role in avian sexual differentiation. Clearly defining the roles of the DMRT1 and PKC1W genes in gonadal development, and ultimately determining whether avian sex is dependent on Z or W, will require transgenic experiments.     

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.     

Genetic mapping of Z chromosome and identification of W chromosome-specific markers in the silkworm, Bombyx mori

In the silkworm, Bombyx mori, the female is the heterogametic (ZW) sex and the male is homogametic (ZZ). The female heterogamety is a typical situation in the insect order Lepidoptera. Although the W chromosome in silkworm is strongly female determining, no W-linked gene for a morphological character has been found on it. The Z chromosome carries important traits of economic value as well as genes for various phenotypic traits, but only 2% of molecular information based on its relative size is known. Studies conducted so far indicate that the Z-linked genes are not dosage compensated. In the present study, we constructed a genetic map of randomly amplified polymorphic DNA fragments (RAPD), simple sequence repeats (SSR), and fluorescent intersimple sequence repeat PCR (FISSR) markers for the Z chromosome using a backcross mapping population. A total of 16 Z-linked markers were identified, characterized, and mapped using od, a recessive trait for translucent skin as an anchor marker yielding a total recombination map of 334.5 cM. The linkage distances obtained suggested that the markers were distributed throughout the Z chromosome. Four RAPD and four SSR markers that were linked to W chromosome were also identified. The proposed mapping approach should be useful to identify and map sex-linked traits in the silkworm. The economic and evolutionary significance of Z- and W-linked genes in silkworm, in particular, and lepidopterans, in general, is discussed.     

Do Avian Mitochondria Recombine?

The dogma of strict maternal inheritance of mitochondria is now being tested with population genetics methods on sequence data from many species. In this study we investigated whether recombination occurs in the mitochondria of the blue tit (Parus caeruleus) by studying polymorphisms in the mitochondrial control region and in a recently identified (A)_n microsatellite on the W chromosome. The female heterogamety of avian sex chromosomes allows a test of whether mitochondrial recombination affects genealogical inference by comparison of mitochondrial and W-linked sequence variation. There is no discrepancy between mitochondrial and W-linked genealogies in blue tits, consistent with no recombination. We also analyzed mitochondrial sequence variation in both blue tits and peregrine falcons (Falco peregrinus) using a coalescent-based approach which accounts for recurrent mutation; in neither bird species did we find evidence of recombination. We conclude that it is unlikely that mitochondrial recombination has large effects on mitochondrial genetic variability in birds.     

Cotton rat<Emphasis Type="Italic"> Sihi-</Emphasis>M3 is a minimally oligomorphic<Emphasis Type="Italic"> Mhc</Emphasis> I-b molecule that binds the chemotactic peptide fMLF under stringent conditions

The leading model for class I-b evolution suggests non-polymorphic I-b genes evolve by gene duplication from polymorphic I-a genes. We recently found N-formyl peptide-specific orthologs of the class I-b gene H2-M3 in the rodent subfamily Sigmodontinae. To test if sigmodont M3 is a I-b gene, we sequenced M3 from wild cotton rats (Sigmodon hispidus) diverse at the class II locus, Sihi-DQA. These haplotypes carry a single allele of M3 that closely resembles H2-M3. However, peptide-binding assays showed that cotton rat M3 bound the chemotactic N-formylpeptide fMLF better than did rat or mouse M3. The Ala116Lys substitution in cotton rat M3 might enhance binding of fMLF and is one of eight residues of M3 that interact with ligand residues P3 and P4 and that are positively selected, with a d_N/d_S ratio of 1.8. Thus, M3 is a class I-b gene in both sigmodontine and murine murids, but positive selection operates on a small subset of residues in the traditionally defined antigen recognition site.Electronic Supplementary Material Supplementary material is available in the online version of this article at     

Nonrandom Representation of Sex-Biased Genes on Chicken Z Chromosome

Several lines of evidence suggest that the X chromosome of various animal species has an unusual complement of genes with sex-biased or sex-specific expression. However, the study of the X chromosome gene content in different organisms provided conflicting results. The most striking contrast concerns the male-biased genes, which were reported to be almost depleted from the X chromosome in Drosophila but overrepresented on the X chromosome in mammals. To elucidate the reason for these discrepancies, we analysed the gene content of the Z chromosome in chicken. Our analysis of the publicly available expressed sequence tags (EST) data and genome draft sequence revealed a significant underrepresentation of ovary-specific genes on the chicken Z chromosome. For the brain-expressed genes, we found a significant enrichment of male-biased genes but an indication of underrepresentation of female-biased genes on the Z chromosome. This is the first report on the nonrandom gene content in a homogametic sex chromosome of a species with heterogametic female individuals. Further comparison of gene contents of the independently evolved X and Z sex chromosomes may offer new insight into the evolutionary processes leading to the nonrandom genomic distribution of sex-biased and sex-specific genes. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users. [Reviewing Editor: Dr. Manyuan Long]