Accumulation of Spock and Worf,two novel non-LTR retrotransposons,on the neo-Y chromosome of Drosophila miranda

Transposable elements constitute a major fraction of eukaryotic genomes. Here, I characterize two novel non-LTR retrotransposons, cloned from the neo-Y chromosome of Drosophila miranda. Worf is 4.1 kb in size and shows homology to the T1-2 non-LTR transposon characterized in Anopheles. Spock is 4.9 kb in size and shows similarity to the Doc element of D. melanogaster. Southern blot analysis of both elements yielded stronger signals for male DNA. In situ hybridization to polytene chromosomes revealed that both elements are accumulating on the neo-Y chromosome of D. miranda. PCR analysis was conducted to investigate the frequency of spock and worf and of the previously identified transposons, TRIM and TRAM, at individual chromosomal sites among 12 strains of D. miranda. Contrary to the observation that element frequencies are usually kept low at individual sites in Drosophila, the four transposons investigated are fixed at their genomic locations on the neo-Y chromosome. These results support the hypothesis that transposons accumulate in nonrecombining regions and may be one cause of the heteromorphism of sex chromosomes.     

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.     

A Duplication Including the Y Allele of Lcp2 and the Trim Retrotransposon at the Lcp Locus on the Degenerating Neo-Y Chromosome of Drosophila Miranda: Molecular Structure and Mechanisms by Which It May Have Arisen

Evolutionary changes during the process of sex chromosome differentiation in Drosophila miranda are associated with massive DNA rearrangements. Comparing the DNA structure of the larval cuticle protein (Lcp) region from the X2 and neo-Y chromosome pair, we observed insertions, deletions and a large duplication at the neo-Y chromosomal locus. The duplication encompasses a complete copy of the neo-Y allele of Lcp2, and the ISY3 and the ISY4 insertion sequences. The latter was identified as a retrotransposon, termed TRIM. ISY3 shows DNA sequence similarity to P element homologs identified in the Drosophila obscura species group. We were interested in mechanistic aspects generating the duplication. We cannot exclude unequivocally that unequal sister-chromatid exchange could give rise to the observed duplication; however, recombination is a rare event in Drosophila males. Location and sequence of the retrotransposon TRIM served as molecular markers allowing us to reconstruct two intrachromosomal transposition events that could lead to the observed duplication.     

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.     

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.     

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.     

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.     

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.     

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.     

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.     

Dosage compensation and dietary glucose repression of larval amylase activity inDrosophila miranda

The functional locus for alpha-amylase (Amy) in Drosophila miranda is in the evolutionarily new X2 chromosome. X2 evolved from an autosome in response to an ancestral autosome-Y translocation that gave rise to the "neo-Y" chromosome of this species. Y-linked Amy, if still present in the ancestrally translocated element, is unexpressed. Dosage compensation for amylase activity was examined in larvae of the S 204 strain. Since dietary glucose is known to repress Amy expression in Drosophila melanogaster, dosage compensation of amylase activity in male larvae of D. miranda was tested by rearing larvae of both sexes on yeast diets with or without a glucose supplement. The WT 10 strain of Drosophila persimilis, a sibling species in which Amy is autosomally linked, was used as a reference for tests of amylase activity differences between the sexes. On the diet with glucose, Amy expression was repressed in both WT 10 and S 204 larvae and male larvae of S 204 displayed dosage compensation for amylase activity. On the nonrepressing diet consisting of yeast alone, S 204 continued to display dosage compensation.     

New transposable elements identified as insertions in rice transposon Tnr1

Tnr1 (235 bp long) is a transposable element in rice. Polymerase chain reactions (PCRs) done with a primer(s) that hybridizes to terminal inverted repeat sequences (TIRs) of Tnr1 detected new Tnr1 members with one or two insertions in rice genomes. Six identified insertion sequences (Tnr4, Tnr5, Tnr11, Tnr12, Tnr13 and RIRE9) did not have extensive homology to known transposable elements, rather they had structural features characteristic of transposable elements. Tnr4 (1767 bp long) had imperfect 64-bp TIRs and appeared to generate duplication of a 9-bp sequence at the target site. However, the TIR sequences were not homologous to those of known transposable elements, indicative that Tnr4 is a new transposable element. Tnr5 (209 bp long) had imperfect 46-bp TIRs and appeared to generate duplication of sequence TTA like that of some elements of the Tourist family. Tnr11 (811 bp long) had 73-bp TIRs with significant homology to those of Tnr1 and Stowaway and appeared to generate duplication of sequence TA, indicative that Tnr11 is a transposable element of the Tnr1/Stowaway family. Tnr12 (2426 bp long) carried perfect 9-bp TIRs, which began with 5'-CACTA- -3' from both ends and appeared to generate duplication of a 3-bp target sequence, indicative that Tnr12 is a transposable element of the En/Spm family. Tnr13 (347 bp long) had 31-bp TIRs and appeared to generate duplication of an 8-bp target sequence. Two sequences, one the transposon-like element Crackle, had partial homology in the Tnr13 ends. All five insertions appear to be defective elements derived from autonomous ones encoding the transposase gene. All had characteristic tandem repeat sequences which may be recognized by transposase. The sixth insertion sequence, named RIRE9 (3852 bp long), which begins with 5'-TG- -3' and ends with 5'- -CA-3', appeared to generate duplication of a 5-bp target sequence. These and other structural features indicate that this insertion is a solo LTR (long terminal repeat) of a retrotransposon. The transposable elements described above could be identified as insertions into Tnr1, which do not deleteriously affect the growth of rice cells.     

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.     

Genetic organization of insertion element IS2 based on a revised nucleotide sequence

We identified a transposable element resident in the chromosome of Escherichia coli K-12 strain HB101. This is an approx. 4400-bp-long transposon flanked by two copies of insertion sequence (IS) 1 element in direct orientation. One of the IS1 elements was found to be integrated into an IS2 element between IS2 bp 139 and bp 140 with the large moiety of IS2 within the transposon. The sequence of this part of IS2 differs from the published sequence of galOP-308::IS2 at a number of positions. Restriction analysis of the published allele, however, indicated that both alleles may in fact be identical. Since six of the eight differences found alter open reading frames, the revised sequence presents a new outlook for the potential genetic organization of IS2.     

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.