Evolution of amino-acid sequences and codon usage on the Drosophila miranda neo-sex chromosomes

We have studied patterns of DNA sequence variation and evolution for 22 genes located on the neo-X and neo-Y chromosomes of Drosophila miranda. As found previously, nucleotide site diversity is greatly reduced on the neo-Y chromosome, with a severely distorted frequency spectrum. There is also an accelerated rate of amino-acid sequence evolution on the neo-Y chromosome. Comparisons of nonsynonymous and silent variation and divergence suggest that amino-acid sequences on the neo-X chromosome are subject to purifying selection, whereas this is much weaker on the neo-Y. The same applies to synonymous variants affecting codon usage. There is also an indication of a recent relaxation of selection on synonymous mutations for genes on other chromosomes. Genes that are weakly expressed on the neo-Y chromosome appear to have a faster rate of accumulation of both nonsynonymous and unpreferred synonymous mutations than genes with high levels of expression, although the rate of accumulation when both types of mutation are pooled is higher for the neo-Y chromosome than for the neo-X chromosome even for highly expressed genes.     

Contrasting patterns of molecular evolution of the genes on the new and old sex chromosomes of Drosophila miranda

In organisms with chromosomal sex determination, sex is determined by a set of dimorphic sex chromosomes that are thought to have evolved from a set of originally homologous chromosomes. The chromosome inherited only through the heterogametic sex (the Y chromosome in the case of male heterogamety) often exhibits loss of genetic activity for most of the genes carried on its homolog and is hence referred to as degenerate. The process by which the proto-Y chromosome loses its genetic activity has long been the subject of much speculation. We present a DNA sequence variation analysis of marker genes on the evolving sex chromosomes (neo-sex chromosomes) of Drosophila miranda. Due to its relatively recent origin, the neo-Y chromosome of this species is presumed to be still experiencing the forces responsible for the loss of its genetic activity. Indeed, several previous studies have confirmed the presence of some active loci on this chromosome. The genes on the neo-Y chromosome surveyed in the current study show generally lower levels of variation compared with their counterparts on the neo-X chromosome or an X-linked gene. This is in accord with a reduced effective population size of the neo-Y chromosome. Interestingly, the rate of replacement nucleotide substitutions for the neo-Y linked genes is significantly higher than that for the neo-X linked genes. This is not expected under a model where the faster evolution of the X chromosome is postulated to be the main force driving the degeneration of the Y chromosome.     

Reduced levels of microsatellite variability on the neo-Y chromosome of Drosophila miranda

BACKGROUND: In many species, sex is determined by a system involving X and Y chromosomes, the latter having lost much of their genetic activity. Sex chromosomes have evolved independently many times, and several different mechanisms responsible for the degeneration of the Y chromosome have been proposed. Here, we have taken advantage of the secondary sex chromosome pair in Drosophila miranda to test for the effects of evolutionary forces involved in the early stages of Y-chromosome degeneration. Because of a fusion of one of the autosomes to the Y chromosome, a neo-Y chromosome and a neo-X chromosome have been formed, resulting in the transmission of formerly autosomal genes in association with the sex chromosomes. RESULTS: We found a 25-fold lower level of variation at microsatellites located on the neo-Y chromosome compared with homologous loci on the neo-X chromosome, or with autosomal and X-linked microsatellites. Sequence analyses of the region flanking the microsatellites suggested that the neo-sex chromosomes originated about 1 million years ago. CONCLUSIONS: Variability of the neo-Y chromosome of D. miranda is substantially reduced below expectations at mutation-drift equilibrium. Such a reduction is predicted by theories of the degeneration of the Y chromosome. Another possibility is that there is little or no mutation at microsatellite loci on a non-recombining chromosome such as the neo-Y, but this seems inconsistent with other data.     

Sex Chromosome Turnover Contributes to Genomic Divergence between Incipient Stickleback Species

