Characterizing gene movements between chromosomes in Drosophila

Genes occasionally change their location in the genome through inter-chromosomal duplication and loss. These changes happen as mistakes during recombination or through retrotransposition. In Han and Hahn 2011,(1) we surveyed the genomes of ten Drosophila species, to identify and characterize the gene transposition events in the history of these species. In the paper, we showed that the rate of gene transposition in Drosophila is higher than previously appreciated. To understand the process of gene transposition, we examined the sequences, locations, and functions of the transposed genes. Based on the elevated rate of sequence evolution in transposed genes and the frequent movements near the centromeres and telomeres, we could not reject the hypothesis that these are mutations fixed through relaxed selection. But, by examining the functions of transposed genes more carefully, we found that genes with male-specific functions and genes with female-specific functions move in opposite directions involving the X chromosome. We also found an over-representation of chromosome related functions among the transposed genes. These observations suggest the possibility of particular selection pressures contributing to gene transpositions in Drosophila.     

THE CONTRIBUTION OF GENE MOVEMENT TO THE “TWO RULES OF SPECIATION”

The two “rules of speciation”—the Large X‐effect and Haldane's rule—hold throughout the animal kingdom, but the underlying genetic mechanisms that cause them are still unclear. Two predominant explanations—the “dominance theory” and faster male evolution—both have some empirical support, suggesting that the genetic basis of these rules is likely multifarious. We revisit one historical explanation for these rules, based on dysfunctional genetic interactions involving genes recently moved between chromosomes. We suggest that gene movement specifically off or onto the X chromosome is another mechanism that could contribute to the two rules, especially as X chromosome movements can be subject to unique sex‐specific and sex chromosome specific consequences in hybrids. Our hypothesis is supported by patterns emerging from comparative genomic data, including a strong bias in interchromosomal gene movements involving the X and an overrepresentation of male reproductive functions among chromosomally relocated genes. In addition, our model indicates that the contribution of gene movement to the two rules in any specific group will depend upon key developmental and reproductive parameters that are taxon specific. We provide several testable predictions that can be used to assess the importance of gene movement as a contributor to these rules in the future.     

Sex chromosome-to-autosome transposition events counter Y-chromosome gene loss in mammals

BackgroundAlthough the mammalian X and Y chromosomes evolved from a single pair of autosomes, they are highly differentiated: the Y chromosome is dramatically smaller than the X and has lost most of its genes. The surviving genes are a specialized set with extraordinary evolutionary longevity. Most mammalian lineages have experienced delayed, or relatively recent, loss of at least one conserved Y-linked gene. An extreme example of this phenomenon is in the Japanese spiny rat, where the Y chromosome has disappeared altogether. In this species, many Y-linked genes were rescued by transposition to new genomic locations, but until our work presented here, this has been considered an isolated case.ResultsWe describe eight cases of genes that have relocated to autosomes in mammalian lineages where the corresponding Y-linked gene has been lost. These gene transpositions originated from either the X or Y chromosomes, and are observed in diverse mammalian lineages: occurring at least once in marsupials, apes, and cattle, and at least twice in rodents and marmoset. For two genes - EIF1AX/Y and RPS4X/Y - transposition to autosomes occurred independently in three distinct lineages.ConclusionsRescue of Y-linked gene loss through transposition to autosomes has previously been reported for a single isolated rodent species. However, our findings indicate that this compensatory mechanism is widespread among mammalian species. Thus, Y-linked gene loss emerges as an additional driver of gene transposition from the sex chromosomes, a phenomenon thought to be driven primarily by meiotic sex chromosome inactivation.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0667-4) contains supplementary material, which is available to authorized users.     

Retrogenes Moved Out of the Z Chromosome in the Silkworm

