Molecular evolution between Drosophila melanogaster and D. simulans: reduced codon bias, faster rates of amino acid substitution, and larger proteins in D. melanogaster

Elevated evolutionary rates in the laboratory strain of Saccharomyces cerevisiae

By using the maximum likelihood method, we made a genome-wide comparison of the evolutionary rates in the lineages leading to the laboratory strain (S288c) and a wild strain (YJM789) of Saccharomyces cerevisiae and found that genes in the laboratory strain tend to evolve faster than in the wild strain. The pattern of elevated evolution suggests that relaxation of selection intensity is the dominant underlying reason, which is consistent with recurrent bottlenecks in the S. cerevisiae laboratory strain population. Supporting this conclusion are the following observations: (i) the increases in nonsynonymous evolutionary rate occur for genes in all functional categories; (ii) most of the synonymous evolutionary rate increases in S288c occur in genes with strong codon usage bias; (iii) genes under stronger negative selection have a larger increase in nonsynonymous evolutionary rate; and (iv) more genes with adaptive evolution were detected in the laboratory strain, but they do not account for the majority of the increased evolution. The present discoveries suggest that experimental and possible industrial manipulations of the laboratory strain of yeast could have had a strong effect on the genetic makeup of this model organism. Furthermore, they imply an evolution of laboratory model organisms away from their wild counterparts, questioning the relevancy of the models especially when extensive laboratory cultivation has occurred. In addition, these results shed light on the evolution of livestock and crop species that have been under human domestication for years.

Patterns of selection on synonymous and nonsynonymous variants in Drosophila miranda

We have investigated patterns of within-species polymorphism and between-species divergence for synonymous and nonsynonymous variants at a set of autosomal and X-linked loci of Drosophila miranda. D. pseudoobscura and D. affinis were used for the between-species comparisons. The results suggest the action of purifying selection on nonsynonymous, polymorphic variants. Among synonymous polymorphisms, there is a significant excess of synonymous mutations from preferred to unpreferred codons and of GC to AT mutations. There was no excess of GC to AT mutations among polymorphisms at noncoding sites. This suggests that selection is acting to maintain the use of preferred codons. Indirect evidence suggests that biased gene conversion in favor of GC base pairs may also be operating. The joint intensity of selection and biased gene conversion, in terms of the product of effective population size and the sum of the selection and conversion coefficients, was estimated to be approximately 0.65.

Genomic heterogeneity of background substitutional patterns in Drosophila melanogaster

Mutation is the underlying force that provides the variation upon which evolutionary forces can act. It is important to understand how mutation rates vary within genomes and how the probabilities of fixation of new mutations vary as well. If substitutional processes across the genome are heterogeneous, then examining patterns of coding sequence evolution without taking these underlying variations into account may be misleading. Here we present the first rigorous test of substitution rate heterogeneity in the Drosophila melanogaster genome using almost 1500 nonfunctional fragments of the transposable element DNAREP1_DM. Not only do our analyses suggest that substitutional patterns in heterochromatic and euchromatic sequences are different, but also they provide support in favor of a recombination-associated substitutional bias toward G and C in this species. The magnitude of this bias is entirely sufficient to explain recombination-associated patterns of codon usage on the autosomes of the D. melanogaster genome. We also document a bias toward lower GC content in the pattern of small insertions and deletions (indels). In addition, the GC content of noncoding DNA in Drosophila is higher than would be predicted on the basis of the pattern of nucleotide substitutions and small indels. However, we argue that the fast turnover of noncoding sequences in Drosophila makes it difficult to assess the importance of the GC biases in nucleotide substitutions and small indels in shaping the base composition of noncoding sequences.

Linking phylogenetics with population genetics to reconstruct the geographic origin of a species

Reconstructing ancestral geographic origins is critical for understanding the long-term evolution of a species. Bayesian methods have been proposed to test biogeographic hypotheses while accommodating uncertainty in phylogenetic reconstruction. However, the problem that certain taxa may have a disproportionate influence on conclusions has not been addressed. Here, we infer the geographic origin of Drosophila simulans using 2,014 bp of the period locus from 63 lines collected from 18 countries. We also analyze two previously published datasets, alcohol dehydrogenase related and NADH:ubiquinone reductase 75 kDa subunit precursor. Phylogenetic inferences of all three loci support Madagascar as the geographic origin of D. simulans. Our phylogenetic conclusions are robust to taxon resampling and to the potentially confounding effects of recombination. To test our phylogenetically derived hypothesis we develop a randomization test of the population genetics prediction that sequences from the geographic origin should contain more genetic polymorphism than those from derived populations. We find that the Madagascar population has elevated genetic polymorphism relative to non-Madagascar sequences. These data are corroborated by mitochondrial DNA sequence data.

Experimental reduction of codon bias in the Drosophila alcohol dehydrogenase gene results in decreased ethanol tolerance of adult flies

The ethanol tolerance of adult transgenic flies of Drosophila containing between zero and ten unpreferred synonymous mutations that reduced codon bias in the alcohol dehydrogenase (Adh) gene was assayed. As the amino acid sequences of the ADH protein were identical in the four genotypes assayed, differences in ethanol tolerance were due to differences in the abundance of ADH protein, presumably driven by the effects of codon bias on translational efficiency. The ethanol tolerance of genotypes decreased with the number of unpreferred synonymous mutations, and a positive correlation between ADH protein abundance and ethanol tolerance was observed. This work confirms that the fitness effects of unpreferred synonymous mutations that reduce codon bias in a highly expressed gene are experimentally measurable in Drosophila melanogaster.

Molecular phylogeny of the Drosophila melanogaster species subgroup

Although molecular and phenotypic evolution have been studied extensively in Drosophila melanogaster and its close relatives, phylogenetic relationships within the D. melanogaster species subgroup remain unresolved. In particular, recent molecular studies have not converged on the branching orders of the D. yakuba-D. teissieri and D. erecta-D. orena species pairs relative to the D. melanogaster-D. simulans-D. mauritiana-D. sechellia species complex. Here, we reconstruct the phylogeny of the melanogaster species subgroup using DNA sequence data from four nuclear genes. We have employed "vectorette PCR" to obtain sequence data for orthologous regions of the Alcohol dehydrogenase (Adh), Alcohol dehydrogenase related (Adhr), Glucose dehydrogenase (Gld), and rosy (ry) genes (totaling 7164 bp) from six melanogaster subgroup species (D. melanogaster, D. simulans, D. teissieri, D. yakuba, D. erecta, and D. orena) and three species from subgroups outside the melanogaster species subgroup [D. eugracilis (eugracilis subgroup), D. mimetica (suzukii subgroup), and D. lutescens (takahashii subgroup)]. Relationships within the D. simulans complex are not addressed. Phylogenetic analyses employing maximum parsimony, neighbor-joining, and maximum likelihood methods strongly support a D. yakuba-D. teissieri and D. erecta-D. orena clade within the melanogaster species subgroup. D. eugracilis is grouped closer to the melanogaster subgroup than a D. mimetica-D. lutescens clade. This tree topology is supported by reconstructions employing simple (single parameter) and more complex (nonreversible) substitution models.

