Lack of nucleotide polymorphism in the Y-linked sperm flagellar dynein gene Dhc-Yh3 of Drosophila melanogaster and D. simulans

Transposable element accumulation drives size differences among polymorphic Y Chromosomes in Drosophila

Y Chromosomes of many species are gene poor and show low levels of nucleotide variation, yet often display high amounts of structural diversity. Dobzhansky cataloged several morphologically distinct Y Chromosomes in Drosophila pseudoobscura that differ in size and shape, but the molecular causes of their dramatic size differences are unclear. Here we use cytogenetics and long-read sequencing to study the sequence content of polymorphic Y Chromosomes in D. pseudoobscura We show that Y Chromosomes differ almost 2-fold in size, ranging from 30 to 60 Mb. Most of this size difference is caused by a handful of active transposable elements (TEs) that have recently expanded on the largest Y Chromosome, with different elements being responsible for Y expansion on differently sized D. pseudoobscura Y's. We show that Y Chromosomes differ in their heterochromatin enrichment, expression of Y-enriched TEs, and also influence expression of dozens of autosomal and X-linked genes. The same helitron element that showed the most drastic amplification on the largest Y in D. pseudoobscura independently amplified on a polymorphic large Y Chromosome in D. affinis, suggesting that some TEs are inherently more prone to become deregulated on Y Chromosomes.

Sex-Ratio Meiotic Drive Shapes the Evolution of the Y Chromosome in Drosophila simulans

The recent emergence and spread of X-linked segregation distorters-called "Paris" system-in the worldwide species Drosophila simulans has elicited the selection of drive-resistant Y chromosomes. Here, we investigate the evolutionary history of 386 Y chromosomes originating from 29 population samples collected over a period of 20 years, showing a wide continuum of phenotypes when tested against the Paris distorters, from high sensitivity to complete resistance (males sire ∼95% to ∼40% female progeny). Analyzing around 13 kb of Y-linked gene sequences in a representative subset of nine Y chromosomes, we identified only three polymorphic sites resulting in three haplotypes. Remarkably, one of the haplotypes is associated with resistance. This haplotype is fixed in all samples from Sub-Saharan Africa, the region of origin of the drivers. Exceptionally, with the spread of the drivers in Egypt and Morocco, we were able to record the replacement of the sensitive lineage by the resistant haplotype in real time, within only a few years. In addition, we performed in situ hybridization, using satellite DNA probes, on a subset of 21 Y chromosomes from six locations. In contrast to the low molecular polymorphism, this revealed extensive structural variation suggestive of rapid evolution, either neutral or adaptive. Moreover, our results show that intragenomic conflicts can drive astonishingly rapid replacement of Y chromosomes and suggest that the emergence of Paris segregation distorters in East Africa occurred less than half a century ago.

The Nix locus on the male-specific homologue of chromosome 1 in Aedes albopictus is a strong candidate for a male-determining factor

Background

Global concern over the rapid expansion of the Asian tiger mosquito, Aedes albopictus, and its vector competence has highlighted an urgent need to improve currently available population control methods, like the Sterile Insect Technique. Knowledge of the sex determination cascade is a prerequisite for the development of early-stage sexing systems. To this end, we have characterised the putative sex determination gene, Nix, in this species. In Aedes species the chromosome complement consists of three pairs of chromosomes. The sex determination alleles are linked to the smallest homomorphic chromosome.

Results

We identified the male-specific chromosome 1 of Ae. albopictus that carries the putative male-determining gene Nix. We have also characterised the complete genomic sequence of the Nix gene which is composed of two exons and a short intron. The gene displays different levels of intron retention during development. Comparison of DNA sequences covering most of the Nix gene from individuals across the species range revealed no polymorphism.

Conclusions

Our characterisation of the Nix gene in Ae. albopictus represents an initial step in the analysis of the sex determination cascade in this species. We found evidence of intron retention (IR) in Nix. IR might play a role in regulating the expression of Nix during development. Our results provide the basis for the development of new genetic control strategies.

Electronic supplementary material

The online version of this article (10.1186/s13071-018-3215-8) contains supplementary material, which is available to authorized users.

Heterochromatin-Enriched Assemblies Reveal the Sequence and Organization of the Drosophila melanogaster Y Chromosome

Heterochromatic regions of the genome are repeat-rich and poor in protein coding genes, and are therefore underrepresented in even the best genome assemblies. One of the most difficult regions of the genome to assemble are sex-limited chromosomes. The Drosophila melanogaster Y chromosome is entirely heterochromatic, yet has wide-ranging effects on male fertility, fitness, and genome-wide gene expression. The genetic basis of this phenotypic variation is difficult to study, in part because we do not know the detailed organization of the Y chromosome. To study Y chromosome organization in D. melanogaster, we develop an assembly strategy involving the in silico enrichment of heterochromatic long single-molecule reads and use these reads to create targeted de novo assemblies of heterochromatic sequences. We assigned contigs to the Y chromosome using Illumina reads to identify male-specific sequences. Our pipeline extends the D. melanogaster reference genome by 11.9 Mb, closes 43.8% of the gaps, and improves overall contiguity. The addition of 10.6 MB of Y-linked sequence permitted us to study the organization of repeats and genes along the Y chromosome. We detected a high rate of duplication to the pericentric regions of the Y chromosome from other regions in the genome. Most of these duplicated genes exist in multiple copies. We detail the evolutionary history of one sex-linked gene family, crystal-Stellate While the Y chromosome does not undergo crossing over, we observed high gene conversion rates within and between members of the crystal-Stellate gene family, Su(Ste), and PCKR, compared to genome-wide estimates. Our results suggest that gene conversion and gene duplication play an important role in the evolution of Y-linked genes.

Genetic Diversity on the Sex Chromosomes

Levels and patterns of genetic diversity can provide insights into a population’s history. In species with sex chromosomes, differences between genomic regions with unique inheritance patterns can be used to distinguish between different sets of possible demographic and selective events. This review introduces the differences in population history for sex chromosomes and autosomes, provides the expectations for genetic diversity across the genome under different evolutionary scenarios, and gives an introductory description for how deviations in these expectations are calculated and can be interpreted. Predominantly, diversity on the sex chromosomes has been used to explore and address three research areas: 1) Mating patterns and sex-biased variance in reproductive success, 2) signatures of selection, and 3) evidence for modes of speciation and introgression. After introducing the theory, this review catalogs recent studies of genetic diversity on the sex chromosomes across species within the major research areas that sex chromosomes are typically applied to, arguing that there are broad similarities not only between male-heterogametic (XX/XY) and female-heterogametic (ZZ/ZW) sex determination systems but also any mating system with reduced recombination in a sex-determining region. Further, general patterns of reduced diversity in nonrecombining regions are shared across plants and animals. There are unique patterns across populations with vastly different patterns of mating and speciation, but these do not tend to cluster by taxa or sex determination system.