Sex chromosomes turn over rapidly in some taxonomic groups, where closely related species have different sex chromosomes. Although there are many examples of sex chromosome turnover, we know little about the functional roles of sex chromosome turnover in phenotypic diversification and genomic evolution. The sympatric pair of Japanese threespine stickleback (Gasterosteus aculeatus) provides an excellent system to address these questions: the Japan Sea species has a neo-sex chromosome system resulting from a fusion between an ancestral Y chromosome and an autosome, while the sympatric Pacific Ocean species has a simple XY sex chromosome system. Furthermore, previous quantitative trait locus (QTL) mapping demonstrated that the Japan Sea neo-X chromosome contributes to phenotypic divergence and reproductive isolation between these sympatric species. To investigate the genomic basis for the accumulation of genes important for speciation on the neo-X chromosome, we conducted whole genome sequencing of males and females of both the Japan Sea and the Pacific Ocean species. No substantial degeneration has yet occurred on the neo-Y chromosome, but the nucleotide sequence of the neo-X and the neo-Y has started to diverge, particularly at regions near the fusion. The neo-sex chromosomes also harbor an excess of genes with sex-biased expression. Furthermore, genes on the neo-X chromosome showed higher non-synonymous substitution rates than autosomal genes in the Japan Sea lineage. Genomic regions of higher sequence divergence between species, genes with divergent expression between species, and QTL for inter-species phenotypic differences were found not only at the regions near the fusion site, but also at other regions along the neo-X chromosome. Neo-sex chromosomes can therefore accumulate substitutions causing species differences even in the absence of substantial neo-Y degeneration.     

Genomic degradation of a young Y chromosome in Drosophila miranda

Background

Y chromosomes are derived from ordinary autosomes and degenerate because of a lack of recombination. Well-studied Y chromosomes only have few of their original genes left and contain little information about their evolutionary origin. Here, we take advantage of the recently formed neo-Y chromosome of Drosophila miranda to study the processes involved in Y degeneration on a genomic scale.

Results

We obtained sequence information from 14 homologous bacterial artificial chromosome (BAC) clones from the neo-X and neo-Y chromosome of D. miranda, encompassing over 2.5 Mb of neo-sex-linked DNA. A large fraction of neo-Y DNA is composed of repetitive and transposable-element-derived DNA (20% of total DNA) relative to their homologous neo-X linked regions (1%). The overlapping regions of the neo-sex linked BAC clones contain 118 gene pairs, half of which are pseudogenized on the neo-Y. Pseudogenes evolve significantly faster on the neo-Y than functional genes, and both functional and non-functional genes show higher rates of protein evolution on the neo-Y relative to their neo-X homologs. No heterogeneity in levels of degeneration was detected among the regions investigated. Functional genes on the neo-Y are under stronger evolutionary constraint on the neo-X, but genes were found to degenerate randomly on the neo-Y with regards to their function or sex-biased expression patterns.

Conclusion

Patterns of genome evolution in D. miranda demonstrate that degeneration of a recently formed Y chromosome can proceed very rapidly, by both an accumulation of repetitive DNA and degeneration of protein-coding genes. Our data support a random model of Y inactivation, with little heterogeneity in degeneration among genomic regions, or between functional classes of genes or genes with sex-biased expression patterns.     

Accelerated Adaptive Evolution on a Newly Formed X Chromosome

Sex chromosomes originated from ordinary autosomes, and their evolution is characterized by continuous gene loss from the ancestral Y chromosome. Here, we document a new feature of sex chromosome evolution: bursts of adaptive fixations on a newly formed X chromosome. Taking advantage of the recently formed neo-X chromosome of Drosophila miranda, we compare patterns of DNA sequence variation at genes located on the neo-X to genes on the ancestral X chromosome. This contrast allows us to draw inferences of selection on a newly formed X chromosome relative to background levels of adaptation in the genome while controlling for demographic effects. Chromosome-wide synonymous diversity on the neo-X is reduced 2-fold relative to the ancestral X, as expected under recent and recurrent directional selection. Several statistical tests employing various features of the data consistently identify 10%–15% of neo-X genes as targets of recent adaptive evolution but only 1%–3% of genes on the ancestral X. In addition, both the rate of adaptation and the fitness effects of adaptive substitutions are estimated to be roughly an order of magnitude higher for neo-X genes relative to genes on the ancestral X. Thus, newly formed X chromosomes are not passive players in the evolutionary process of sex chromosome differentiation, but respond adaptively to both their sex-biased transmission and to Y chromosome degeneration, possibly through demasculinization of their gene content and the evolution of dosage compensation.     

