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

Pool-GWAS on reproductive dormancy in Drosophila simulans suggests a polygenic architecture

The genetic basis of adaptation to different environments has been of long-standing interest to evolutionary biologists. Dormancy is a well-studied adaptation to facilitate overwintering. In Drosophila melanogaster, a moderate number of genes with large effects have been described, which suggests a simple genetic basis of dormancy. On the other hand, genome-wide scans for dormancy suggest a polygenic architecture in insects. In D. melanogaster, the analysis of the genetic architecture of dormancy is complicated by the presence of cosmopolitan inversions. Here, we performed a genome-wide scan to characterize the genetic basis of this ecologically extremely important trait in the sibling species of D. melanogaster, D. simulans that lacks cosmopolitan inversions. We performed Pool-GWAS in a South African D. simulans population for dormancy incidence at 2 temperature regimes (10 and 12°C, LD 10:14). We identified several genes with SNPs that showed a significant association with dormancy (P-value < 1e-13), but the overall modest response suggests that dormancy is a polygenic trait with many loci of small effect. Our results shed light on controversies on reproductive dormancy in Drosophila and have important implications for the characterization of the genetic basis of this trait.

Light dependent courtship behavior in Drosophila simulans and D. melanogaster

Differences in courtship signals and perception are well-known among Drosophila species. One such described difference is the dependency on light, and thus presumably vision, for copulation success. Many studies have described a difference in light-dependent copulation success between D. melanogaster and D. simulans, identifying D. simulans as a light-dependent species, and D. melanogaster as a light-independent one. However, many of these studies use assays of varying design and few strains to represent the entire species. Here, we attempt to better characterize this purported difference using 11 strains of each species, paired by collection location, in behavioral assays conducted at two different exposure times. We show that, while there is a species-wide difference in magnitude of light-dependent copulation success, D. melanogaster copulation success is, on average, still impaired in the dark at both exposure times we measured. Additionally, there is significant variation in strain-specific ability to copulate in the dark in both species across two different exposure times. We find that this variation correlates strongly with longitude in D. melanogaster, but not in D. simulans. We hypothesize that differences in species history and demography may explain behavioral variation. Finally, we use courtship assays to show that light-dependent copulation success in one D. simulans strain is driven in part by both males and females. We discuss potential differences in courtship signals and/or signal importance between these species and potential for further comparative studies for functional characterization.

Rapid and Predictable Evolution of Admixed Populations Between Two Drosophila Species Pairs

The consequences of hybridization are varied, ranging from the origin of new lineages, introgression of some genes between species, to the extinction of one of the hybridizing species. We generated replicate admixed populations between two pairs of sister species of Drosophila: D. simulans and D. mauritiana; and D. yakuba and D. santomea Each pair consisted of a continental species and an island endemic. The admixed populations were maintained by random mating in discrete generations for over 20 generations. We assessed morphological, behavioral, and fitness-related traits from each replicate population periodically, and sequenced genomic DNA from the populations at generation 20. For both pairs of species, species-specific traits and their genomes regressed to those of the continental species. A few alleles from the island species persisted, but they tended to be proportionally rare among all sites in the genome and were rarely fixed within the populations. This paucity of alleles from the island species was particularly pronounced on the X-chromosome. These results indicate that nearly all foreign genes were quickly eliminated after hybridization and that selection against the minor species genome might be similar across experimental replicates.

piRNA silencing contributes to interspecies hybrid sterility and reproductive isolation in Drosophila melanogaster

The piRNA pathway is an adaptive mechanism that maintains genome stability by repression of selfish genomic elements. In the male germline of Drosophila melanogaster repression of Stellate genes by piRNAs generated from Supressor of Stellate (Su(Ste)) locus is required for male fertility, but both Su(Ste) piRNAs and their targets are absent in other Drosophila species. We found that D. melanogaster genome contains multiple X-linked non-coding genomic repeats that have sequence similarity to the protein-coding host gene vasa. In the male germline, these vasa-related AT-chX repeats produce abundant piRNAs that are antisense to vasa; however, vasa mRNA escapes silencing due to imperfect complementarity to AT-chX piRNAs. Unexpectedly, we discovered AT-chX piRNAs target vasa of Drosophila mauritiana in the testes of interspecies hybrids. In the majority of hybrid flies, the testes were strongly reduced in size and germline content. A minority of hybrids maintained wild-type array of premeiotic germ cells in the testes, but in them harmful Stellate genes were derepressed due to the absence of Su(Ste) piRNAs, and meiotic failures were observed. Thus, the piRNA pathway contributes to reproductive isolation between D. melanogaster and closely related species, causing hybrid male sterility via misregulation of two different host protein factors.

Loss of cytoplasmic incompatibility and minimal fecundity effects explain relatively low Wolbachia frequencies in Drosophila mauritiana

Maternally transmitted Wolbachia bacteria infect about half of all insect species. Many Wolbachia cause cytoplasmic incompatibility (CI) and reduced egg hatch when uninfected females mate with infected males. Although CI produces a frequency-dependent fitness advantage that leads to high equilibrium Wolbachia frequencies, it does not aid Wolbachia spread from low frequencies. Indeed, the fitness advantages that produce initial Wolbachia spread and maintain non-CI Wolbachia remain elusive. wMau Wolbachia infecting Drosophila mauritiana do not cause CI, despite being very similar to CI-causing wNo from Drosophila simulans (0.068% sequence divergence over 682,494 bp), suggesting recent CI loss. Using draft wMau genomes, we identify a deletion in a CI-associated gene, consistent with theory predicting that selection within host lineages does not act to increase or maintain CI. In the laboratory, wMau shows near-perfect maternal transmission; but we find no significant effect on host fecundity, in contrast to published data. Intermediate wMau frequencies on the island of Mauritius are consistent with a balance between unidentified small, positive fitness effects and imperfect maternal transmission. Our phylogenomic analyses suggest that group-B Wolbachia, including wMau and wPip, diverged from group-A Wolbachia, such as wMel and wRi, 6-46 million years ago, more recently than previously estimated.

Gene flow mediates the role of sex chromosome meiotic drive during complex speciation

During speciation, sex chromosomes often accumulate interspecific genetic incompatibilities faster than the rest of the genome. The drive theory posits that sex chromosomes are susceptible to recurrent bouts of meiotic drive and suppression, causing the evolutionary build-up of divergent cryptic sex-linked drive systems and, incidentally, genetic incompatibilities. To assess the role of drive during speciation, we combine high-resolution genetic mapping of X-linked hybrid male sterility with population genomics analyses of divergence and recent gene flow between the fruitfly species, Drosophila mauritiana and D. simulans. Our findings reveal a high density of genetic incompatibilities and a corresponding dearth of gene flow on the X chromosome. Surprisingly, we find that a known drive element recently migrated between species and, rather than contributing to interspecific divergence, caused a strong reduction in local sequence divergence, undermining the evolution of hybrid sterility. Gene flow can therefore mediate the effects of selfish genetic elements during speciation.