The Y Chromosome Modulates Splicing and Sex-Biased Intron Retention Rates in Drosophila

The Drosophila Y chromosome is a 40-Mb segment of mostly repetitive DNA; it harbors a handful of protein-coding genes and a disproportionate amount of satellite repeats, transposable elements, and multicopy DNA arrays. Intron retention (IR) is a type of alternative splicing (AS) event by which one or more introns remain within the mature transcript. IR recently emerged as a deliberate cellular mechanism to modulate gene expression levels and has been implicated in multiple biological processes. However, the extent of sex differences in IR and the contribution of the Y chromosome to the modulation of AS and IR rates has not been addressed. Here we showed pervasive IR in the fruit fly Drosophila melanogaster with thousands of novel IR events, hundreds of which displayed extensive sex bias. The data also revealed an unsuspected role for the Y chromosome in the modulation of AS and IR. The majority of sex-biased IR events introduced premature termination codons and the magnitude of sex bias was associated with gene expression differences between the sexes. Surprisingly, an extra Y chromosome in males (X^YY genotype) or the presence of a Y chromosome in females (X^XY genotype) significantly modulated IR and recapitulated natural differences in IR between the sexes. Our results highlight the significance of sex-biased IR in tuning sex differences and the role of the Y chromosome as a source of variable IR rates between the sexes. Modulation of splicing and IR rates across the genome represent new and unexpected outcomes of the Drosophila Y chromosome.

Interspecific Y chromosome variation is sufficient to rescue hybrid male sterility and is influenced by the grandparental origin of the chromosomes

Y chromosomes display population variation within and between species. Co-evolution within populations is expected to produce adaptive interactions between Y chromosomes and the rest of the genome. One consequence is that Y chromosomes from disparate populations could disrupt harmonious interactions between co-evolved genetic elements and result in reduced male fertility, sterility or inviability. Here we address the contribution of 'heterospecific Y chromosomes' to fertility in hybrid males carrying a homozygous region of Drosophila mauritiana introgressed in the Drosophila simulans background. In order to detect Y chromosome-autosome interactions, which may go unnoticed in a single-species background of autosomes, we constructed hybrid genotypes involving three sister species: Drosophila simulans, D. mauritiana, and D. sechellia. These engineered strains varied due to: (i) species origin of the Y chromosome (D. simulans or D. sechellia); (ii) location of the introgressed D. mauritiana segment on the D. simulans third chromosome, and (iii) grandparental genomic background (three genotypes of D. simulans). We find complex interactions between the species origin of the Y chromosome, the identity of the D. mauritiana segment and the grandparental genetic background donating the chromosomes. Unexpectedly, the interaction of the Y chromosome and one segment of D. mauritiana drastically reduced fertility in the presence of Ysim, whereas the fertility is partially rescued by the Y chromosome of D. sechellia when it descends from a specific grandparental genotype. The restoration of fertility occurs in spite of an autosomal and X-linked genome that is mostly of D. simulans origin. These results illustrate the multifactorial basis of genetic interactions involving the Y chromosome. Our study supports the hypothesis that the Y chromosome can contribute significantly to the evolution of reproductive isolation and highlights the conditional manifestation of infertility in specific genotypic combinations.

Positive and purifying selection on the Drosophila Y chromosome

Y chromosomes, with their reduced effective population size, lack of recombination, and male-limited transmission, present a unique collection of constraints for the operation of natural selection. Male-limited transmission may greatly increase the efficacy of selection for male-beneficial mutations, but the reduced effective size also inflates the role of random genetic drift. Together, these defining features of the Y chromosome are expected to influence rates and patterns of molecular evolution on the Y as compared with X-linked or autosomal loci. Here, we use sequence data from 11 genes in 9 Drosophila species to gain insight into the efficacy of natural selection on the Drosophila Y relative to the rest of the genome. Drosophila is an ideal system for assessing the consequences of Y-linkage for molecular evolution in part because the gene content of Drosophila Y chromosomes is highly dynamic, with orthologous genes being Y-linked in some species whereas autosomal in others. Our results confirm the expectation that the efficacy of natural selection at weakly selected sites is reduced on the Y chromosome. In contrast, purifying selection on the Y chromosome for strongly deleterious mutations does not appear to be compromised. Finally, we find evidence of recurrent positive selection for 4 of the 11 genes studied here. Our results thus highlight the variable nature of the mode and impact of natural selection on the Drosophila Y chromosome.

Introgression in the Drosophila subobscura--D. Madeirensis sister species: evidence of gene flow in nuclear genes despite mitochondrial differentiation

Species hybridization, and thus the potential for gene flow, was once viewed as reproductive mistake. However, recent analysis based on large datasets and newly developed models suggest that gene exchange is not as rare as originally suspected. To investigate the history and speciation of the closely related species Drosophila subobscura, D. madeirensis, and D. guanche, we obtained polymorphism and divergence data for 26 regions throughout the genome, including the Y chromosome and mitochondrial DNA. We found that the D. subobscura X/autosome ratio of silent nucleotide diversity is significantly smaller than the 0.75 expected under neutrality. This pattern, if held genomewide, may reflect a faster accumulation of beneficial mutations on the X chromosome than on autosomes. We also detected evidence of gene flow in autosomal regions, while sex chromosomes remain distinct. This is consistent with the large X effect on hybrid male sterility seen in this system and the presence of two X chromosome inversions fixed between species. Overall, our data conform to chromosomal speciation models in which rearrangements are proposed to serve as gene flow barriers. Contrary to other observations in Drosophila, the mitochondrial genome appears resilient to gene flow in the presence of nuclear exchange.

Meta-analysis reveals that genes regulated by the Y chromosome in Drosophila melanogaster are preferentially localized to repressive chromatin

