Search

Chromosome inversions are of unique importance in the evolution of genomes and species because when heterozygous with a standard arrangement chromosome, they suppress meiotic crossovers within the inversion. In Drosophila species, heterozygous inversions also cause the interchromosomal effect, whereby the presence of a heterozygous inversion induces a dramatic increase in crossover frequencies in the remainder of the genome within a single meiosis. To date, the interchromosomal effect has been studied exclusively in species that also have high frequencies of inversions in wild populations. We took advantage of a recently developed approach for generating inversions in Drosophila simulans , a species that does not have inversions in wild populations, to ask if there is an interchromosomal effect. We used the existing chromosome 3R balancer and generated a new chromosome 2L balancer to assay for the interchromosomal effect genetically and cytologically. We found no evidence of an interchromosomal effect in D. simulans. To gain insights into the underlying mechanistic reasons, we qualitatively analyzed the relationship between meiotic double-stranded break (DSB) formation and synaptonemal complex (SC) assembly. We found that the SC is assembled prior to DSB formation as in D. melanogaster ; however, we show that the SC is assembled prior to localization of the oocyte determination factor Orb, whereas in D. melanogaster , SC formation does not begin until the Orb is localized. Together, our data show no evidence that heterozygous inversions in D. simulans induce an interchromosomal effect and that there are differences in the developmental programming of the early stages of meiosis. [ABSTRACT FROM AUTHOR]

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.

We develop analytical and simulation tools for evolve-and-resequencing experiments and apply them to a new study of rapid evolution in Drosophila simulans. Likelihood test statistics applied to pooled population sequencing data suggest parallel evolution of 138 SNPs across the genome. This number is reduced by orders of magnitude from previous studies (thousands or tens of thousands), owing to differences in both experimental design and statistical analysis. Whole genome simulations calibrated from Drosophila genetic data sets indicate that major features of the genome-wide response could be explained by as few as 30 loci under strong directional selection with a corresponding hitchhiking effect. Smaller effect loci are likely also responding, but are below the detection limit of the experiment. Finally, SNPs showing strong parallel evolution in the experiment are intermediate in frequency in the natural population (usually 30-70%) indicative of balancing selection in nature. These loci also exhibit elevated differentiation among natural populations of D. simulans, suggesting environmental heterogeneity as a potential balancing mechanism. [ABSTRACT FROM AUTHOR]

Crosses between Drosophila melanogaster females and Drosophila simulans males produce hybrid sons that die at the larval stage. This hybrid lethality is suppressed by loss-of-function mutations in the D. melanogaster Hybrid male rescue (Hmr) or in the D. simulans Lethal hybrid rescue (Lhr) genes. Previous studies have shown that Hmr and Lhr interact with heterochromatin proteins and suppress expression of transposable elements within D. melanogaster. It also has been proposed that Hmr and Lhr function at the centromere. We examined mitotic divisions in larval brains from Hmr and Lhr single mutants and Hmr; Lhr double mutants in D. melanogaster. In none of the mutants did we observe defects in metaphase chromosome alignment or hyperploid cells, which are hallmarks of centromere or kinetochore dysfunction. In addition, we found that Hmr-HA and Lhr-HA do not colocalize with centromeres either during interphase or mitotic division. However, all mutants displayed anaphase bridges and chromosome aberrations resulting from the breakage of these bridges, predominantly at the euchromatin–heterochromatin junction. The few dividing cells present in hybrid males showed fuzzy and irregularly condensed chromosomes with unresolved sister chromatids. Despite this defect in condensation, chromosomes in hybrids managed to align on the metaphase plate and undergo anaphase. We conclude that there is no evidence for a centromeric function of Hmr and Lhr within D. melanogaster nor for a centromere defect causing hybrid lethality. Instead, we find that Hmr and Lhr are required in D. melanogaster for detachment of sister chromatids during anaphase. [ABSTRACT FROM AUTHOR]