Redefining reproductive dormancy in Drosophila as a general stress response to cold temperatures

Organisms regularly encounter unfavorable conditions and the genetic adaptations facilitating survival have been of long-standing interest to evolutionary biologists. Winter is one particularly stressful condition for insects, during which they encounter low temperatures and scarcity of food. Despite dormancy being a well-studied adaptation to facilitate overwintering, there is still considerable controversy about the distribution of dormancy among natural populations and between species in Drosophila. The current definition of dormancy as developmental arrest of oogenesis at the previtellogenic stage (stage 7) distinguishes dormancy from general stress related block of oogenesis at early vitellogenic stages (stages 8 - 9). In an attempt to resolve this, we scrutinized reproductive dormancy in D. melanogaster and D. simulans. We show that dormancy shows the same hallmarks of arrest of oogenesis at stage 9, as described for other stressors and propose a new classification for dormancy. Applying this modified classification, we show that both species express dormancy in cosmopolitan and African populations, further supporting that dormancy uses an ancestral pathway induced by environmental stress. While we found significant differences between individuals and the two Drosophila species in their sensitivity to cold temperature stress, we also noted that extreme temperature stress (8 °C) resulted in very strong dormancy incidence, which strongly reduced the differences seen at less extreme temperatures. We conclude that dormancy in Drosophila should not be considered a special trait, but is better understood as a generic stress response occurring at low temperatures.

Variation in the Intensity of Selection on Codon Bias over Time Causes Contrasting Patterns of Base Composition Evolution in Drosophila

Four-fold degenerate coding sites form a major component of the genome, and are often used to make inferences about selection and demography, so that understanding their evolution is important. Despite previous efforts, many questions regarding the causes of base composition changes at these sites in Drosophila remain unanswered. To shed further light on this issue, we obtained a new whole-genome polymorphism data set from D. simulans. We analyzed samples from the putatively ancestral range of D. simulans, as well as an existing polymorphism data set from an African population of D. melanogaster. By using D. yakuba as an outgroup, we found clear evidence for selection on 4-fold sites along both lineages over a substantial period, with the intensity of selection increasing with GC content. Based on an explicit model of base composition evolution, we suggest that the observed AT-biased substitution pattern in both lineages is probably due to an ancestral reduction in selection intensity, and is unlikely to be the result of an increase in mutational bias towards AT alone. By using two polymorphism-based methods for estimating selection coefficients over different timescales, we show that the selection intensity on codon usage has been rather stable in D. simulans in the recent past, but the long-term estimates in D. melanogaster are much higher than the short-term ones, indicating a continuing decline in selection intensity, to such an extent that the short-term estimates suggest that selection is only active in the most GC-rich parts of the genome. Finally, we provide evidence for complex evolutionary patterns in the putatively neutral short introns, which cannot be explained by the standard GC-biased gene conversion model. These results reveal a dynamic picture of base composition evolution.

Genomics of parallel adaptation at two timescales in Drosophila

Two interesting unanswered questions are the extent to which both the broad patterns and genetic details of adaptive divergence are repeatable across species, and the timescales over which parallel adaptation may be observed. Drosophila melanogaster is a key model system for population and evolutionary genomics. Findings from genetics and genomics suggest that recent adaptation to latitudinal environmental variation (on the timescale of hundreds or thousands of years) associated with Out-of-Africa colonization plays an important role in maintaining biological variation in the species. Additionally, studies of interspecific differences between D. melanogaster and its sister species D. simulans have revealed that a substantial proportion of proteins and amino acid residues exhibit adaptive divergence on a roughly few million years long timescale. Here we use population genomic approaches to attack the problem of parallelism between D. melanogaster and a highly diverged conger, D. hydei, on two timescales. D. hydei, a member of the repleta group of Drosophila, is similar to D. melanogaster, in that it too appears to be a recently cosmopolitan species and recent colonizer of high latitude environments. We observed parallelism both for genes exhibiting latitudinal allele frequency differentiation within species and for genes exhibiting recurrent adaptive protein divergence between species. Greater parallelism was observed for long-term adaptive protein evolution and this parallelism includes not only the specific genes/proteins that exhibit adaptive evolution, but extends even to the magnitudes of the selective effects on interspecific protein differences. Thus, despite the roughly 50 million years of time separating D. melanogaster and D. hydei, and despite their considerably divergent biology, they exhibit substantial parallelism, suggesting the existence of a fundamental predictability of adaptive evolution in the genus.

Both local adaptation on short timescales and the long-term accumulation of adaptive differences between species have recently been investigated using comparative genomic and population genomic approaches in several species. However, the repeatability of adaptive evolution at the genetic level is poorly understood. Here we attack this problem by comparing patterns of long and short-term adaptation in Drosophila melanogaster to patterns of adaptation on two timescales in a highly diverged congener, Drosophila hydei. We found, despite the fact that these species diverged from a common ancestor roughly 50 million years ago, the population genomics of latitudinal allele frequency differentiation shows that there is a substantial shared set of genes likely playing a role in the short term adaptive divergence of populations in both species. Analyses of longer-term adaptive protein divergence for the D. hydei-D. mojavensis and D. melanogaster-D. simulans clades reveal a striking level of parallel adaptation. This parallelism includes not only the specific genes/proteins that exhibit adaptive evolution, but extends even to the magnitudes of the selective effects on interspecific protein differences.

The Drosophila simulans Y chromosome interacts with the autosomes to influence male fitness

The Y chromosome should degenerate because it cannot recombine. However, male-limited transmission increases selection efficiency for male-benefit alleles on the Y, and therefore, Y chromosomes should contribute significantly to variation in male fitness. This means that although the Drosophila Y chromosome is small and gene-poor, Y-linked genes are vital for male fertility in Drosophila melanogaster and the Y chromosome has large male fitness effects. It is unclear whether the same pattern is seen in the closely related Drosophila simulans. We backcrossed Y chromosomes from three geographic locations into five genetic backgrounds and found strong Y and genetic background effects on male fertility. There was a significant Y-background interaction, indicating substantial epistasis between the Y and autosomal genes affecting male fertility. This supports accumulating evidence that interactions between the Y chromosome and the autosomes are key determinants of male fitness.

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.

Parallel Gene Expression Differences between Low and High Latitude Populations of Drosophila melanogaster and D. simulans

Gene expression variation within species is relatively common, however, the role of natural selection in the maintenance of this variation is poorly understood. Here we investigate low and high latitude populations of Drosophila melanogaster and its sister species, D. simulans, to determine whether the two species show similar patterns of population differentiation, consistent with a role for spatially varying selection in maintaining gene expression variation. We compared at two temperatures the whole male transcriptome of D. melanogaster and D. simulans sampled from Panama City (Panama) and Maine (USA). We observed a significant excess of genes exhibiting differential expression in both species, consistent with parallel adaptation to heterogeneous environments. Moreover, the majority of genes showing parallel expression differentiation showed the same direction of differential expression in the two species and the magnitudes of expression differences between high and low latitude populations were correlated across species, further bolstering the conclusion that parallelism for expression phenotypes results from spatially varying selection. However, the species also exhibited important differences in expression phenotypes. For example, the genomic extent of genotype × environment interaction was much more common in D. melanogaster. Highly differentiated SNPs between low and high latitudes were enriched in the 3’ UTRs and CDS of the geographically differently expressed genes in both species, consistent with an important role for cis-acting variants in driving local adaptation for expression-related phenotypes.