A selective sweep associated with a recent gene transposition in Drosophila miranda

In Drosophila miranda, a chromosome fusion between the Y chromosome and the autosome corresponding to Muller's element C has created a new sex chromosome system. The chromosome attached to the ancestral Y chromosome is transmitted paternally and hence is not exposed to crossing over. This chromosome, conventionally called the neo-Y, and the homologous neo-X chromosome display many properties of evolving sex chromosomes. We report here the transposition of the exuperantia1 (exu1) locus from a neo-sex chromosome to the ancestral X chromosome of D. miranda. Exu1 is known to have several critical developmental functions, including a male-specific role in spermatogenesis. The ancestral location of exu1 is conserved in the sibling species of D. miranda, as well as in a more distantly related species. The transposition of exu1 can be interpreted as an adaptive fixation, driven by a selective advantage conferred by its effect on dosage compensation. This explanation is supported by the pattern of within-species sequence variation at exu1 and the nearby exu2 locus. The implications of this phenomenon for genome evolution are discussed.     

The Epigenome of Evolving Drosophila Neo-Sex Chromosomes: Dosage Compensation and Heterochromatin Formation

Sex chromosomes originated from autosomes but have evolved a highly specialized chromatin structure. Drosophila Y chromosomes are composed entirely of silent heterochromatin, while male X chromosomes have highly accessible chromatin and are hypertranscribed as a result of dosage compensation. Here, we dissect the molecular mechanisms and functional pressures driving heterochromatin formation and dosage compensation of the recently formed neo-sex chromosomes of Drosophila miranda. We show that the onset of heterochromatin formation on the neo-Y is triggered by an accumulation of repetitive DNA. The neo-X has evolved partial dosage compensation and we find that diverse mutational paths have been utilized to establish several dozen novel binding consensus motifs for the dosage compensation complex on the neo-X, including simple point mutations at pre-binding sites, insertion and deletion mutations, microsatellite expansions, or tandem amplification of weak binding sites. Spreading of these silencing or activating chromatin modifications to adjacent regions results in massive mis-expression of neo-sex linked genes, and little correspondence between functionality of genes and their silencing on the neo-Y or dosage compensation on the neo-X. Intriguingly, the genomic regions being targeted by the dosage compensation complex on the neo-X and those becoming heterochromatic on the neo-Y show little overlap, possibly reflecting different propensities along the ancestral chromosome that formed the sex chromosome to adopt active or repressive chromatin configurations. Our findings have broad implications for current models of sex chromosome evolution, and demonstrate how mechanistic constraints can limit evolutionary adaptations. Our study also highlights how evolution can follow predictable genetic trajectories, by repeatedly acquiring the same 21-bp consensus motif for recruitment of the dosage compensation complex, yet utilizing a diverse array of random mutational changes to attain the same phenotypic outcome.     

Reduced sequence variability on the Neo-Y chromosome of Drosophila americana americana.

Sex chromosomes are generally morphologically and functionally distinct, but the evolutionary forces that cause this differentiation are poorly understood. Drosophila americana americana was used in this study to examine one aspect of sex chromosome evolution, the degeneration of nonrecombining Y chromosomes. The primary X chromosome of D. a. americana is fused with a chromosomal element that was ancestrally an autosome, causing this homologous chromosomal pair to segregate with the sex chromosomes. Sequence variation at the Alcohol Dehydrogenase (Adh) gene was used to determine the pattern of nucleotide variation on the neo-sex chromosomes in natural populations. Sequences of Adh were obtained for neo-X and neo-Y chromosomes of D. a. americana, and for Adh of D. a. texana, in which it is autosomal. No significant sequence differentiation is present between the neo-X and neo-Y chromosomes of D. a. americana or the autosomes of D. a. texana. There is a significantly lower level of sequence diversity on the neo-Y chromosome relative to the neo-X in D. a. americana. This reduction in variability on the neo-Y does not appear to have resulted from a selective sweep. Coalescent simulations of the evolutionary transition of an autosome into a Y chromosome indicate there may be a low level of recombination between the neo-X and neo-Y alleles of Adh and that the effective population size of this chromosome may have been reduced below the expected value of 25% of the autosomal effective size, possibly because of the effects of background selection or sexual selection.     