Patterns of evolutionary constraints in intronic and intergenic DNA of Drosophila

We develop methods to infer levels of evolutionary constraints in the genome by comparing rates of nucleotide substitution in noncoding DNA with rates predicted from rates of synonymous site evolution in adjacent genes or other putatively neutrally evolving sites, while accounting for differences in base composition. We apply the methods to estimate levels of constraint in noncoding DNA of Drosophila. In introns, constraint (the estimated fraction of mutations that are selectively eliminated) is absolute at the 5' and 3' splice junction dinucleotides, and averages 72% in base pairs 3-6 at the 5'-end. Constraint at the 5' base pairs 3-6 is significantly lower in the lineage leading to Drosophila melanogaster than in Drosophila simulans, a finding that agrees with other features of genome evolution in Drosophila and indicates that the effect of selection on intron function has been weaker in the melanogaster lineage. Elsewhere in intron sequences, the rate of nucleotide substitution is significantly higher than at synonymous sites. By using intronic sites outside splice control regions as a putative neutrally evolving standard, constraint in the 500 bp of intergenic DNA upstream and downstream regions of protein-coding genes averages approximately 44%. Although the estimated level of constraint in intergenic regions close to genes is only about one-half of that of amino acid sites, selection against single-nucleotide mutations in intergenic DNA makes a substantial contribution to the mutation load in Drosophila.

Selection on codon usage in Drosophila americana

Synonymous codons are not used at random, significantly influencing the base composition of the genome. The selection-mutation-drift model proposes that this bias reflects natural selection in favor of a subset of preferred codons. Previous estimates in Drosophila of the intensity of selective forces involved seem too large to be reconciled with theoretical predictions of the level of codon bias. This probably results from confounding effects of the demographic histories of the species concerned. We have studied three species of the virilis group of Drosophila, which are more likely to satisfy the assumptions of the evolutionary models. We analyzed the patterns of polymorphism and divergence in a sample of 18 genes and applied a new method for estimating the intensity of selection on synonymous mutations based on the frequencies of unpreferred mutations among polymorphic sites. This yielded estimates of selection intensities (N(e)s) of the order of 0.65, which is more compatible with the observed levels of codon bias. Our results support the action of both selection and mutational bias on codon usage bias and suggest that codon usage and genome base composition in the D. americana lineage are in approximate equilibrium. Biased gene conversion may also contribute to the observed patterns.

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.

Estimating the distribution of fitness effects from DNA sequence data: implications for the molecular clock

We present a method for estimating the distribution of fitness effects of new amino acid mutations when those mutations can be assumed to be slightly advantageous, slightly deleterious, or strongly deleterious. We apply the method to mitochondrial data from several different species. In the majority of the data sets, the shape of the distribution is approximately exponential. Our results provide an estimate of the distribution of fitness effects of weakly selected mutations and provide a possible explanation for why the molecular clock is fairly constant across taxa and time.

Translational selection and yeast proteome evolution

The primary structures of peptides may be adapted for efficient synthesis as well as proper function. Here, the Saccharomyces cerevisiae genome sequence, DNA microarray expression data, tRNA gene numbers, and functional categorizations of proteins are employed to determine whether the amino acid composition of peptides reflects natural selection to optimize the speed and accuracy of translation. Strong relationships between synonymous codon usage bias and estimates of transcript abundance suggest that DNA array data serve as adequate predictors of translation rates. Amino acid usage also shows striking relationships with expression levels. Stronger correlations between tRNA concentrations and amino acid abundances among highly expressed proteins than among less abundant proteins support adaptation of both tRNA abundances and amino acid usage to enhance the speed and accuracy of protein synthesis. Natural selection for efficient synthesis appears to also favor shorter proteins as a function of their expression levels. Comparisons restricted to proteins within functional classes are employed to control for differences in amino acid composition and protein size that reflect differences in the functional requirements of proteins expressed at different levels.

Natural selection drives Drosophila immune system evolution

Evidence from disparate sources suggests that natural selection may often play a role in the evolution of host immune system proteins. However, there have been few attempts to make general population genetic inferences on the basis of analysis of several immune-system-related genes from a single species. Here we present DNA polymorphism and divergence data from 34 genes thought to function in the innate immune system of Drosophila simulans and compare these data to those from 28 nonimmunity genes sequenced from the same lines. Several statistics, including average K(A)/K(S) ratio, average silent heterozygosity, and average haplotype diversity, significantly differ between the immunity and nonimmunity genes, suggesting an important role for directional selection in immune system protein evolution. In contrast to data from mammalian immunoglobulins and other proteins, we find no strong evidence for the selective maintenance of protein diversity in Drosophila immune system proteins. This may be a consequence of Drosophila's generalized innate immune response.

In vivo introduction of unpreferred synonymous codons into the Drosophila Adh gene results in reduced levels of ADH protein

The evolution of codon bias, the unequal usage of synonymous codons, is thought to be due to natural selection for the use of preferred codons that match the most abundant species of isoaccepting tRNA, resulting in increased translational efficiency and accuracy. We examined this hypothesis by introducing 1, 6, and 10 unpreferred codons into the Drosophila alcohol dehydrogenase gene (Adh). We observed a significant decrease in ADH protein production with number of unpreferred codons, confirming the importance of natural selection as a mechanism leading to codon bias. We then used this empirical relationship to estimate the selection coefficient (s) against unpreferred synonymous mutations and found the value (s >or= 10(-5)) to be approximately one order of magnitude greater than previous estimates from population genetics theory. The observed differences in protein production appear to be too large to be consistent with current estimates of the strength of selection on synonymous sites in D. melanogaster.

Quantifying the slightly deleterious mutation model of molecular evolution

We have attempted to quantify the frequency and effects of slightly deleterious mutations (SDMs), those that have selective effects close to the reciprocal of the effective population size of a species, by comparing the level of selective constraint in protein-coding genes of related species that have different present-day effective population sizes. In our two comparisons, the species with the smaller effective population size showed lower constraint, implying that SDMs had become fixed. The fixation of SDMs was supported by the observation of a higher fraction of radical to conservative amino acid substitutions in species with smaller effective population sizes. The fraction of strongly deleterious mutations (which rarely become fixed) is >70% in most species. Only approximately 10% or fewer of mutations seem to behave as SDMs, but SDMs could comprise a substantial fraction of mutations in protein-coding genes that have a chance of becoming fixed between species.

