Molecular population genetics and evolution of Drosophila meiosis genes

While many functional elements of the meiotic process are well characterized in model organisms, the genetic basis of most of the natural phenotypic variation observed in meiotic pathways has not been determined. To begin to address this issue, we characterized patterns of polymorphism and divergence in the protein-coding regions of 33 genes across 31 lines of Drosophila melanogaster and 6 lines of Drosophila simulans. We sequenced genes known to be involved in chromosome segregation, recombination, DNA repair, and related heterochromatin binding. As expected, we found several of the genes to be highly conserved, consistent with purifying selection. However, a subset of genes showed patterns of polymorphism and divergence typical of other types of natural selection. Moreover, several intriguing differences between the two Drosophila lineages were evident: along the D. simulans lineage we consistently found evidence of adaptive protein evolution, whereas along the D. melanogaster lineage several loci exhibited patterns consistent with the maintenance of protein variation.     

Molecular evolution of the ligands of the insulin-signaling pathway: dilp genes in the genus Drosophila

Drosophila melanogaster, unlike mammals, has seven insulin-like peptides (DILPS). In Drosophila, all seven genes (dilp1-7) are single copy in the 12 species studied, except for D. grimshawi with two tandem copies of dilp2. Our comparative analysis revealed that genes dilp1-dilp7 exhibit differential functional constraint, which is indicative of some functional divergence. Species of the subgenera Sophophora and Drosophila differ in some traits likely affected by the insulin-signaling pathway, such as adult body size. It is in the branch connecting the two subgenera that we found the footprint left by positive selection driving nonsynonymous changes at some dilp1 codons to fixation. Finally, the similar rate at which the two dilp2 copies of D. grimshawi have evolved since their duplication and the presence of a putative regulatory region highly conserved between the two paralogs would suggest that both copies were preserved either because of subfunctionalization or dose dependency rather than by the neofunctionalization of one of the two copies.     

Population genomics: whole-genome analysis of polymorphism and divergence in Drosophila simulans

The population genetic perspective is that the processes shaping genomic variation can be revealed only through simultaneous investigation of sequence polymorphism and divergence within and between closely related species. Here we present a population genetic analysis of Drosophila simulans based on whole-genome shotgun sequencing of multiple inbred lines and comparison of the resulting data to genome assemblies of the closely related species, D. melanogaster and D. yakuba. We discovered previously unknown, large-scale fluctuations of polymorphism and divergence along chromosome arms, and significantly less polymorphism and faster divergence on the X chromosome. We generated a comprehensive list of functional elements in the D. simulans genome influenced by adaptive evolution. Finally, we characterized genomic patterns of base composition for coding and noncoding sequence. These results suggest several new hypotheses regarding the genetic and biological mechanisms controlling polymorphism and divergence across the Drosophila genome, and provide a rich resource for the investigation of adaptive evolution and functional variation in D. simulans.     

Odorant receptor (Or) genes: polymorphism and divergence in the D. melanogaster and D. pseudoobscura lineages

Background

In insects, like in most invertebrates, olfaction is the principal sensory modality, which provides animals with essential information for survival and reproduction. Odorant receptors are involved in this response, mediating interactions between an individual and its environment, as well as between individuals of the same or different species. The adaptive importance of odorant receptors renders them good candidates for having their variation shaped by natural selection.

Methodology/Principal Findings

We analyzed nucleotide variation in a subset of eight Or genes located on the 3L chromosomal arm of Drosophila melanogaster in a derived population of this species and also in a population of Drosophila pseudoobscura. Some heterogeneity in the silent polymorphism to divergence ratio was detected in the D. melanogaster/D. simulans comparison, with a single gene (Or67b) contributing ∼37% to the test statistic. However, no other signals of a very recent selective event were detected at this gene. In contrast, at the speciation timescale, the MK test uncovered the footprint of positive selection driving the evolution of two of the encoded proteins in both D. melanogaster —OR65c and OR67a —and D. pseudoobscura —OR65b1 and OR67c.