Two genes encoding protein components of the nuclear pore complex Nup160 and Nup96 cause lethality in F2-like hybrid genotypes between Drosophila simulans and Drosophila melanogaster. In particular, D. simulans Nup160 and Nup96 each cause inviability when hemizygous or homozygous in species hybrids that are also hemizygous (or homozygous) for the D. melanogaster X chromosome. The hybrid lethality of Nup160, however, is genetically complex, depending on one or more unknown additional factors in the autosomal background. Here we study the genetics and evolution of Nup160-mediated hybrid lethality in three ways. First, we test for variability in Nup160-mediated hybrid lethality within and among the three species of the D. simulans clade— D. simulans, D. sechellia, and D. mauritiana. We show that the hybrid lethality of Nup160 is fixed in D. simulans and D. sechellia but absent in D. mauritiana. Second, we explore how the hybrid lethality of Nup160 depends on other loci in the autosomal background. We find that D. simulans Nup160-mediated hybrid lethality does not depend on the presence of D. melanogaster Nup96, and we find that D. simulans and D. mauritiana are functionally differentiated at Nup160 as well as at other autosomal factor(s). Finally, we use population genetics data to show that Nup160 has experienced histories of recurrent positive selection both before and after the split of the three D. simulans clade species ~240,000 years ago. Our genetic results suggest that a hybrid lethal Nup160 allele evolved before the split of the three D. simulans clade species, whereas the other autosomal factor(s) evolved more recently. [ABSTRACT FROM AUTHOR]

Male sexual characters are often among the first traits to diverge between closely related species and identifying the genetic basis of such changes can contribute to our understanding of their evolutionary history. However, little is known about the genetic architecture or the specific genes underlying the evolution of male genitalia. The morphology of the claspers, posterior lobes, and anal plates exhibit striking differences between Drosophila mauritiana and D. simulans. Using QTL and introgression-based high-resolution mapping, we identified several small regions on chromosome arms 3L and 3R that contribute to differences in these traits. However, we found that the loci underlying the evolution of clasper differences between these two species are independent from those that contribute to posterior lobe and anal plate divergence. Furthermore, while most of the loci affect each trait in the same direction and act additively, we also found evidence for epistasis between loci for clasper bristle number. In addition, we conducted an RNAi screen in D. melanogasterto investigate if positional and expression candidate genes located on chromosome 3L, are also involved in genital development. We found that six of these genes, including components of Wnt signaling and male-specific lethal 3 (msl3), regulate the development of genital traits consistent with the effects of the introgressed regions where they are located and that thus represent promising candidate genes for the evolution these traits. [ABSTRACT FROM AUTHOR]

Hybrid incompatibility (HI) genes are frequently observed to be rapidly evolving under selection. This observation has led to the attractive conjecture that selection-derived protein-sequence divergence is culpable for incompatibilities in hybrids. The Drosophila simulans HI gene Lethal hybrid rescue (Lhr) is an intriguing case, because despite having experienced rapid sequence evolution, its HI properties are a shared function inherited from the ancestral state. Using an unusual D. simulans Lhr hybrid rescue allele, Lhr², we here identify a conserved stretch of 10 amino acids in the C terminus of LHR that is critical for causing hybrid incompatibility. Altering these 10 amino acids weakens or abolishes the ability of Lhr to suppress the hybrid rescue alleles Lhr¹ or Hmr¹, respectively. Besides single- amino-acid substitutions, Lhr orthologs differ by a 16-aa indel polymorphism, with the ancestral deletion state fixed in D. melanogaster and the derived insertion state at very high frequency in D. simulans. Lhr² is a rare D. simulans allele that has the ancestral deletion state of the 16-aa polymorphism. Through a series of transgenic constructs we demonstrate that the ancestral deletion state contributes to the rescue activity of Lhr². This indel is thus a polymorphism that can affect the HI function of Lhr. [ABSTRACT FROM AUTHOR]

Finding the genes underlying complex traits is difficult. We show that new sequencing technology combined with traditional genetic techniques can efficiently identify genetic regions underlying a complex and quantitative behavioral trait. As a proof of concept we used phenotype-based introgression to backcross loci that control innate food preference in Drosophila simulans into the genomic background of D. sechellia, which expresses the opposite preference. We successfully mapped D. simulans introgression regions in a small mapping population (30 flies) with whole-genome resequencing using light coverage (~1x). We found six loci contributing to D. simulans food preference, one of which overlaps a previously discovered allele. This approach is applicable to many systems, does not rely on laborious marker development or genotyping, does not require existing high quality reference genomes, and needs only small mapping populations. Because introgression is used, researchers can scale mapping population size, replication, and number of backcross generations to their needs. Finally, in contrast to more widely used mapping techniques like F2 bulk-segregant analysis, our method produces near-isogenic lines that can be kept and reused indefinitely. [ABSTRACT FROM AUTHOR]

