Protein variation in Drosophila simulans, and comparison of genes from centromeric versus noncentromeric regions of chromosome 3

Molecular population genetics of female-expressed mating-induced serine proteases in Drosophila melanogaster

Population genetic analyses have shown that directional selection causes amino acid substitution in several seminal fluid proteins (Acps) and that in general, Acps tend to diverge rapidly. If rapid, adaptive divergence of such male reproduction-related genes is driven by sexual conflict, we might also expect to observe rapid, adaptive evolution in female reproduction-related genes, especially those mediating conflicts between the sexes. Female expressed genes differentially expressed shortly after mating were recently identified using whole genome expression micro-arrays. Such genes may play roles in storing sperm and mediating effects of seminal fluid proteins. Here, we report the results of a molecular population genetic survey from five female reproductive tract expressed serine proteases that show increased transcription shortly after mating. These genes are evolving rapidly, in some cases under directional selection, consistent with models of conflict.

Novel genes derived from noncoding DNA in Drosophila melanogaster are frequently X-linked and exhibit testis-biased expression

Descriptions of recently evolved genes suggest several mechanisms of origin including exon shuffling, gene fission/fusion, retrotransposition, duplication-divergence, and lateral gene transfer, all of which involve recruitment of preexisting genes or genetic elements into new function. The importance of noncoding DNA in the origin of novel genes remains an open question. We used the well annotated genome of the genetic model system Drosophila melanogaster and genome sequences of related species to carry out a whole-genome search for new D. melanogaster genes that are derived from noncoding DNA. Here, we describe five such genes, four of which are X-linked. Our RT-PCR experiments show that all five putative novel genes are expressed predominantly in testes. These data support the idea that these novel genes are derived from ancestral noncoding sequence and that new, favored genes are likely to invade populations under selective pressures relating to male reproduction.

Estimating the genomewide rate of adaptive protein evolution in Drosophila

When polymorphism and divergence data are available for multiple loci, extended forms of the McDonald-Kreitman test can be used to estimate the average proportion of the amino acid divergence due to adaptive evolution--a statistic denoted alpha. But such tests are subject to many biases. Most serious is the possibility that high estimates of alpha reflect demographic changes rather than adaptive substitution. Testing for between-locus variation in alpha is one possible way of distinguishing between demography and selection. However, such tests have yielded contradictory results, and their efficacy is unclear. Estimates of alpha from the same model organisms have also varied widely. This study clarifies the reasons for these discrepancies, identifying several method-specific biases in widely used estimators and assessing the power of the methods. As part of this process, a new maximum-likelihood estimator is introduced. This estimator is applied to a newly compiled data set of 115 genes from Drosophila simulans, each with each orthologs from D. melanogaster and D. yakuba. In this way, it is estimated that alpha approximately 0.4+/-0.1, a value that does not vary substantially between different loci or over different periods of divergence. The implications of these results are discussed.

Recently evolved genes identified from Drosophila yakuba and D. erecta accessory gland expressed sequence tags

