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

Genome-wide analysis of retrogene polymorphisms in Drosophila melanogaster

Gene duplication via retrotransposition has been shown to be an important mechanism in evolution, affecting gene dosage and allowing for the acquisition of new gene functions. Although fixed retrotransposed genes have been found in a variety of species, very little effort has been made to identify retrogene polymorphisms. Here, we examine 37 Illumina-sequenced North American Drosophila melanogaster inbred lines and present the first ever data set and analysis of polymorphic retrogenes in Drosophila. We show that this type of polymorphism is quite common, with any two gametes in the North American population differing in the presence or absence of six retrogenes, accounting for ~13% of gene copy-number heterozygosity. These retrogenes were identified by a straightforward method that can be applied using any type of DNA sequencing data. We also use a variant of this method to conduct a genome-wide scan for intron presence/absence polymorphisms, and show that any two chromosomes in the population likely differ in the presence of multiple introns. We show that these polymorphisms are all in fact deletions rather than intron gain events present in the reference genome. Finally, by leveraging the known location of the parental genes that give rise to the retrogene polymorphisms, we provide direct evidence that natural selection is responsible for the excess of fixations of retrogenes moving off of the X chromosome in Drosophila. Further efforts to identify retrogene and intron presence/absence polymorphisms will undoubtedly improve our understanding of the evolution of gene copy number and gene structure.

Distinguishing among evolutionary models for the maintenance of gene duplicates

Determining the evolutionary forces responsible for the maintenance of gene duplicates is key to understanding the processes leading to evolutionary adaptation and novelty. In his highly prescient book, Susumu Ohno recognized that duplicate genes are fixed and maintained within a population with 3 distinct outcomes: neofunctionalization, subfunctionalization, and conservation of function. Subsequent researchers have proposed a multitude of population genetic models that lead to these outcomes, each differing largely in the role played by adaptive natural selection. In this paper, I present a nonmathematical review of these models, their predictions, and the evidence collected in support of each of them. Though the various outcomes of gene duplication are often strictly associated with the presence or absence of adaptive natural selection, I argue that determining the outcome of duplication is orthogonal to determining whether natural selection has acted. Despite an ever-growing field of research into the fate of gene duplicates, there is not yet clear evidence for the preponderance of one outcome over the others, much less evidence for the importance of adaptive or nonadaptive forces in maintaining these duplicates.

Many or most genes in Arabidopsis transposed after the origin of the order Brassicales

Previous to this work, typical genes were thought to move from one position to another infrequently. On the contrary, we now estimate that between one-fourth and three-fourths of the genes in Arabidopsis transposed in the Brassicales. We used the CoGe comparative genomics system to perform and visualize multiple orthologous chromosomal alignments. Using this tool, we found large differences between different categories of genes. Ten of the gene families examined, including genes in most transcription factor families, exhibited a median frequency of 5% transposed genes. In contrast, other gene families were composed largely of transposed genes: NB-LRR disease-resistance genes, genes encoding MADS-box and B3 transcription factors, and genes encoding F-box proteins. A unique method involving transposition-rich regions of genome allowed us to obtain an indirect estimate of the positional stability of the average gene. The observed differences between gene families raise important questions concerning the causes and consequences of gene transposition.

Two independent duplications forming the Cyp307a genes in Drosophila

The conserved relationship between orthologs of many cytochrome P450 genes involved in ecdysone synthesis is not reflected in the evolution of the Drosophila Cyp307a genes. In Drosophila melanogaster Cyp307a1 (spook) and Cyp307a2 (spookier) both play essential roles in ecdysone synthesis and may possess biochemically redundant functions. Using phylogenetic analyses we show that the Drosophila Cyp307a genes were formed from two independent duplication events depicting a complicated evolutionary scenario. An initial duplication, from a Cyp307a2 ancestral gene produced the Cyp307a1 gene that has been maintained only in the Sophophoran subgenus. A second duplication in the Drosophila subgenus formed an additional paralog, Cyp307a3. Microsynteny is conserved for Cyp307a2 throughout the Drosophila species, but is not conserved between Cyp307a1 and Cyp307a3. These are located in different genomic positions in the Sophophora and Drosophila subgenera, respectively. Cyp307a3 appears to encode a functional gene product and is expressed in a different spatial and temporal manner to Cyp307a1. This suggests some level of functional divergence between the Cyp307a paralogs in different Drosophila species.

A New Retroposed Gene in Drosophila Heterochromatin Detected by Microarray-Based Comparative Genomic Hybridization