Analysis of genomic data requires an efficient way to calculate likelihoods across very large numbers of loci. We describe a general method for finding the distribution of genealogies: we allow migration between demes, splitting of demes [as in the isolation-with-migration (IM) model], and recombination between linked loci. These processes are described by a set of linear recursions for the generating function of branch lengths. Under the infinite-sites model, the probability of any configuration of mutations can be found by differentiating this generating function. Such calculations are feasible for small numbers of sampled genomes: as an example, we show how the generating function can be derived explicitly for three genes under the two-deme IM model. This derivation is done automatically, using Mathematica. Given data from a large number of unlinked and nonrecombining blocks of sequence, these results can be used to find maximum-likelihood estimates of model parameters by tabulating the probabilities of all relevant mutational configurations and then multiplying across loci. The feasibility of the method is demonstrated by applying it to simulated data and to a data set previously analyzed by Wang and Hey (2010) consisting of 26,141 loci sampled from Drosophila simulans and D. melanogaster. Our results suggest that such likelihood calculations are scalable to genomic data as long as the numbers of sampled individuals and mutations per sequence block are small. [ABSTRACT FROM AUTHOR]

Several recent studies have suggested that F1 hybrid male lethality in crosses between Drosophila melanogaster and D. simulans is due to a failure in dosage compensation, caused by incompatibilities between D. simulans dosage compensation proteins and the D. melanogaster X chromosome. Contrary to the predictions of this hypothesis, mutations in four essential D. melanogaster dosage compensation genes are shown here to moderately increase rather than decrease hybrid male viability. [ABSTRACT FROM AUTHOR]

Selfish genes, such as meiotic drive elements, propagate themselves through a population without increasing the fitness of host organisms. X-linked (or Y-linked) meiotic drive elements reduce the transmission of the Y (X) chromosome and skew progeny and population sex ratios, leading to intense conflict among genomic compartments. Drosophila simulans is unusual in having a least three distinct systems of X chromosome meiotic drive. Here, we characterize naturally occurring genetic variation at the Winters sex-ratio driver (Distorter on the X or Dox), its progenitor gene (Mother of Dox or MDox), and its suppressor gene (Not Much Yang or Nmy), which have been previously mapped and characterized. We survey three North American populations as well as 13 globally distributed strains and present molecular polymorphism data at the three loci. We find that all three genes show signatures of selection in North America, judging from levels of polymorphism and skews in the site-frequency spectrum. These signatures likely result from the biased transmission of the driver and selection on the suppressor for the maintenance of equal sex ratios. Coalescent modeling indicates that the timing of selection is more recent than the age of the alleles, suggesting that the driver and suppressor are coevolving under an evolutionary ''arms race.'' None of the Winters sex-ratio genes are fixed in D. simulans, and at all loci we find ancestral alleles, which lack the gene insertions and exhibit high levels of nucleotide polymorphism compared to the derived alleles. In addition, we find several ''null'' alleles that have mutations on the derived Dox background, which result in loss of drive function. We discuss the possible causes of the maintenance of presence-absence polymorphism in the Winters sex-ratio genes. [ABSTRACT FROM AUTHOR]

A large-effect QTL for divergence in sex-comb tooth number between Drosophila simulans and D. mauritiana was previously mapped to 73A-84AB. Here we identify genes that are likely contributors to this divergence. We first improved the mapping resolution in the 73A-84AB region using 12 introgression lines and 62 recombinant nearly isogenic lines. To further narrow the list of candidate genes, we assayed leg-specific expression and identified genes with transcript-level evolution consistent with a potential role in sex-comb divergence. Sex combs are formed on the prothoracic (front) legs, but not on the mesothoracic (middle) legs of Drosophila males. We extracted RNA from the prothoracic and mesothoracic pupal legs of two species to determine which of the genes expressed differently between leg types were also divergent for gene expression. Two good functional candidate genes, Scr and dsx, are located in one of our line-scale QTL regions. In addition, three previously uncharacterized genes (CG15186, CG2016, and CG2791) emerged as new candidates. These genes are located in regions strongly associated with sex-comb tooth number differences and are expressed differently between leg tissues and between species. Further supporting the potential involvement of these genes in sex-comb divergence, we found a significant difference in sex-comb tooth number between co-isogenic D. melanogaster lines with and without P-element insertions at CG2791 . [ABSTRACT FROM AUTHOR]