Previous studies on organisms with well-differentiated X and Y chromosomes, such as Drosophila and mammals, consistently detected an excess of genes moving out of the X chromosome and gaining testis-biased expression. Several selective evolutionary mechanisms were shown to be associated with this nonrandom gene traffic, which contributed to the evolution of the X chromosome and autosomes. If selection drives gene traffic, such traffic should also exist in species with Z and W chromosomes, where the females are the heterogametic sex. However, no previous studies on gene traffic in species with female heterogamety have found any nonrandom chromosomal gene movement. Here, we report an excess of retrogenes moving out of the Z chromosome in an organism with the ZW sex determination system, Bombyx mori. In addition, we showed that those "out of Z" retrogenes tended to have ovary-biased expression, which is consistent with the pattern of non-retrogene traffic recently reported in birds and symmetrical to the retrogene movement in mammals and fruit flies out of the X chromosome evolving testis functions. These properties of gene traffic in the ZW system suggest a general role for the heterogamety of sex chromosomes in determining the chromosomal locations and the evolution of sex-biased genes.     

Gene content evolution on the X chromosome

Compared with autosomes, the X chromosome shows different patterns of evolution as a result of its hemizygosity in males. Additionally, inactivation of the X during spermatogenesis can make the X chromosome an unfavorable location for male-specific genes. These factors can help to explain why in many species gene content of the X chromosome differs from that of autosomes. Indeed, the X chromosome in mouse is enriched for male-specific genes while they are depleted on the X in Drosophila but show neither of these trends in mosquito. Here, we will discuss recent findings on the ancestral and neo-X chromosomes in Drosophila that support sexual antagonism as a force shaping gene content evolution of sex chromosomes and suggest that selection could be driving male-biased genes off the X.     

On the genomic location of the exuperantia1 gene in Drosophila miranda: the limits of in situ hybridization experiments

In situ hybridization to Drosophila polytene chromosomes is a powerful tool for determining the chromosomal location of genes. Using in situ hybridization experiments, Yi and Charlesworth recently reported the transposition of the exuperantia1 gene (exu1) from a neo-sex chromosome to the ancestral X chromosome of Drosophila miranda, close to exuperantia2 (exu2). By characterizing sequences flanking exu1, however, we found the position of exu1 to be conserved on the neo-sex chromosome. Further, the exu2 gene was found to be tandemly duplicated on the X chromosome of D. miranda. The misleading hybridization signal of exu1 may be caused by multiple copies of exu2, which interfere with the hybridization of the exu1 probe to its genomic position on the neo-X chromosome. This suggests that flanking DNA should be used to confirm the positions of members of gene families.     

Sex chromosomes evolved from independent ancestral linkage groups in winged insects

The evolution of a pair of chromosomes that differ in appearance between males and females (heteromorphic sex chromosomes) has occurred repeatedly across plants and animals. Recent work has shown that the male heterogametic (XY) and female heterogametic (ZW) sex chromosomes evolved independently from different pairs of homomorphic autosomes in the common ancestor of birds and mammals but also that X and Z chromosomes share many convergent molecular features. However, little is known about how often heteromorphic sex chromosomes have either evolved convergently from different autosomes or in parallel from the same pair of autosomes and how universal patterns of molecular evolution on sex chromosomes really are. Among winged insects with sequenced genomes, there are male heterogametic species in both the Diptera (e.g., Drosophila melanogaster) and the Coleoptera (Tribolium castaneum), female heterogametic species in the Lepidoptera (Bombyx mori), and haplodiploid species in the Hymenoptera (e.g., Nasonia vitripennis). By determining orthologous relationships among genes on the X and Z chromosomes of insects with sequenced genomes, we are able to show that these chromosomes are not homologous to one another but are homologous to autosomes in each of the other species. These results strongly imply that heteromorphic sex chromosomes have evolved independently from different pairs of ancestral chromosomes in each of the insect orders studied. We also find that the convergently evolved X chromosomes of Diptera and Coleoptera share genomic features with each other and with vertebrate X chromosomes, including excess gene movement from the X to the autosomes. However, other patterns of molecular evolution--such as increased codon bias, decreased gene density, and the paucity of male-biased genes on the X--differ among the insect X and Z chromosomes. Our results provide evidence for both differences and nearly universal similarities in patterns of evolution among independently derived sex chromosomes.     