The Drosophila Y chromosome is a degenerated, heterochromatic chromosome with few functional genes. Despite this, natural variation on the Y chromosome in D. melanogaster has substantial trans-acting effects on the regulation of X-linked and autosomal genes. It is not clear, however, whether these genes simply represent a random subset of the genome or whether specific functional properties are associated with susceptibility to regulation by Y-linked variation. Here, we present a meta-analysis of four previously published microarray studies of Y-linked regulatory variation (YRV) in D. melanogaster. We show that YRV genes are far from a random subset of the genome: They are more likely to be in repressive chromatin contexts, be expressed tissue specifically, and vary in expression within and between species than non-YRV genes. Furthermore, YRV genes are more likely to be associated with the nuclear lamina than non-YRV genes and are generally more likely to be close to each other in the nucleus (although not along chromosomes). Taken together, these results suggest that variation on the Y chromosome plays a role in modifying how the genome is distributed across chromatin compartments, either via changes in the distribution of DNA-binding proteins or via changes in the spatial arrangement of the genome in the nucleus.

Chromatin-associated proteins HP1 and Mod(mdg4) modify Y-linked regulatory variation in the drosophila testis

Chromatin remodeling is crucial for gene regulation. Remodeling is often mediated through chemical modifications of the DNA template, DNA-associated proteins, and RNA-mediated processes. Y-linked regulatory variation (YRV) refers to the quantitative effects that polymorphic tracts of Y-linked chromatin exert on gene expression of X-linked and autosomal genes. Here we show that naturally occurring polymorphisms in the Drosophila melanogaster Y chromosome contribute disproportionally to gene expression variation in the testis. The variation is dependent on wild-type expression levels of mod(mdg4) as well as Su(var)205; the latter gene codes for heterochromatin protein 1 (HP1) in Drosophila. Testis-specific YRV is abolished in genotypes with heterozygous loss-of-function mutations for mod(mdg4) and Su(var)205 but not in similar experiments with JIL-1. Furthermore, the Y chromosome differentially regulates several ubiquitously expressed genes. The results highlight the requirement for wild-type dosage of Su(var)205 and mod(mdg4) in enabling naturally occurring Y-linked regulatory variation in the testis. The phenotypes that emerge in the context of wild-type levels of the HP1 and Mod(mdg4) proteins might be part of an adaptive response to the environment.

Natural variation of the Y chromosome suppresses sex ratio distortion and modulates testis-specific gene expression in Drosophila simulans

X-linked sex-ratio distorters that disrupt spermatogenesis can cause a deficiency in functional Y-bearing sperm and a female-biased sex ratio. Y-linked modifiers that restore a normal sex ratio might be abundant and favored when a X-linked distorter is present. Here we investigated natural variation of Y-linked suppressors of sex-ratio in the Winters systems and the ability of these chromosomes to modulate gene expression in Drosophila simulans. Seventy-eight Y chromosomes of worldwide origin were assayed for their resistance to the X-linked sex-ratio distorter gene Dox. Y chromosome diversity caused males to sire ∼63% to ∼98% female progeny. Genome-wide gene expression analysis revealed hundreds of genes differentially expressed between isogenic males with sensitive (high sex ratio) and resistant (low sex ratio) Y chromosomes from the same population. Although the expression of about 75% of all testis-specific genes remained unchanged across Y chromosomes, a subset of post-meiotic genes was upregulated by resistant Y chromosomes. Conversely, a set of accessory gland-specific genes and mitochondrial genes were downregulated in males with resistant Y chromosomes. The D. simulans Y chromosome also modulated gene expression in XXY females in which the Y-linked protein-coding genes are not transcribed. The data suggest that the Y chromosome might exert its regulatory functions through epigenetic mechanisms that do not require the expression of protein-coding genes. The gene network that modulates sex ratio distortion by the Y chromosome is poorly understood, other than that it might include interactions with mitochondria and enriched for genes expressed in post-meiotic stages of spermatogenesis.

Surprising differences in the variability of Y chromosomes in African and cosmopolitan populations of Drosophila melanogaster

The nonrecombining Drosophila melanogaster Y chromosome is heterochromatic and has few genes. Despite these limitations, there remains ample opportunity for natural selection to act on the genes that are vital for male fertility and on Y factors that modulate gene expression elsewhere in the genome. Y chromosomes of many organisms have low levels of nucleotide variability, but a formal survey of D. melanogaster Y chromosome variation had yet to be performed. Here we surveyed Y-linked variation in six populations of D. melanogaster spread across the globe. We find surprisingly low levels of variability in African relative to Cosmopolitan (i.e., non-African) populations. While the low levels of Cosmopolitan Y chromosome polymorphism can be explained by the demographic histories of these populations, the staggeringly low polymorphism of African Y chromosomes cannot be explained by demographic history. An explanation that is entirely consistent with the data is that the Y chromosomes of Zimbabwe and Uganda populations have experienced recent selective sweeps. Interestingly, the Zimbabwe and Uganda Y chromosomes differ: in Zimbabwe, a European Y chromosome appears to have swept through the population.

Y chromosome mediates ribosomal DNA silencing and modulates the chromatin state in Drosophila

Although the Drosophila Y chromosome is degenerated, heterochromatic, and contains few genes, increasing evidence suggests that it plays an important role in regulating the expression of numerous autosomal and X-linked genes. Here we use 15 Y chromosomes originating from a single founder 550 generations ago to study the role of the Y chromosome in regulating rRNA gene transcription, position-effect variegation (PEV), and the link among rDNA copy number, global gene expression, and chromatin regulation. Based on patterns of rRNA gene transcription indicated by transcription of the retrotransposon R2 that specifically inserts into the 28S rRNA gene, we show that X-linked rDNA is silenced in males. The silencing of X-linked rDNA expression by the Y chromosome is consistent across populations and independent of genetic background. These Y chromosomes also vary more than threefold in rDNA locus size and cause dramatically different levels of PEV suppression. The degree of suppression is negatively associated with the number and fraction of rDNA units without transposon insertions, but not with total rDNA locus size. Gene expression profiling revealed hundreds of differentially expressed genes among these Y chromosome introgression lines, as well as a divergent global gene expression pattern between the low-PEV and high-PEV flies. Our findings suggest that the Y chromosome is involved in diverse phenomena related to transcriptional regulation including X-linked rDNA silencing and suppression of PEV phenotype. These results further expand our understanding of the role of the Y chromosome in modulating global gene expression, and suggest a link with modifications of the chromatin state.

How closely does genetic diversity in finite populations conform to predictions of neutral theory? Large deficits in regions of low recombination