The cross Drosophila melanogaster females to D. simulans males typically produces lethal F1 hybrid males. F1 male lethality is suppressed when the D. simulans Lhr¹ hybrid rescue strain is used. Viability of these F1 males carrying Lhr¹ is in turn substantially reduced when the hybrids are heterozygous for some mutant alleles of the D. melanogaster Nup96 gene. I show here that similar patterns of Nup96-dependent lethality occur when other hybrid rescue mutations are used to create F1 males, demonstrating that Nup96 does not reduce hybrid viability by suppressing the Lhr¹ rescue effect. The penetrance of this Nup96-dependent lethality dos not correlate with the penetrance of the F1 hybrid rescue, arguing that these two phenomena reflect genetically independent processes. D. simulans, together with two additional sister species, forms a clade that speciated after the divergence of their common ancestor from D. melanogaster. I report there that Nup96- reduces F1 viability in D. melanogaster hybrids with one of these sister species, D. sechellia, but not with the other, D. mauritiana. These results suggest that Nup96-dependent lethality evolved after the speciation of D. melanogaster from the common ancestor of the simulans clade and is caused by an interaction among Nup96, unknown gene(s) on the D. melanogaster X chromosome, and unknown autosomal gene(s), at least some of which have diverged in D. simulans and D. sechellia but not in D. mauritiana. The genetic properties of Nup96 are also discussed relative to other hybrid lethal genes. [ABSTRACT FROM AUTHOR]

Nucleotide variation at the genomic region encompassing the odorant-binding protein genes OS-E and OS-F (OS region) was surveyed in two populations of Drosophila simulans, one from Europe and the other from Africa. We found that the European population shows an atypical and large haplotype structure, which extends throughout the ~5-kb surveyed genomic region. This structure is depicted by two major haplotype groups segregating at intermediate frequency in the sample, one haplogroup with nearly no variation, and the other at levels more typical for this species. This pattern of variation was in- compatible with neutral predictions for a population at a stationary equilibrium. Nevertheless, neutrality tests contrasting polymorphism and divergence data fail to detect any departure from the standard neutral model in this species, whereas they confirm the non-neutral behavior previously observed at the OS-E gene in D. melanogaster. Although positive Darwinian selection may have been responsible for the observed unusual nucleotide variation structure, coalescent simulation results do not allow rejecting the hypothesis that the pattern was generated by a recent bottleneck in the history of European populations of B. simulans. [ABSTRACT FROM AUTHOR]

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 α¯. But such tests are subject to many biases. Most serious is the possibility that high estimates of α¯ reflect demographic changes rather than adaptive substitution. Testing for between-locus variation in α is one possible way of distinguishing between demography and selection. However, such tests have yielded contradictory results, and their efficacy is unclear. Estimates of α¯ 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 α¯ 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. [ABSTRACT FROM AUTHOR]

African populations of Drosophila simulans are thought to be ancestral in this model species and are increasingly used for testing general hypotheses in evolutionary genetics. It is often assumed that African populations are more likely to be at a neutral mutation drift equilibrium than other populations. Here we examine population structuring and the demographic profile in nine populations of D. simulans. We surveyed sequence variation in four X-linked genes (runt, sevenless, Sex-lethal, and vermilion) that have been used in a parallel study in the closely related species D. melanogaster. We found that an eastern group of populations from continental Africa and Indian Ocean islands (Kenya, Tanzania, Madagascar, and Mayotte Island) is widespread, shows little differentiation, and has probably undergone demographic expansion. The other two African populations surveyed (Cameroon and Zimbabwe) show no evidence of population expansion and are markedly differentiated from each other as well as from the populations from the eastern group. Two other populations, Europe and Antilles, are probably recent invaders to these areas. The Antilles population is probably derived from Europe through a substantial bottleneck. The history of these populations should be taken into account when drawing general conclusions from variation patterns. [ABSTRACT FROM AUTHOR]