Conclusions

The powerful polymorphism/divergence approach provided evidence for adaptive evolution at a rather high proportion of the Or genes studied after relatively recent speciation events. It did not provide, however, clear evidence for very recent selective events in either D. melanogaster or D. pseudoobscura.     

Adaptive gene expression divergence inferred from population genomics

Detailed studies of individual genes have shown that gene expression divergence often results from adaptive evolution of regulatory sequence. Genome-wide analyses, however, have yet to unite patterns of gene expression with polymorphism and divergence to infer population genetic mechanisms underlying expression evolution. Here, we combined genomic expression data—analyzed in a phylogenetic context—with whole genome light-shotgun sequence data from six Drosophila simulans lines and reference sequences from D. melanogaster and D. yakuba. These data allowed us to use molecular population genetics to test for neutral versus adaptive gene expression divergence on a genomic scale. We identified recent and recurrent adaptive evolution along the D. simulans lineage by contrasting sequence polymorphism within D. simulans to divergence from D. melanogaster and D. yakuba. Genes that evolved higher levels of expression in D. simulans have experienced adaptive evolution of the associated 3′ flanking and amino acid sequence. Concomitantly, these genes are also decelerating in their rates of protein evolution, which is in agreement with the finding that highly expressed genes evolve slowly. Interestingly, adaptive evolution in 5′ cis-regulatory regions did not correspond strongly with expression evolution. Our results provide a genomic view of the intimate link between selection acting on a phenotype and associated genic evolution.     

Comparative population genomics of latitudinal variation in Drosophila simulans and Drosophila melanogaster

Examples of clinal variation in phenotypes and genotypes across latitudinal transects have served as important models for understanding how spatially varying selection and demographic forces shape variation within species. Here, we examine the selective and demographic contributions to latitudinal variation through the largest comparative genomic study to date of Drosophila simulans and Drosophila melanogaster, with genomic sequence data from 382 individual fruit flies, collected across a spatial transect of 19 degrees latitude and at multiple time points over 2 years. Consistent with phenotypic studies, we find less clinal variation in D. simulans than D. melanogaster, particularly for the autosomes. Moreover, we find that clinally varying loci in D. simulans are less stable over multiple years than comparable clines in D. melanogaster. D. simulans shows a significantly weaker pattern of isolation by distance than D. melanogaster and we find evidence for a stronger contribution of migration to D. simulans population genetic structure. While population bottlenecks and migration can plausibly explain the differences in stability of clinal variation between the two species, we also observe a significant enrichment of shared clinal genes, suggesting that the selective forces associated with climate are acting on the same genes and phenotypes in D. simulans and D. melanogaster.     

Molecular evolution of sex-biased genes in the Drosophila ananassae subgroup

Background  

Genes with sex-biased expression often show rapid molecular evolution between species. Previous population genetic and comparative genomic studies of Drosophila melanogaster and D. simulans revealed that male-biased genes have especially high rates of adaptive evolution. To test if this is also the case for other lineages within the melanogaster group, we investigated gene expression in D. ananassae, a species that occurs in structured populations in tropical and subtropical regions. We used custom-made microarrays and published microarray data to characterize the sex-biased expression of 129 D. ananassae genes whose D. melanogaster orthologs had been classified previously as male-biased, female-biased, or unbiased in their expression and had been studied extensively at the population-genetic level. For 43 of these genes we surveyed DNA sequence polymorphism in a natural population of D. ananassae and determined divergence to the sister species D. atripex and D. phaeopleura.     

Rapid Concerted Evolution via Gene Conversion at the <Emphasis Type="BoldItalic">Drosophila hsp70</Emphasis> Genes