A survey of chromosomal and nucleotide sequence variation in Drosophila miranda

There have recently been several studies of the evolution of Y chromosome degeneration and dosage compensation using the neo-sex chromosomes of Drosophila miranda as a model system. To understand these evolutionary processes more fully, it is necessary to document the general pattern of genetic variation in this species. Here we report a survey of chromosomal variation, as well as polymorphism and divergence data, for 12 nuclear genes of D. miranda. These genes exhibit varying levels of DNA sequence polymorphism. Compared to its well-studied sibling species D. pseudoobscura, D. miranda has much less nucleotide sequence variation, and the effective population size of this species is inferred to be several-fold lower. Nevertheless, it harbors a few inversion polymorphisms, one of which involves the neo-X chromosome. There is no convincing evidence for a recent population expansion in D. miranda, in contrast to D. pseudoobscura. The pattern of population subdivision previously observed for the X-linked gene period is not seen for the other loci, suggesting that there is no general population subdivision in D. miranda. However, data on an additional region of period confirm population subdivision for this gene, suggesting that local selection is operating at or near period to promote differentiation between populations.     

Retroelements: tools for sex chromosome evolution

Many eukaryotic taxa inherit a heteromorphic sex chromosome pair. It is a generally accepted hypothesis that the sex chromosome pair is derived from a pair of homologous autosomes that has developed after the occurrence of a sex differentiator in an evolutionary process into two structurally and functionally different partners. In most of the analyzed systems the occurrence of the dominant sex differentiator is paralleled by the suppression of recombination within and close by that region. The recombinational isolation can spread in an evolutionary selection process from neighboring regions finally over the whole chromosome. Suppression of recombination strongly biases the distribution of retrotransposons in the genome. Our results and that from others indicate that the major force driving the evolution of Y chromosomes are retrotransposons, remodeling euchromatic chromosome structures into heterochromatic ones. In our model, intact or already eroded retrotransposons become trapped due to their inherent transposition mechanisms in non-recombining regions. The massive accumulation of retrotransposons interferes strongly with the activity of genes. We hypothesize that Y chromosome degeneration is a stepwise evolutionary process: (1) Massive accumulation of retrotransposons occurs in the non-recombining regions. (2) Heterochromatic nucleation centers are formed as a consequence of genomic defense against invasive parasitic elements; the established nucleation centers become epigenetically inherited. (3) Spreading of heterochromatin from the nucleation centers into flanking regions induces in an adaptive process gene silencing of neighbored genes that could either be still intact or in an already eroded condition, e.g., showing point mutations, deletions, insertions; the retroelements should be subjects to the same forces of deterioration as the genes themselves. (4) Constitutive silenced genes are not committed to the same genetic selection pressure as active genes and therefore more exposed to the decay process. (5) Gene dosage balance is reestablished by the parallel evolution of dosage compensation mechanisms. The evolving secondary sex chromosomes, neo-X and neo-Y, of Drosophila miranda are revealed to be a unique and potent model system to catch the evolutionary Y deterioration process in progress.     

Lack of Degeneration of Loci on the Neo-Y Chromosome of Drosophila Americana Americana