Cytoplasmic incompatibility (CI) is a reproductive incompatibility induced by maternally transmitted bacteria of the genera Wolbachia and Cardinium. In the simplest form of CI, offspring from infected males and uninfected females suffer from increased mortality. However, it has been noted that crosses between males and females carrying different strains of infection are often also incompatible. The evolutionary processes leading to the emergence of new CI-compatibility types are still not resolved. Here, we develop a model that extends previous theoretical approaches by including segregation of bacterial strains during transmission as well as a continuum of breeding systems ranging from inbreeding (complete sib mating) to outbreeding (complete sib-mating avoidance). Our results demonstrate that (1) with segregation of strains, evolution is unlikely to lead to new CI types that co-occur as a double infection with the preexisting one, (2) inbreeding substantially hampers the evolution of new CI types, and (3) outbreeding facilitates the evolution of new CI types. Our model also provides a hypothesis on the evolutionary origin of CI. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

The study of base composition evolution in Drosophila has been achieved mostly through the analysis of coding sequences. Third codon position GC content, however, is influenced by both neutral forces (e.g., mutation bias) and natural selection for codon usage optimization. In this article, large data sets of noncoding DNA sequence polymorphism in D. melanogaster and D. simulans were gathered from public databases to try to disentangle these two factors—noncoding sequences are not affected by selection for codon usage. Allele frequency analyses revealed an asymmetric pattern of AT vs. GC noncoding polymorphisms: AT ↵ GC mutations are less numerous, and tend to segregate at a higher frequency, than GC ↵ AT ones, especially at GC-rich loci. This is indicative of nonstationary evolution of base composition and/or of GC-biased allele transmission. Fitting population genetics models to the allele frequency spectra confirmed this result and favored the hypothesis of a biased transmission. These results, together with previous reports, suggest that GC-biased gene conversion has influenced base composition evolution in Drosophila and explain the correlation between intron and exon GC content. [ABSTRACT FROM AUTHOR]

The rRNA gene (rDNA) loci of all arthropod lineages contain non-LTR retrotransposable elements that have evolved to specifically insert into the 28S rRNA genes. Extensive in vitro experiments have been conducted to investigate the mechanism of R2 retrotransposition but little is known of the insertion frequency or cellular factors that might regulate R2 activity. In this article, isofemale lines obtained from a population of Drosophila simulans were surveyed for recent R2 insertions. Within most lines, all individuals showed the same collection of R2 insertions, providing no evidence for recent R2 activity. However, in a few of the isofemale lines, virtually all individuals differed in their R2 insertion profiles. The descendants of individual pairs of flies from these ‘active lines’ rapidly accumulated new insertions. The frequent insertion of new R2 elements was associated with the elimination of old t12 elements from the rDNA locus. The existence of lines in which R2 retrotransposes frequently and lines in which the elements appear dormant suggests that cellular mechanisms that can regulate the activity of R2 exist. Retrotransposition activity was correlated with the number of full-length R2 elements but not with the size of the rDNA locus or the number of uninserted units. [ABSTRACT FROM AUTHOR]

Immune system genes in a California population sample of Drosophila simulans were shown to bear several hallmarks of the effects of past directional selection. One potential effect of directional selection is an increase in linkage disequilibrium among the polymorphic sites that are linked to the site under selection. In this study, we focus on three D. simulans immunity loci, Hmu, Sr-CI/Sr-CIII, and Tehao, for which the polymorphic sites are in nearly perfect linkage disequilibrium, an unusual finding even with respect to other immunity genes sampled from the same lines. The most likely explanation for this finding is that, at each locus, two divergent alleles have been selected to intermediate frequencies in the recent past. The extent to which the linkage disequilibrium extends to the flanks of each of the immunity genes is minimal, suggesting that the favored mutations actually occurred within the immunity genes themselves. Furthermore, the excess linkage disequilibrium found in the California population is not found in an African D. simulans population sample and may be a result of novel pathogen-mediated selection pressures encountered during establishment of non-African populations. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Understanding how genetic variation is maintained begins with a comprehensive description of what types of genetic variation exist, the extent and magnitude of the variation, and patterns discernable in that variation. However, such studies have focused primarily on DNA sequence data and have ignored genetic variation at other hierarchical levels of genetic information. Microarray technology permits an examination of genetic variation at the level of mRNA abundance. Utilizing a round-robin design, we present a quantitative description of variation in mRNA abundance in terms of GCA (general combining ability or additive variance). We test whether genes significant for GGA are randomly distributed across chromosomes and use a nonparametric approach to demonstrate that the magnitude of the variation is not random for GCA. We find that there is a paucity of genes significant for GGA on the X relative to the autosomes. The overall magnitude of the effects for GCA on the K tends to be lower than that on the autosomes and is contributed by rare alleles of larger effect, Due to male hemizygosity, GCA for K-linked phenotypes must be due to trans-acting factors, while GCA for autosomal phenotypes may be due to cis- or trans-acting factors. The contrast in the amount of variation between the X and the autosomes suggests that both cis and trans factors contribute to variation for expression in P. simulans with the preponderance of effects being trans. This nonrandom patterning of genetic variation in gene expression data with respect to chromosomal context may be due to hemizygosity in the male. [ABSTRACT FROM AUTHOR]