We analyzed nucleotide variation in the hsp70 genes of Drosophila melanogaster (five genes) and D. simulans (four genes) to characterize the homogenizing and diversifying roles of gene conversion in their evolution. Gene conversion within and between the 87A7 and 87C1 gene clusters homogenize the hsp70 coding regions; in both D. melanogaster and D. simulans, same-cluster paralogues are virtually identical, and large intercluster conversion tracts diminish 87A7/87C1 divergence. Same-cluster paralogues share many polymorphisms, consistent with frequent intracluster conversion. Shared polymorphism is highly biased toward silent variation; homogenizing conversion interacts with purifying selection. In contrast to the coding regions, some hsp70 flanking regions show conversion-mediated diversification. Strong reductions of nucleotide variability and linkage disequilibria among conversion-mediated sites in hsp70Ab and hsp70Bb alleles sampled from a single natural population are consistent with a selective sweep. Comparison of the D. melanogaster and D. simulans hsp70 genes reveals whole-family fixed differences, consistent with rapid propagation of novel mutations among duplicate genes. These results suggest that the homogenizing and diversifying roles of conversion interact to drive dynamic concerted evolution of the hsp70 genes. Received: 25 June 2001 / Accepted: 10 October 2001     

Population Genomics: Whole-Genome Analysis of Polymorphism and Divergence in Drosophila simulans

The population genetic perspective is that the processes shaping genomic variation can be revealed only through simultaneous investigation of sequence polymorphism and divergence within and between closely related species. Here we present a population genetic analysis of Drosophila simulans based on whole-genome shotgun sequencing of multiple inbred lines and comparison of the resulting data to genome assemblies of the closely related species, D. melanogaster and D. yakuba. We discovered previously unknown, large-scale fluctuations of polymorphism and divergence along chromosome arms, and significantly less polymorphism and faster divergence on the X chromosome. We generated a comprehensive list of functional elements in the D. simulans genome influenced by adaptive evolution. Finally, we characterized genomic patterns of base composition for coding and noncoding sequence. These results suggest several new hypotheses regarding the genetic and biological mechanisms controlling polymorphism and divergence across the Drosophila genome, and provide a rich resource for the investigation of adaptive evolution and functional variation in D. simulans.     

Drosophila Melanogaster, Drosophila Simulans: so Similar yet so Different

During the last two decades, the two cosmopolitan species Drosophila melanogaster and Drosophila simulans have been compared with regard to numerous characteristics, ranging from their geographic distribution and ecology to their DNA polymorphism. Various traits have been compared, including morphology, physiology, sexual behavior, allozymes and other proteins, chromosomal inversions, mitochondrial and nuclear DNA, transposable elements, wolbachia etc. Such comparisons reveal similarities and differences between the two species, depending on the trait considered. In most cases, the between-population variability of D. simulans is lower than that of D. melanogaster, but the two species exhibit similar levels of within-population variability. One of the main exceptions is the nucleotide polymorphism of several nuclear regions. Although several hypotheses have been proposed to explain these observations, the evolutionary dynamics of these two species are far from being understood. How have two species sharing a common ancestor in the recent past accumulated so many differences? A brief history of comparisons of the two species, from the first in 1919 by A.H. Sturtevant, and a summary of the hypotheses proposed to explain the similarities and the differences between these species are presented and discussed.     

Genomic Patterns of Geographic Differentiation in Drosophila simulans

Geographic patterns of genetic differentiation have long been used to understand population history and to learn about the biological mechanisms of adaptation. Here we present an examination of genomic patterns of differentiation between northern and southern populations of Australian and North American Drosophila simulans, with an emphasis on characterizing signals of parallel differentiation. We report on the genomic scale of differentiation and functional enrichment of outlier SNPs. While, overall, signals of shared differentiation are modest, we find the strongest support for parallel differentiation in genomic regions that are associated with regulation. Comparisons to Drosophila melanogaster yield potential candidate genes involved in local adaptation in both species, providing insight into common selective pressures and responses. In contrast to D. melanogaster, in D. simulans we observe patterns of variation that are inconsistent with a model of temperate adaptation out of a tropical ancestral range, highlighting potential differences in demographic and colonization histories of this cosmopolitan species pair.     

Absence of Protein Polymorphism in the Ras Genes of Drosophila melanogaster