Molecular population genetics of Xdh and the evolution of base composition in Drosophila

Few loci have been measured for DNA polymorphism and divergence in several species. Here we report such data from the protein-coding region of xanthine dehydrogenase (Xdh) in 22 species of Drosophila. Many of our samples were from closely related species, allowing us to confidently assign substitutions to individual lineages. Surprisingly, Xdh appears to be fixing more A/T mutations than G/C mutations in most lineages, leading to evolution of higher A/T content in the recent past. We found no compelling evidence for selection on protein variation, though some aspects of the data support the notion that a significant fraction of amino acid polymorphisms are slightly deleterious. Finally, we found no convincing evidence that levels of silent heterozygosity are associated with rates of protein evolution.

Changing effective population size and the McDonald-Kreitman test

Artifactual evidence of adaptive amino acid substitution can be generated within a McDonald-Kreitman test if some amino acid mutations are slightly deleterious and there has been an increase in effective population size. Here I investigate the conditions under which this occurs. I show that fairly small increases in effective population size can generate artifactual evidence of positive selection if there is no selection upon synonymous codon use. This problem is exacerbated by the removal of low-frequency polymorphisms. However, selection on synonymous codon use restricts the conditions under which artifactual evidence of adaptive evolution is produced.

Rates and patterns of molecular evolution in inbred and outbred Arabidopsis

The evolution of self-fertilization is associated with a large reduction in the effective rate of recombination and a corresponding decline in effective population size. If many spontaneous mutations are slightly deleterious, this shift in the breeding system is expected to lead to a reduced efficacy of natural selection and genome-wide changes in the rates of molecular evolution. Here, we investigate the effects of the breeding system on molecular evolution in the highly self-fertilizing plant Arabidopsis thaliana by comparing its coding and noncoding genomic regions with those of its close outcrossing relative, the self-incompatible A. lyrata. More distantly related species in the Brassicaceae are used as outgroups to polarize the substitutions along each lineage. In contrast to expectations, no significant difference in the rates of protein evolution is observed between selfing and outcrossing Arabidopsis species. Similarly, no consistent overall difference in codon bias is observed between the species, although for low-biased genes A. lyrata shows significantly higher major codon usage. There is also evidence of intron size evolution in A. thaliana, which has consistently smaller introns than its outcrossing congener, potentially reflecting directional selection on intron size. The results are discussed in the context of heterogeneity in selection coefficients across loci and the effects of life history and population structure on rates of molecular evolution. Using estimates of substitution rates in coding regions and approximate estimates of divergence and generation times, the genomic deleterious mutation rate (U) for amino acid substitutions in Arabidopsis is estimated to be approximately 0.2-0.6 per generation.

Codon bias differentiates between the duplicated amylase loci following gene duplication in Drosophila

We examined the pattern of synonymous substitutions in the duplicated Amylase (Amy) genes (called the Amy1- and Amy3-type genes, respectively) in the Drosophila montium species subgroup. The GC content at the third synonymous codon sites of the Amy1-type genes was higher than that of the Amy3-type genes, while the GC content in the 5'-flanking region was the same in both genes. This suggests that the difference in the GC content at third synonymous sites between the duplicated genes is not due to the temporal or regional changes in mutation bias. We inferred the direction of synonymous substitutions along branches of a phylogeny. In most lineages, there were more synonymous substitutions from G/C (G or C) to A/T (A or T) than from A/T to G/C. However, in one lineage leading to the Amy1-type genes, which is immediately after gene duplication but before speciation of the montium species, synonymous substitutions from A/T to G/C were predominant. According to a simple model of synonymous DNA evolution in which major codons are selectively advantageous within each codon family, we estimated the selection intensity for specific lineages in a phylogeny on the basis of inferred patterns of synonymous substitutions. Our result suggested that the difference in GC content at synonymous sites between the two Amy-type genes was due to the change of selection intensity immediately after gene duplication but before speciation of the montium species.

On the abundance and distribution of transposable elements in the genome of Drosophila melanogaster

The abundance and distribution of transposable elements (TEs) in a representative part of the euchromatic genome of Drosophila melanogaster were studied by analyzing the sizes and locations of TEs of all known families in the genomic sequences of chromosomes 2R, X, and 4. TEs contribute to up to 2% of the sequenced DNA, which corresponds roughly to the euchromatin of these chromosomes. This estimate is lower than that previously available from in situ data and suggests that TEs accumulate in the heterochromatin more intensively than was previously thought. We have also found that TEs are not distributed at random in the chromosomes and that their abundance is more strongly associated with local recombination rates, rather than with gene density. The results are compatible with the ectopic exchange model, which proposes that selection against deleterious effects of chromosomal rearrangements is a major force opposing element spread in the genome of this species. Selection against insertional mutations also influences the observed patterns, such as an absence of insertions in coding regions. The results of the analyses are discussed in the light of recent findings on the distribution of TEs in other species.

Population, evolutionary and genomic consequences of interference selection

Weakly selected mutations are most likely to be physically clustered across genomes and, when sufficiently linked, they alter each others' fixation probability, a process we call interference selection (IS). Here we study population genetics and evolutionary consequences of IS on the selected mutations themselves and on adjacent selectively neutral variation. We show that IS reduces levels of polymorphism and increases low-frequency variants and linkage disequilibrium, in both selected and adjacent neutral mutations. IS can account for several well-documented patterns of variation and composition in genomic regions with low rates of crossing over in Drosophila. IS cannot be described simply as a reduction in the efficacy of selection and effective population size in standard models of selection and drift. Rather, IS can be better understood with models that incorporate a constant "traffic" of competing alleles. Our simulations also allow us to make genome-wide predictions that are specific to IS. We show that IS will be more severe at sites in the center of a region containing weakly selected mutations than at sites located close to the edge of the region. Drosophila melanogaster genomic data strongly support this prediction, with genes without introns showing significantly reduced codon bias in the center of coding regions. As expected, if introns relieve IS, genes with centrally located introns do not show reduced codon bias in the center of the coding region. We also show that reasonably small differences in the length of intermediate "neutral" sequences embedded in a region under selection increase the effectiveness of selection on the adjacent selected sequences. Hence, the presence and length of sequences such as introns or intergenic regions can be a trait subject to selection in recombining genomes. In support of this prediction, intron presence is positively correlated with a gene's codon bias in D. melanogaster. Finally, the study of temporal dynamics of IS after a change of recombination rate shows that nonequilibrium codon usage may be the norm rather than the exception.