While gene expression variation in natural populations is common, the population genetic processes responsible for the maintenance of this variation remain obscure. Here we study geographic differences in gene expression in recently established low and high latitude populations of two closely related species of Drosophila. We observe substantial parallelism in expression differences and expression plasticity between populations, which supports the idea that spatially varying selection correlated with latitude contributes to the maintenance of gene expression variation in these species. Comparison of inter-population sequence differentiation and expression differentiation suggests that cis-acting variants play a role in geographic expression differentiation.

Drosophila suzukii: the genetic footprint of a recent, worldwide invasion

Native to Asia, the soft-skinned fruit pest Drosophila suzukii has recently invaded the United States and Europe. The eastern United States represents the most recent expansion of their range, and presents an opportunity to test alternative models of colonization history. Here, we investigate the genetic population structure of this invasive fruit fly, with a focus on the eastern United States. We sequenced six X-linked gene fragments from 246 individuals collected from a total of 12 populations. We examine patterns of genetic diversity within and between populations and explore alternative colonization scenarios using approximate Bayesian computation. Our results indicate high levels of nucleotide diversity in this species and suggest that the recent invasions of Europe and the continental United States are independent demographic events. More broadly speaking, our results highlight the importance of integrating population structure into demographic models, particularly when attempting to reconstruct invasion histories. Finally, our simulation results illustrate the general challenge in reconstructing invasion histories using genetic data and suggest that genome-level data are often required to distinguish among alternative demographic scenarios.

Enzymatic characterization of recombinant α-amylase in the Drosophila melanogaster species subgroup: is there an effect of specialization on digestive enzyme?

We performed a comparative study on the enzymological features of purified recombinant α-amylase of three species belonging to the Drosophila melanogaster species subgroup: D. melanogaster, D. erecta and D. sechellia. D. erecta and D. sechellia are specialist species, with host plant Pandanus candelabrum (Pandanaceae) and Morinda citrifolia (Rubiaceae), respectively. The temperature optima were around 57-60℃ for the three species. The pH optima were 7.2 for D. melanogaster, 8.2 for D. erecta and 8.5 for D. sechellia. The kcat and Km were also estimated for each species with different substrates. The specialist species D. erecta and D. sechellia display a higher affinity for starch than D. melanogaster. α-Amylase activity is higher on starch than on glycogen in all species. α-Amylases of D. erecta and D. sechellia have a higher activity on maltooligosaccharides (G6 and G7) than on starch, contrary to D. melanogaster. Such differences in the enzymological features between the species might reflect adaptation to different ecological niches and feeding habits.

Genome diversity and divergence in Drosophila mauritiana: multiple signatures of faster X evolution

Drosophila mauritiana is an Indian Ocean island endemic species that diverged from its two sister species, Drosophila simulans and Drosophila sechellia, approximately 240,000 years ago. Multiple forms of incomplete reproductive isolation have evolved among these species, including sexual, gametic, ecological, and intrinsic postzygotic barriers, with crosses among all three species conforming to Haldane’s rule: F₁ hybrid males are sterile and F₁ hybrid females are fertile. Extensive genetic resources and the fertility of hybrid females have made D. mauritiana, in particular, an important model for speciation genetics. Analyses between D. mauritiana and both of its siblings have shown that the X chromosome makes a disproportionate contribution to hybrid male sterility. But why the X plays a special role in the evolution of hybrid sterility in these, and other, species remains an unsolved problem. To complement functional genetic analyses, we have investigated the population genomics of D. mauritiana, giving special attention to differences between the X and the autosomes. We present a de novo genome assembly of D. mauritiana annotated with RNAseq data and a whole-genome analysis of polymorphism and divergence from ten individuals. Our analyses show that, relative to the autosomes, the X chromosome has reduced nucleotide diversity but elevated nucleotide divergence; an excess of recurrent adaptive evolution at its protein-coding genes; an excess of recent, strong selective sweeps; and a large excess of satellite DNA. Interestingly, one of two centimorgan-scale selective sweeps on the D. mauritiana X chromosome spans a region containing two sex-ratio meiotic drive elements and a high concentration of satellite DNA. Furthermore, genes with roles in reproduction and chromosome biology are enriched among genes that have histories of recurrent adaptive protein evolution. Together, these genome-wide analyses suggest that genetic conflict and frequent positive natural selection on the X chromosome have shaped the molecular evolutionary history of D. mauritiana, refining our understanding of the possible causes of the large X-effect in speciation.

Distinguishing migration from isolation using genes with intragenic recombination: detecting introgression in the Drosophila simulans species complex

BACKGROUND

Determining the presence or absence of gene flow between populations is the target of some statistical methods in population genetics. Until recently, these methods either avoided the use of recombining genes, or treated recombination as a nuisance parameter. However, genes with recombination contribute additional information for the detection of gene flow (i.e. through linkage disequilibrium).

METHODS

We present three summary statistics based on the spatial arrangement of fixed differences, and shared and exclusive polymorphisms that are sensitive to the presence and direction of gene flow. Power and false positive rate for tests based on these statistics are studied by simulation.

RESULTS

The application of these tests to populations from the Drosophila simulans species complex yielded results consistent with migration between D. simulans and its two endemic sister species D. mauritiana and D. sechellia, and between populations D. mauritiana on the islands of the Mauritius and Rodrigues.

CONCLUSIONS

We demonstrate the sensitivity of the developed statistics to the presence and direction of gene flow, and characterize their power as a function of differentiation level and recombination rate. The properties of these statistics make them especially suitable for analyzing high-throughput sequencing data or for their integration within the approximate Bayesian computation framework.

Population history and pathways of spread of the plant pathogen Phytophthora plurivora

Human activity has been shown to considerably affect the spread of dangerous pests and pathogens worldwide. Therefore, strict regulations of international trade exist for particularly harmful pathogenic organisms. Phytophthora plurivora, which is not subject to regulations, is a plant pathogen frequently found on a broad range of host species, both in natural and artificial environments. It is supposed to be native to Europe while resident populations are also present in the US. We characterized a hierarchical sample of isolates from Europe and the US and conducted coalescent-, migration, and population genetic analysis of sequence and microsatellite data, to determine the pathways of spread and the demographic history of this pathogen. We found P. plurivora populations to be moderately diverse but not geographically structured. High levels of gene flow were observed within Europe and unidirectional from Europe to the US. Coalescent analyses revealed a signal of a recent expansion of the global P. plurivora population. Our study shows that P. plurivora has most likely been spread around the world by nursery trade of diseased plant material. In particular, P. plurivora was introduced into the US from Europe. International trade has allowed the pathogen to colonize new environments and/or hosts, resulting in population growth.