The fraction of the genome associated with male reproduction in Drosophila may be unusually dynamic. For example, male reproduction-related genes show higher-than-average rates of protein divergence and gene expression evolution compared to most Drosophila genes. Drosophila male reproduction may also be enriched for novel genetic functions. Our earlier work, based on accessory gland protein genes (Acp's) in D. simulans and D. melanogaster, suggested that the melanogaster subgroup Acp's may be lost and/or gained on a relatively rapid timescale. Here we investigate this possibility more thoroughly through description of the accessory gland transcriptome in two melanogaster subgroup species, D. yakuba and D. erecta. A genomic analysis of previously unknown genes isolated from cDNA libraries of these species revealed several cases of genes present in one or both species, yet absent from ingroup and outgroup species. We found no evidence that these novel genes are attributable primarily to duplication and divergence, which suggests the possibility that Acp's or other genes coding for small proteins may originate from ancestrally noncoding DNA.

Molecular evolution in the Drosophila melanogaster species subgroup: frequent parameter fluctuations on the timescale of molecular divergence

Although mutation, genetic drift, and natural selection are well established as determinants of genome evolution, the importance (frequency and magnitude) of parameter fluctuations in molecular evolution is less understood. DNA sequence comparisons among closely related species allow specific substitutions to be assigned to lineages on a phylogenetic tree. In this study, we compare patterns of codon usage and protein evolution in 22 genes (>11,000 codons) among Drosophila melanogaster and five relatives within the D. melanogaster subgroup. We assign changes to eight lineages using a maximum-likelihood approach to infer ancestral states. Uncertainty in ancestral reconstructions is taken into account, at least to some extent, by weighting reconstructions by their posterior probabilities. Four of the eight lineages show potentially genomewide departures from equilibrium synonymous codon usage; three are decreasing and one is increasing in major codon usage. Several of these departures are consistent with lineage-specific changes in selection intensity (selection coefficients scaled to effective population size) at silent sites. Intron base composition and rates and patterns of protein evolution are also heterogeneous among these lineages. The magnitude of forces governing silent, intron, and protein evolution appears to have varied frequently, and in a lineage-specific manner, within the D. melanogaster subgroup.

Rapid evolution and gene-specific patterns of selection for three genes of spermatogenesis in Drosophila

Hybrid males resulting from crosses between closely related species of Drosophila are sterile. The F1 hybrid sterility phenotype is mainly due to defects occurring during late stages of development that relate to sperm individualization, and so genes controlling sperm development may have been subjected to selective diversification between species. It is also possible that genes of spermatogenesis experience selective constraints given their role in a developmental pathway. We analyzed the molecular evolution of three genes playing a role during the sperm developmental pathway in Drosophila at an early (bam), a mid (aly), and a late (dj) stage. The complete coding region of these genes was sequenced in different strains of Drosophila melanogaster and Drosophila simulans. All three genes showed rapid divergence between species, with larger numbers of nonsynonymous to synonymous differences between species than polymorphisms. Although this could be interpreted as evidence for positive selection at all three genes, formal tests of selection do not support such a conclusion. Departures from neutrality were detected only for dj and bam but not aly. The role played by selection is unique and determined by gene-specific characteristics rather than site of expression. In dj, the departure was due to a high proportion of neutral synonymous polymorphisms in D. simulans, and there was evidence of purifying selection maintaining a high lysine amino acid protein content that is characteristic of other DNA-binding proteins. The earliest spermatogenesis gene surveyed, which plays a role in both male and female gametogenesis, was bam, and its significant departure from neutrality was due to an excess of nonsynonymous substitutions between species. Bam is degraded at the end of mitosis, and rapid evolutionary changes among species might be a characteristic shared with other degradable transient proteins. However, the large number of nonsynonymous changes between D. melanogaster and D. simulans and a phylogenetic comparative analysis among species confirms evidence of positive selection driving the evolution of Bam and suggests an yet unknown germ cell line developmental adaptive change between these two species.

Elevated polymorphism and divergence in the class C scavenger receptors of Drosophila melanogaster and D. simulans

Scavenger receptor proteins are involved in the cellular internalization of a broad variety of foreign material, including pathogenic bacteria during phagocytosis. I find here that nonsynonymous divergence in three class C scavenger receptors (Sr-C's) between Drosophila melanogaster and D. simulans and between each of these species and D. yakuba is approximately four times the typical genome average. These genes also exhibit unusually high levels of segregating nonsynonymous polymorphism in D. melanogaster and D. simulans populations. A fourth Sr-C is comparatively conserved. McDonald-Kreitman tests reveal a significant excess of replacement fixations between D. melanogaster and D. simulans in the Sr-C's, but tests of polymorphic site frequency spectra do not support models of directional selection. It is possible that the molecular functions of SR-C proteins are sufficiently robust to allow exceptionally high amino acid substitution rates without compromising organismal fitness. Alternatively, SR-Cs may evolve under diversifying selection, perhaps as a result of pressure from pathogens. Interestingly, Sr-CIII and Sr-CIV are polymorphic for premature stop codons. Sr-CIV is also polymorphic for an in-frame 101-codon deletion and for the absence of one intron.

Patterns of Microsatellite Variability Among X Chromosomes and Autosomes Indicate a High Frequency of Beneficial Mutations in Non-African D. simulans

We analyzed microsatellite variability at 42 X-linked and 39 autosomal loci from African and European populations of Drosophila simulans. The African D. simulans harbored significantly more microsatellite variability than the European flies. In the European population, X-linked polymorphism was more reduced than autosomal variation, whereas there was no significant difference between chromosomes in the African population. Previous studies also observed a similar pattern but failed to distinguish between a demographic event and a selection scenario. We performed extensive computer simulations using a wide range of demographic scenarios to distinguish between the two hypotheses. Approximate summary likelihood estimates differed dramatically among X chromosomes and autosomes. Furthermore, our experimental data showed a surplus of X-linked microsatellites with a significantly reduced variability in non-African D. simulans. We conclude that our data are not compatible with a neutral scenario. Thus, the reduced variability at X-linked loci is most likely caused by selective sweeps associated with the out-of-Africa habitat expansion of D. simulans.

Functional divergence caused by ancient positive selection of a Drosophila hybrid incompatibility locus

Interspecific hybrid lethality and sterility are a consequence of divergent evolution between species and serve to maintain the discrete identities of species. The evolution of hybrid incompatibilities has been described in widely accepted models by Dobzhansky and Muller where lineage-specific functional divergence is the essential characteristic of hybrid incompatibility genes. Experimentally tractable models are required to identify and test candidate hybrid incompatibility genes. Several Drosophila melanogaster genes involved in hybrid incompatibility have been identified but none has yet been shown to have functionally diverged in accordance with the Dobzhansky-Muller model. By introducing transgenic copies of the X-linked Hybrid male rescue (Hmr) gene into D. melanogaster from its sibling species D. simulans and D. mauritiana, we demonstrate that Hmr has functionally diverged to cause F1 hybrid incompatibility between these species. Consistent with the Dobzhansky-Muller model, we find that Hmr has diverged extensively in the D. melanogaster lineage, but we also find extensive divergence in the sibling-species lineage. Together, these findings implicate over 13% of the amino acids encoded by Hmr as candidates for causing hybrid incompatibility. The exceptional level of divergence at Hmr cannot be explained by neutral processes because we use phylogenetic methods and population genetic analyses to show that the elevated amino-acid divergence in both lineages is due to positive selection in the distant past—at least one million generations ago. Our findings suggest that multiple substitutions driven by natural selection may be a general phenomenon required to generate hybrid incompatibility alleles.

Transgenic experiments show that the HMR gene has functionally diverged in Drosophila melanogaster and its sibling species and causes the death of hybrid offspring in interspecific crosses

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.

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.

A rapidly evolving MYB-related protein causes species isolation in Drosophila

Matings among different species of animals or plants often result in sterile or lethal hybrids. Identifying the evolutionary forces that create hybrid incompatibility alleles is fundamental to understanding the process of speciation, but very few such alleles have been identified, particularly in model organisms that are amenable to experimental manipulation. We report here the cloning of the first, to our knowledge, Drosophila melanogaster gene involved in hybrid incompatibilities, Hybrid male rescue (Hmr). Hmr causes lethality and female sterility in hybrids among D. melanogaster and its sibling species. We have found that Hmr encodes a protein with homology to a family of MYB-related DNA-binding transcriptional regulators. The HMR protein has evolved both amino acid substitutions and insertions and deletions at an extraordinarily high rate between D. melanogaster and its sibling species, including in its predicted DNA-binding domain. Our results suggest that hybrid lethality may result from disruptions in gene regulation, and we also propose that rapid evolution may be a hallmark of speciation genes in general.