A Genomic Comparison of Faster-Sex, Faster-X, and Faster-Male Evolution Between Drosophila melanogaster and Drosophila pseudoobscura

A genomic comparison of Drosophila melanogaster and Drosophila pseudoobscura provides a unique opportunity to investigate factors involved in sequence divergence. The chromosomal arrangements of these species include an autosomal segment in D. melanogaster which is homologous to part of the X chromosome in D. pseudoobscura. Using orthologues to calculate rates of nonsynonymous (d_N) substitutions, we found genes on the X chromosome to be significantly more diverged than those on the autosomes, but it is not true for segment 3L-XR which is autosomal in D. melanogaster (3L) and X-linked in D. pseudoobscura (XR). We also found that the median d_N values for genes having reproductive functions in either the male, the female, or both sexes are higher than those for sequences without reproductive function and even higher for sequences involved in male-specific function. These estimates of divergence for male sex-related sequences are most likely underestimates, as the very rapidly evolving reproductive genes would tend to lose homology sooner and thus not be included in the comparison of orthologues. We also noticed a high proportion of male reproductive genes among the othologous genes with the highest rates of d_N. Reproductive genes with and without an orthologue in D. pseudoobscura were compared among D. melanogaster, D. simulans, and D. yakuba and it was found that there were in fact higher rates of divergence in the group without a D. pseudoobscura orthologue. These results, from widely separated taxa, bolster the thesis that sexual system genes experience accelerated rates of change in comparison to nonsexual genes in evolution and speciation. [Reviewing Editor: Dr. Willie J. Swanson]     

Historical effective size and the level of genetic diversity in Drosophila melanogaster and Drosophila pseudoobscura

We report the results of a sequential gel electrophoretic study of protein variation in Drosophila melanogaster and its comparison with D. pseudoobscura. The number of alleles and mean heterozygosity were lower in D. melanogaster than in D. pseudoobscura. On the other hand, geographical populations of Drosophila melanogaster have been shown to be much more differentiated than those of D. pseudoobscura. The results suggest that in D. melanogaster low-frequency alleles have been lost during the colonization process and that major alleles have become differentiated among populations. Population bottlenecks, due to various causes, appear to have played a significant role in the shaping of genetic variation in natural populations of many species. It is proposed that a comparison of genetic variation at homologous gene loci between related species can bring out effects of historical bottlenecks and provide an alternative approach for analyzing causes of genetic variation in natural populations.We thank the Natural Science and Engineering Research Council of Canada for financial support (Grant A0235 to R.S.S.).     

Strong Purifying Selection on the Odysseus Gene in Two Clades of Sibling Species of the Drosophila montium Species Subgroup

The Odysseus (OdsH) gene was duplicated from its ancestral neuron-expressed gene, unc-4, and then evolved very rapidly under strong positive Darwinian selection as a speciation gene causing hybrid-male sterility between closely related species of the Drosophila simulans clade. Has OdsH also experienced similar positive selection between Drosophila sibling species other than those of the simulans clade? We cloned and sequenced OdsH and unc-4 from two clades of the Drosophila montium species subgroup, the Drosophila lini and the Drosophila kikkawai clades. The ratios of Ka/Ks for OdsH were remarkably low between sibling species of these two clades, suggesting that OdsH has been subjected to strong purifying selection in these two clades.     

Evidence of Positive Darwinian Selection in Putative Meningococcal Vaccine Antigens