Sequence analysis of 27 alleles of each of the three Ras-related genes in Drosophila melanogaster indicates that they all have low levels of polymorphism but may experience slightly different evolutionary pressures. No amino acid replacement substitutions were indicated in any of the sequences, or in the sibling species D. simulans and D. mauritiana. The Dras1 gene, which is the major ras homologue in Drosophila, has less within-species variation in D. melanogaster relative to the amount of divergence from the sibling species than does Dras2, although the contrast was not significant by the HKA test. Dras2 appears to be maintaining two classes of haplotype in D. melanogaster, one of which is closer to the alleles observed in the sibling species, suggesting that this is not likely to be a pseudogene despite the absence of a mutant phenotype. Although differences in level of expression may affect the function of the genes, it is concluded that genetic variation in the Ras signal transduction pathways cannot be attributed to catalytic variation in the Ras proteins. Received: 5 November 1998 / Accepted: 26 March 1999     

Structural and genetic studies of the proliferation disrupter genes of Drosophila simulans and D. melanogaster

To examine whether structural and functional differences exist in the proliferation disrupter (prod) genes between Drosophila simulans and D. melanogaster, we analyzed and compared both genes. The exon–intron structure of the genes was found to be the same – three exons were interrupted by two introns, although a previous report suggested that only one intron existed in D. melanogaster. The prod genes of D. simulans and D. melanogaster both turn out to encode 346 amino acids, not 301 as previously reported for D. melanogaster. The numbers of nucleotide substitutions in the prod genes was 0.0747 ±  per synonymous site and 0.0116 ± 0.0039 per replacement site, both comparable to those previously known for homologous genes between D. simulans and D. melanogaster. Genetic analysis demonstrated that D. simulans PROD can compensate for a deficiency of D. melanogaster PROD in hybrids. The PRODs of D. simulans and D. melanogaster presumably share the same function and a conserved working mechanism. The prod gene showed no significant interaction with the lethality of the male hybrid between these species. This revised version was published online in July 2006 with corrections to the Cover Date.     

Effects of female behavior on wing display and courtship of male Drosophila simulans and Drosophila melanogaster

Corrolations between female rejection behaviors and male wing display were calculated for both Drosophila simulans and Drosophila melanogaster intraspicific pair-matings. No significant correlations were found for D. melanogaster, but in D. simulans flicking by the female appeared to be associated with a shift in male wing display pattern resulting in higher levels of vibration. Flicking did not appear to discourage courtship by males in either species.     

Genomic Evidence of Rapid and Stable Adaptive Oscillations over Seasonal Time Scales in Drosophila

In many species, genomic data have revealed pervasive adaptive evolution indicated by the fixation of beneficial alleles. However, when selection pressures are highly variable along a species'' range or through time adaptive alleles may persist at intermediate frequencies for long periods. So called “balanced polymorphisms” have long been understood to be an important component of standing genetic variation, yet direct evidence of the strength of balancing selection and the stability and prevalence of balanced polymorphisms has remained elusive. We hypothesized that environmental fluctuations among seasons in a North American orchard would impose temporally variable selection on Drosophila melanogaster that would drive repeatable adaptive oscillations at balanced polymorphisms. We identified hundreds of polymorphisms whose frequency oscillates among seasons and argue that these loci are subject to strong, temporally variable selection. We show that these polymorphisms respond to acute and persistent changes in climate and are associated in predictable ways with seasonally variable phenotypes. In addition, our results suggest that adaptively oscillating polymorphisms are likely millions of years old, with some possibly predating the divergence between D. melanogaster and D. simulans. Taken together, our results are consistent with a model of balancing selection wherein rapid temporal fluctuations in climate over generational time promotes adaptive genetic diversity at loci underlying polygenic variation in fitness related phenotypes.     

Normal Segregation of a Foreign-Species Chromosome During Drosophila Female Meiosis Despite Extensive Heterochromatin Divergence