Variation of gene expression associated with colonisation of an anthropized environment: comparison between African and European populations of Drosophila simulans

The comparison of transcriptome profiles among populations is a powerful tool for investigating the role of gene expression change in adaptation to new environments. In this study, we use massively parallel sequencing of 3' cDNAs obtained from large samples of adult males, to compare a population of Drosophila simulans from a natural reserve within its ancestral range (eastern Africa) with a derived population collected in the strongly anthropized Rhône valley (France). The goal was to scan for adaptation linked to the invasion of new environments by the species. Among 15,090 genes retained for the analysis, 794 were found to be differentially expressed between the two populations. We observed an increase in expression of reproduction-related genes in eastern Africa, and an even stronger increase in expression of Cytochrome P450, Glutathione transferase and Glucuronosyl transferase genes in the derived population. These three gene families are involved in detoxification processes, which suggests that pesticides are a major environmental pressure for the species in this area. The survey of the Cyp6g1 upstream region revealed the insertion of a transposable element, Juan, in the regulatory sequence that is almost fixed in the Rhône Valley, but barely present in Mayotte. This shows that Cyp6g1 has undergone parallel evolution in derived populations of D. simulans as previously shown for D. melanogaster. The increasing amount of data produced by comparative population genomics and transcriptomics should permit the identification of additional genes associated with functional divergence among those differentially expressed.

A selective sweep across species boundaries in Drosophila

Adaptive mutations that accumulate during species divergence are likely to contribute to reproductive incompatibilities and hinder gene flow; however, there may also be a class of mutations that are generally advantageous and can spread across species boundaries. In this study, we characterize a 15 kb region on chromosome 3R that has introgressed from the cosmopolitan generalist species Drosophila simulans into the island endemic D. sechellia, which is an ecological specialist. The introgressed haplotype is fixed in D. sechellia over almost the entirety of the resequenced region, whereas a core region of the introgressed haplotype occurs at high frequency in D. simulans. The observed patterns of nucleotide variation and linkage disequilibrium are consistent with a recently completed selective sweep in D. sechellia and an incomplete sweep in D. simulans. Independent estimates of both the time to the introgression and sweep events are all close to 10,000 years before the present. Interestingly, the most likely target of selection is a highly occupied transcription factor binding region. This work confirms that it is possible for mutations to be globally advantageous, despite their occurrence in divergent genomic and ecological backgrounds.

Pupariation site preference within and between Drosophila sibling species

Holometabolous insects pass through a sedentary pupal stage and often choose a location for pupation that is different from the site of larval feeding. We have characterized a difference in pupariation site choice within and between sibling species of Drosophila. We found that, in nature, Drosophila sechellia pupariate within their host fruit, Morinda citrifolia, and that they perform this behavior in laboratory assays. In contrast, in the laboratory, geographically diverse strains of Drosophila simulans vary in their pupariation site preference; D. simulans lines from the ancestral range in southeast Africa pupariate on fruit, or a fruit substitute, whereas populations from Europe or the New World select sites off of fruit. We explored the genetic basis for the evolved preference in puariation site preference by performing quantitative trait locus mapping within and between species. We found that the interspecific difference is controlled largely by loci on chromosomes X and II. In contrast, variation between two strains of D. simulans appears to be highly polygenic, with the majority of phenotypic effects due to loci on chromosome III. These data address the genetic basis of how new traits arise as species diverge and populations disperse.

Genomic variation in natural populations of Drosophila melanogaster

This report of independent genome sequences of two natural populations of Drosophila melanogaster (37 from North America and 6 from Africa) provides unique insight into forces shaping genomic polymorphism and divergence. Evidence of interactions between natural selection and genetic linkage is abundant not only in centromere- and telomere-proximal regions, but also throughout the euchromatic arms. Linkage disequilibrium, which decays within 1 kbp, exhibits a strong bias toward coupling of the more frequent alleles and provides a high-resolution map of recombination rate. The juxtaposition of population genetics statistics in small genomic windows with gene structures and chromatin states yields a rich, high-resolution annotation, including the following: (1) 5'- and 3'-UTRs are enriched for regions of reduced polymorphism relative to lineage-specific divergence; (2) exons overlap with windows of excess relative polymorphism; (3) epigenetic marks associated with active transcription initiation sites overlap with regions of reduced relative polymorphism and relatively reduced estimates of the rate of recombination; (4) the rate of adaptive nonsynonymous fixation increases with the rate of crossing over per base pair; and (5) both duplications and deletions are enriched near origins of replication and their density correlates negatively with the rate of crossing over. Available demographic models of X and autosome descent cannot account for the increased divergence on the X and loss of diversity associated with the out-of-Africa migration. Comparison of the variation among these genomes to variation among genomes from D. simulans suggests that many targets of directional selection are shared between these species.

Abundant genetic variability in Drosophila simulans for hybrid female lethality in interspecific crosses to Drosophila melanogaster

Intrinsic postzygotic reproductive isolation is thought to result from the substitution of multiple harmless or beneficial genetic differences between species that are incidentally deleterious when combined in species hybrids, causing hybrid sterility or inviability. Genetic variability for hybrid sterility or inviability phenotypes is, however, rarely assessed in natural populations. Here, we assess variation for Drosophila simulans-encoded maternal factor(s) that cause lethality in D. simulans-Drosophila melanogaster F(1) hybrid females. First, we survey genetic variability in the strength of D. simulans-mediated maternal effect hybrid lethality among 37 geographic and laboratory isolates. We find abundant variability in the strength of maternal effect hybrid lethality, ranging from complete lethality to none. Second, we assess maternal effect hybrid lethality for a subset of wild isolates made heterozygous with two so-called hybrid rescue strains. The results suggest that the D. simulans maternal effect hybrid lethality involves a diversity of alleles and/or multiple loci.

Genome sequencing reveals complex speciation in the Drosophila simulans clade

The three species of the Drosophila simulans clade—the cosmopolitan species, D. simulans, and the two island endemic species, D. mauritiana and D. sechellia—are important models in speciation genetics, but some details of their phylogenetic and speciation history remain unresolved. The order and timing of speciation are disputed, and the existence, magnitude, and timing of gene flow among the three species remain unclear. Here we report on the analysis of a whole-genome four-species sequence alignment that includes all three D. simulans clade species as well as the D. melanogaster reference sequence. The alignment comprises novel, paired short-read sequence data from a single highly inbred line each from D. simulans, D. mauritiana, and D. sechellia. We are unable to reject a species phylogeny with a basal polytomy; the estimated age of the polytomy is 242,000 yr before the present. However, we also find that up to 4.6% of autosomal and 2.2% of X-linked regions have evolutionary histories consistent with recent gene flow between the mainland species (D. simulans) and the two island endemic species (D. mauritiana and D. sechellia). Our findings thus show that gene flow has occurred throughout the genomes of the D. simulans clade species despite considerable geographic, ecological, and intrinsic reproductive isolation. Last, our analysis of lineage-specific changes confirms that the D. sechellia genome has experienced a significant excess of slightly deleterious changes and a dearth of presumed favorable changes. The relatively reduced efficacy of natural selection in D. sechellia is consistent with its derived, persistently reduced historical effective population size.