Maintenance of transposable element copy number in natural populations of Drosophila melanogaster and D. simulans

To investigate the main forces controlling the containment of transposable elements (TE) in natural populations, we analyzed the copia, mdg1, and 412 elements in various populations of Drosophila melanogaster and D. simulans. A lower proportion of insertion sites on the X chromosome in comparison with the autosomes suggests that selection against the detrimental effects of TE insertions is the major force containing TE copies in populations of Drosophila. This selection effect hypothesis is strengthened by the absence of the negative correlation between recombination rate and TE copy number along the chromosomes, which was expected under the alternative ectopic exchange model (selection against the deleterious rearrangements promoted by recombination between TE insertions). A cline in 412 copy number in relation to latitude was observed among the natural populations of D. simulans, with very high numbers existing in some local populations (around 60 copies in a sample from Canberra, Australia). An apparent absence of selection effects in this Canberra sample and a value of transposition rate equal to 1–2 × 10^-3 whatever the population and its copy number agree with the idea of recent but temporarily drastic TE movements in local populations. The high values of transposition rate in D. simulans clearly disfavor the hypothesis that the low amount of transposable elements in this species could result from a low transposition rate. This revised version was published online in August 2006 with corrections to the Cover Date.     

Codon Bias and Noncoding GC Content Correlate Negatively with Recombination Rate on the Drosophila X Chromosome

The patterns and processes of molecular evolution may differ between the X chromosome and the autosomes in Drosophila melanogaster. This may in part be due to differences in the effective population size between the two chromosome sets and in part to the hemizygosity of the X chromosome in Drosophila males. These and other factors may lead to differences both in the gene complements of the X and the autosomes and in the properties of the genes residing on those chromosomes. Here we show that codon bias and recombination rate are correlated strongly and negatively on the X chromosome, and that this correlation cannot be explained by indirect relationships with other known determinants of codon bias. This is in dramatic contrast to the weak positive correlation found on the autosomes. We explored possible explanations for these patterns, which required a comprehensive analysis of the relationships among multiple genetic properties such as protein length and expression level. This analysis highlights conserved features of coding sequence evolution on the X and the autosomes and illuminates interesting differences between these two chromosome sets.[Reviewing editor: Dr. Richard Kliman]     

Structural features of transposed human VK genes and implications for the mechanism of their transpositions.

The genes encoding the variable, joining and constant regions of human immunoglobulin light chains have been localized to the short arm of chromosome 2. However, several VK genes lie outside of the locus: a single copy cluster of five VK genes is located on chromosome 22; an isolated but amplified VkI gene is found on chromosome 1; and several isolated VkI genes are on as-yet-unidentified chromosomes other than chromosome 2. Vk genes not contained within the kappa locus are termed orphons. We have attempted to gain insight into the mechanism of transposition of both the chromosome 22 cluster and the several amplified VkI genes by searching in the kappa locus for a parent copy of the former, and by analyzing the junctions between transposed VKI-containing segments and adjacent non-amplified regions. The chromosome 22 orphon cluster must have been non-duplicatively transposed. Sequence features at the junctions of this and other orphon regions are direct and inverted repeats, and, in one case, an Alu repeat. These unusual features may have predisposed the orphon regions to transposition by serving as target sites for enzymes involved in recombination.     

Duplicative and conservative transpositions of larval serum protein 1 genes in the genus Drosophila