Interactions between natural selection, recombination and gene density in the genes of Drosophila

In Drosophila, as in many organisms, natural selection leads to high levels of codon bias in genes that are highly expressed. Thus codon bias is an indicator of the intensity of one kind of selection that is experienced by genes and can be used to assess the impact of other genomic factors on natural selection. Among 13,000 genes in the Drosophila genome, codon bias has a slight positive, and strongly significant, association with recombination--as expected if recombination allows natural selection to act more efficiently when multiple linked sites segregate functional variation. The same reasoning leads to the expectation that the efficiency of selection, and thus average codon bias, should decline with gene density. However, this prediction is not confirmed. Levels of codon bias and gene expression are highest for those genes in an intermediate range of gene density, a pattern that may be the result of a tradeoff between the advantages for gene expression of close gene spacing and disadvantages arising from regulatory conflicts among tightly packed genes. These factors appear to overlay the more subtle effect of linkage among selected sites that gives rise to the association between recombination rate and codon bias.

Patterns of genetic variation at a chromosome 4 locus of Drosophila melanogaster and D. simulans

DNA sequence surveys of Drosophila melanogaster populations show a strong positive correlation between the recombination rate experienced by a locus and its level of nucleotide polymorphism. In particular, surveys of the fourth chromosome gene ci(D) show greatly reduced levels of nucleotide variation; this observation was originally interpreted in terms of selective sweeps occurring on the nonrecombining fourth chromosome. Subsequent theoretical work has, however, uncovered several other selective processes that can reduce variation. In this study, we revisit the Drosophila fourth chromosome, investigating variation in 5-6 kb of the gene ankyrin in D. melanogaster and D. simulans. Silent nucleotide site diversity is approximately 5 x 10(-4) for both species, consistent with the previous observations of low variation at ci(D). Given the observed frequency spectra at ankyrin, coalescent simulations indicate that reduced diversity in the region is unlikely to be due to a selective sweep alone. We find evidence for recombinational exchange at this locus, and both species appear to be fixed for an insertion of the transposable element HB in an intron of ankyrin.

A maximum likelihood method for analyzing pseudogene evolution: implications for silent site evolution in humans and rodents

We present a new likelihood method for detecting constrained evolution at synonymous sites and other forms of nonneutral evolution in putative pseudogenes. The model is applicable whenever the DNA sequence is available from a protein-coding functional gene, a pseudogene derived from the protein-coding gene, and an orthologous functional copy of the gene. Two nested likelihood ratio tests are developed to test the hypotheses that (1) the putative pseudogene has equal rates of silent and replacement substitutions; and (2) the rate of synonymous substitution in the functional gene equals the rate of substitution in the pseudogene. The method is applied to a data set containing 74 human processed-pseudogene loci, 25 mouse processed-pseudogene loci, and 22 rat processed-pseudogene loci. Using the informatics resources of the Human Genome Project, we localized 67 of the human-pseudogene pairs in the genome and estimated the GC content of a large surrounding genomic region for each. We find that, for pseudogenes deposited in GC regions similar to those of their paralogs, the assumption of equal rates of silent and replacement site evolution in the pseudogene is upheld; in these cases, the rate of silent site evolution in the functional genes is approximately 70% the rate of evolution in the pseudogene. On the other hand, for pseudogenes located in genomic regions of much lower GC than their functional gene, we see a sharp increase in the rate of silent site substitutions, leading to a large rate of rejection for the pseudogene equality likelihood ratio test.

Gene expression and molecular evolution

The combination of complete genome sequence information and estimates of mRNA abundances have begun to reveal causes of both silent and protein sequence evolution. Translational selection appears to explain patterns of synonymous codon usage in many prokaryotes as well as a number of eukaryotic model organisms (with the notable exception of vertebrates). Relationships between gene length and codon usage bias, however, remain unexplained. Intriguing correlations between expression patterns and protein divergence suggest some general mechanisms underlying protein evolution.

The evolutionary analysis of "orphans" from the Drosophila genome identifies rapidly diverging and incorrectly annotated genes

In genome projects of eukaryotic model organisms, a large number of novel genes of unknown function and evolutionary history ("orphans") are being identified. Since many orphans have no known homologs in distant species, it is unclear whether they are restricted to certain taxa or evolve rapidly, either because of a lack of constraints or positive Darwinian selection. Here we use three criteria for the selection of putatively rapidly evolving genes from a single sequence of Drosophila melanogaster. Thirteen candidate genes were chosen from the Adh region on the second chromosome and 1 from the tip of the X chromosome. We succeeded in obtaining sequence from 6 of these in the closely related species D. simulans and D. yakuba. Only 1 of the 6 genes showed a large number of amino acid replacements and in-frame insertions/deletions. A population survey of this gene suggests that its rapid evolution is due to the fixation of many neutral or nearly neutral mutations. Two other genes showed "normal" levels of divergence between species. Four genes had insertions/deletions that destroy the putative reading frame within exons, suggesting that these exons have been incorrectly annotated. The evolutionary analysis of orphan genes in closely related species is useful for the identification of both rapidly evolving and incorrectly annotated genes.

The frequency distribution of nucleotide variation in Drosophila simulans

Patterns of codon bias in Drosophila suggest that silent mutations can be classified into two types: unpreferred (slightly deleterious) and preferred (slightly beneficial). Results of previous analyses of polymorphism and divergence in Drosophila simulans were interpreted as supporting a mutation-selection-drift model in which slightly deleterious, silent mutants make significantly greater contributions to polymorphism than to divergence. Frequencies of unpreferred polymorphisms were inferred to be lower than frequencies of other silent polymorphisms. Here, I analyzed additional D. simulans data to reevaluate the support for these ideas. I found that D. simulans has fixed more unpreferred than preferred mutations, suggesting that this lineage has not been at mutation-selection-drift equilibrium at silent sites. Frequencies of polarized unpreferred polymorphisms are not skewed toward rare alleles. However, frequencies of unpolarized unpreferred codons are lower in high-bias genes than in low-bias genes. This supports the idea that unpreferred codons are borderline deleterious mutations. Purifying selection on silent sites appears to be stronger at twofold-degenerate codons than at fourfold-degenerate codons. Finally, I found that X-linked polymorphisms occur at a higher average frequency than polymorphisms on chromosome arm 3R, even though an average X-linked site is significantly less likely to be polymorphic than an average site on 3R. This result supports a previous analysis of D. simulans indicating different population genetics of X-linked versus autosomal mutations.