The recent action of positive selection is expected to influence patterns of intraspecific DNA sequence variation in chromosomal regions linked to the selected locus. These effects include decreased polymorphism, increased linkage disequilibrium, and an increased frequency of derived variants. These effects are all expected to dissipate with distance from the selected locus due to recombination. Therefore, in regions of high recombination, it should be possible to localize a target of selection to a relatively small interval. Previously described patterns of intraspecific variation in three tandemly arranged, testes-expressed genes (janusA, janusB, and ocnus) in Drosophila simulans included all three of these features. Here we expand the original sample and also survey nucleotide polymorphism at three neighboring loci. On the basis of recombination events between derived and ancestral alleles, we localize the target of selection to a 1.5-kb region surrounding janusB. A composite-likelihood-ratio test based on the spatial distribution and frequency of derived polymorphic variants corroborates this result and provides an estimate of the strength of selection. However, the data are difficult to reconcile with the simplest model of positive selection, whereas a new composite-likelihood method suggests that the data are better described by a model in which the selected allele has not yet gone to fixation. [ABSTRACT FROM AUTHOR]

Accessory gland proteins are a major component of Drosophila seminal fluid. These proteins have a variety of functions and may be subject to sexual selection and/or antagonistic evolution between the sexes. Most population genetic data from these proteins are frown D. melanogaster and D. simulans. Here, we extend the population genetic analysis of Acp genes to the other simulans complex species, D. mauritiana and D. sechellia. We sequenced population samples of seven Acp's from D. mauritiana, D. sechellia, and D. simulans. We investigated the population genetics of these genes on individual simulans complex lineages and compared Acp polymorphism and divergence to polymorphism and divergence from a set of nonAcp loci in the same species. Polymorphism and divergence data frown the simulans complex revealed little evidence for adaptive protein evolution at individual loci. However, we observed a dramatically inflated index of dispersion for amino acid substitutions in the simulans complex at Acp genes, but not at non-Acp genes. This pattern of episodic bursts of protein evolution in Acp's provides the strongest evidence to date that the population genetic mechanisms driving Acp divergence are different from the mechanisms driving evolution at most Drosophila genes. [ABSTRACT FROM AUTHOR]

The genetic architecture of traits influencing sexual isolation can give insight into the evolution of reproductive isolation and hence speciation. Here we report a quantitative trait loci (QTL) analysis of the difference in mean interpulse interval (IPI), an important component of the male courtship song, between Drosophila simulans and D. sechellia. Using a backcross analysis, we find six QTL that explain a total of 40.7% of the phenotypic variance. Three candidate genes are located in the intervals bounded by two of the QTL and there are no significant QTL on the X chromosome. The values of mean IPI for hybrid individuals imply the presence of dominant alleles or epistasis. Because unisexual hybrid sterility prevents an F2 analysis, we cannot distinguish dominant from additive genetic effects at the scale of QTL. A comparison with a study of QTL for intraspecific variation in D. melanogaster shows that, for these strains, the QTL we have identified for interspecific variation cannot be those that contribute to intraspecific variation. We find that the QTL have bidirectional effects, which indicates that the genetic architecture is compatible with divergence due to genetic drift, although other possibilities are discussed. [ABSTRACT FROM AUTHOR]