Meningococcal meningitidis is a life-threatening disease. In Europe and the United States the majority of cases are caused by virulent meningococcal strains belonging to serogroup B. Presently there is no effective vaccine against serogroup B strains, as traditional vaccine antigens such as polysaccharide capsules are unusable as they lead to autoimmunity. The year 2000 saw the publication of the complete genome of Neisseria meningitidis MC58, a virulent serogroup B bacterium. Working in conjunction with the sequencing project, researchers endeavored to locate highly conserved membrane-associated proteins that elicit an immune response. It is hoped that these proteins will provide a basis for novel vaccines against serogroup B strains. A number of potential vaccine antigens have been located and are presently in phase I clinical trials. Recently many reports pertaining to the evidence of positive Darwinian selection in membrane proteins of pathogens have been reported. This study utilized in silico methods to test for evidence of historical positive Darwinian selection in seven such vaccine candidates. We found that two of these proteins show signatures of adaptive evolution, while the remaining proteins show evidence of strong purifying selection. This has significant implications for the design of a vaccine against serogroup B strains, as it has been shown that vaccines that target epitopes that are under strong purifying selection are better than those that target variable epitopes.[Reviewing Editor: Rasmus Nielsen]     

Molecular Evolution and Positive Selection of the Symbiotic Gene NORK in Medicago truncatula

Understanding the selective constraints of partner specificity in mutually beneficial symbiosis is a significant, yet largely unexplored, prospect of evolutionary biology. These selective constraints can be explored through the study of nucleotide polymorphism at loci controlling specificity. The membrane-anchored receptor NORK (nodulation receptor kinase) of the legume Medicago truncatula controls early steps of root infection by two symbiotic microorganisms: nitrogen-fixing bacteria (rhizobia) and endomycorrhizal fungi (Glomales). We analyzed the diversity of the gene NORK by sequencing 4 kilobases in 28 inbred lines sampled from natural populations. We detected 33 polymorphic sites with only one nonsynonymous change. Analysis based on Tajima’s D and Fay and Wu’s H summary statistics revealed no departure from the neutral model. We analyzed divergence using sequences from the closely related species M. coerulea. The McDonald-Kreitman test indicated a significant excess of nonsynonymous changes contributing to this divergence. Furthermore, maximum-likelihood analysis of a molecular phylogeny of a few legume species indicated that a number of amino acid sites, likely located in the receptor domain of the protein, evolved under the regime of positive selection. Further research should focus on the rate and direction of molecular coevolution between microorganisms’ signaling molecules and legumes’ receptors. [Reviewing Editor: Dr. Deborah Charlesworth] Sequence data were deposited in the GenBank database under accession nos. AY676428 to AY676457 and AJ884582.     

Sequence and evolution of the Drosophila pseudoobscura glycerol-3-phosphate dehydrogenase locus

The Gpdh genomic region has been cloned and sequenced in Drosophila pseudoobscura. A total of 6.8 kb of sequence was obtained, encompassing all eight exons of the gene. The exons have been aligned with the sequence from D. melanogaster, and the rates of synonymous and nonsynonymous substitution have been compared to those of other genes sequenced in these two species. Gpdh has the lowest rate of nonsynonymous substitution yet seen in genes sequenced in both D. pseudoobscura and D. melanogaster. No insertion/deletion events were observed, and the overall architecture of the gene (i.e., intron sites, etc.) is conserved. An interesting amino acid reversal was noted between the D. melanogaster Fast allele and the D. pseudoobscura gene.     

A New Test for Selection Applied to Codon Usage in Drosophila simulans and D. mauritiana

In many organisms, synonymous codon usage is biased by a history of natural selection. However, codon bias, itself, does not indicate that selection is ongoing; it may be a vestige of past selection. Simple statistical tests have been devised to infer ongoing selection on codon usage by comparing the derived state frequency spectra at polymorphic sites segregating either derived preferred codons or derived unpreferred codons; if selection is effective, the frequency of derived states should be higher in the former. We propose a new test that uses the inferred degree of preference, essentially calculating the correlation of derived state frequency and the difference in preference between the derived and the ancestral states; the correlation should be positive if selection is effective. When implementing the test, derived and ancestral states can be assigned by parsimony or on the basis of relative probability. In either case, statistical significance is estimated by a simple permutation test. We explored the statistical power of the test by sampling polymorphism data from 14 loci in 16 strains of D. simulans, finding that the test retains 80% power even when quite a few of the data are discarded. The power of the test likely reflects better use of multiple features of the data, combining population frequencies of polymorphic variants and quantitative estimates of codon preferences. We also applied this novel test to 14 newly sequenced loci in five strains of D. mauritiana, showing for the first time ongoing selection on codon usage in this species.     