Genome evolution at the genus level: comparison of three complete genomes of hyperthermophilic archaea

We have compared three complete genomes of closely related hyperthermophilic species of Archaea belonging to the Pyrococcus genus: Pyrococcus abyssi, Pyrococcus horikoshii, and Pyrococcus furiosus. At the genomic level, the comparison reveals a differential conservation among four regions of the Pyrococcus chromosomes correlated with the location of genetic elements mediating DNA reorganization. This discloses the relative contribution of the major mechanisms that promote genomic plasticity in these Archaea, namely rearrangements linked to the replication terminus, insertion sequence-mediated recombinations, and DNA integration within tRNA genes. The combination of these mechanisms leads to a high level of genomic plasticity in these hyperthermophilic Archaea, at least comparable to the plasticity observed between closely related bacteria. At the proteomic level, the comparison of the three Pyrococcus species sheds light on specific selection pressures acting both on their coding capacities and evolutionary rates. Indeed, thanks to two independent methods, the "reciprocal best hits" approach and a new distance ratio analysis, we detect the false orthology relationships within the Pyrococcus lineage. This reveals a high amount of differential gains and losses of genes since the divergence of the three closely related species. The resulting polymorphism is probably linked to an adaptation of these free-living organisms to differential environmental constraints. As a corollary, we delineate the set of orthologous genes shared by the three species, that is, the genes that may characterize the Pyrococcus genus. In this conserved core, the amino acid substitution rate is equal between P. abyssi and P. horikoshii for most of their shared proteins, even for fast-evolving ones. In contrast, strong discrepancies exist among the substitution rates observed in P. furiosus relative to the two other species, which is in disagreement with the molecular clock hypothesis.

Genome Evolution at the Genus Level: Comparison of Three Complete Genomes of Hyperthermophilic Archaea

We have compared three complete genomes of closely related hyperthermophilic species of Archaea belonging to thePyrococcus genus: Pyrococcus abyssi, Pyrococcus horikoshii, and Pyrococcus furiosus. At the genomic level, the comparison reveals a differential conservation among four regions of the Pyrococcus chromosomes correlated with the location of genetic elements mediating DNA reorganization. This discloses the relative contribution of the major mechanisms that promote genomic plasticity in these Archaea, namely rearrangements linked to the replication terminus, insertion sequence-mediated recombinations, and DNA integration within tRNA genes. The combination of these mechanisms leads to a high level of genomic plasticity in these hyperthermophilic Archaea, at least comparable to the plasticity observed between closely related bacteria. At the proteomic level, the comparison of the threePyrococcus species sheds light on specific selection pressures acting both on their coding capacities and evolutionary rates. Indeed, thanks to two independent methods, the “reciprocal best hits“ approach and a new distance ratio analysis, we detect the false orthology relationships within the Pyrococcus lineage. This reveals a high amount of differential gains and losses of genes since the divergence of the three closely related species. The resulting polymorphism is probably linked to an adaptation of these free-living organisms to differential environmental constraints. As a corollary, we delineate the set of orthologous genes shared by the three species, that is, the genes that may characterize the Pyrococcus genus. In this conserved core, the amino acid substitution rate is equal between P. abyssi and P. horikoshii for most of their shared proteins, even for fast-evolving ones. In contrast, strong discrepancies exist among the substitution rates observed in P. furiosus relative to the two other species, which is in disagreement with the molecular clock hypothesis.

Local changes in GC/AT substitution biases and in crossover frequencies on Drosophila chromosomes

I present here evidence of remarkable local changes in GC/AT substitution biases and in crossover frequencies on Drosophila chromosomes. The substitution pattern at 10 loci in the telomeric region of the X chromosome was studied for four species of the Drosophila melanogaster species subgroup. Drosophila orena and Drosophila erecta are clearly the most closely related species pair (the erecta complex) among the four species studied; however, the overall data at the 10 loci revealed a clear dichotomy in the silent substitution patterns between the AT-biased- substitution melanogaster and erecta lineages and the GC-biased-substitution yakuba and orena lineages, suggesting two or more independent changes in GC/AT substitution biases. More importantly, the results indicated a between- loci heterogeneity in GC/AT substitution bias in this small region independently in the yakuba and orena lineages. Indeed, silent substitutions in the orena lineage were significantly biased toward G and C at the consecutive yellow, lethal of scute, and asense loci, but they were significantly biased toward A and T at sta. The substitution bias toward G and C was centered in different areas in yakuba (significantly biased at EG:165H7.3, EG:171D11.2, and suppressor of sable). The similar silent substitution patterns in coding and noncoding regions, furthermore, suggested mutational biases as a cause of the substitution biases. On the other hand, previous study reveals that Drosophila yakuba has about 20-fold higher crossover frequencies in the telomeric region of the X chromosome than does D. melanogaster; this study revealed that the total genetic map length of the yakuba X chromosome was only about 1.5 times as large as that of melanogaster and that the map length of the X-telomeric y-sta region did not differ between Drosophila yakuba and D. erecta. Taken together, the data strongly suggested that an approximately 20- fold reduction in the X-telomeric crossover frequencies occurred in the ancestral population of D. melanogaster after the melanogaster-yakuba divergence but before the melanogaster-simulans divergence.

Intraspecific variation in symbiont genomes: bottlenecks and the aphid-buchnera association

Buchnera are maternally transmitted bacterial endosymbionts that synthesize amino acids that are limiting in the diet of their aphid hosts. Previous studies demonstrated accelerated sequence evolution in Buchnera compared to free-living bacteria, especially for nonsynonymous substitutions. Two mechanisms may explain this acceleration: relaxed purifying selection and increased fixation of slightly deleterious alleles under drift. Here, we test the divergent predictions of these hypotheses for intraspecific polymorphism using Buchnera associated with natural populations of the ragweed aphid, Uroleucon ambrosiae. Contrary to expectations under relaxed selection, U. ambrosiae from across the United States yielded strikingly low sequence diversity at three Buchnera loci (dnaN, trpBC, trpEG), revealing polymorphism three orders of magnitude lower than in enteric bacteria. An excess of nonsynonymous polymorphism and of rare alleles was also observed. Local sampling of additional dnaN sequences revealed similar patterns of polymorphism and no evidence of food plant-associated genetic structure. Aphid mitochondrial sequences further suggested that host bottlenecks and large-scale dispersal may contribute to genetic homogenization of aphids and symbionts. Together, our results support reduced N(e) as a primary cause of accelerated sequence evolution in Buchnera. However, our study cannot rule out the possibility that mechanisms other than bottlenecks also contribute to reduced N(e) at aphid and endosymbiont loci.