Genetic diversity, population structure and Wolbachia infection status in a worldwide sample of Drosophila melanogaster and D. simulans populations

Drosophila melanogaster and its close relatives have been extremely important model species in the development of population genetic models that serve to explain patterns of diversity in natural populations, a major goal of evolutionary biology. A detailed picture of the evolutionary history of these species is beginning to emerge, as the relative importance of forces including demographic changes and natural selection is established. A continuing aim is to characterise levels of genetic diversity in a large number of populations of these species, covering a wide geographic area. We have used collections from five previously un-sampled wild populations of D. melanogaster and two of D. simulans, across three continents. We estimated levels of genetic diversity within, and divergence between, these populations, and looked for evidence of genetic structure both between ancestral and derived populations, and amongst derived populations. We also investigated the prevalence of infection with the bacterial endosymbiont Wolbachia. We found that D. melanogaster populations from Sub-Saharan Africa are the most diverse, and that divergence is highest between these and non-Sub-Saharan populations. There is strong evidence for structuring of populations between Sub-Saharan Africa and the rest of the world, and some evidence for weak structure amongst derived populations. Populations from Sub-Saharan Africa also differ in the prevalence of Wolbachia infection, with very low levels of infection compared to populations from the rest of the world.

Inter-island divergence within Drosophila mauritiana, a species of the D. simulans complex: Past history and/or speciation in progress?

Speciation with gene flow may be more common than generally thought, which makes detailed understanding of the extent and pattern of genetic divergence between geographically isolated populations useful. Species of the Drosophila simulans complex provide a good model for speciation and evolutionary studies, and hence understanding their population genetic structure will increase our understanding of the context in which speciation has occurred. Here, we describe genetic diversity and genetic differentiation of two distant populations of D. mauritiana (Mauritius and Rodrigues Islands) at mitochondrial and nuclear loci. We surveyed the two populations for their mitochondrial haplotypes, eight nuclear genes and 18 microsatellite loci. A new mitochondrial type is fixed in the Rodrigues population of D. mauritiana. The two populations are highly differentiated, their divergence appears relatively ancient (100,000 years) compared to the origin of the species, around 0.25MYA, and they exhibit very limited gene flow. However, they have similar levels of divergence from their sibling, D. simulans. Both nuclear genes and microsatellites revealed contrasting demographic histories between the two populations, expansion for the Mauritius population and stable population size for the Rodrigues Island population. The discovery of pronounced geographic structure within D. mauritiana combined to genetic structuring and low gene flow between the two island populations illuminates the evolutionary history of the species and clearly merits further attention in the broad context of speciation.

Rapid rise and fall of selfish sex-ratio X chromosomes in Drosophila simulans: spatiotemporal analysis of phenotypic and molecular data

Sex-ratio drive, which has been documented in several Drosophila species, is induced by X-linked segregation distorters. Contrary to Mendel's law of independent assortment, the sex-ratio chromosome (X(SR)) is inherited by more than half the offspring of carrier males, resulting in a female-biased sex ratio. This segregation advantage allows X(SR) to spread in populations, even if it is not beneficial for the carriers. In the cosmopolitan species D. simulans, the Paris sex-ratio is caused by recently emerged selfish X(SR) chromosomes. These chromosomes have triggered an intragenomic conflict, and their propagation has been halted over a large area by the evolution of complete drive suppression. Previous molecular population genetics analyses revealed a selective sweep indicating that the invasion of X(SR) chromosomes was very recent in Madagascar (likely less than 100 years ago). Here, we show that X(SR) chromosomes are now declining at this location as well as in Mayotte and Kenya. Drive suppression is complete in the three populations, which display little genetic differentiation and share swept haplotypes, attesting to a common and very recent ancestry of the X(SR) chromosomes. Patterns of DNA sequence variation also indicate a fitness cost of the segmental duplication involved in drive. The data suggest that X(SR) chromosomes started declining first on the African continent, then in Mayotte, and finally in Madagascar and strongly support a scenario of rapid cycling of X chromosomes. Once drive suppression has evolved, standard X(ST) chromosomes locally replace costly X(SR) chromosomes in a few decades.

Population transcriptomics: insights from Drosophila simulans, Drosophila sechellia and their hybrids

Sequence differentiation has been widely studied between populations and species, whereas interest in expression divergence is relatively recent. Using microarrays, we compared four geographically distinct populations of Drosophila simulans and a population of Drosophila sechellia, and interspecific hybrids. We observed few differences between populations, suggesting a slight population structure in D. simulans. This structure was observed in direct population comparisons, as well as in interspecific comparisons (hybrids vs. parents, D. sechellia vs. D. simulans). Expression variance is higher in the French and Zimbabwean populations than in the populations from the ancestral range of D. simulans (Kenya and Seychelles). This suggests a large scale phenomenon of decanalization following the invasion of a new environment. Comparing D. simulans and D. sechellia, we revealed 304 consistently differentially expressed genes, with striking overrepresentation of genes of the cytochrome P450 family, which could be related to their role in detoxification as well as in hormone regulation. We also revealed differences in genes involved in Juvenile hormone and Dopamine differentiation. We finally observed very few differentially expressed genes between hybrids and parental populations, with an overrepresentation of X-linked genes.

Recent and recurrent selective sweeps of the antiviral RNAi gene Argonaute-2 in three species of Drosophila

Antagonistic host–parasite interactions can drive rapid adaptive evolution in genes of the immune system, and such arms races may be an important force shaping polymorphism in the genome. The RNA interference pathway gene Argonaute-2 (AGO2) is a key component of antiviral defense in Drosophila, and we have previously shown that genes in this pathway experience unusually high rates of adaptive substitution. Here we study patterns of genetic variation in a 100-kbp region around AGO2 in three different species of Drosophila. Our data suggest that recent independent selective sweeps in AGO2 have reduced genetic variation across a region of more than 50 kbp in Drosophila melanogaster, D. simulans, and D. yakuba, and we estimate that selection has fixed adaptive substitutions in this gene every 30–100 thousand years. The strongest signal of recent selection is evident in D. simulans, where we estimate that the most recent selective sweep involved an allele with a selective advantage of the order of 0.5–1% and occurred roughly 13–60 Kya. To evaluate the potential consequences of the recent substitutions on the structure and function of AGO2, we used fold-recognition and homology-based modeling to derive a structural model for the Drosophila protein, and this suggests that recent substitutions in D. simulans are overrepresented at the protein surface. In summary, our results show that selection by parasites can consistently target the same genes in multiple species, resulting in areas of the genome that have markedly reduced genetic diversity.

Populations, hybrids and the systematic concepts of species and subspecies in Chagas disease triatomine vectors inferred from nuclear ribosomal and mitochondrial DNA