Esterase Patterns and Phylogenetic Relationships of Drosophila Species in the Saltans Subgroup (Saltans Group)

The esterase patterns of sixteen strains from four species in the saltans subgroup were analyzed using polyacrylamide gel electrophoresis. Thirty-four esterase bands were detected. By using and naphthyl acetates as substrates, they were classified in 18 -esterases (they hydrolyse the -naphtyl substrate), 15 -esterases (they hydrolyse the -naphtyl substrate) and 1 /-esterase (it hydrolyses the and -naphtyl substrates). Among the -esterases, three were detected exclusively in males. Malathion, Eserine and pCMB were used as inhibitors in order to characterize biochemically the esterases. The results indicated the presence of cholinesterases, carboxylesterases and acetylesterases. The degree of mobility of the bands in the gels, their specificity to and naphthyl acetates and the results of the inhibition tests allowed us to recognize tentatively nine genetic loci. Phylogenetic relationships among species inferred on the basis of the esterase patterns by PAUP 4.0b8, with neighbor-joining search and a bootstrap analysis showed that, although the four species are closely related, D. septentriosaltans, D. saltans and D. austrosaltans are closer to each other than to D. prosaltans. These results showed to be consistent with phylogenetic relationships previously inferred from inversion polymorphism.     

Variability on the dot chromosome in the Drosophila simulans clade

A recent study suggested that recent nuclear gene introgression between Drosophila simulans and D. mauritiana may have obscured efforts to estimate the phylogeny of the species of the D. simulans clade, which includes these two species and D. sechellia. Here, we report sequence variation of an intron of the eyeless gene in this species group. This gene should introgress freely between these species because it is not linked to any known barriers to gene exchange. We have also reevaluated levels of sequence divergence among species in this clade, noting differences between loci in regions of low recombination (as in all chromosome 4 loci) relative to other loci. Overall, none of the data analyzed were consistent with recent introgression exclusively between D. simulans and D. mauritiana.     

Evolutionary experiments on mate recognition in the Drosophila serrata species complex

It is becoming increasingly apparent that at least some aspects of the evolution of mate recognition may be amenable to manipulation in evolutionary experiments. Quantitative genetic analyses that focus on the genetic consequences of evolutionary processes that result in mate recognition evolution may eventually provide an understanding of the genetic basis of the process of speciation. We review a series of experiments that have attempted to determine the genetic basis of the response to natural and sexual selection on mate recognition in the Drosophila serrata species complex. The genetic basis of mate recognition has been investigated at three levels: (1) between the species of D. serrata and D. birchii using interspecific hybrids, (2) between populations of D. serrata that are sympatric and allopatric with respect to D. birchii, and (3) within populations of D. serrata. These experiments suggest that it may be possible to use evolutionary experiments to observe important events such as the reinforcement of mate recognition, or the generation of the genetic associations that are central to many sexual selection models.     

Reduced Selection for Codon Usage Bias in <Emphasis Type="Italic">Drosophila miranda</Emphasis>

Biased codon usage in many species results from a balance among mutation, weak selection, and genetic drift. Here I show that selection to maintain biased codon usage is reduced in Drosophila miranda relative to its ancestor. Analyses of mutation patterns in noncoding DNA suggest that the extent of this reduction cannot be explained by changes in mutation bias or by biased gene conversion. Low levels of variability in D. miranda relative to its sibling species, D. pseudoobscura, suggest that it has a much smaller effective population size. Reduced codon usage bias in D. miranda may thus result from the reduced efficacy of selection against newly arising mutations to unpreferred codons. [Reviewing Editor: Dr. Richard Kliman]     

Positive Selection on HIV Accessory Proteins and the Analysis of Molecular Adaptation After Interspecies Transmission