Population genomics: genome-wide sampling of insect populations

Modern population genetics underwent a major paradigm shift during the last decade of the 20th century with the discovery that thousands of genes of known function and position in a genome can be analyzed simultaneously in a single individual. The impact of this technology on insect population genetics is potentially profound. Sampling distributions of genetic statistics can now be derived from many individual loci or among many segregating sites within a gene. Inferences regarding random mating, gene flow, effective population sizes, disequilibrium, and relatedness among populations can now be based on patterns of variation at many loci. More importantly, genome-wide sampling enables population geneticists to distinguish effects that act on the whole genome from those that act on individual loci or nucleotides. We introduce the term "population genomics" to describe the process of simultaneous sampling of numerous variable loci within a genome and the inference of locus-specific effects from the sample distributions. The four critical assumptions implicit in the population genomics approach are explained in detail. Studies adopting this paradigm are reviewed, and the steps necessary to complete a population genomics study are outlined.

Mitochondrial genomics flies high

A new pair of papers reports the complete coding sequences for 28 mitochondrial DNAs in Drosophila. By examining the patterns of polymorphism and divergence among functionally distinct classes of synonymous and nonsynonymous nucleotide sites, Bill Ballard provides a comprehensive whole-genome picture of how mtDNA sequence evolution can depart from the strictly neutral and nearly neutral models of molecular evolution.

Inferring parameters of mutation, selection and demography from patterns of synonymous site evolution in Drosophila

Selection acting on codon usage can cause patterns of synonymous evolution to deviate considerably from those expected under neutrality. To investigate the quantitative relationship between parameters of mutation, selection, and demography, and patterns of synonymous site divergence, we have developed a novel combination of population genetic models and likelihood methods of phylogenetic sequence analysis. Comparing 50 orthologous gene pairs from Drosophila melanogaster and D. virilis and 27 from D. melanogaster and D. simulans, we show considerable variation between amino acids and genes in the strength of selection acting on codon usage and find evidence for both long-term and short-term changes in the strength of selection between species. Remarkably, D. melanogaster shows no evidence of current selection on codon usage, while its sister species D. simulans experiences only half the selection pressure for codon usage of their common ancestor. We also find evidence for considerable base asymmetries in the rate of mutation, such that the average synonymous mutation rate is 20-30% higher than in noncoding regions. A Bayesian approach is adopted to investigate how accounting for selection on codon usage influences estimates of the parameters of mutation.

Neutrality tests of conservative-radical amino acid changes in nuclear- and mitochondrially-encoded proteins

The neutralist-selectionist debate should not be viewed as a dichotomy but as a continuum. While the strictly neutral model suggests a neutralist-selectionist dichotomy, the nearly neutral model is a continuous model spanning strict neutrality through weak selection (Ns approximately 1) to deterministic selection (Ns>3). We illustrate these points with polymorphism and divergence data from a sample of 73 genes (31 mitochondrial, 36 nuclear genes from Drosophila, and six Arabidopsis data sets). In an earlier study we used the McDonald-Kreitman (MK) test to show that amino acid replacement polymorphism in animal mitochondrial genes and Arabidopsis genes show a consistent trend toward negative selection, whereas nuclear genes from Drosophila span a range from negative selection, through neutrality, to positive selection. Here we analyze a subset of these genes (13 Drosophila nuclear, ten mitochondrial, and six Arabidopsis nuclear) for polymorphism and divergence of conservative and radical amino acid replacements (a protein-based conservative-radical MK, or pMK, test). The distinct patterns of selection between the different genomes is not apparent with the pMK test. Different definitions of conservative and radical (based on amino acid polarity, volume or charge) give inconsistent results across genes. We suggest that segregating fitness difference between silent and replacement mutations are more visible to selection than are segregating fitness differences between conservative and radical amino acid mutations. New data on the variation among genes with different opportunities for positive and negative selection are as important to the continuum view of the neutralist-selectionist debate as is the distribution of selection coefficients within individual genes.

Post-genomics and the neutral theory: variation and conservation in the tumor necrosis factor-alpha promoter

In the post-genomics era, molecular evolutionary geneticists have come to possess the molecular, statistical, and computational tools for estimating the relative importance of selection and random genetic drift in virtually any gene in almost any organism. We have examined single-nucleotide polymorphisms (SNPs) and nucleotide divergence across a region of approximately 1 kb in the promoter of the human tumor necrosis factor alpha (TNF-alpha) gene. TNF-alpha, which plays an important role in lymphocyte biology and in the pathogenesis of infectious and autoimmune diseases, is tightly regulated at the level of transcription through sequence-specific binding of transcription factors to cognate binding sites in a relatively small region of the 5' non-coding region of the gene. Analysis of the promoter region in 207 human chromosomes revealed nine SNPs, none of which were located in regions known to be important in transcriptional activation. Comparison with one promoter sequence in each of seven species of primates revealed 162 nucleotide sites occupied by a monomorphic nucleotide in the human sample but occupied by a different nucleotide in at least one of the primate sequences (a 'fixed human difference'). The fixed human differences were found outside the regions known to be important in transcriptional activation, and their large number suggests that they might be effectively neutral (Ns<<1). With regard to the human SNPs, although the hypothesis Ns approximately 0 cannot be rejected, the sample configurations suggest that the substitutions might be mildly deleterious. We emphasize the analytical insight to be gained from interspecific comparisons: through the interspecific comparisons, 3.1% of the total sequence information yielded 94.7% of the variable nucleotides. This combined approach, using interspecific comparisons and human polymorphism together with data from functional analyses, provides valuable insights into the evolutionary history and regulation of a key gene in the human immune response.

Extensive amino acid polymorphism at the pgm locus is consistent with adaptive protein evolution in Drosophila melanogaster

PGM plays a central role in the glycolytic pathway at the branch point leading to glycogen metabolism and is highly polymorphic in allozyme studies of many species. We have characterized the nucleotide diversity across the Pgm gene in Drosophila melanogaster and D. simulans to investigate the role that protein polymorphism plays at this crucial metabolic branch point shared with several other enzymes. Although D. melanogaster and D. simulans share common allozyme mobility alleles, we find these allozymes are the result of many different amino acid changes at the nucleotide level. In addition, specific allozyme classes within species contain several amino acid changes, which may explain the absence of latitudinal clines for PGM allozyme alleles, the lack of association of PGM allozymes with the cosmopolitan In(3L)P inversion, and the failure to detect differences between PGM allozymes in functional studies. We find a significant excess of amino acid polymorphisms within D. melanogaster when compared to the complete absence of fixed replacements with D. simulans. There is also strong linkage disequilibrium across the 2354 bp of the Pgm locus, which may be explained by a specific amino acid haplotype that is high in frequency yet contains an excess of singleton polymorphisms. Like G6pd, Pgm shows strong evidence for a branch point enzyme that exhibits adaptive protein evolution.