In Chagas disease, triatomine vectors are the main target for control measures because of the absence of effective drugs. The broad usefulness of nuclear rDNA and mtDNA sequences explains why triatomine studies using these markers have increased so pronouncedly in recent years. This indicates the appropriateness of an updated review about these molecular markers, concentrating on aspects useful for research on Chagas disease vectors. A comparative analysis is presented on the efficiency, weight of their different characteristics, limitations and problems of each of the different DNA markers in the light of the results obtained in studies on populations, hybrids, subspecies and species of the subfamily Triatominae. The use of a standardized composite haplotype code nomenclature for both nuclear rDNA and mtDNA markers is strongly encouraged to avoid difficulties in comparative studies. Triatomine aspects related to concerted evolution, microsatellites, minisatellites and insertions/deletions in nuclear rDNA and silent/non-silent mutations, pseudogenes and weaknesses of partial sequences in mtDNA are analysed. Introgression and hybrids, nuclear and mitochondrial DNA strengths, and compared evolutionary rates of nuclear rDNA and mtDNA in triatomines are discussed. Many conclusions are obtained thanks to the availability, for the first time in triatomines, of a complete sequence of a protein-coding mtDNA gene as ND1 from very numerous triatomine species covering from different populations of a species up to members belonging to different tribes. The evolutionary rates of each nuclear rDNA marker and mtDNA marker are analysed by comparison at subspecies level (intrapopulational, interpopulational, between morphs, and between subspecies) and species level (close and distant species of the same genus, species of different genera, and species of different tribes). Weaknesses of mtDNA for systematic-taxonomic purposes detected recently and newly in insects and triatomines, respectively, are discussed in detail. Emphasis is given to taxonomic units and biological entities presenting well-known problematics, both from the systematic-taxonomic and/or epidemiological-control points of view, as well as to molecular situations which can give rise to erroneous conclusions. All these aspects constitute the background on which the key question about the systematic concepts of species and subspecies in triatomines is focused. The global purpose is to facilitate future work on triatomines by highlighting present gaps, how better choice the appropriate markers, and marker aspects which should be taken into account. Key characteristics as alpha, CI and transformation rate matrices ought to be obtained and noted to get appropriate results and allow correct interpretations. The main aim is to offer a baseline for future fundamental research on triatomines and applied research on transmission, epidemiology and control measures related to Chagas disease vectors.

Introgression as a likely cause of mtDNA paraphyly in two allopatric skippers (Lepidoptera: Hesperiidae)

Gene transfer between species during interspecific hybridization is a widely accepted reality in plants but is considered a relatively rare phenomenon among animals. Here we describe a unique case of mitochondrial DNA (mtDNA) paraphyly in the skipper genus, Erynnis, that involves well-diverged allopatric species. Using molecular evidence from both mitochondrial and nuclear genomes, we found high levels of intraspecific divergence in the mitochondrial genome within E. propertius (over 4% pair-wise sequence divergence) but no such differentiation in the nuclear genome. Sequence comparisons with related Erynnis suggest that past, but recent and infrequent introgression between E. propertius and E. horatius is the most reasonable explanation for the observed pattern of mtDNA paraphyly. This example of putative introgression highlights the complexity of mtDNA evolution and suggests that similar processes could be operating in other taxa that have not been extensively sampled. Our observations reinforce the importance of involving multiple genes with different modes of inheritance in the analysis of population history of congeneric taxa.

Non-neutral processes drive the nucleotide composition of non-coding sequences in Drosophila

The nature of the forces affecting base composition is a key question in genome evolution. There is uncertainty as to whether differences in the GC contents of non-coding sequences reflect differences in mutational bias, or in the intensity of selection or biased gene conversion. We have used a polymorphism dataset for non-coding sequences on the X chromosome of Drosophila simulans to examine this question. The proportion of GC→AT versus AT→GC polymorphic mutations in a locus is correlated with its GC content. This implies the action of forces that favour GC over AT base pairs, which are apparently strongest in GC-rich sequences.

Positive selection at the binding sites of the male-specific lethal complex involved in dosage compensation in Drosophila

In many taxa, males and females differ with respect to their sex chromosomes, and dosage compensation mechanisms have evolved to equalize X-linked gene transcription. In Drosophila, the male-specific lethal (MSL) complex binds to hundreds of sites along the male X chromosome and mediates twofold hypertranscription of the single male X. Two recent studies found evidence for lineage-specific adaptive evolution in all five core protein-coding genes of the MSL complex in Drosophila melanogaster. In particular, dramatic positive selection was detected in domains shown to be responsible for their specific targeting to the X chromosome. Here I use population genetics to show that three previously characterized MSL-binding DNA segments on the X themselves underwent adaptive evolution in D. melanogaster, but not in its close relatives D. simulans and D. yakuba. MSL components have been shown to not correctly target the D. melanogaster X chromosome in hybrids between D. melanogaster and D. simulans. My finding supports the idea of selection-driven coevolution among DNA-protein interactions of the dosage compensation machinery and suggests that misregulated dosage compensation could contribute to male hybrid inviability in Drosophila.

Positive and negative selection on noncoding DNA in Drosophila simulans

There is now a wealth of evidence that some of the most important regions of the genome are found outside those that encode proteins, and noncoding regions of the genome have been shown to be subject to substantial levels of selective constraint, particularly in Drosophila. Recent work has suggested that these regions may also have been subject to the action of positive selection, with large fractions of noncoding divergence having been driven to fixation by adaptive evolution. However, this work has focused on Drosophila melanogaster, which is thought to have experienced a reduction in effective population size (N(e)), and thus a reduction in the efficacy of selection, compared with its closest relative Drosophila simulans. Here, we examine patterns of evolution at several classes of noncoding DNA in D. simulans and find that all noncoding DNA is subject to the action of negative selection, indicated by reduced levels of polymorphism and divergence and a skew in the frequency spectrum toward rare variants. We find that the signature of negative selection on noncoding DNA and nonsynonymous sites is obscured to some extent by purifying selection acting on preferred to unpreferred synonymous codon mutations. We investigate the extent to which divergence in noncoding DNA is inferred to be the product of positive selection and to what extent these inferences depend on selection on synonymous sites and demography. Based on patterns of polymorphism and divergence for different classes of synonymous substitution, we find the divergence excess inferred in noncoding DNA and nonsynonymous sites in the D. simulans lineage difficult to reconcile with demographic explanations.

African Drosophila melanogaster and D. simulans populations have similar levels of sequence variability, suggesting comparable effective population sizes

Drosophila melanogaster and D. simulans are two closely related species with a similar distribution range. Many studies suggested that D. melanogaster has a smaller effective population size than D. simulans. As most evidence was derived from non-African populations, we readdressed this question by sequencing 10 X-linked loci in five African D. simulans and six African D. melanogaster populations. Contrary to previous results, we found no evidence for higher variability, and thus larger effective population size, in D. simulans. Our observation of similar levels of variability of both species will have important implications for the interpretation of patterns of molecular evolution.

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.