The extent of genetic degeneration of the neo-Y chromosome of Drosophila americana americana has been investigated. Three loci, coding for the enzymes enolase, phosphoglycerate kinase and alcohol dehydrogenase, have been localized to chromosome 4 of D. a. americana, which forms the neo-Y and neo-X chromosomes. Crosses between D. a. americana and D. virilis or D. montana showed that the loci coding for these enzymes carry active alleles on the neo-Y chromosome in all wild-derived strains of americana that were tested. Intercrosses between a genetically marked stock of virilis and strains of americana were carried out, creating F(3) males that were homozygous for sections of the neo-Y chromosome. The sex ratios in the F(3) generation of the intercrosses showed that no lethal alleles have accumulated on any of the neo-Y chromosomes tested. There was evidence for more minor reductions in fitness, but this seems to be mainly caused by deleterious alleles that are specific to each strain. A similar picture was provided by examination of the segregation ratios of two marker genes among the F(3) progeny. Overall, the data suggest that the neo-Y chromosome has undergone very little degeneration, certainly not to the extent of having lost the functions of vital genes. This is consistent with the recent origin of the neo-Y and neo-X chromosomes, and the slow rates at which the forces that cause Y chromosome degeneration are likely to work.     

The Amylase gene cluster on the evolving sex chromosomes of Drosophila miranda.

On the basis of chromosomal homology, the Amylase gene cluster in Drosophila miranda must be located on the secondary sex chromosome pair, neo-X (X2) and neo-Y, but is autosomally inherited in all other Drosophila species. Genetic evidence indicates no active amylase on the neo-Y chromosome and the X2-chromosomal locus already shows dosage compensation. Several lines of evidence strongly suggest that the Amy gene cluster has been lost already from the evolving neo-Y chromosome. This finding shows that a relatively new neo-Y chromosome can start to lose genes and hence gradually lose homology with the neo-X. The X2-chromosomal Amy1 is intact and Amy2 contains a complete coding sequence, but has a deletion in the 3''-flanking region. Amy3 is structurally eroded and hampered by missing regulatory motifs. Functional analysis of the X2-chromosomal Amy1 and Amy2 regions from D. miranda in transgenic D. melanogaster flies reveals ectopic AMY1 expression. AMY1 shows the same electrophoretic mobility as the single amylase band in D. miranda, while ectopic AMY2 expression is characterized by a different mobility. Therefore, only the Amy1 gene of the resident Amy cluster remains functional and hence Amy1 is the dosage compensated gene.     

Evidence for degeneration of the Y chromosome in the dioecious plant Silene latifolia

The human Y--probably because of its nonrecombining nature--has lost 97% of its genes since X and Y chromosomes started to diverge [1, 2]. There are clear signs of degeneration in the Drosophila miranda neoY chromosome (an autosome fused to the Y chromosome), with neoY genes showing faster protein evolution [3-6], accumulation of unpreferred codons [6], more insertions of transposable elements [5, 7], and lower levels of expression [8] than neoX genes. In the many other taxa with sex chromosomes, Y degeneration has hardly been studied. In plants, many genes are expressed in pollen [9], and strong pollen selection may oppose the degeneration of plant Y chromosomes [10]. Silene latifolia is a dioecious plant with young heteromorphic sex chromosomes [11, 12]. Here we test whether the S. latifolia Y chromosome is undergoing genetic degeneration by analyzing seven sex-linked genes. S. latifolia Y-linked genes tend to evolve faster at the protein level than their X-linked homologs, and they have lower expression levels. Several Y gene introns have increased in length, with evidence for transposable-element accumulation. We detect signs of degeneration in most of the Y-linked gene sequences analyzed, similar to those of animal Y-linked and neo-Y chromosome genes.     

Sequence differentiation associated with an inversion on the neo-X chromosome of Drosophila americana