Interspecific comparative molecular analyses of transposed genes and their flanking regions can help to elucidate the time, direction, and mechanism of gene transposition. In the Drosophila melanogaster genome, three Larval serum protein 1 (Lsp1) genes (alpha, beta and gamma) are present and each of them is located on a different chromosome, suggesting multiple transposition events. We have characterized the molecular organization of Lsp1 genes in D. buzzatii, a species of the Drosophila subgenus and in D. pseudoobscura, a species of the Sophophora subgenus. Our results show that only two Lsp1 genes (beta and gamma) exist in these two species. The same chromosomal localization and genomic organization, different from that of D. melanogaster, is found in both species for the Lsp1beta and Lsp1gamma genes. Overall, at least two duplicative and two conservative transpositions are necessary to explain the present chromosomal distribution of Lsp1 genes in the three Drosophila species. Clear evidence for implication of snRNA genes in the transposition of Lsp1beta in Drosophila has been found. We suggest that an ectopic exchange between highly similar snRNA sequences was responsible for the transposition of this gene. We have also identified the putative cis-acting regulatory regions of these genes, which seemingly transposed along with the coding sequences.     

Population Structure of Morphological Traits in Clarkia Dudleyana I. Comparison of F(st) between Allozymes and Morphological Traits

The X-linked white gene when transposed to autosomes retains only partial dosage compensation. One copy of the gene in males expresses more than one copy but less than two copies in females. When inserted in ectopic X chromosome sites, the mini-white gene of the CaspeR vector can be fully dosage compensated and can even achieve hyperdosage compensation, meaning that one copy in males gives more expression than two copies in females. As sequences are removed gradually from the 5' end of the gene, we observe a progressive transition from hyperdosage compensation to full dosage compensation to partial dosage compensation. When the deletion reaches -17, the gene can no longer dosage compensate fully even on the X chromosome. A deletion reaching +173, 4 bp preceeding the AUG initiation codon, further reduces dosage compensation both on the X chromosome and on autosomes. This truncated gene can still partially dosage compensate on autosomes, indicating the presence of dosage compensation determinants in the protein coding region. We conclude that full dosage compensation requires an X chromosome environment and that the white gene contains multiple dosage-compensation determinants, some near the promoter and some in the coding region.     

Plant sex chromosomes: molecular structure and function

Recent molecular and genomic studies carried out in a number of model dioecious plant species, including Asparagus officinalis, Carica papaya, Silene latifolia, Rumex acetosa and Marchantia polymorpha, have shed light on the molecular structure of both homomorphic and heteromorphic sex chromosomes, and also on the gene functions they have maintained since their evolution from a pair of autosomes. The molecular structure of sex chromosomes in species from different plant families represents the evolutionary pathway followed by sex chromosomes during their evolution. The degree of Y chromosome degeneration that accompanies the suppression of recombination between the Xs and Ys differs among species. The primitive Ys of A. officinalis and C. papaya have only diverged from their homomorphic Xs in a short male-specific and non-recombining region (MSY), while the heteromorphic Ys of S. latifolia, R. acetosa and M. polymorpha have diverged from their respective Xs. As in the Y chromosomes of mammals and Drosophila, the accumulation of repetitive DNA, including both transposable elements and satellite DNA, has played an important role in the divergence and size enlargement of plant Ys, and consequently in reducing gene density. Nevertheless, the degeneration process in plants does not appear to have reached the Y-linked genes. Although a low gene density has been found in the sequenced Y chromosome of M. polymorpha, most of its genes are essential and are expressed in the vegetative and reproductive organs in both male and females. Similarly, most of the Y-linked genes that have been isolated and characterized up to now in S. latifolia are housekeeping genes that have X-linked homologues, and are therefore expressed in both males and females. Only one of them seems to be degenerate with respect to its homologous region in the X. Sequence analysis of larger regions in the homomorphic X and Y chromosomes of papaya and asparagus, and also in the heteromorphic sex chromosomes of S. latifolia and R. acetosa, will reveal the degenerative changes that the Y-linked gene functions have experienced during sex chromosome evolution.     

Evolution on the X chromosome: unusual patterns and processes