Studies examining positive selection on accessory proteins of HIV are rare, although these proteins play an important role in pathogenesis in vivo. Moreover, despite the biological relevance of analyses of molecular adaptation after viral transmission between species, the issue is still poorly studied. Here we present evidence that accessory proteins are subjected to positive selective forces exclusively in HIV. This scenario suggests that accessory protein genes are under adaptive evolution in HIV clades, while in SIVcpz such a phenomenon could not be detected. As a result, we show that comparative studies are critical to carry out functional investigation of positively selected protein sites, as they might help to achieve a better comprehension of the biology of HIV pathogenesis.     

Clonal analysis in hybrids between Drosophila melanogaster and Drosophila simulans

We have analysed the viability of cellular clones induced by mitotic recombination in Drosophila melanogaster/D. simulans hybrid females during larval growth. These clones contain a portion of either melanogaster or simulans genomes in homozygosity. Analysis has been carried out for the X and the second chromosomes, as well as for the 3L chromosome arm. Clones were not found in certain structures, and in others they appeared in a very low frequency. Only in abdominal tergites was a significant number of clones observed, although their frequency was lower than in melanogaster abdomens. The bigger the portion of the genome that is homozygous, the less viable is the recombinant melano-gaster/simulans hybrid clone. The few clones that appeared may represent cases in which mitotic recombination took place in distal chromosome intervals, so that the clones contained a small portion of either melanogaster or simulans chromosomes in homozygosity. Moreover, Lhr, a gene of D. simulans that suppresses the lethality of male and female melanogaster/simulans hybrids, does not suppress the lethality of the recombinant melanogaster/simulans clones. Thus, it appears that there is not just a single gene, but at least one per tested chromosome arm (and maybe more) that cause hybrid lethality. Therefore, the two species, D. melanogaster and D. simulans, have diverged to such a degree that the absence of part of the genome of one species cannot be substituted by the corresponding part of the genome of the other, probably due to the formation of co-adapted gene complexes in both species following their divergent evolution after speciation. The disruption of those coadapted gene complexes would cause the lethality of the recombinant hybrid clones.     

Evidence for Positive Selection on a Sexual Reproduction Gene in the Diatom Genus Thalassiosira (Bacillariophyta)

Single likelihood ancestor counting (SLAC), fixed effects likelihood (FEL), and several random effects likelihood (REL) methods were utilized to identify positively and negatively selected sites in sexually induced gene 1 (Sig1) of four different Thalassiosira species. The SLAC analysis did not find any sites affected by positive selection but suggested 13 sites influenced by negative selection. The SLAC approach may be too conservative because of low sequence divergence. The FEL and REL analyses revealed over 60 negatively selected sites and two positively selected sites that were unique to each method. The REL method may not be able to reliably identify individual sites under selection when applied to short sequences with low divergence. Instead, we proposed a new alignment-wide test for adaptive evolution based on codon models with variation in synonymous and nonsynonymous substitution rates among sites and found evidence for diversifying evolution without relying on site-by-site testing. The performance of the FEL and REL approaches was evaluated by subjecting the tests to a type I error rate simulation analysis, using the specific characteristics of the Sig1 data set. Simulation results indicated that the FEL test had reasonable Type I errors, while REL might have been too liberal, suggesting that the two positively selected sites identified by FEL (codons 94 and 174) are not likely to be false positives. The evolution of these codon sites, one of which is located in functional domain II, appears to be associated with divergence among the three major Thalassiosira lineages. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users. [Reviewing Editor: Dr. Martin Kreitman]     

Reproductive fitness and artificial selection in animal breeding: culling on fitness prevents a decline in reproductive fitness in lines of Drosophila melanogaster selected for increased inebriation time