The abundance and composition of heterochromatin changes rapidly between species and contributes to hybrid incompatibility and reproductive isolation. Heterochromatin differences may also destabilize chromosome segregation and cause meiotic drive, the non-Mendelian segregation of homologous chromosomes. Here we use a range of genetic and cytological assays to examine the meiotic properties of a Drosophila simulans chromosome 4 (sim-IV) introgressed into D. melanogaster. These two species differ by ∼12–13% at synonymous sites and several genes essential for chromosome segregation have experienced recurrent adaptive evolution since their divergence. Furthermore, their chromosome 4s are visibly different due to heterochromatin divergence, including in the AATAT pericentromeric satellite DNA. We find a visible imbalance in the positioning of the two chromosome 4s in sim-IV/mel-IV heterozygote and also replicate this finding with a D. melanogaster 4 containing a heterochromatic deletion. These results demonstrate that heterochromatin abundance can have a visible effect on chromosome positioning during meiosis. Despite this effect, however, we find that sim-IV segregates normally in both diplo and triplo 4 D. melanogaster females and does not experience elevated nondisjunction. We conclude that segregation abnormalities and a high level of meiotic drive are not inevitable byproducts of extensive heterochromatin divergence. Animal chromosomes typically contain large amounts of noncoding repetitive DNA that nevertheless varies widely between species. This variation may potentially induce non-Mendelian transmission of chromosomes. We have examined the meiotic properties and transmission of a highly diverged chromosome 4 from a foreign species within the fruitfly Drosophila melanogaster. This chromosome has substantially less of a simple sequence repeat than does D. melanogaster 4, and we find that this difference results in altered positioning when chromosomes align during meiosis. Yet this foreign chromosome segregates at normal frequencies, demonstrating that chromosome segregation can be robust to major differences in repetitive DNA abundance.     

A Comprehensive Study of Genic Variation in Natural Populations of Drosophila melanogaster. III. Variations in Genetic Structure and Their Causes between Drosophila melanogaster and Its Sibling Species Drosophila simulans

The natural populations of Drosophila melanogaster and Drosophila simulans were compared for their genetic structure. A total of 114 gene-protein loci were studied in four mainland (from Europe and Africa) and an island (Seychelle) populations of D. simulans and the results were compared with those obtained on the same set of homologous loci in fifteen worldwide populations of D. melanogaster. The main results are as follows: (1) D. melanogaster shows a significantly higher proportion of loci polymorphic than D. simulans (52% vs. 39%, P<0.05), (2) both species have similar mean heterozygosity and mean number of alleles per locus, (3) the two species share some highly polymorphic loci but they do not share loci that show high geographic differentiation, and (4) D. simulans shows significantly less geographic differentiation than D. melanogaster. The differences in genetic differentiation between the two species are limited to loci located on the X and second chromosomes only; loci on the third chromosome show similar level of geographic differentiation in both species. These two species have previously been shown to differ in their pattern of variation for chromosomal polymorphisms, quantitative and physiological characters, two-dimensional electrophoretic (2DE) proteins, middle repetitive DNA and mitochondrial DNA. Variation in niche-widths and/or genetic "strategies" of adaptation appear to be the main causes of differences in the genetic structure of these two species.     

Drosophila melanogaster males respond differently at the behavioural and genome‐wide levels to Drosophila melanogaster and Drosophila simulans females

Drosophila melanogaster are found in sympatry with Drosophila simulans, and matings between the species produce nonfertile hybrid offspring at low frequency. Evolutionary theory predicts that females choose mates, so males should alter their behaviour in response to female cues. We show that D. melanogaster males quickly decrease courtship towards D. simulans females. Courtship levels are reduced within 5 min of exposure to a heterospecific female, and overall courtship is significantly lower than courtship towards conspecific females. To understand changes at the molecular level during mate choice, we performed microarray analysis on D. melanogaster males that courted heterospecific D. simulans females and found nine genes have altered expression compared with controls. In contrast, males that court conspecific females alter expression of at least 35 loci. The changes elicited by conspecific courtship likely modulate nervous system function to reinforce positive conspecific signals and dampen the response to heterospecific signals.