A genomic pattern of new gene origination is often dependent on a genomic method that can efficiently identify a statistically adequate number of recently originated genes. The heterochromatic regions have often been viewed as genomic deserts with low coding potential and thus a low flux of new genes. However, increasing reports revealed unexpected roles of heterochromatic regions in the evolution of genes and genomes. We identified recently retroposed genes that originated in heterochromatic regions in Drosophila, by developing microarray-based comparative genomic hybridization (CGH) with multiple species. This new gene family, named Ifc-2h, originated in the common ancestor of the clade of D. simulans, D. mauritiana, and D. sechellia. The sequence features and phylogenetic distribution indicated that Ifc-2h resulted from the retroposition from its parental gene, Infertile crescent (Ifc), and integrated into heterochromatic region of common ancestor of the three sibling species 2 million years ago. Expression analysis revealed that Ifc-2h had developed a new expression pattern by recruiting a putative regulatory element from its target sequence. The distribution of indel variation in Ifc-2h of D. simulans and D. mauritiana revealed a significant sequence constraint, suggesting that the Ifc-2h gene may be functional. These analyses cast fresh insight into the evolution of heterochromatin and the origin of its coding regions.

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.

Rates and patterns of chromosomal evolution in Drosophila pseudoobscura and D. miranda

Comparisons of gene orders between species permit estimation of the rate of chromosomal evolution since their divergence from a common ancestor. We have compared gene orders on three chromosomes of Drosophila pseudoobscura with its close relative, D. miranda, and the distant outgroup species, D. melanogaster, by using the public genome sequences of D. pseudoobscura and D. melanogaster and approximately 50 in situ hybridizations of gene probes in D. miranda. We find no evidence for extensive transfer of genes among chromosomes in D. miranda. The rates of chromosomal rearrangements between D. miranda and D. pseudoobscura are far higher than those found before in Drosophila and approach those for nematodes, the fastest rates among higher eukaryotes. In addition, we find that the D. pseudoobscura chromosome with the highest level of inversion polymorphism (Muller's element C) does not show an unusually fast rate of evolution with respect to chromosome structure, suggesting that this classic case of inversion polymorphism reflects selection rather than mutational processes. On the basis of our results, we propose possible ancestral arrangements for the D. pseudoobscura C chromosome, which are different from those in the current literature. We also describe a new method for correcting for rearrangements that are not detected with a limited set of markers.

Possible effects of (trans)gene flow from crops on the genetic diversity from landraces and wild relatives

Gene flow is a potential concern associated with the use of transgenic crops because it could affect genetic diversity of related landraces and wild relatives. This concern has taken on added importance with the looming introduction of transgenic crops in centers of crop domestication (Mexico, China) and those producing pharmaceutical compounds. For gene flow to take place among cultivars and their wild relatives, several steps have to be fulfilled, including the presence of cultivars or wild relatives within pollen or seed dispersal range, the ability to produce viable and fertile hybrids, at least partial overlap in flowering time, actual gene flow by pollen or seed, and the establishment of crop genes in the domesticated or wild recipient populations. In contrast with domestication genes, which often make crops less adapted to natural ecosystems, transgenes frequently represent gains of function, which might release wild relatives from constraints that limit their fitness. In most sexually reproducing organisms, the chromosomal region affected by selection of a single gene amounts to a small percentage of the total genome size. Because of gene flow, the level of genetic diversity present in the domesticated gene pool becomes a crucial factor affecting the genetic diversity of the wild gene pool. For some crops, such as cotton and maize, the introduction of transgenic technologies has led to a consolidation of the seed industry and a reduction in the diversity of the elite crop gene pool. Thus, diversity in improved varieties grown by farmers needs to be monitored. Several areas deserve further study, such as the actual magnitude of gene flow and its determinants in different agroecosystems, the long-term effects of gene flow on genetic diversity both across gene pools and within genomes, the expression of transgenes in new genetic backgrounds, and the effects of socio-economic factors on genetic diversity.

Variation after a selective sweep in a subdivided population

The effect of genetic hitchhiking on neutral variation is analyzed in subdivided populations with differentiated demes. After fixation of a favorable mutation, the consequences on particular subpopulations can be radically different. In the subpopulation where the mutation first appeared by mutation, variation at linked neutral loci is expected to be reduced, as predicted by the classical theory. However, the effect in the other subpopulations, where the mutation is introduced by migration, can be the opposite. This effect depends on the level of genetic differentiation of the subpopulations, the selective advantage of the mutation, the recombination frequency, and the population size, as stated by analytical derivations and computer simulations. The characteristic outcomes of the effect are three. First, the genomic region of reduced variation around the selected locus is smaller than that predicted in a panmictic population. Second, for more distant neutral loci, the amount of variation increases over the level they had before the hitchhiking event. Third, for these loci, the spectrum of gene frequencies is dominated by an excess of alleles at intermediate frequencies when compared with the neutral theory. At these loci, hitchhiking works like a system that takes variation from the between-subpopulation component and introduces it into the subpopulations. The mechanism can also operate in other systems in which the genetic variation is distributed in clusters with limited exchange of variation, such as chromosome arrangements or genomic regions closely linked to targets of balancing selection.

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.

Dntf-2r, a young Drosophila retroposed gene with specific male expression under positive Darwinian selection

A direct approach to investigating new gene origination is to examine recently evolved genes. We report a new gene in the Drosophila melanogaster subgroup, Drosophila nuclear transport factor-2-related (Dntf-2r). Its sequence features and phylogenetic distribution indicate that Dntf-2r is a retroposed functional gene and originated in the common ancestor of D. melanogaster, D. simulans, D. sechellia, and D. mauritiana, within the past 3-12 million years (MY). Dntf-2r evolved more rapidly than the parental gene, under positive Darwinian selection as revealed by the McDonald-Kreitman test and other evolutionary analyses. Comparative expression analysis shows that Dntf-2r is male specific whereas the parental gene, Dntf-2, is widely expressed in D. melanogaster. In agreement with its new expression pattern, the Dntf-2r putative promoter sequence is similar to the late testis promoter of beta2-tubulin. We discuss the possibility that the action of positive selection in Dntf-2r is related to its putative male-specific functions.

Low occurrence of gene transposition events during the evolution of the genus Drosophila

The role played by gene transpositions during the evolution of eukaryotic genomes is still poorly understood and indeed has been analyzed in detail only in nematodes. In Drosophila, a limited number of transpositions have been detected by comparing the chromosomal location of genes between different species. The relative importance of gene transposition versus other types of chromosomal rearrangements, for example, inversions, has not yet been evaluated. Here, we use physical mapping to perform an extensive search for long-distance gene transpositions and assess their impact during the evolution of the Drosophila genome. We compare the relative order of 297 molecular markers that cover 60% of the euchromatic fraction of the genome between two related Drosophila species and conclude that the frequency of gene transpositions is very low, namely one order of magnitude lower than that of nematodes. In addition, gene transpositions seem to be events almost exclusively associated with genes of repetitive nature such as the Histone gene complex (HIS-C).

Sexual antagonism and X inactivation--the SAXI hypothesis

X inactivation has evolved in the soma of mammalian females so that both sexes have the same ratio of X:autosomal gene expression. The X chromosome in the germ cells of XY males is also precociously inactivated for reasons that remain unclear. Unlike X inactivation in the soma, this germline X inactivation is not restricted to mammals but has evolved independently in several animal phyla. Thus, germline X inactivation might have been the precursor of somatic X inactivation in mammals. We now propose a hypothesis for the evolution of germline X inactivation. The hypothesis predicts a redistribution of late spermatogenic genes from the X chromosome to the autosomes, leading eventually to germline X inactivation as the X chromosome becomes 'demasculinized'. Sexual antagonism could be the mechanism driving this redistribution. Recent expression and genetic studies in mammals, nematodes and Drosophila support this hypothesis, and expression data on taxa that have not evolved germline X inactivation, such as birds and butterflies, should shed further light on it.

Chromosomal elements evolve at different rates in the Drosophila genome

Recent results indicate that the rate of chromosomal rearrangement in the genus Drosophila is the highest found so far in any eukaryote. This conclusion is based chiefly on the comparative mapping analysis of a single chromosomal element (Muller's element E) in two species, D. melanogaster and D. repleta, representing the two farthest lineages within the genus (the Sophophora and Drosophila subgenera, respectively). We have extended the analysis to two other chromosomal elements (Muller's elements A and D) and tested for differences in rate of evolution among chromosomes. With this purpose, detailed physical maps of chromosomes X and 4 of D. repleta were constructed by in situ hybridization of 145 DNA probes (gene clones, cosmids, and P1 phages) and their gene arrangements compared with those of the homologous chromosomes X and 3L of D. melanogaster. Both chromosomal elements have been extensively reshuffled over their entire length. The number of paracentric inversions fixed has been estimated as 118 +/- 17 for element A and 56 +/- 8 for element D. Comparison with previous data for elements E and B shows that there are fourfold differences in evolution rate among chromosomal elements, with chromosome X exhibiting the highest rate of rearrangement. Combining all results, we estimated that 393 paracentric inversions have been fixed in the whole genome since the divergence between D. repleta and D. melanogaster. This amounts to an average rate of 0.053 disruptions/Mb/myr, corroborating the high rate of rearrangement in the genus Drosophila.

Adaptive hitchhiking effects on genome variability

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

Regions of lower crossing over harbor more rare variants in African populations of Drosophila melanogaster

A correlation between diversity levels and rates of recombination is predicted both by models of positive selection, such as hitchhiking associated with the rapid fixation of advantageous mutations, and by models of purifying selection against strongly deleterious mutations (commonly referred to as "background selection"). With parameter values appropriate for Drosophila populations, only the first class of models predicts a marked skew in the frequency spectrum of linked neutral variants, relative to a neutral model. Here, we consider 29 loci scattered throughout the Drosophila melanogaster genome. We show that, in African populations, a summary of the frequency spectrum of polymorphic mutations is positively correlated with the meiotic rate of crossing over. This pattern is demonstrated to be unlikely under a model of background selection. Models of weakly deleterious selection are not expected to produce both the observed correlation and the extent to which nucleotide diversity is reduced in regions of low (but nonzero) recombination. Thus, of existing models, hitchhiking due to the recurrent fixation of advantageous variants is the most plausible explanation for the data.