Levels of genetic diversity in finite populations are crucial in conservation and evolutionary biology. Genetic diversity is required for populations to evolve and its loss is related to inbreeding in random mating populations, and thus to reduced population fitness and increased extinction risk. Neutral theory is widely used to predict levels of genetic diversity. I review levels of genetic diversity in finite populations in relation to predictions of neutral theory. Positive associations between genetic diversity and population size, as predicted by neutral theory, are observed for microsatellites, allozymes, quantitative genetic variation and usually for mitochondrial DNA (mtDNA). However, there are frequently significant deviations from neutral theory owing to indirect selection at linked loci caused by balancing selection, selective sweeps and background selection. Substantially lower genetic diversity than predicted under neutrality was found for chromosomes with low recombination rates and high linkage disequilibrium (compared with 'normally' recombining chromosomes within species and adjusted for different copy numbers and mutation rates), including W (median 100% lower) and Y (89% lower) chromosomes, dot fourth chromosomes in Drosophila (94% lower) and mtDNA (67% lower). Further, microsatellite genetic and allelic diversity were lost at 12 and 33% faster rates than expected in populations adapting to captivity, owing to widespread selective sweeps. Overall, neither neutral theory nor most versions of the genetic draft hypothesis are compatible with all empirical results.

Interspecific Y chromosome introgressions disrupt testis-specific gene expression and male reproductive phenotypes in Drosophila

The Drosophila Y chromosome is a degenerated, heterochromatic chromosome with few functional genes. Nonetheless, natural variation on the Y chromosome in Drosophila melanogaster has substantial trans-acting effects on the regulation of X-linked and autosomal genes. However, the contribution of Y chromosome divergence to gene expression divergence between species is unknown. In this study, we constructed a series of Y chromosome introgression lines, in which Y chromosomes from either Drosophila sechellia or Drosophila simulans are introgressed into a common D. simulans genetic background. Using these lines, we compared genome-wide gene expression and male reproductive phenotypes between heterospecific and conspecific Y chromosomes. We find significant differences in expression for 122 genes, or 2.84% of all genes analyzed. Genes down-regulated in males with heterospecific Y chromosomes are significantly biased toward testis-specific expression patterns. These same lines show reduced fecundity and sperm competitive ability. Taken together, these results imply a significant role for Y/X and Y/autosome interactions in maintaining proper expression of male-specific genes, either directly or via indirect effects on male reproductive tissue development or function.

Evolution of sex chromosomes in insects

Sex chromosomes have many unusual features relative to autosomes. Y (or W) chromosomes lack genetic recombination, are male- (female-) limited, and show an abundance of genetically inert heterochromatic DNA but contain few functional genes. X (or Z) chromosomes also show sex-biased transmission (i.e., X chromosomes show female-biased and Z-chromosomes show male-biased inheritance) and are hemizygous in the heterogametic sex. Their unusual ploidy level and pattern of inheritance imply that sex chromosomes play a unique role in many biological processes and phenomena, including sex determination, epigenetic chromosome-wide regulation of gene expression, the distribution of genes in the genome, genomic conflict, local adaptation, and speciation. The vast diversity of sex chromosome systems in insects--ranging from the classical male heterogametic XY system in Drosophila to ZW systems in Lepidoptera or mobile genes determining sex as found in house flies--implies that insects can serve as unique model systems to study various functional and evolutionary aspects of these different processes.

Molecular evolution of a Y chromosome to autosome gene duplication in Drosophila

In contrast to the rest of the genome, the Y chromosome is restricted to males and lacks recombination. As a result, Y chromosomes are unable to respond efficiently to selection, and newly formed Y chromosomes degenerate until few genes remain. The rapid loss of genes from newly formed Y chromosomes has been well studied, but gene loss from highly degenerate Y chromosomes has only recently received attention. Here, we identify and characterize a Y to autosome duplication of the male fertility gene kl-5 that occurred during the evolution of the testacea group species of Drosophila. The duplication was likely DNA based, as other Y-linked genes remain on the Y chromosome, the locations of introns are conserved, and expression analyses suggest that regulatory elements remain linked. Genetic mapping reveals that the autosomal copy of kl-5 resides on the dot chromosome, a tiny autosome with strongly suppressed recombination. Molecular evolutionary analyses show that autosomal copies of kl-5 have reduced polymorphism and little recombination. Importantly, the rate of protein evolution of kl-5 has increased significantly in lineages where it is on the dot versus Y linked. Further analyses suggest this pattern is a consequence of relaxed purifying selection, rather than adaptive evolution. Thus, although the initial fixation of the kl-5 duplication may have been advantageous, slightly deleterious mutations have accumulated in the dot-linked copies of kl-5 faster than in the Y-linked copies. Because the dot chromosome contains seven times more genes than the Y and is exposed to selection in both males and females, these results suggest that the dot suffers the deleterious effects of genetic linkage to more selective targets compared with the Y chromosome. Thus, a highly degenerate Y chromosome may not be the worst environment in the genome, as is generally thought, but may in fact be protected from the accumulation of deleterious mutations relative to other nonrecombining regions that contain more genes.

Epigenetic effects of polymorphic Y chromosomes modulate chromatin components, immune response, and sexual conflict

Genetic conflicts between sexes and generations provide a foundation for understanding the functional evolution of sex chromosomes and sexually dimorphic phenotypes. Y chromosomes of Drosophila contain multi-megabase stretches of satellite DNA repeats and a handful of protein-coding genes that are monomorphic within species. Nevertheless, polymorphic variation in heterochromatic Y chromosomes of Drosophila result in genome-wide gene expression variation. Here we show that such naturally occurring Y-linked regulatory variation (YRV) can be detected in somatic tissues and contributes to the epigenetic balance of heterochromatin/euchromatin at three distinct loci showing position-effect variegation (PEV). Moreover, polymorphic Y chromosomes differentially affect the expression of thousands of genes in XXY female genotypes in which Y-linked protein-coding genes are not transcribed. The data show a disproportionate influence of YRV on the variable expression of genes whose protein products localize to the nucleus, have nucleic-acid binding activity, and are involved in transcription, chromosome organization, and chromatin assembly. These include key components such as HP1, Trithorax-like (GAGA factor), Su(var)3-9, Brahma, MCM2, ORC2, and inner centromere protein. Furthermore, mitochondria-related genes, immune response genes, and transposable elements are also disproportionally affected by Y chromosome polymorphism. These functional clusterings may arise as a consequence of the involvement of Y-linked heterochromatin in the origin and resolution of genetic conflicts between males and females. Taken together, our results indicate that Y chromosome heterochromatin serves as a major source of epigenetic variation in natural populations that interacts with chromatin components to modulate the expression of biologically relevant phenotypic variation.

Y not a dead end: epistatic interactions between Y-linked regulatory polymorphisms and genetic background affect global gene expression in Drosophila melanogaster