Population genomics, the study of genome-wide patterns of sequence variation within and between closely related species, can provide a comprehensive view of the relative importance of mutation, recombination, natural selection, and genetic drift in evolution. It can also provide fundamental insights into the biological attributes of organisms that are specifically shaped by adaptive evolution. One approach for generating population genomic datasets is to align DNA sequences from whole-genome shotgun projects to a standard reference sequence. We used this approach to carry out whole-genome analysis of polymorphism and divergence in Drosophila simulans, a close relative of the model system, D. melanogaster. We find that polymorphism and divergence fluctuate on a large scale across the genome and that these fluctuations are probably explained by natural selection rather than by variation in mutation rates. Our analysis suggests that adaptive protein evolution is common and is often related to biological processes that may be associated with gene expression, chromosome biology, and reproduction. The approaches presented here will have broad applicability to future analysis of population genomic variation in other systems, including humans.

Low-coverage genome sequences from multiple Drosophila simulans strains provide the first comprehensive view of polymorphism and divergence in the fruit fly.

Mitochondrial DNA variation is associated with measurable differences in life-history traits and mitochondrial metabolism in Drosophila simulans

Recent studies have used a variety of theoretical arguments to show that mitochondrial (mt) DNA rarely evolves as a strictly neutral marker and that selection operates on the mtDNA of many species. However, the vast majority of researchers are not convinced by these arguments because data linking mtDNA variation with phenotypic differences are limited. We investigated sequence variation in the three mtDNA and nine nuclear genes (including all isoforms) that encode the 12 subunits of cytochrome c oxidase of the electron transport chain in Drosophila. We then studied cytochrome c oxidase activity as a key aspect of mitochondrial bioenergetics and four life-history traits. In Drosophila simulans, sequence data from the three mtDNA encoded cytochrome c oxidase genes show that there are 76 synonymous and two nonsynonymous fixed differences among flies harboring siII compared with siIII mtDNA. In contrast, 13 nuclear encoded genes show no evidence of genetic subdivision associated with the mtDNA. Flies with siIII mtDNA had higher cytochrome c oxidase activity and were more starvation resistant. Flies harboring siII mtDNA had greater egg size and fecundity, and recovered faster from cold coma. These data are consistent with a causative role for mtDNA variation in these phenotypic differences, but we cannot completely rule out the involvement of nuclear genes. The results of this study have significant implications for the use of mtDNA as an assumed neutral marker and show that evolutionary shifts can involve changes in mtDNA despite the small number of genes encoded in the organelle genome.

Sympatric Drosophila simulans flies with distinct mtDNA show age related differences in mitochondrial metabolism

The primary causes of age-related changes in mitochondrial metabolism are not known. The goal of this study is to document the influence of naturally occurring mtDNA variation on age dependent changes in mitochondrial respiration, hydrogen peroxide (H(2)O(2)) generation and antioxidant defenses in the fly Drosophila simulans. Possible changes include an increase in rates of reactive oxygen species production with age and/or an age dependent decrease in antioxidant response. For this study we have used flies harboring distinct siII and siIII mtDNA types. Previously we have shown that males harboring siII mtDNA had higher rates of mitochondrial H(2)O(2) production from complex III at 11d compared to males with the siIII mtDNA type. Here, we corroborate those results and show that Drosophila harboring the siII and siIII mtDNA types exhibit significantly different patterns of pro-oxidant and antioxidant activities as they age. Flies harboring siII mtDNA had higher rates of mitochondrial H(2)O(2) production and manganese superoxide dismutase activity at 11 and 18d of age than siIII mtDNA harboring flies. Copper-zinc superoxide dismutase activity increased from 11 to 25d in siII flies while the accumulation of oxidized glutathione did not change between 11 and 25d. In contrast, siIII harboring flies showed an age dependent increase in H(2)O(2) production, reaching higher production rates on day 25 than that observed in siII flies. Copper-zinc superoxide dismutase activities did not change between 11 and 25d while the oxidized glutathione accumulation increased with age. The results show antioxidant levels correlate with pro-oxidant levels in siII but not siIII flies. These results demonstrate our ability to correlate mtDNA variation with differences in whole mitochondrial physiology and individual complex biochemistry.

A screen for immunity genes evolving under positive selection in Drosophila

Genes involved in the immune system tend to have higher rates of adaptive evolution than other genes in the genome, probably because they are coevolving with pathogens. We have screened a sample of Drosophila genes to identify those evolving under positive selection. First, we identified rapidly evolving immunity genes by comparing 140 loci in Drosophila erecta and D. yakuba. Secondly, we resequenced 23 of the fastest evolving genes from the independent species pair D. melanogaster and D. simulans, and identified those under positive selection using a McDonald-Kreitman test. There was strong evidence of adaptive evolution in two serine proteases (persephone and spirit) and a homolog of the Anopheles serpin SRPN6, and weaker evidence in another serine protease and the death domain protein dFADD. These results add to mounting evidence that immune signalling pathway molecules often evolve rapidly, possibly because they are sites of host-parasite coevolution.

Sympatric Drosophila simulans flies with distinct mtDNA show difference in mitochondrial respiration and electron transport

The role of mitochondrial DNA (mtDNA) in mitochondrial metabolism is understudied yet humans harboring specific mtDNA types age at dissimilar rates, are unequally susceptible to various diseases, and differentially adapt to various environmental conditions. This study compares mitochondrial respiration, proton leak and electron transport of Drosophila simulans males with distinct mtDNA haplogroups (siII and -III) that were collected in sympatry in Kenya. Despite the large divergence among haplogroups there is very low intrahaplogroup variation and no correlated variation in the nuclear genome has been detected. We show that repeatable bioenergetic differences exist between 11d old males harboring siII and siIII mtDNA. Males with siIII mtDNA showed higher (i) state 3 respiration rates from isolated mitochondria for both complex I and complex III based substrates, and (ii) complex IV (cytochrome c oxidase) activity. Males harboring siIII mtDNA had lower (i) hydrogen peroxide formation by both complexes I and III, (ii) proton leak from isolated mitochondria, (iii) mitochondrial ATPase activity, and (iv) mitochondrial cytochrome content. In combination, the results suggest that mitochondria isolated from siIII mtDNA harboring males have more efficient metabolism than siII mtDNA harboring males.

An Andean origin of Phytophthora infestans inferred from mitochondrial and nuclear gene genealogies

Phytophthora infestans (Mont.) de Bary caused the 19th century Irish Potato Famine. We assessed the genealogical history of P. infestans using sequences from portions of two nuclear genes (beta-tubulin and Ras) and several mitochondrial loci P3, (rpl14, rpl5, tRNA) and P4 (Cox1) from 94 isolates from South, Central, and North America, as well as Ireland. Summary statistics, migration analyses and the genealogy of current populations of P. infestans for both nuclear and mitochondrial loci are consistent with an "out of South America" origin for P. infestans. Mexican populations of P. infestans from the putative center of origin in Toluca Mexico harbored less nucleotide and haplotype diversity than Andean populations. Coalescent-based genealogies of all loci were congruent and demonstrate the existence of two lineages leading to present day haplotypes of P. infestans on potatoes. The oldest lineage associated with isolates from the section Anarrhichomenun including Solanum tetrapetalum from Ecuador was identified as Phytophthora andina and evolved from a common ancestor of P. infestans. Nuclear and mitochondrial haplotypes found in Toluca Mexico were derived from only one of the two lineages, whereas haplotypes from Andean populations in Peru and Ecuador were derived from both lineages. Haplotypes found in populations from the U.S. and Ireland was derived from both ancestral lineages that occur in South America suggesting a common ancestry among these populations. The geographic distribution of mutations on the rooted gene genealogies demonstrate that the oldest mutations in P. infestans originated in South America and are consistent with a South American origin.

