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.     

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.     

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.     

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.     

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.     

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.     

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.     

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.     

Searching for guidance cues: Follow the Sidestep trail

Recombination restriction between evolving sex chromosomes leads to the degeneration of the chromosome that is present only in the heterogametic sex (the Y chromosome in XY species). The evolutionary forces driving Y chromosome degeneration, however, are still under debate and include positive and negative selection models. In a recent study, we showed that the rate of accumulation of loss-of-function mutations on the neo-Y chromosome of Drosophila miranda is compatible with the process of Muller's ratchet, the stochastic loss of the best mutational class of individuals from a small asexual population. Purifying selection at amino acid sites can accelerate the ratchet, and the speed of degeneration depends on the number of genes still present on the evolving Y chromosome. Our study shows that Y chromosome degeneration does not require the action of selective sweeps at linked sites, and can take place under realistic parameters of purifying selection only.     

The Enigma of Y Chromosome Degeneration: Tram, a Novel Retrotransposon Is Preferentially Located on the Neo-Y Chromosome of Drosophila Miranda

We have cloned a novel transposable element from the neo-Y chromosome of Drosophila miranda. The size of the element, designated as TRAM, is 3.452 bp, including on both sides long terminal direct repeats (LTRs) of 372 bp, respectively. The element is flanked by a 5-bp target site duplication, ATATG. The putative primer binding site (PBS) for minus-strand priming is complementary to 18 nucleotides of the 3'-end of tRNA(Trp). Data base screens for DNA sequence identities were negative, apart from the sequence motif of the PBS. The deduced amino acid sequence from the large ORF does not reveal identities described for other transposons. In situ hybridizations with TRAM subclones show a biased distribution in the genome, with a massive accumulation of TRAM in the neo-Y chromosome, while the former homologue, the X2 chromosome is devoid of TRAM sites. The enriched occurrence of the TRAM element at the evolving neo-Y chromosome of D.miranda adds compelling evidence in favor of the view that Y chromosome degeneration is driven by the accumulation of transposable elements.     

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.     

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 degeneration of Y chromosomes

Y chromosomes are genetically degenerate, having lost most of the active genes that were present in their ancestors. The causes of this degeneration have attracted much attention from evolutionary theorists. Four major theories are reviewed here: Muller's ratchet, background selection, the Hill Robertson effect with weak selection, and the 'hitchhiking' of deleterious alleles by favourable mutations. All of these involve a reduction in effective population size as a result of selective events occurring in a non-recombining genome, and the consequent weakening of the efficacy of selection. We review the consequences of these processes for patterns of molecular evolution and variation at loci on Y chromosomes, and discuss the results of empirical studies of these patterns for some evolving Y-chromosome and neo-Y-chromosome systems. These results suggest that the effective population sizes of evolving Y or neo-Y chromosomes are severely reduced, as expected if some or all of the hypothesized processes leading to degeneration are operative. It is, however, currently unclear which of the various processes is most important; some directions for future work to help to resolve this question are discussed.     

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.     

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.     

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.