The Y chromosome, inherited without meiotic recombination from father to son, carries relatively few genes in most species. This is consistent with predictions from evolutionary theory that nonrecombining chromosomes lack variation and degenerate rapidly. However, recent work has suggested a dynamic role for the Y chromosome in gene regulation, a finding with important implications for spermatogenesis and male fitness. We studied Y chromosomes from two populations of Drosophila melanogaster that had previously been shown to have major effects on the thermal tolerance of spermatogenesis. We show that these Y chromosomes differentially modify the expression of hundreds of autosomal and X-linked genes. Genes showing Y-linked regulatory variation (YRV) also show an association with immune response and pheromone detection. Indeed, genes located proximal to the euchromatin-heterochromatin boundary of the X chromosome appear particularly responsive to Y-linked variation, including a substantial number of odorant-binding genes. Furthermore, the data show significant regulatory interactions between the Y chromosome and the genetic background of autosomes and X chromosome. Altogether, our findings support the view that interpopulation, Y-linked regulatory polymorphisms can differentially modulate the expression of many genes important to male fitness, and they also point to complex interactions between the Y chromosome and genetic background affecting global gene expression.

Functional copies of the Mst77F gene on the Y chromosome of Drosophila melanogaster

The Y chromosome of Drosophila melanogaster has <20 protein-coding genes. These genes originated from the duplication of autosomal genes and have male-related functions. In 1993, Russell and Kaiser found three Y-linked pseudogenes of the Mst77F gene, which is a testis-expressed autosomal gene that is essential for male fertility. We did a thorough search using experimental and computational methods and found 18 Y-linked copies of this gene (named Mst77Y-1-Mst77Y-18). Ten Mst77Y genes encode defective proteins and the other eight are potentially functional. These eight genes produce approximately 20% of the functional Mst77F-like mRNA, and molecular evolutionary analysis shows that they evolved under purifying selection. Hence several Mst77Y genes have all the features of functional genes. Mst77Y genes are present only in D. melanogaster, and phylogenetic analysis confirmed that the duplication is a recent event. The identification of functional Mst77Y genes reinforces the previous finding that gene gains play a prominent role in the evolution of the Drosophila Y chromosome.

Pervasive natural selection in the Drosophila genome?

Over the past four decades, the predominant view of molecular evolution saw little connection between natural selection and genome evolution, assuming that the functionally constrained fraction of the genome is relatively small and that adaptation is sufficiently infrequent to play little role in shaping patterns of variation within and even between species. Recent evidence from Drosophila, reviewed here, suggests that this view may be invalid. Analyses of genetic variation within and between species reveal that much of the Drosophila genome is under purifying selection, and thus of functional importance, and that a large fraction of coding and noncoding differences between species are adaptive. The findings further indicate that, in Drosophila, adaptations may be both common and strong enough that the fate of neutral mutations depends on their chance linkage to adaptive mutations as much as on the vagaries of genetic drift. The emerging evidence has implications for a wide variety of fields, from conservation genetics to bioinformatics, and presents challenges to modelers and experimentalists alike.

The temporal dynamics of processes underlying Y chromosome degeneration