Summary The maintenance of reproductive fitness in lines subjected to artificial selection is one of the major problems in animal breeding. The decline in reproductive performance has neither been predictable from heritabilities and genetic correlations, nor have conventional selection indices been adequate to avoid the problem. Gowe (1983) has suggested that the heritabilities of reproductive traits are non-linear, with heritabilities being higher on the lower fitness side. Consequently, he has predicted that culling on reproductive fitness in artificial selection lines will be effective in preventing the usual declines in fitness. An experimental evaluation of Gowe's prediction has been carried out by comparing fitnesses of replicated lines of three treatments: selection for increased inebriation time without culling on fitness (HO), selection for inebriation time with culling of 20% (4/20) of selected females on reproductive fitness (HS), and unselected controls (C). Response to selection for inebriation time in the two selection treatments was similar. After 25 generations, the competitive index, a measure of reproductive fitness, was significantly lower in the HO treatment than the HS treatment, while the HS treatment did not differ from the control lines or the base population. These results demonstrate for the first time that culling on reproductive fitness in selection lines can be used to prevent the usual decline in reproductive performance.     

Selection against homozygotes in the ADH polymorphic system of Drosophila melanogaster associated with the lethal factor l(2) Stm

In the natural populations ⁺Tüb, ⁺Prov, and ⁺Rov, similar Adh ^F allele frequencies occur (q _F=0.11, 0.18, and 0.08, respectively). However, there is a discrepancy in that the Adh ^F allele in ⁺Tüb is closely linked to the lethal factor 1(2)Stm, which reduces relative fitness of the F phenotype to zero. In spite of this, polymorphism is maintained also in ⁺Tüb, because the heterozygotes are superior to the homozygous S type (relative fitness=0.88). Under laboratory culture conditions, in ⁺Tüb the relative fitness of the S genotype further decreases to 0.6. After outcrossing the lethal factor, relative fitnesses for S, FS, and F become 0.6, 1, and 0.48, respectively, implying that fitness for S remains the same. Relative values for S, FS, and F in ⁺Prov, not affected by the lethal factor, are calculated by the maximum average fitness method to be 1, 1.2, and 0.2 under the assumption that heterozygous FS are similarly superior to S as in the natural ⁺Tüb population and all allele frequencies found are stable equilibrium values.     

Genomic impacts of chromosomal inversions in parapatric Drosophila species

Chromosomal inversions impact genetic variation and facilitate speciation in part by reducing recombination in heterokaryotypes. We generated multiple whole-genome shotgun sequences of the parapatric species pair Drosophila pseudoobscura and Drosophila persimilis and their sympatric outgroup (Drosophila miranda) and compared the average pairwise differences for neutral sites within, just outside and far outside of the three large inversions. Divergence between D. pseudoobscura and D. persimilis is high inside the inversions and in the suppressed recombination regions extending 2.5 Mb outside of inversions, but significantly lower in collinear regions further from the inversions. We observe little evidence of decreased divergence predicted to exist in the centre of inversions, suggesting that gene flow through double crossovers or gene conversion is limited within the inversion, or selection is acting within the inversion to maintain divergence in the face of gene flow. In combination with past studies, we provide evidence that inversions in this system maintain areas of high divergence in the face of hybridization, and have done so for a substantial period of time. The left arm of the X chromosome and chromosome 2 inversions appear to have arisen in the lineage leading to D. persimilis approximately 2 Ma, near the time of the split of D. persimilis-D. pseudoobscura-D. miranda, but likely fixed within D. persimilis much more recently, as diversity within D. persimilis is substantially reduced inside and near these two inversions. We also hypothesize that the inversions in D. persimilis may provide an empirical example of the 'mixed geographical mode' theory of inversion origin and fixation, whereby allopatry and secondary contact both play a role.     

Indirect thermal selection in Drosophila melanogaster and adaptive consequences

Summary Short-term indirect selection in Drosophila melanogaster for heat-sensitivity and heat resistance resulted in two strains, one heat sensitive and another heat resistant, and correlated responses were found for the rate of heat shock protein synthesis, behavioral patterns (asymmetrical sexual isolation) and fitness components (fecundity, fertility, viability, developmental time), as well as for several enzyme activities (MDH, G-6-PDH, ADH, ACHE). These responses associated with temperature selection may reflect the effects of differential inbreeding depression caused by homozygosity of temperature sensitive mutations with different pleiotropic effects. Selection even of a very short duration can induce significant adaptive and evolutionary changes.