Sex chromosomes originate from pairs of autosomes that acquire controlling genes in the sex-determining cascade. Universal mechanisms apparently influence the evolution of sex chromosomes, because this chromosomal pair is characteristically heteromorphic in a broad range of organisms. To examine the pattern of initial differentiation between sex chromosomes, sequence analyses were performed on a pair of newly formed sex chromosomes in Drosophila americana. This species has neo-sex chromosomes as a result of a centromeric fusion between the X chromosome and an autosome. Sequences were analyzed from the Alcohol dehydrogenase (Adh), big brain (bib), and timeless (tim) gene regions, which represent separate positions along this pair of neo-sex chromosomes. In the northwestern range of the species, the bib and Adh regions exhibit significant sequence differentiation for neo-X chromosomes relative to neo-Y chromosomes from the same geographic region and other chromosomal populations of D. americana. Furthermore, a nucleotide site defining a common haplotype in bib is shown to be associated with a paracentric inversion [In(4)ab] on the neo-X chromosome, and this inversion suppresses recombination between neo-X and neo-Y chromosomes. These observations are consistent with the inversion acting as a recombination modifier that suppresses exchange between these neo-sex chromosomes, as predicted by models of sex chromosome evolution.     

Biased distribution of repetitive elements: a landmark for neo-Y chromosome evolution in Drosophila miranda

It is generally assumed that the sex chromosomes developed from a pair of homologs. Over evolution, the proto-Y chromosome, with a very short differential segment, matured in its final stage into a heterochromatic and, for the most part, genetically eroded Y chromosome. The constraints on the evolution of the proto-Y chromosome have been speculated upon since the sex chromosomes were discovered. Several models have been suggested. Drosophila miranda has proved to be a unique and potent model system to study Y-chromosome evolution. We use selected test genes distributed along the neo-Y chromosome as entry gates to analyze the molecular mechanisms involved in the process of Y-chromosome evolution. Here, we report our findings on the Krüppel gene (Kr), which is located distally on the neo-sex chromosome pair.     

Evidence for genetic drift in endosymbionts (Buchnera): analyses of protein-coding genes.

Buchnera, the bacterial endosymbionts of aphids, undergo severe population bottlenecks during maternal transmission through their hosts. Previous studies suggest an increased effect of drift within these strictly asexual, small populations, resulting in an increased fixation of slightly deleterious mutations. This study further explores sequence evolution in Buchnera using three approaches. First, patterns of codon usage were compared across several homologous Escherichia coli and Buchnera loci, in order to test the prediction that selection for the use of optimal codons is less effective in small populations. A chi 2-based measure of codon bias was developed to adjust for the overall A + T richness of silent positions in the endosymbionts. In contrast to E. coli homologues, adaptive codon bias across Buchnera loci is markedly low, and patterns of codon usage lack a strong relationship with gene expression level. These data suggest that codon usage in Buchnera has been shaped largely by mutational pressure and drift rather than by selection for translational efficiency. One exception to the overall lack of bias is groEL, which is known to be constitutively overexpressed in Buchnera and other endosymbionts. Second, relative-rate tests show elevated rates of sequence evolution of numerous protein-coding loci across Buchnera, compared to E. coli. Finally, consistently higher ratios of nonsynonymous to synonymous substitutions in Buchnera loci relative to the enteric bacteria strongly suggest the accumulation of nonsynonymous substitutions in endosymbiont lineages. Combined, these results suggest a decreased effectiveness of purifying selection in purging endosymbiont populations of slightly deleterious mutations, particularly those affecting codon usage and amino acid identity.     

X-linked genes evolve higher codon bias in Drosophila and Caenorhabditis

Comparing patterns of molecular evolution between autosomes and sex chromosomes (such as X and W chromosomes) can provide insight into the forces underlying genome evolution. Here we investigate patterns of codon bias evolution on the X chromosome and autosomes in Drosophila and Caenorhabditis. We demonstrate that X-linked genes have significantly higher codon bias compared to autosomal genes in both Drosophila and Caenorhabditis. Furthermore, genes that become X-linked evolve higher codon bias gradually, over tens of millions of years. We provide several lines of evidence that this elevation in codon bias is due exclusively to their chromosomal location and not to any other property of X-linked genes. We present two possible explanations for these observations. One possibility is that natural selection is more efficient on the X chromosome due to effective haploidy of the X chromosomes in males and persistently low effective numbers of reproducing males compared to that of females. Alternatively, X-linked genes might experience stronger natural selection for higher codon bias as a result of maladaptive reduction of their dosage engendered by the loss of the Y-linked homologs.