Although the X chromosome is usually similar to the autosomes in size and cytogenetic appearance, theoretical models predict that its hemizygosity in males may cause unusual patterns of evolution. The sequencing of several genomes has indeed revealed differences between the X chromosome and the autosomes in the rates of gene divergence, patterns of gene expression and rates of gene movement between chromosomes. A better understanding of these patterns should provide valuable information on the evolution of genes located on the X chromosome. It could also suggest solutions to more general problems in molecular evolution, such as detecting selection and estimating mutational effects on fitness.     

Adaptive protein evolution of X-linked and autosomal genes in Drosophila: implications for faster-X hypotheses

Patterns of sex chromosome and autosome evolution can be used to elucidate the underlying genetic basis of adaptative change. Evolutionary theory predicts that X-linked genes will adapt more rapidly than autosomes if adaptation is limited by the availability of beneficial mutations and if such mutations are recessive. In Drosophila, rates of molecular divergence between species appear to be equivalent between autosomes and the X chromosome. However, molecular divergence contrasts are difficult to interpret because they reflect a composite of adaptive and nonadaptive substitutions between species. Predictions based on faster-X theory also assume that selection is equally effective on the X and autosomes; this might not be true because the effective population sizes of X-linked and autosomal genes systematically differ. Here, population genetic and divergence data from Drosophila melanogaster, Drosophila simulans, and Drosophila yakuba are used to estimate the proportion of adaptive amino acid substitutions occurring in the D. melanogaster lineage. After gene composition and effective population size differences between chromosomes are controlled, X-linked and autosomal genes are shown to have equivalent rates of adaptive divergence with approximately 30% of amino acid substitutions driven by positive selection. The results suggest that adaptation is either unconstrained by a lack of beneficial genetic variation or that beneficial mutations are not recessive and are thus highly visible to natural selection whether on sex chromosomes or on autosomes.     

Accelerated evolution of sex chromosomes in aphids, an x0 system

Sex chromosomes play a role in many important biological processes, including sex determination, genomic conflicts, imprinting, and speciation. In particular, they exhibit several unusual properties such as inheritance pattern, hemizygosity, and reduced recombination, which influence their response to evolutionary factors (e.g., drift, selection, and demography). Here, we examine the evolutionary forces driving X chromosome evolution in aphids, an XO system where females are homozygous (XX) and males are hemizygous (X0) at sex chromosomes. We show by simulations that the unusual mode of transmission of the X chromosome in aphids, coupled with cyclical parthenogenesis, results in similar effective population sizes and predicted levels of genetic diversity for X chromosomes and autosomes under neutral evolution. These results contrast with expectations from standard XX/XY or XX/X0 systems (where the effective population size of the X is three-fourths that of autosomes) and have deep consequences for aphid X chromosome evolution. We then localized 52 microsatellite markers on the X and 351 on autosomes. We genotyped 167 individuals with 356 of these loci and found similar levels of allelic richness on the X and on the autosomes, as predicted by our simulations. In contrast, we detected higher dN and dN/dS ratio for X-linked genes compared with autosomal genes, a pattern compatible with either positive or relaxed selection. Given that both types of chromosomes have similar effective population sizes and that the single copy of the X chromosome of male aphids exposes its recessive genes to selection, some degree of positive selection seems to best explain the higher rates of evolution of X-linked genes. Overall, this study highlights the particular relevance of aphids to study the evolutionary factors driving sex chromosomes and genome evolution.     

Contrasting the efficacy of selection on the X and autosomes in Drosophila