Microsatellite variation and differentiation in African and non-African populations of Drosophila simulans

Drosophila simulans originated in sub-Saharan Africa or Madagascar and colonized the rest of the world after the last glaciation about 10 000 years ago. Consistent with this demographic history, sub-Saharan African populations have been shown to harbour higher levels of microsatellite and sequence variation than cosmopolitan populations. Nevertheless, only limited information is available on the population structure of D. simulans. Here, we analysed X-linked and autosomal microsatellite loci in four sub-Saharan African, one North African, one Israeli, and two European D. simulans populations. Bayesian clustering algorithms combined the North African, Israeli, and European populations into a single cosmopolitan group. The four sub-Saharan populations were split into two separate groups. Pairwise F(ST) analysis, however, indicated significant population differentiation between all eight populations surveyed. A significant signal for population reduction in cosmopolitan populations was found only for X-linked loci.

Pervasive adaptive evolution among interactors of the Drosophila hybrid inviability gene, Nup96

Nup96 is involved in a lethal hybrid incompatibility between 2 fruit fly species, Drosophila melanogaster and Drosophila simulans. Recurrent adaptive evolution drove the rapid functional divergence of Nup96 in both the D. melanogaster and the D. simulans lineages. Functional divergence of Nup96 between these 2 species is unexpected as Nup96 encodes part of the Nup107 subcomplex, an architectural component of nuclear pore complexes, the macromolecular channels in nuclear envelopes that mediate nucleocytoplasmic traffic in all eukaryotes. Here we study the evolutionary histories of 5 of Nup96's protein interactors--3 stable Nup107 subcomplex proteins (Nup75, Nup107, and Nup133) and 2 mobile nucleoporins (Nup98 and Nup153)--and show that all 5 have experienced recurrent adaptive evolution. These results are consistent with selection-driven coevolution among molecular interactors within species causing the incidental evolution of incompatible interactions seen in hybrids between species. We suggest that genetic conflict-driven processes may have contributed to the rapid molecular evolution of Nup107 subcomplex genes.

History and structure of sub-Saharan populations of Drosophila melanogaster

Drosophila melanogaster is an important model organism in evolutionary genetics, yet little is known about the population structure and the demographic history of this species within sub-Saharan Africa, which is thought to contain its ancestral range. We surveyed nucleotide variation at four 1-kb fragments in 240 individual lines representing 21 sub-Saharan and 4 Palearctic population samples of D. melanogaster. In agreement with recent studies, we find a small but significant level of genetic differentiation within sub-Saharan Africa. A clear geographic pattern is observed, with eastern and western African populations composing two genetically distinct groups. This pattern may have resulted from a relatively recent establishment of D. melanogaster in western Africa. Eastern populations show greater evidence for long-term stability, consistent with the hypothesis that eastern Africa contains the ancestral range of the species. Three sub-Saharan populations show evidence for cosmopolitan introgression. Apart from those cases, the closest relationships between Palearctic and sub-Saharan populations involve a sample from the rift zone (Uganda), suggesting that the progenitors of Palearctic D. melanogaster might have come from this region. Finally, we find a large excess of singleton polymorphisms in the full data set, which is best explained by a combination of population growth and purifying selection.

A Wolbachia-associated fitness benefit depends on genetic background in Drosophila simulans

The alpha-proteobacteria Wolbachia infect a number of insect species and influence host reproduction to favour the spread of infected females through a population. The fitness effect of this infection is important in understanding the spread and maintenance of Wolbachia within and among host populations. However, a full elucidation of fitness effect requires careful control of host genetic background. Here, I transferred a single clone of Wolbachia (the wHa strain) into three genetically distinct isofemale lines of the fly Drosophila simulans using microinjection methodology. These lines carried one of the three described mitochondrial haplogroups (siI, siII or siIII) and differ in nuclear genome as well. Population cage assays showed that wHa-infected siIII flies enjoyed a dramatic fitness benefit compared to uninfected siIII. In contrast, wHa did not affect the fitness of siI or siII flies. This study points to the importance of host-by-symbiont interaction terms that may play an important role in organismal-fitness.

Origination of an X-linked testes chimeric gene by illegitimate recombination in Drosophila

The formation of chimeric gene structures provides important routes by which novel proteins and functions are introduced into genomes. Signatures of these events have been identified in organisms from wide phylogenic distributions. However, the ability to characterize the early phases of these evolutionary processes has been difficult due to the ancient age of the genes or to the limitations of strictly computational approaches. While examples involving retrotransposition exist, our understanding of chimeric genes originating via illegitimate recombination is limited to speculations based on ancient genes or transfection experiments. Here we report a case of a young chimeric gene that has originated by illegitimate recombination in Drosophila. This gene was created within the last 2–3 million years, prior to the speciation of Drosophila simulans, Drosophila sechellia, and Drosophila mauritiana. The duplication, which involved the Bällchen gene on Chromosome 3R, was partial, removing substantial 3′ coding sequence. Subsequent to the duplication onto the X chromosome, intergenic sequence was recruited into the protein-coding region creating a chimeric peptide with ~ 33 new amino acid residues. In addition, a novel intron-containing 5′ UTR and novel 3′ UTR evolved. We further found that this new X-linked gene has evolved testes-specific expression. Following speciation of the D. simulans complex, this novel gene evolved lineage-specifically with evidence for positive selection acting along the D. simulans branch.

Illegitimate recombination, the non-homologous recombination that occurs between DNA sequences with few or no identical nucleotides, is a general phenomenon that has been known to cause many medically important deleterious changes. However, little is known about the positive side of such a process. For example, little is known about its relative role in the origin of new gene functions that confer increased fitness to species. This work contributes to the understanding of the significance of this process. Here the authors report on a young chimeric gene that has originated by illegitimate recombination in Drosophila. The term “chimeric gene” refers to gene structures—both coding and noncoding—which have been generated from distinct parental loci. This chimeric gene was created within the last 2–3 million years, prior to the speciation of Drosophila simulans, Drosophila sechellia, and Drosophila mauritiana. A gene on Chromosome 3R was duplicated onto the X chromosome and recruited intergenic sequence, creating a chimeric peptide. It was found that this new X-linked gene has evolved testes-specific expression. Following speciation of the D. simulans complex, this novel gene evolved lineage-specifically under positive Darwinian selection.

Contrasted polymorphism patterns in a large sample of populations from the evolutionary genetics model Drosophila simulans

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.