Y chromosomes originate from ordinary autosomes and degenerate by accumulating deleterious mutations. This accumulation results from a lack of recombination on the Y and is driven by interference among deleterious mutations (Muller's ratchet and background selection) and the fixation of beneficial alleles (genetic hitchhiking). Here I show that the relative importance of these processes is expected to vary over the course of Y chromosome evolution due to changes in the number of active genes. The dominant mode of degeneration on a newly formed gene-rich Y chromosome is expected to be Muller's ratchet and/or background selection due to the large numbers of deleterious mutations arising in active genes. However, the relative importance of these modes of degeneration declines rapidly as active genes are lost. In contrast, the rate of degeneration due to hitchhiking is predicted to be highest on Y chromosomes containing an intermediate number of active genes. The temporal dynamics of these processes imply that a gradual restriction of recombination, as inferred in mammals, will increase the importance of genetic hitchhiking relative to Muller's ratchet and background selection.

Active miniature transposons from a plant genome and its nonrecombining Y chromosome

Mechanisms involved in eroding fitness of evolving Y chromosomes have been the focus of much theoretical and empirical work. Evolving Y chromosomes are expected to accumulate transposable elements (TEs), but it is not known whether such accumulation contributes to their genetic degeneration. Among TEs, miniature inverted-repeat transposable elements are nonautonomous DNA transposons, often inserted in introns and untranslated regions of genes. Thus, if they invade Y-linked genes and selection against their insertion is ineffective, they could contribute to genetic degeneration of evolving Y chromosomes. Here, we examine the population dynamics of active MITEs in the young Y chromosomes of the plant Silene latifolia and compare their distribution with those in recombining genomic regions. To isolate active MITEs, we developed a straightforward approach on the basis of the assumption that recent transposon insertions or excisions create singleton or low-frequency size polymorphisms that can be detected in alleles from natural populations. Transposon display was then used to infer the distribution of MITE insertion frequencies. The overall frequency spectrum showed an excess of singleton and low-frequency insertions, which suggests that these elements are readily removed from recombining chromosomes. In contrast, insertions on the Y chromosomes were present at high frequencies. Their potential contribution to Y degeneration is discussed.

Polymorphic Y chromosomes harbor cryptic variation with manifold functional consequences

The paucity of polymorphisms in single-copy genes on the Y chromosome of Drosophila contrasts with data indicating that this chromosome has polymorphic phenotypic effects on sex ratio, temperature sensitivity, behavior, and fitness. We show that the Y chromosome of D. melanogaster harbors substantial genetic diversity in the form of polymorphisms for genetic elements that differentially affect the expression of hundreds of X-linked and autosomal genes. The affected genes are more highly expressed in males, more meagerly expressed in females, and more highly divergent between species. Functionally, they affect microtubule stability, lipid and mitochondrial metabolism, and the thermal sensitivity of spermatogenesis. Our findings provide a mechanism for adaptive phenotypic variation associated with the Y chromosome.

Hitchhiking effects of recurrent beneficial amino acid substitutions in the Drosophila melanogaster genome

Several recent studies have estimated that a large fraction of amino acid divergence between species of Drosophila was fixed by positive selection, using statistical approaches based on the McDonald-Kreitman test. However, little is known about associated selection coefficients of beneficial amino acid mutations. Recurrent selective sweeps associated with adaptive substitutions should leave a characteristic signature in genome variability data that contains information about the frequency and strength of selection. Here, I document a significant negative correlation between the level and the frequency of synonymous site polymorphism and the rate of protein evolution in highly recombining regions of the X chromosome of D. melanogaster. This pattern is predicted by recurrent adaptive protein evolution and suggests that adaptation is an important determinant of patterns of neutral variation genome-wide. Using a maximum likelihood approach, I estimate the product of the rate and strength of selection under a recurrent genetic hitchhiking model, lambda2N(e)s approximately 3 x 10(-8). Using an approach based on the McDonald-Kreitman test, I estimate that approximately 50% of divergent amino acids were driven to fixation by positive selection, implying that beneficial amino acid substitutions are of weak effect on average, on the order of 10(-5) (i.e., 2N(e)s approximately 40). Two implications of these results are that most adaptive substitutions will be difficult to detect in genome scans of selection and that population size (and genetic drift) may be an important determinant of the evolutionary dynamics of protein adaptation.

The Hill–Robertson effect: evolutionary consequences of weak selection and linkage in finite populations

The 'Hill-Robertson (HR) effect' describes that linkage between sites under selection will reduce the overall effectiveness of selection in finite populations. Here we discuss the major concepts associated with the HR effect and present results of computer simulations focusing on the linkage effects generated by multiple sites under weak selection. Most models of linkage and selection forecast differences in effectiveness of selection between chromosomes or chromosomal regions involving a number of genes. The abundance and physical clustering of weakly selected mutations across genomes, however, justify the investigation of HR effects at a very local level and we pay particular attention to linkage effects among selected sites of the same gene. Overall, HR effects caused by weakly selected mutations predict differences in effectiveness of selection between genes that differ in exon-intron structures and across genes. Under this scenario, introns might play an advantageous role reducing intragenic HR effects. Finally, we summarize observations that are consistent with local HR effects in Drosophila, discuss potential consequences on population genetic studies and suggest future lines of research.

Unusual pattern of nucleotide sequence variation at the OS-E and OS-F genomic regions of Drosophila simulans

Nucleotide variation at the genomic region encompassing the odorant-binding protein genes OS-E and OS-F (OS region) was surveyed in two populations of Drosophila simulans, one from Europe and the other from Africa. We found that the European population shows an atypical and large haplotype structure, which extends throughout the approximately 5-kb surveyed genomic region. This structure is depicted by two major haplotype groups segregating at intermediate frequency in the sample, one haplogroup with nearly no variation, and the other at levels more typical for this species. This pattern of variation was incompatible with neutral predictions for a population at a stationary equilibrium. Nevertheless, neutrality tests contrasting polymorphism and divergence data fail to detect any departure from the standard neutral model in this species, whereas they confirm the non-neutral behavior previously observed at the OS-E gene in D. melanogaster. Although positive Darwinian selection may have been responsible for the observed unusual nucleotide variation structure, coalescent simulation results do not allow rejecting the hypothesis that the pattern was generated by a recent bottleneck in the history of European populations of D. simulans.

Increased nucleotide diversity with transient Y linkage in Drosophila americana

Recombination shapes nucleotide variation within genomes. Patterns are thought to arise from the local recombination landscape, influencing the degree to which neutral variation experiences hitchhiking with selected variation. This study examines DNA polymorphism along Chromosome 4 (element B) of Drosophila americana to identify effects of hitchhiking arising as a consequence of Y-linked transmission. A centromeric fusion between the X and 4(th) chromosomes segregates in natural populations of D. americana. Frequency of the X-4 fusion exhibits a strong positive correlation with latitude, which has explicit consequences for unfused 4(th) chromosomes. Unfused Chromosome 4 exists as a non-recombining Y chromosome or as an autosome proportional to the frequency of the X-4 fusion. Furthermore, Y linkage along the unfused 4 is disrupted as a function of the rate of recombination with the centromere. Inter-population and intra-chromosomal patterns of nucleotide diversity were assayed using six regions distributed along unfused 4(th) chromosomes derived from populations with different frequencies of the X-4 fusion. No difference in overall level of nucleotide diversity was detected among populations, yet variation along the chromosome exhibits a distinct pattern in relation to the X-4 fusion. Sequence diversity is inflated at loci experiencing the strongest Y linkage. These findings are inconsistent with the expected reduction in nucleotide diversity resulting from hitchhiking due to background selection or selective sweeps. In contrast, excessive polymorphism is accruing in association with transient Y linkage, and furthermore, hitchhiking with sexually antagonistic alleles is potentially responsible.

Interspecific divergence, intrachromosomal recombination, and phylogenetic utility of Y-chromosomal genes in Drosophila

Reconstruction of phylogenetic relationships among recently diverged species is complicated by three general problems: segregation of polymorphisms that pre-date species divergence, gene flow during and after speciation, and intra-locus recombination. In light of these difficulties, the Y chromosome offers several important advantages over other genomic regions as a source of phylogenetic information. These advantages include the absence of recombination, rapid coalescence, and reduced opportunity for interspecific introgression due to hybrid male sterility. In this report, we test the phylogenetic utility of Y-chromosomal sequences in two groups of closely related and partially inter-fertile Drosophila species. In the D. bipectinata species complex, Y-chromosomal loci unambiguously recover the phylogeny most consistent with previous multi-locus analysis and with reproductive relationships, and show no evidence of either post-speciation gene flow or persisting ancestral polymorphisms. In the D. simulans species complex, the situation is complicated by the duplication of at least one Y-linked gene region, followed by intrachromosomal recombination between the duplicate genes that scrambles their genealogy. We suggest that Y-chromosomal sequences are a useful tool for resolving phylogenetic relationships among recently diverged species, especially in male-heterogametic organisms that conform to Haldane's rule. However, duplication of Y-linked genes may not be uncommon, and special care should be taken to distinguish between orthologous and paralogous sequences.

Male sterility at extreme temperatures: a significant but neglected phenomenon for understanding Drosophila climatic adaptations

The thermal range for viability is quite variable among Drosophila species and it has long been known that these variations are correlated with geographic distribution: temperate species are on average more cold tolerant but more heat sensitive than tropical species. At both ends of their viability range, sterile males have been observed in all species investigated so far. This symmetrical phenomenon restricts the temperature limits within which permanent cultures can be kept in the laboratory. Thermal heat sterility thresholds are very variable across species from 23 degrees C in heat sensitive species up to 31 degrees C in heat tolerant species. In Drosophila melanogaster, genetic variations are observed among geographic populations. Tropical populations are more tolerant to heat induced sterility and recover more rapidly than temperate ones. A genetic analysis revealed that about 50% of the difference observed between natural populations was due to the Y chromosome. Natural populations have not reached a selection limit, however: thermal tolerance was still increased by keeping strains at a high temperature, close to the sterility threshold. On the low temperature side, a symmetrical reverse phenomenon seems to exist: temperate populations are more tolerant to cold than tropical ones. Compared to Mammals, drosophilids exhibit two major differences: first, male sterility occurs not only at high temperature, but also at a low temperature; second, sterility thresholds are not evolutionarily constrained, but highly variable. Altogether, significant and sometimes major genetic variations have been observed between species, between geographic races of the same species, and even between strains kept in the laboratory under different thermal regimes. In each case, it is easily argued that the observed variations correspond to adaptations to climatic conditions, and that male sterility is a significant component of fitness and a target of natural selection.

Historical biogeography of Drosophila simulans based on Y-chromosomal sequences

Y-chromosomal sequences have been used for phylogeographic studies in humans and other mammals, but so far have been ignored as a source of historical information in Drosophila and other insects with X/Y sex determination. Here, we present the first phylogeographic study of Drosophila simulans based on the Y chromosome. Geographic distribution of Y-chromosomal haplotypes suggests a high degree of population subdivision within Africa, as well as between the African and cosmopolitan groups of populations. Consistent with earlier studies based on autosomal and X-linked loci, our results suggest that D. simulans originated in Madagascar or East Africa, and that the South and West African populations of this species are derived.

Multilocus analysis of introgression between two sympatric sister species of Drosophila: Drosophila yakuba and D. santomea

Drosophila yakuba is widely distributed in sub-Saharan Africa, while D. santomea is endemic to the volcanic island of São Tomé in the Atlantic Ocean, 280 km west of Gabon. On São Tomé, D. yakuba is found mainly in open lowland forests, and D. santomea is restricted to the wet misty forests at higher elevations. At intermediate elevations, the species form a hybrid zone where hybrids occur at a frequency of approximately 1%. To determine the extent of gene flow between these species we studied polymorphism and divergence patterns in 29 regions distributed throughout the genome, including mtDNA and three genes on the Y chromosome. This multilocus approach, together with the comparison to the two allopatric species D. mauritiana and D. sechellia, allowed us to distinguish between forces that should affect all genes and forces that should act on some genes (e.g., introgression). Our results show that D. yakuba mtDNA has replaced that of D. santomea and that there is also significant introgression for two nuclear genes, yellow and salr. The majority of genes, however, has remained distinct. These two species therefore do not form a "hybrid swarm" in which much of the genome shows substantial introgression while disruptive selection maintains distinctness for only a few traits (e.g., pigmentation and male genitalia).

Chicken W: a genetically uniform chromosome in a highly variable genome

The Y chromosome of organisms with male heterogamety is expected to show reduced levels of genetic diversity, because the effective population size is one-fourth that of autosomes. However, studies in mammals, flies, and plants show that Y chromosome diversity is lower than expected even when differences in effective population size are taken into account. This may be explained by skewed reproductive success among males, leading to low male effective population size, or by a strong role of selection in shaping levels of nucleotide diversity in nonrecombining chromosomes. We tested these hypotheses in a system with female heterogamety by estimating nucleotide diversity in the female-specific W chromosome of the domestic chicken by resequencing of 7,643 base pairs in 47 birds from 10 highly divergent breeds. The screening revealed only one single segregating site, which is in sharp contrast to our previous observation, using a similar panel of birds of, on average, one segregating site every 39 base pairs in autosomal sequence. When taking sex-specific mutation rates and differences in effective population size into account, the observed degree of W chromosome polymorphism is 28-fold lower than expected for the frequency of segregating sites and 13-fold lower than expected for estimates of nucleotide diversity (autosomes, 6.5 x 10(-3); W, 7.0 x 10(-5)). We note that selection is the only factor that can explain the reduced diversity in the sex-limited chromosome irrespective of mode of reproduction or whether there is male or female heterogamety. Reduced variability in female-specific W chromosomes is not easily explained by sexual selection.

Substitution rates in a new Silene latifolia sex-linked gene, SlssX/Y

Dioecious white campion Silene latifolia has sex chromosomal sex determination, with homogametic (XX) females and heterogametic (XY) males. This species has become popular in studies of sex chromosome evolution. However, the lack of genes isolated from the X and Y chromosomes of this species is a major obstacle for such studies. Here, I report the isolation of a new sex-linked gene, Slss, with strong homology to spermidine synthase genes of other species. The new gene has homologous intact copies on the X and Y chromosomes (SlssX and SlssY, respectively). Synonymous divergence between the SlssX and SlssY genes is 4.7%, and nonsynonymous divergence is 1.4%. Isolation of a homologous gene from nondioecious S. vulgaris provided a root to the gene tree and allowed the estimation of the silent and replacement substitution rates along the SlssX and SlssY lineages. Interestingly, the Y-linked gene has higher synonymous and nonsynonymous substitution rates. The elevated synonymous rate in the SlssY gene, compared with SlssX, confirms our previous suggestion that the S. latifolia Y chromosome has a higher mutation rate, compared with the X chromosome. When differences in silent substitution rate are taken into account, the Y-linked gene still demonstrates significantly faster accumulation of nonsynonymous substitutions, which is consistent with the theoretical prediction of relaxed purifying selection in Y-linked genes, leading to the accumulation of nonsynonymous substitutions and genetic degeneration of the Y-linked genes.

Linking phylogenetics with population genetics to reconstruct the geographic origin of a species

Reconstructing ancestral geographic origins is critical for understanding the long-term evolution of a species. Bayesian methods have been proposed to test biogeographic hypotheses while accommodating uncertainty in phylogenetic reconstruction. However, the problem that certain taxa may have a disproportionate influence on conclusions has not been addressed. Here, we infer the geographic origin of Drosophila simulans using 2,014 bp of the period locus from 63 lines collected from 18 countries. We also analyze two previously published datasets, alcohol dehydrogenase related and NADH:ubiquinone reductase 75 kDa subunit precursor. Phylogenetic inferences of all three loci support Madagascar as the geographic origin of D. simulans. Our phylogenetic conclusions are robust to taxon resampling and to the potentially confounding effects of recombination. To test our phylogenetically derived hypothesis we develop a randomization test of the population genetics prediction that sequences from the geographic origin should contain more genetic polymorphism than those from derived populations. We find that the Madagascar population has elevated genetic polymorphism relative to non-Madagascar sequences. These data are corroborated by mitochondrial DNA sequence data.

Recent radiation of endemic Caribbean Drosophila of the dunni subgroup inferred from multilocus DNA sequence variation

Studies of island endemism provide a unique opportunity to elucidate fundamental mechanisms of speciation. Here we examine intra- and interspecific DNA sequence variation at four unlinked genetic loci among populations of the Drosophila dunni subgroup to provide a detailed genealogical portrait of the process of speciation among these island endemic species. Our data indicate two major rounds of diversification that have shaped the D. dunni subgroup. The first occurred 1.6-2.6 million years ago and separated three major lineages, one in Puerto Rico and the Virgin Islands, a second in the northern Lesser Antilles and Barbados, and a third in St. Vincent and Grenada. A second round of diversification occurred in the last 96,000 years in the northern Lesser Antilles and Barbados. The four distinct species that resulted from this recent round of diversification maintain relatively high amounts of genetic variation, similar to that of a closely related mainland species, and share extensive ancestral polymorphism. These data suggest a minimal role for population bottlenecks linked to founder events in the history of the D. dunni subgroup. Further, the recent divergence of these island populations highlights the extremely rapid development of reproductive isolation and distinct patterns of abdominal pigmentation that has occurred in these species.

Protein Evolution and Codon Usage Bias on the Neo-Sex Chromosomes of Drosophila miranda

The neo-sex chromosomes of Drosophila miranda constitute an ideal system to study the effects of recombination on patterns of genome evolution. Due to a fusion of an autosome with the Y chromosome, one homolog is transmitted clonally. Here, I compare patterns of molecular evolution of 18 protein-coding genes located on the recombining neo-X and their homologs on the nonrecombining neo-Y chromosome. The rate of protein evolution has significantly increased on the neo-Y lineage since its formation. Amino acid substitutions are accumulating uniformly among neo-Y-linked genes, as expected if all loci on the neo-Y chromosome suffer from a reduced effectiveness of natural selection. In contrast, there is significant heterogeneity in the rate of protein evolution among neo-X-linked genes, with most loci being under strong purifying selection and two genes showing evidence for adaptive evolution. This observation agrees with theory predicting that linkage limits adaptive protein evolution. Both the neo-X and the neo-Y chromosome show an excess of unpreferred codon substitutions over preferred ones and no difference in this pattern was observed between the chromosomes. This suggests that there has been little or no selection maintaining codon bias in the D. miranda lineage. A change in mutational bias toward AT substitutions also contributes to the decline in codon bias. The contrast in patterns of molecular evolution between amino acid mutations and synonymous mutations on the neo-sex-linked genes can be understood in terms of chromosome-specific differences in effective population size and the distribution of selective effects of mutations.

Effective Population Size and Population Subdivision in Demographically Structured Populations

A fast-timescale approximation is applied to the coalescent process in a single population, which is demographically structured by sex and/or age. This provides a general expression for the probability that a pair of alleles sampled from the population coalesce in the previous time interval. The effective population size is defined as the reciprocal of twice the product of generation time and the coalescence probability. Biologically explicit formulas for effective population size with discrete generations and separate sexes are derived for a variety of different modes of inheritance. The method is also applied to a nuclear gene in a population of partially self-fertilizing hermaphrodites. The effects of population subdivision on a demographically structured population are analyzed, using a matrix of net rates of movement of genes between different local populations. This involves weighting the migration probabilities of individuals of a given age/sex class by the contribution of this class to the leading left eigenvector of the matrix describing the movements of genes between age/sex classes. The effects of sex-specific migration and nonrandom distributions of offspring number on levels of genetic variability and among-population differentiation are described for different modes of inheritance in an island model. Data on DNA sequence variability in human and plant populations are discussed in the light of the results.

Adaptive hitchhiking effects on genome variability

The continuing deluge of nucleotide polymorphism data is providing insights into the role of adaptation in shaping genome-wide patterns of variability and molecular evolution. Population genetic models in which linkage and selection interact (i.e. hitchhiking) predict that selection can leave 'footprints' in closely linked genomic regions. New analytical approaches show promise for distinguishing the signature of adaptation from that of several non-adaptive alternatives. Accounting for the effects of population structure and history poses a challenge for future investigations.

Identification of five new genes on the Y chromosome of Drosophila melanogaster

The heterochromatic state of the Drosophila Y chromosome has made the cloning and identification of Y-linked genes a challenging process. Here, we report application of a procedure to identify Y-linked gene fragments from the unmapped residue of the whole genome sequencing effort. Previously identified Y-linked genes appear in sequenced scaffolds as individual exons, apparently because many introns have become heterochromatic, growing to enormous size and becoming virtually unclonable. A TBLASTN search using all known proteins as query sequences, tested against a blastable database of the unmapped fragments, produced a number of matches consistent with this scenario. Reverse transcription-PCR and genetic methods were used to confirm those that are expressed, Y-linked genes. The five genes reported here include three protein phosphatases (Pp1-Y1, Pp1-Y2, and PPr-Y), an occludin-related gene (ORY), and a coiled-coils gene (CCY). This brings the total to nine protein-coding genes identified on the Drosophila Y chromosome. ORY and CCY may correspond, respectively, to the fertility factors ks-1 and ks-2, whereas the three protein phosphatases represent novel genes. There remains a strong functional coherence to male function among the genes on the Drosophila Y chromosome.

Contrasting levels of nucleotide diversity on the avian Z and W sex chromosomes

Sex chromosomes may provide a context for studying the local effects of mutation rate on molecular evolution, since the two types of sex chromosomes are generally exposed to different mutational environments in male and female germ lines. Importantly, recent studies of some vertebrates have provided evidence for a higher mutation rate among males than among females. Thus, in birds, the Z chromosome, which spends two thirds of its time in the male germ line, is exposed to more mutations than the female-specific W chromosome. We show here that levels of nucleotide diversity are drastically higher on the avian Z chromosome than in paralogous sequences on the W chromosome. In fact, no intraspecific polymorphism whatsoever was seen in about 3.4 kb of CHD1W intron sequence from a total of >150 W chromosome copies of seven different bird species. In contrast, the amount of genetic variability in paralogous sequences on the Z chromosome was significant, with an average pairwise nucleotide diversity (d) of 0.0020 between CHD1Z introns and with 37 segregating sites in a total of 3.8 kb of Z sequence. The contrasting levels of genetic variability on the avian sex chromosomes are thus in a direction predicted from a male-biased mutation rate. However, although a low gene number, as well as some other factors, argues against background selection and/or selective sweeps shaping the genetic variability of the avian W chromosome, we cannot completely exclude selection as a contributor to the low levels of variation on the W chromosome.