Molecular population genetics of male accessory gland proteins in Drosophila

Drosophila seminal proteins have an unusually high rate of molecular sequence evolution, suggesting either a high rate of neutral substitution or rapid adaptive evolution. To further quantify patterns of polymorphism and divergence in genes encoding seminal proteins, also called accessory gland proteins (Acp's), we conducted a sequencing survey of 10 Acp genes in samples of Drosophila melanogaster and D. simulans (Acp29AB, Acp32CD, Acp33A, Acp36DE, Acp53Ea, Acp62F, Acp63F, Acp76A, Acp95EF, and Acp98AB). Mean heterozygosity at replacement sites in D. simulans was 0.0074 for Acp genes and 0.0013 for a set of 19 non-Acp genes, and mean melanogaster-simulans divergence at replacement sites was 0.0497 for Acp genes and 0.0107 at non-Acp genes. The elevated divergence of Acp genes is thus accompanied by elevated within-species polymorphism. In addition to the already-reported departures of Acp26A, Acp29AB, and Acp70A from neutrality, our data reject neutrality at Acp29AB and Acp36DE in the direction of excess replacements in interspecific comparisons.

The correlation between intron length and recombination in drosophila. Dynamic equilibrium between mutational and selective forces

Intron length is negatively correlated with recombination in both Drosophila melanogaster and humans. This correlation is not likely to be the result of mutational processes alone: evolutionary analysis of intron length polymorphism in D. melanogaster reveals equivalent ratios of deletion to insertion in regions of high and low recombination. The polymorphism data do reveal, however, an excess of deletions relative to insertions (i.e., a deletion bias), with an overall deletion-to-insertion events ratio of 1.35. We propose two types of selection favoring longer intron lengths. First, the natural mutational bias toward deletion must be opposed by strong selection in very short introns to maintain the minimum intron length needed for the intron splicing reaction. Second, selection will favor insertions in introns that increase recombination between mutations under the influence of selection in adjacent exons. Mutations that increase recombination, even slightly, will be selectively favored because they reduce interference among selected mutations. Interference selection acting on intron length mutations must be very weak, as indicated by frequency spectrum analysis of Drosophila intron length polymorphism, making the equilibrium for intron length sensitive to changes in the recombinational environment and population size. One consequence of this sensitivity is that the advantage of longer introns is expected to decrease inversely with the rate of recombination, thus leading to a negative correlation between intron length and recombination rate. Also in accord with this model, intron length differs between closely related Drosophila species, with the longest variant present more often in D. melanogaster than in D. simulans. We suggest that the study of the proposed dynamic model, taking into account interference among selected sites, might shed light on many aspects of the comparative biology of genome sizes including the C value paradox.

Relationship of codon bias to mRNA concentration and protein length in Saccharomyces cerevisiae

In 1982, Ikemura reported a strikingly unequal usage of different synonymous codons, in five Saccharomyces cerevisiae nuclear genes having high protein levels. To study this trend in detail, we examined data from three independent studies that used oligonucleotide arrays or SAGE to estimate mRNA concentrations for nearly all genes in the genome. Correlation coefficients were calculated for the relationship of mRNA concentration to four commonly used measures of synonymous codon usage bias: the codon adaptation index (CAI), the codon bias index (CBI), the frequency of optimal codons (F(op)), and the effective number of codons (N(c)). mRNA concentration was best approximated as an exponential function of each of these four measures. Of the four, the CAI was the most strongly correlated with mRNA concentration (r(s)=0.62+/-0.01, n=2525, p<10(-17)). When we controlled for CAI, mRNA concentration and protein length were negatively correlated (partial r(s)=-0.23+/-0.01, n=4765, p<10(-17)). This may result from selection to reduce the size of abundant proteins to minimize transcriptional and translational costs. When we controlled for mRNA concentration, protein length and CAI were positively correlated (partial r(s)=0.16+/-0.01, n=4765, p<10(-17)). This may reflect more effective selection in longer genes against missense errors during translation. The correlation coefficients between the mRNA levels of individual genes, as measured by different investigators and methods, were low, in the range r(s)=0.39-0.68.

Molecular evolution of the histone 3 multigene family in the Drosophila melanogaster species subgroup

Molecular evolution of the histone multigene family was studied by cloning and sequencing regions of the histone 3 gene in the Drosophila melanogaster species subgroup. Analysis of the nucleotide substitution pattern showed that in the coding region synonymous changes occurred more frequently to A or T in contrast to the GC-rich base composition, while in the 3' region the nucleotide substitutions were most likely in equilibrium. These results suggested that the base composition at the third codon position of the H3 gene, i.e., codon usage, has been changing to A or T in the Drosophila melanogaster species subgroup.

Contrasting patterns of nonneutral evolution in proteins encoded in nuclear and mitochondrial genomes

We report that patterns of nonneutral DNA sequence evolution among published nuclear and mitochondrially encoded protein-coding loci differ significantly in animals. Whereas an apparent excess of amino acid polymorphism is seen in most (25/31) mitochondrial genes, this pattern is seen in fewer than half (15/36) of the nuclear data sets. This differentiation is even greater among data sets with significant departures from neutrality (14/15 vs. 1/6). Using forward simulations, we examined patterns of nonneutral evolution using parameters chosen to mimic the differences between mitochondrial and nuclear genetics (we varied recombination rate, population size, mutation rate, selective dominance, and intensity of germ line bottleneck). Patterns of evolution were correlated only with effective population size and strength of selection, and no single genetic factor explains the empirical contrast in patterns. We further report that in Arabidopsis thaliana, a highly self-fertilizing plant with effectively low recombination, five of six published nuclear data sets also exhibit an excess of amino acid polymorphism. We suggest that the contrast between nuclear and mitochondrial nonneutrality in animals stems from differences in rates of recombination in conjunction with a distribution of selective effects. If the majority of mutations segregating in populations are deleterious, high linkage may hinder the spread of the occasional beneficial mutation.

Nucleotide polymorphism at the RpII215 gene in Drosophila subobscura. Weak selection on synonymous mutations

Nucleotide variation in an 8.1-kb fragment encompassing the RpII215 gene, which encodes the largest subunit of the RNA polymerase II complex, is analyzed in a sample of 11 chromosomes from a natural population of Drosophila subobscura. No amino acid polymorphism was detected among the 157 segregating sites. The observed numbers of preferred and unpreferred derived synonymous mutations can be explained by neutral mutational processes. In contrast, preferred mutations segregate at significantly higher frequency than unpreferred mutations, suggesting the action of natural selection. The polymorphism to divergence ratio is different for preferred and unpreferred changes, in agreement with their beneficial and deleterious effects on fitness, respectively. Preferred and unpreferred codons are nonrandomly distributed in the RpII215 gene, leading to a heterogeneous distribution of polymorphic to fixed synonymous differences across this coding region. This intragenic variation of the polymorphism/divergence ratio cannot be explained by different patterns of gene expression, mutation, or recombination rates, and therefore it indicates that selection coefficients for synonymous mutations can vary extensively across a coding region. The application of nucleotide composition stationarity tests in coding and flanking noncoding regions, assumed to behave neutrally, allows the detection of the action of natural selection when stationarity holds in the noncoding region.

Effect of inversion polymorphism on the neutral nucleotide variability of linked chromosomal regions in Drosophila

Recombination is a main factor determining nucleotide variability in different regions of the genome. Chromosomal inversions, which are ubiquitous in the genus Drosophila, are known to reduce and redistribute recombination, and thus their specific effect on nucleotide variation may be of major importance as an explanatory factor for levels of DNA variation. Here, we use the coalescent approach to study this effect. First, we develop analytical expressions to predict nucleotide variability in old inversion polymorphisms that have reached mutation-drift-flux equilibrium. The effects on nucleotide variability of a new arrangement appearing in the population and reaching a stable polymorphism are then studied by computer simulation. We show that inversions modulate nucleotide variability in a complex way. The establishment of an inversion polymorphism involves a partial selective sweep that eliminates part of the variability in the population. This is followed by a slow convergence to the equilibrium values. During this convergence, regions close to the breakpoints exhibit much lower variability than central regions. However, at equilibrium, regions close to the breakpoints have higher levels of variability and differentiation between arrangements than regions in the middle of the inverted segment. The implications of these findings for overall variability levels during the evolution of Drosophila species are discussed.

Genetic analysis of bristle loss in hybrids between Drosophila melanogaster and D. simulans provides evidence for divergence of cis-regulatory sequences in the achaete-scute gene complex

The two closely related species of Drosophila, D. melanogaster and D. simulans, display an identical bristle pattern on the notum, but hybrids between the two are lacking a variable number of bristles. We show that the loss is temperature-dependent and provide evidence for two periods of temperature sensitivity. A first period of heat sensitivity occurs during larval development and corresponds to the time when the prepattern of expression of genes whose products activate achaete-scute in the proneural clusters preceding bristle precursor formation is established. A second period of cold sensitivity corresponds to the time of emergence of the bristle precursor cells and the maintenance of their neural fate, a process requiring high levels of Achaete-Scute. Expression of achaete-scute at these two critical periods depends on cis-regulatory elements of the achaete-scute complex (AS-C). The differences between males, which have only one copy of the X-linked AS-C from D. simulans, and females, which have copies from both parental species, are compared, together with the effects of crossing in different rearrangements of the D. melanogaster AS-C that delete regulatory and/or coding sequences. We provide evidence that bristle loss in the hybrids may result from a decrease in the level of transcription at the AS-C and argue that interaction between trans-acting factors and cis-regulatory elements within the AS-C has diverged between the two species.

Adaptive evolution of relish, a Drosophila NF-kappaB/IkappaB protein

NF-kappaB and IkappaB proteins have central roles in regulation of inflammation and innate immunity in mammals. Homologues of these proteins also play an important role in regulation of the Drosophila immune response. Here we present a molecular population genetic analysis of Relish, a Drosophila NF-kappaB/IkappaB protein, in Drosophila simulans and D. melanogaster. We find strong evidence for adaptive protein evolution in D. simulans, but not in D. melanogaster. The adaptive evolution appears to be restricted to the IkappaB domain. A possible explanation for these results is that Relish is a site of evolutionary conflict between flies and their microbial pathogens.

Large number of replacement polymorphisms in rapidly evolving genes of Drosophila. Implications for genome-wide surveys of DNA polymorphism

We present a survey of nucleotide polymorphism of three novel, rapidly evolving genes in populations of Drosophila melanogaster and D. simulans. Levels of silent polymorphism are comparable to other loci, but the number of replacement polymorphisms is higher than that in most other genes surveyed in D. melanogaster and D. simulans. Tests of neutrality fail to reject neutral evolution with one exception. This concerns a gene located in a region of high recombination rate in D. simulans and in a region of low recombination rate in D. melanogaster, due to an inversion. In the latter case it shows a very low number of polymorphisms, presumably due to selective sweeps in the region. Patterns of nucleotide polymorphism suggest that most substitutions are neutral or nearly neutral and that weak (positive and purifying) selection plays a significant role in the evolution of these genes. At all three loci, purifying selection of slightly deleterious replacement mutations appears to be more efficient in D. simulans than in D. melanogaster, presumably due to different effective population sizes. Our analysis suggests that current knowledge about genome-wide patterns of nucleotide polymorphism is far from complete with respect to the types and range of nucleotide substitutions and that further analysis of differences between local populations will be required to understand the forces more completely. We note that rapidly diverging and nearly neutrally evolving genes cannot be expected only in the genome of Drosophila, but are likely to occur in large numbers also in other organisms and that their function and evolution are little understood so far.

Local recombination and mutation effects on molecular evolution in Drosophila

I studied the cause of the significant difference in the synonymous-substitution pattern found in the achaete-scute complex genes in two Drosophila lineages, higher codon bias in Drosophila yakuba, and lower bias in D. melanogaster. Besides these genes, the functionally unrelated yellow gene showed the same substitution pattern, suggesting a region-dependent phenomenon in the X-chromosome telomere. Because the numbers of A/T --> G/C substitutions were not significantly different from those of G/C --> A/T in the yellow noncoding regions of these species, a AT/GC mutational bias could not completely account for the synonymous-substitution biases. In contrast, we did find an approximately 14-fold difference in recombination rates in the X-chromosome telomere regions between the two species, suggesting that the reduction of recombination rates in this region resulted in the reduction of the efficacy of selection in D. melanogaster. In addition, the D. orena yellow showed a 5% increase in the G + C content at silent sites in the coding and noncoding regions since the divergence from D. erecta. This pattern was significantly different from those at the orena Adh and Amy loci. These results suggest that local changes in recombination rates and mutational pressures are contributing to the irregular synonymous-substitution patterns in Drosophila.