To investigate the relative efficacy of both positive and purifying natural selection on the X chromosome and the autosomes in Drosophila, we compared rates and patterns of molecular evolution between these chromosome sets using the newly available alignments of orthologous genes from 12 species. Parameters that may influence the relative X versus autosomal substitution rates include the relative effective population sizes, the male and female germline mutation rates, the distribution of allelic effects on fitness, and the degree of dominance of novel mutations. Our analysis reveals that codon usage bias is consistently greater for X-linked genes, suggesting that purifying selection consistently has greater efficacy on the X chromosome than on the autosomes across the Drosophila phylogeny. However, our results are less consistent with respect to the efficacy of positive selection, with only some lineages showing a higher substitution rate on the X chromosome. This suggests that either the distribution of selective effects of mutations or other relevant parameters are sufficiently variable across species to tip the balance in different ways in individual lineages. These data suggest that rates of substitution are not solely governed by adaptive evolution. This genome-wide analysis provides a clear picture that the efficacy of selection varies intragenomically and that this effect is markedly more consistent across the phylogeny in the case of purifying selection. Our results also suggest that simple models that predict systematic differences in rates of evolution between the X and the autosomes can only be made to be compatible with these Drosophila data if the relevant population genetic parameters that drive substitution rates differ among species and chromosomal contexts.     

The distribution of transposable elements within and between chromosomes in a population of Drosophila melanogaster. II. Inferences on the nature of selection against elements.

Data were collected on the distribution of nine families of transposable elements among a sample of autosomes isolated from a natural population of Drosophila melanogaster, by means of in situ hybridization of biotinylated probes to polytene chromosomes. There is no general tendency for elements to accumulate at the tips of chromosomes. Elements tend to be present in excess of random expectation in the euchromatin proximal to the centromeres of the major autosomes, and on chromosome four. There is considerable heterogeneity between different families in the extent of this excess. The overall abundance of element families is inversely related to the extent to which they accumulate proximally. The level of proximal accumulation for the major autosomes is similar to that on the fourth chromosome, but less than that for the X chromosome. There is an overall deficiency of elements in the mid-section of the X compared with the mid-sections of the major autosomes, with considerable heterogeneity between families. The magnitude of this deficiency is positively related to the extent to which elements accumulate proximally. No such deficiency is seen if the proximal regions of the X and autosomes are compared. There is a small and non-significant excess of elements in third chromosomes carrying inversions. There is some between-year heterogeneity in element abundance. The implications of these findings are discussed, and it is concluded that they generally support the hypothesis that transposable element abundance is regulated primarily by the deleterious fitness consequences of meiotic ectopic exchange between elements. If this is the case, such exchange must be very infrequent in the proximal euchromatin, and the elements detected in population surveys of this kind must be inserted into sites where they have negligible mutational effects on fitness.     

A new resource of locally transposed Dissociation elements for screening gene-knockout lines in silico on the Arabidopsis genome

We transposed Dissociation (Ds) elements from three start loci on chromosome 5 in Arabidopsis (Nossen ecotype) by using a local transposition system. We determined partial genomic sequences flanking the Ds elements and mapped the elements' insertion sites in 1,173 transposed lines by comparison with the published genomic sequence. Most of the lines contained a single copy of the Ds element. One-half of the lines contained Ds on chromosome 5; in particular, insertion "hot spots" near the three start loci were clearly observed. In the other lines, the Ds elements were transposed across chromosomes. We found other insertion hot spots at the tops of chromosomes 2 and 4, near nucleolus organizer regions 2 and 4, respectively. Another characteristic feature was that the Ds elements tended to transpose near the chromosome ends and rarely transposed near centromeres. The distribution patterns differed among the three start loci, even though they possessed the same Ds construct. More than one-half of the Ds elements were inserted irregularly into the genome; that is, they did not retain the perfect inverted repeat sequence of Ds nor leave perfect target site duplications. This precise analysis of distribution patterns will contribute to a comprehensive understanding of the transposing mechanism. From these Ds insertion sites, we have constructed a database for screening gene-knockout mutants in silico. In 583 of the 1,173 lines, the Ds elements were inserted into protein-coding genes, which suggests that these lines are gene-knockout mutants. The database and individual lines will be available freely for academic use from the RIKEN Bio-Resource Center (http://www.brc.riken.go.jp/Eng/index.html).