A selective sweep across species boundaries in Drosophila

IntroUNET: Identifying introgressed alleles via semantic segmentation

A growing body of evidence suggests that gene flow between closely related species is a widespread phenomenon. Alleles that introgress from one species into a close relative are typically neutral or deleterious, but sometimes confer a significant fitness advantage. Given the potential relevance to speciation and adaptation, numerous methods have therefore been devised to identify regions of the genome that have experienced introgression. Recently, supervised machine learning approaches have been shown to be highly effective for detecting introgression. One especially promising approach is to treat population genetic inference as an image classification problem, and feed an image representation of a population genetic alignment as input to a deep neural network that distinguishes among evolutionary models (i.e. introgression or no introgression). However, if we wish to investigate the full extent and fitness effects of introgression, merely identifying genomic regions in a population genetic alignment that harbor introgressed loci is insufficient-ideally we would be able to infer precisely which individuals have introgressed material and at which positions in the genome. Here we adapt a deep learning algorithm for semantic segmentation, the task of correctly identifying the type of object to which each individual pixel in an image belongs, to the task of identifying introgressed alleles. Our trained neural network is thus able to infer, for each individual in a two-population alignment, which of those individual's alleles were introgressed from the other population. We use simulated data to show that this approach is highly accurate, and that it can be readily extended to identify alleles that are introgressed from an unsampled "ghost" population, performing comparably to a supervised learning method tailored specifically to that task. Finally, we apply this method to data from Drosophila, showing that it is able to accurately recover introgressed haplotypes from real data. This analysis reveals that introgressed alleles are typically confined to lower frequencies within genic regions, suggestive of purifying selection, but are found at much higher frequencies in a region previously shown to be affected by adaptive introgression. Our method's success in recovering introgressed haplotypes in challenging real-world scenarios underscores the utility of deep learning approaches for making richer evolutionary inferences from genomic data.

IntroUNET: identifying introgressed alleles via semantic segmentation

A growing body of evidence suggests that gene flow between closely related species is a widespread phenomenon. Alleles that introgress from one species into a close relative are typically neutral or deleterious, but sometimes confer a significant fitness advantage. Given the potential relevance to speciation and adaptation, numerous methods have therefore been devised to identify regions of the genome that have experienced introgression. Recently, supervised machine learning approaches have been shown to be highly effective for detecting introgression. One especially promising approach is to treat population genetic inference as an image classification problem, and feed an image representation of a population genetic alignment as input to a deep neural network that distinguishes among evolutionary models (i.e. introgression or no introgression). However, if we wish to investigate the full extent and fitness effects of introgression, merely identifying genomic regions in a population genetic alignment that harbor introgressed loci is insufficient-ideally we would be able to infer precisely which individuals have introgressed material and at which positions in the genome. Here we adapt a deep learning algorithm for semantic segmentation, the task of correctly identifying the type of object to which each individual pixel in an image belongs, to the task of identifying introgressed alleles. Our trained neural network is thus able to infer, for each individual in a two-population alignment, which of those individual's alleles were introgressed from the other population. We use simulated data to show that this approach is highly accurate, and that it can be readily extended to identify alleles that are introgressed from an unsampled "ghost" population, performing comparably to a supervised learning method tailored specifically to that task. Finally, we apply this method to data from Drosophila, showing that it is able to accurately recover introgressed haplotypes from real data. This analysis reveals that introgressed alleles are typically confined to lower frequencies within genic regions, suggestive of purifying selection, but are found at much higher frequencies in a region previously shown to be affected by adaptive introgression. Our method's success in recovering introgressed haplotypes in challenging real-world scenarios underscores the utility of deep learning approaches for making richer evolutionary inferences from genomic data.

Divergence and introgression among the virilis group of Drosophila

Speciation with gene flow is now widely regarded as common. However, the frequency of introgression between recently diverged species and the evolutionary consequences of gene flow are still poorly understood. The virilis group of Drosophila contains 12 species that are geographically widespread and show varying levels of prezygotic and postzygotic isolation. Here, we use de novo genome assemblies and whole-genome sequencing data to resolve phylogenetic relationships and describe patterns of introgression and divergence across the group. We suggest that the virilis group consists of three, rather than the traditional two, subgroups. Some genes undergoing rapid sequence divergence across the group were involved in chemical communication and desiccation tolerance, and may be related to the evolution of sexual isolation and adaptation. We found evidence of pervasive phylogenetic discordance caused by ancient introgression events between distant lineages within the group, and more recent gene flow between closely related species. When assessing patterns of genome-wide divergence in species pairs across the group, we found no consistent genomic evidence of a disproportionate role for the X chromosome as has been found in other systems. Our results show how ancient and recent introgressions confuse phylogenetic reconstruction, but may play an important role during early radiation of a group.

Dynamics and Impacts of Transposable Element Proliferation in the Drosophila nasuta Species Group Radiation

Transposable element (TE) mobilization is a constant threat to genome integrity. Eukaryotic organisms have evolved robust defensive mechanisms to suppress their activity, yet TEs can escape suppression and proliferate, creating strong selective pressure for host defense to adapt. This genomic conflict fuels a never-ending arms race that drives the rapid evolution of TEs and recurrent positive selection of genes involved in host defense; the latter has been shown to contribute to postzygotic hybrid incompatibility. However, how TE proliferation impacts genome and regulatory divergence remains poorly understood. Here, we report the highly complete and contiguous (N50 = 33.8-38.0 Mb) genome assemblies of seven closely related Drosophila species that belong to the nasuta species group-a poorly studied group of flies that radiated in the last 2 My. We constructed a high-quality de novo TE library and gathered germline RNA-seq data, which allowed us to comprehensively annotate and compare TE insertion patterns between the species, and infer the evolutionary forces controlling their spread. We find a strong negative association between TE insertion frequency and expression of genes nearby; this likely reflects survivor bias from reduced fitness impact of TEs inserting near lowly expressed, nonessential genes, with limited TE-induced epigenetic silencing. Phylogenetic analyses of insertions of 147 TE families reveal that 53% of them show recent amplification in at least one species. The most highly amplified TE is a nonautonomous DNA element (Drosophila INterspersed Element; DINE) which has gone through multiple bouts of expansions with thousands of full-length copies littered throughout each genome. Across all TEs, we find that TEs expansions are significantly associated with high expression in the expanded species consistent with suppression escape. Thus, whereas horizontal transfer followed by the invasion of a naïve genome has been highlighted to explain the long-term survival of TEs, our analysis suggests that evasion of host suppression of resident TEs is a major strategy to persist over evolutionary times. Altogether, our results shed light on the heterogenous and context-dependent nature in which TEs affect gene regulation and the dynamics of rampant TE proliferation amidst a recently radiated species group.

Widespread introgression across a phylogeny of 155 Drosophila genomes

Genome-scale sequence data have invigorated the study of hybridization and introgression, particularly in animals. However, outside of a few notable cases, we lack systematic tests for introgression at a larger phylogenetic scale across entire clades. Here, we leverage 155 genome assemblies from 149 species to generate a fossil-calibrated phylogeny and conduct multilocus tests for introgression across 9 monophyletic radiations within the genus Drosophila. Using complementary phylogenomic approaches, we identify widespread introgression across the evolutionary history of Drosophila. Mapping gene-tree discordance onto the phylogeny revealed that both ancient and recent introgression has occurred across most of the 9 clades that we examined. Our results provide the first evidence of introgression occurring across the evolutionary history of Drosophila and highlight the need to continue to study the evolutionary consequences of hybridization and introgression in this genus and across the tree of life.

Drosophila sechellia: A Genetic Model for Behavioral Evolution and Neuroecology

Defining the mechanisms by which animals adapt to their ecological niche is an important problem bridging evolution, genetics, and neurobiology. We review the establishment of a powerful genetic model for comparative behavioral analysis and neuroecology, Drosophila sechellia. This island-endemic fly species is closely related to several cosmopolitan generalists, including Drosophila melanogaster, but has evolved extreme specialism, feeding and reproducing exclusively on the noni fruit of the tropical shrub Morinda citrifolia. We first describe the development and use of genetic approaches to facilitate genotype/phenotype associations in these drosophilids. Next, we survey the behavioral, physiological, and morphological adaptations of D. sechellia throughout its life cycle and outline our current understanding of the genetic and cellular basis of these traits. Finally, we discuss the principles this knowledge begins to establish in the context of host specialization, speciation, and the neurobiology of behavioral evolution and consider open questions and challenges in the field.

P-elements strengthen reproductive isolation within the Drosophila simulans species complex

Determining mechanisms that underlie reproductive isolation (RI) is key to understanding how species boundaries are maintained in nature. Transposable elements (TEs) are ubiquitous across eukaryotic genomes. However, the role of TEs in modulating the strength of RI between species is poorly understood. Several species of Drosophila have been found to harbor P-elements (PEs), yet only D. simulans is known to be currently polymorphic for their presence in wild populations. PEs can cause RI between PE-containing (P) and PE-lacking (M) lineages of the same species. However, it is unclear whether they also contribute to the magnitude of RI between species. Here, we use the simulans species complex to assess whether differences in PE status between D. simulans and its sister species, which do not harbor PEs, contribute to multiple barriers to gene flow between species. We show that crosses involving a P D. simulans father and an M mother from a sister species exhibit lower F1 female fecundity than crosses involving an M D. simulans father and an M sister-species mother. We also find that another TE, I-element, might play a minor role in determining the frequency of dysgenesis between species. Our results suggest that the presence of PEs in a species can strengthen isolation from its sister species, providing evidence that TEs can play a role in RI.

Hybridization, Behavioral Patterns, and Pre- and Postzygotic Isolation Between Two Recently Derived Species of Walnut-Infesting Rhagoletis Fruit Flies in the Highlands of Mexico

Episodes of isolation and secondary contact among populations of insects of Nearctic origin during Pleistocene glacial/postglacial climatic cycles had a strong evolutionary influence on the diversity of flies in the genus Rhagoletis in mountainous areas of Mexico. As a series of experiments undertaken to gather support for phylogenetic hypotheses on the origin of three walnut-infesting species in the suavis group, we examined pre- and postzygotic isolation between Rhagoletis completa Cresson, 1929 and R. ramosae Hernández-Ortiz, 1985. Despite morphological, biological, and behavioral differences, these two species were found to be capable of hybridization. Mating experiments in large enclosures revealed asymmetric sexual isolation. There were notable differences in male sexual behavior. While R. ramosae males mated exclusively on host fruit, R. completa males used fruit and alternative mating locations. During fruit-guarding and male-male contests, R. completa and R. ramosae males adopted markedly different wing postures. R. completa females were more reluctant to copulate with heterospecific males than R. ramosae females. During no choice crosses in small enclosures, there was a reduction of egg hatch for the hybrid cross of R. completa males × R. ramosae females. Our results and previous studies on reproductive isolation between other species pairs in the suavis group support a clade in which R. ramosae, R. zoqui Bush, 1966, and R. completa are close relatives all still capable of hybridizing.

Hybrid Sterility, Genetic Conflict and Complex Speciation: Lessons From the Drosophila simulans Clade Species

The three fruitfly species of the Drosophila simulans clade- D. simulans, D. mauritiana, and D. sechellia- have served as important models in speciation genetics for over 40 years. These species are reproductively isolated by geography, ecology, sexual signals, postmating-prezygotic interactions, and postzygotic genetic incompatibilities. All pairwise crosses between these species conform to Haldane's rule, producing fertile F1 hybrid females and sterile F1 hybrid males. The close phylogenetic proximity of the D. simulans clade species to the model organism, D. melanogaster, has empowered genetic analyses of their species differences, including reproductive incompatibilities. But perhaps no phenotype has been subject to more continuous and intensive genetic scrutiny than hybrid male sterility. Here we review the history, progress, and current state of our understanding of hybrid male sterility among the D. simulans clade species. Our aim is to integrate the available information from experimental and population genetics analyses bearing on the causes and consequences of hybrid male sterility. We highlight numerous conclusions that have emerged as well as issues that remain unresolved. We focus on the special role of sex chromosomes, the fine-scale genetic architecture of hybrid male sterility, and the history of gene flow between species. The biggest surprises to emerge from this work are that (i) genetic conflicts may be an important general force in the evolution of hybrid incompatibility, (ii) hybrid male sterility is polygenic with contributions of complex epistasis, and (iii) speciation, even among these geographically allopatric taxa, has involved the interplay of gene flow, negative selection, and positive selection. These three conclusions are marked departures from the classical views of speciation that emerged from the modern evolutionary synthesis.

Pure species discriminate against hybrids in the Drosophila melanogaster species subgroup

Introgression, the exchange of alleles between species, is a common event in nature. This transfer of alleles between species must happen through fertile hybrids. Characterizing the traits that cause defects in hybrids illuminates how and when gene flow is expected to occur. Inviability and sterility are extreme examples of fitness reductions but are not the only type of defects in hybrids. Some traits specific to hybrids are more subtle but are important to determine their fitness. In this report, we study whether F1 hybrids between two species pairs of Drosophila are as attractive as the parental species. We find that in both species pairs, the sexual attractiveness of the F1 hybrids is reduced and that pure species discriminate strongly against them. We also find that the cuticular hydrocarbon (CHC) profile of the female hybrids is intermediate between the parental species. Perfuming experiments show that modifying the CHC profile of the female hybrids to resemble pure species improves their chances of mating. Our results show that behavioral discrimination against hybrids might be an important component of the persistence of species that can hybridize.

Bioinformatic Detection of Positive Selection Pressure in Plant Pathogens: The Neutral Theory of Molecular Sequence Evolution in Action

The genomes of plant pathogenic fungi and oomycetes are often exposed to strong positive selection pressure. During speciation, shifts in host range and preference can lead to major adaptive changes. Furthermore, evolution of total host resistance to most isolates can force rapid evolutionary changes in host-specific pathogens. Crop pathogens are subjected to particularly intense selective pressures from monocultures and fungicides. Detection of the footprints of positive selection in plant pathogen genomes is a worthwhile endeavor as it aids understanding of the fundamental biology of these important organisms. There are two main classes of test for detection of positively selected alleles. Tests based on the ratio of non-synonymous to synonymous substitutions per site detect the footprints of multiple fixation events between divergent lineages. Thus, they are well-suited to the study of ancient adaptation events spanning speciations. On the other hand, tests that scan genomes for local fluctuations in allelic diversity within populations are suitable for detection of recent positive selection in populations. In this review, I briefly describe some of the more widely used tests of positive selection and the theory underlying them. I then discuss various examples of their application to plant pathogen genomes, emphasizing the types of genes that are associated with signatures of positive selection. I conclude with a discussion of the practicality of such tests for identification of pathogen genes of interest and the important features of pathogen ecology that must be taken into account for accurate interpretation.

Hybridization fluctuates with rainfall in Darwin’s tree finches

Hybridization in natural populations may be an adaptive response to shifting climatic regimes, but understanding this can be limited by the timing of sampling effort and confident identification of hybrids. On the Galapagos Islands, Darwin’s finches regularly hybridize; the islands also show extreme annual variation in rainfall, but the effect of annual rainfall on the frequency of finch hybridization is little known. Across a 20-year period on Floreana Island, we compare patterns of hybridization in sympatric Darwin’s tree finches (N = 425; Camaryhnchus spp.) and test for an effect of annual rainfall on (1) the frequency of hybrids (C. pauper × C. parvulus) and (2) the percentage of male hybrid birds produced per year (hybrid recruitment). Annual rainfall correlated with recruitment positively for hybrids, negatively for C. parvulus and not at all for C. pauper. Furthermore, the percentage of hybrids (range: 12–56%) and C. parvulus did not change with sampling year, but the critically endangered C. pauper declined. Our findings indicate that hybrid recruitment is recurring and variable according to annual rainfall in Camarhynchus Darwin’s finches.

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.

Wolbachia Acquisition by Drosophila yakuba-Clade Hosts and Transfer of Incompatibility Loci Between Distantly Related Wolbachia

Maternally transmitted Wolbachia infect about half of insect species, yet the predominant mode(s) of Wolbachia acquisition remains uncertain. Species-specific associations could be old, with Wolbachia and hosts codiversifying (i.e., cladogenic acquisition), or relatively young and acquired by horizontal transfer or introgression. The three Drosophila yakuba-clade hosts [(D. santomea, D. yakuba) D. teissieri] diverged ∼3 MYA and currently hybridize on the West African islands Bioko and São Tomé. Each species is polymorphic for nearly identical Wolbachia that cause weak cytoplasmic incompatibility (CI)-reduced egg hatch when uninfected females mate with infected males. D. yakuba-clade Wolbachia are closely related to wMel, globally polymorphic in D. melanogaster We use draft Wolbachia and mitochondrial genomes to demonstrate that D. yakuba-clade phylogenies for Wolbachia and mitochondria tend to follow host nuclear phylogenies. However, roughly half of D. santomea individuals, sampled both inside and outside of the São Tomé hybrid zone, have introgressed D. yakuba mitochondria. Both mitochondria and Wolbachia possess far more recent common ancestors than the bulk of the host nuclear genomes, precluding cladogenic Wolbachia acquisition. General concordance of Wolbachia and mitochondrial phylogenies suggests that horizontal transmission is rare, but varying relative rates of molecular divergence complicate chronogram-based statistical tests. Loci that cause CI in wMel are disrupted in D. yakuba-clade Wolbachia; but a second set of loci predicted to cause CI are located in the same WO prophage region. These alternative CI loci seem to have been acquired horizontally from distantly related Wolbachia, with transfer mediated by flanking Wolbachia-specific ISWpi1 transposons.

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.

The Role of Transposable Elements in Speciation

Understanding the phenotypic and molecular mechanisms that contribute to genetic diversity between and within species is fundamental in studying the evolution of species. In particular, identifying the interspecific differences that lead to the reduction or even cessation of gene flow between nascent species is one of the main goals of speciation genetic research. Transposable elements (TEs) are DNA sequences with the ability to move within genomes. TEs are ubiquitous throughout eukaryotic genomes and have been shown to alter regulatory networks, gene expression, and to rearrange genomes as a result of their transposition. However, no systematic effort has evaluated the role of TEs in speciation. We compiled the evidence for TEs as potential causes of reproductive isolation across a diversity of taxa. We find that TEs are often associated with hybrid defects that might preclude the fusion between species, but that the involvement of TEs in other barriers to gene flow different from postzygotic isolation is still relatively unknown. Finally, we list a series of guides and research avenues to disentangle the effects of TEs on the origin of new species.

Molecular and ecological signatures of an expanding hybrid zone

Many species are currently changing their distributions and subsequently form sympatric zones with hybridization between formerly allopatric species as one possible consequence. The damselfly Ischnura elegans has recently expanded south into the range of its ecologically and morphologically similar sister species Ischnura graellsii. Molecular work shows ongoing introgression between these species, but the extent to which this species mixing is modulated by ecological niche use is not known. Here, we (1) conduct a detailed population genetic analysis based on molecular markers and (2) model the ecological niche use of both species in allopatric and sympatric regions. Population genetic analyses showed chronic introgression between I. elegans and I. graellsii across a wide part of Spain, and admixture analysis corroborated this, showing that the majority of I. elegans from the sympatric zone could not be assigned to either the I. elegans or I. graellsii species cluster. Niche modeling demonstrated that I. elegans has modified its environmental niche following hybridization and genetic introgression with I. graellsii, making niche space of introgressed I. elegans populations more similar to I. graellsii. Taken together, this corroborates the view that adaptive introgression has moved genes from I. graellsii into I. elegans and that this process is enabling Spanish I. elegans to occupy a novel niche, further facilitating its expansion. Our results add to the growing evidence that hybridization can play an important and creative role in the adaptive evolution of animals.

Supervised machine learning reveals introgressed loci in the genomes of Drosophila simulans and D. sechellia

Hybridization and gene flow between species appears to be common. Even though it is clear that hybridization is widespread across all surveyed taxonomic groups, the magnitude and consequences of introgression are still largely unknown. Thus it is crucial to develop the statistical machinery required to uncover which genomic regions have recently acquired haplotypes via introgression from a sister population. We developed a novel machine learning framework, called FILET (Finding Introgressed Loci via Extra-Trees) capable of revealing genomic introgression with far greater power than competing methods. FILET works by combining information from a number of population genetic summary statistics, including several new statistics that we introduce, that capture patterns of variation across two populations. We show that FILET is able to identify loci that have experienced gene flow between related species with high accuracy, and in most situations can correctly infer which population was the donor and which was the recipient. Here we describe a data set of outbred diploid Drosophila sechellia genomes, and combine them with data from D. simulans to examine recent introgression between these species using FILET. Although we find that these populations may have split more recently than previously appreciated, FILET confirms that there has indeed been appreciable recent introgression (some of which might have been adaptive) between these species, and reveals that this gene flow is primarily in the direction of D. simulans to D. sechellia.

Understanding the extent to which species or diverged populations hybridize in nature is crucially important if we are to understand the speciation process. Accordingly numerous research groups have developed methodology for finding the genetic evidence of such introgression. In this report we develop a supervised machine learning approach for uncovering loci which have introgressed across species boundaries. We show that our method, FILET, has greater accuracy and power than competing methods in discovering introgression, and in addition can detect the directionality associated with the gene flow between species. Using whole genome sequences from Drosophila simulans and Drosophila sechellia we show that FILET discovers quite extensive introgression between these species that has occurred mostly from D. simulans to D. sechellia. Our work highlights the complex process of speciation even within a well-studied system and points to the growing importance of supervised machine learning in population genetics.

Hybrids and horizontal transfer: introgression allows adaptive allele discovery

Evolution has devised countless remarkable solutions to diverse challenges. Understanding the mechanistic basis of these solutions provides insights into how biological systems can be subtly tweaked without maladaptive consequences. The knowledge gained from illuminating these mechanisms is equally important to our understanding of fundamental evolutionary mechanisms as it is to our hopes of developing truly rational plant breeding and synthetic biology. In particular, modern population genomic approaches are proving very powerful in the detection of candidate alleles for mediating consequential adaptations that can be tested functionally. Especially striking are signals gained from contexts involving genetic transfers between populations, closely related species, or indeed between kingdoms. Here we discuss two major classes of these scenarios, adaptive introgression and horizontal gene flow, illustrating discoveries made across kingdoms.

Molecular Population Genetics

Molecular population genetics aims to explain genetic variation and molecular evolution from population genetics principles. The field was born 50 years ago with the first measures of genetic variation in allozyme loci, continued with the nucleotide sequencing era, and is currently in the era of population genomics. During this period, molecular population genetics has been revolutionized by progress in data acquisition and theoretical developments. The conceptual elegance of the neutral theory of molecular evolution or the footprint carved by natural selection on the patterns of genetic variation are two examples of the vast number of inspiring findings of population genetics research. Since the inception of the field, Drosophila has been the prominent model species: molecular variation in populations was first described in Drosophila and most of the population genetics hypotheses were tested in Drosophila species. In this review, we describe the main concepts, methods, and landmarks of molecular population genetics, using the Drosophila model as a reference. We describe the different genetic data sets made available by advances in molecular technologies, and the theoretical developments fostered by these data. Finally, we review the results and new insights provided by the population genomics approach, and conclude by enumerating challenges and new lines of inquiry posed by increasingly large population scale sequence data.

The evolution of reproductive isolation in the Drosophila yakuba complex of species

In the Drosophila melanogaster subgroup, the yakuba species complex, D. yakuba, D. santomea and D. teissieri have identical mitochondrial genomes in spite of nuclear differentiation. The first two species can be readily hybridized in the laboratory and produce fertile females and sterile males. They also form hybrids in natural conditions. Nonetheless, the third species, D. teissieri, was thought to be unable to produce hybrids with either D. yakuba or D. santomea. This in turn posed the conundrum of why the three species shared a single mitochondrial genome. In this report, we show that D. teissieri can indeed hybridize with both D. yakuba and D. santomea. The resulting female hybrids from both crosses are fertile, whereas the hybrid males are sterile. We also characterize six isolating mechanisms that might be involved in keeping the three species apart. Our results open the possibility of studying the history of introgression in the yakuba species complex and dissecting the genetic basis of interspecific differences between these three species by genetic mapping.

Powerful methods for detecting introgressed regions from population genomic data

Understanding the types and functions of genes that are able to cross species boundaries-and those that are not-is an important step in understanding the forces maintaining species as largely independent lineages across the remainder of the genome. With large next-generation sequencing data sets we are now able to ask whether introgression has occurred across the genome, and multiple methods have been proposed to detect the signature of such events. Here, we introduce a new summary statistic that can be used to test for introgression, RNDmin , that makes use of the minimum pairwise sequence distance between two population samples relative to divergence to an outgroup. We find that our method offers a modest increase in power over other, related tests, but that all such tests have high power to detect introgressed loci when migration is recent and strong. RNDmin is robust to variation in the mutation rate, and remains reliable even when estimates of the divergence time between sister species are inaccurate. We apply RNDmin to population genomic data from the African mosquitoes Anopheles quadriannulatus and A. arabiensis, identifying three novel candidate regions for introgression. Interestingly, one of the introgressed loci is on the X chromosome, but outside of an inversion separating these two species. Our results suggest that significant, but rare, sharing of alleles is occurring between species that diverged more than 1 million years ago, and that application of these methods to additional systems are likely to reveal similar results.

How Do Species Interactions Affect Evolutionary Dynamics Across Whole Communities?

Theories of how species evolve in changing environments mostly consider single species in isolation or pairs of interacting species. Yet all organisms live in diverse communities containing many hundreds of species. This review discusses how species interactions influence the evolution of constituent species across whole communities. When species interactions are weak or inconsistent, evolutionary dynamics should be predictable by factors identified by single-species theory. Stronger species interactions, however, can alter evolutionary outcomes and either dampen or promote evolution of constituent species depending on the number of species and the distribution of interaction strengths across the interaction network. Genetic interactions, such as horizontal gene transfer, might also affect evolutionary outcomes. These evolutionary mechanisms in turn affect whole-community properties, such as the level of ecosystem functioning. Successful management of both ecosystems and focal species requires new understanding of evolutionary interactions across whole communities.

Genomic sequencing reveals historical, demographic and selective factors associated with the diversification of the fire-associated fungus Neurospora discreta

Delineating microbial populations, discovering ecologically relevant phenotypes and identifying migrants, hybrids or admixed individuals have long proved notoriously difficult, thereby limiting our understanding of the evolutionary forces at play during the diversification of microbial species. However, recent advances in sequencing and computational methods have enabled an unbiased approach whereby incipient species and the genetic correlates of speciation can be identified by examining patterns of genomic variation within and between lineages. We present here a population genomic study of a phylogenetic species in the Neurospora discreta species complex, based on the resequencing of full genomes (~37 Mb) for 52 fungal isolates from nine sites in three continents. Population structure analyses revealed two distinct lineages in South-East Asia, and three lineages in North America/Europe with a broad longitudinal and latitudinal range and limited admixture between lineages. Genome scans for selective sweeps and comparisons of the genomic landscapes of diversity and recombination provided no support for a role of selection at linked sites on genomic heterogeneity in levels of divergence between lineages. However, demographic inference indicated that the observed genomic heterogeneity in divergence was generated by varying rates of gene flow between lineages following a period of isolation. Many putative cases of exchange of genetic material between phylogenetically divergent fungal lineages have been discovered, and our work highlights the quantitative importance of genetic exchanges between more closely related taxa to the evolution of fungal genomes. Our study also supports the role of allopatric isolation as a driver of diversification in saprobic microbes.

Gene flow between Drosophila yakuba and Drosophila santomea in subunit V of cytochrome c oxidase: A potential case of cytonuclear cointrogression

Introgression is the effective exchange of genetic information between species through natural hybridization. Previous genetic analyses of the Drosophila yakuba—D. santomea hybrid zone showed that the mitochondrial genome of D. yakuba had introgressed into D. santomea and completely replaced its native form. Since mitochondrial proteins work intimately with nuclear‐encoded proteins in the oxidative phosphorylation (OXPHOS) pathway, we hypothesized that some nuclear genes in OXPHOS cointrogressed along with the mitochondrial genome. We analyzed nucleotide variation in the 12 nuclear genes that form cytochrome c oxidase (COX) in 33 Drosophila lines. COX is an OXPHOS enzyme composed of both nuclear‐ and mitochondrial‐encoded proteins and shows evidence of cytonuclear coadaptation in some species. Using maximum‐likelihood methods, we detected significant gene flow from D. yakuba to D. santomea for the entire COX complex. Interestingly, the signal of introgression is concentrated in the three nuclear genes composing subunit V, which shows population migration rates significantly greater than the background level of introgression in these species. The detection of introgression in three proteins that work together, interact directly with the mitochondrial‐encoded core, and are critical for early COX assembly suggests this could be a case of cytonuclear cointrogression.

High nucleotide diversity and limited linkage disequilibrium in Helicoverpa armigera facilitates the detection of a selective sweep

Insecticides impose extreme selective pressures on populations of target pests and so insecticide resistance loci of these species may provide the footprints of 'selective sweeps'. To lay the foundation for future genome-wide scans for selective sweeps and inform genome-wide association study designs, we set out to characterize some of the baseline population genomic parameters of one of the most damaging insect pests in agriculture worldwide, Helicoverpa armigera. To this end, we surveyed nine Z-linked loci in three Australian H. armigera populations. We find that estimates of π are in the higher range among other insects and linkage disequilibrium decays over short distances. One of the surveyed loci, a cytochrome P450, shows an unusual haplotype configuration with a divergent allele at high frequency that led us to investigate the possibility of an adaptive introgression around this locus.

Gene-specific selective sweeps in bacteria and archaea caused by negative frequency-dependent selection

BACKGROUND

Fixation of beneficial genes in bacteria and archaea (collectively, prokaryotes) is often believed to erase pre-existing genomic diversity through the hitchhiking effect, a phenomenon known as genome-wide selective sweep. Recent studies, however, indicate that beneficial genes spread through a prokaryotic population via recombination without causing genome-wide selective sweeps. These gene-specific selective sweeps seem to be at odds with the existing estimates of recombination rates in prokaryotes, which appear far too low to explain such phenomena.

RESULTS

We use mathematical modeling to investigate potential solutions to this apparent paradox. Most microbes in nature evolve in heterogeneous, dynamic communities, in which ecological interactions can substantially impact evolution. Here, we focus on the effect of negative frequency-dependent selection (NFDS) such as caused by viral predation (kill-the-winner dynamics). The NFDS maintains multiple genotypes within a population, so that a gene beneficial to every individual would have to spread via recombination, hence a gene-specific selective sweep. However, gene loci affected by NFDS often are located in variable regions of microbial genomes that contain genes involved in the mobility of selfish genetic elements, such as integrases or transposases. Thus, the NFDS-affected loci are likely to experience elevated rates of recombination compared with the other loci. Consequently, these loci might be effectively unlinked from the rest of the genome, so that NFDS would be unable to prevent genome-wide selective sweeps. To address this problem, we analyzed population genetic models of selective sweeps in prokaryotes under NFDS. The results indicate that NFDS can cause gene-specific selective sweeps despite the effect of locally elevated recombination rates, provided NFDS affects more than one locus and the basal rate of recombination is sufficiently low. Although these conditions might seem to contradict the intuition that gene-specific selective sweeps require high recombination rates, they actually decrease the effective rate of recombination at loci affected by NFDS relative to the per-locus basal level, so that NFDS can cause gene-specific selective sweeps.

CONCLUSION

Because many free-living prokaryotes are likely to evolve under NFDS caused by ubiquitous viruses, gene-specific selective sweeps driven by NFDS are expected to be a major, general phenomenon in prokaryotic populations.

A new method to scan genomes for introgression in a secondary contact model

Secondary contact between divergent populations or incipient species may result in the exchange and introgression of genomic material. We develop a simple DNA sequence measure, called G_min, which is designed to identify genomic regions experiencing introgression in a secondary contact model. G_min is defined as the ratio of the minimum between-population number of nucleotide differences in a genomic window to the average number of between-population differences. Although it is conceptually simple, one advantage of G_min is that it is computationally inexpensive relative to model-based methods for detecting gene flow and it scales easily to the level of whole-genome analysis. We compare the sensitivity and specificity of G_min to those of the widely used index of population differentiation, F_ST, and suggest a simple statistical test for identifying genomic outliers. Extensive computer simulations demonstrate that G_min has both greater sensitivity and specificity for detecting recent introgression than does F_ST. Furthermore, we find that the sensitivity of G_min is robust with respect to both the population mutation and recombination rates. Finally, a scan of G_min across the X chromosome of Drosophila melanogaster identifies candidate regions of introgression between sub-Saharan African and cosmopolitan populations that were previously missed by other methods. These results show that G_min is a biologically straightforward, yet powerful, alternative to F_ST, as well as to more computationally intensive model-based methods for detecting gene flow.

A new genome-wide method to track horizontally transferred sequences: application to Drosophila

Because of methodological breakthroughs and the availability of an increasing amount of whole-genome sequence data, horizontal transfers (HTs) in eukaryotes have received much attention recently. Contrary to similar analyses in prokaryotes, most studies in eukaryotes usually investigate particular sequences corresponding to transposable elements (TEs), neglecting the other components of the genome. We present a new methodological framework for the genome-wide detection of all putative horizontally transferred sequences between two species that requires no prior knowledge of the transferred sequences. This method provides a broader picture of HTs in eukaryotes by fully exploiting complete-genome sequence data. In contrast to previous genome-wide approaches, we used a well-defined statistical framework to control for the number of false positives in the results, and we propose two new validation procedures to control for confounding factors. The first validation procedure relies on a comparative analysis with other species of the phylogeny to validate HTs for the nonrepeated sequences detected, whereas the second one built upon the study of the dynamics of the detected TEs. We applied our method to two closely related Drosophila species, Drosophila melanogaster and D. simulans, in which we discovered 10 new HTs in addition to all the HTs previously detected in different studies, which underscores our method’s high sensitivity and specificity. Our results favor the hypothesis of multiple independent HTs of TEs while unraveling a small portion of the network of HTs in the Drosophila phylogeny.

Unexpected consequences of genetic selection in broilers and turkeys: problems and solutions

1. Genetic theory leads to the expectation that unexpected consequences of genetic selection for production traits will inevitably occur and that these changes are likely to be undesirable. 2. Both artificial selection for production efficiency and "natural" selection for adaptation to the production environment result in selection sweeps that increase the frequencies of rare recessive alleles that have a negative effect on fitness. 3. Fitness is broadly defined as any trait that affects the ability to survive, reproduce and contribute to the next generation, such as musculoskeletal disease in growing broiler chickens and multiple ovulation in adult broiler parents. 4. Welfare concerns about the negative effects of genetic selection on bird welfare are sometimes exaggerated but are nevertheless real. Breeders have paid increasing attention to these traits over several decades and have demonstrated improvement in pedigree flocks. There is an urgent need to monitor changes in commercial flocks to ensure that genetic change is accompanied by improvements in that target population. 5. New technologies for trait measurement, whole genome selection and targeted genetic modification hold out the promise of efficient and rapid improvement of welfare traits in future breeding of broiler chickens and turkeys. The potential of targeted genetic modification for enhancing welfare traits is considerable, but the goal of achieving public acceptability for the progeny of transgenic poultry will be politically challenging.

Testing evolutionary hypotheses for DNA barcoding failure in willows

The goal of DNA barcoding is to enable the rapid identification of taxa from short diagnostic DNA sequence profiles. But how feasible is this objective when many evolutionary processes, such as hybridization and selective sweeps, cause alleles to be shared among related taxa? In this issue of Molecular Ecology, Percy et al. (2014) test the full suite of seven candidate plant barcoding loci in a broad geographic sample of willow species. They show exceptional plastid haplotype sharing between species across continents, with most taxa not possessing a unique barcode sequence. Using population genetic and molecular dating analyses, they implicate hybridization and selective sweeps, but not incomplete lineage sorting, as the historical processes causing widespread haplotype sharing among willow taxa. This study represents an exceptional case of how poorly barcoding can perform, and highlights methodological issues using universal organellar regions for species identification.

Sequential adaptive introgression of the mitochondrial genome in Drosophila yakuba and Drosophila santomea

Interspecific hybridization provides the unique opportunity for species to tap into genetic variation present in a closely related species and potentially take advantage of beneficial alleles. It has become increasingly clear that when hybridization occurs, mitochondrial DNA (mtDNA) often crosses species boundaries, raising the possibility that it could serve as a recurrent target of natural selection and source of species' adaptations. Here we report the sequences of 46 complete mitochondrial genomes of Drosophila yakuba and Drosophila santomea, two sister species known to produce hybrids in nature (∼3%). At least two independent events of mtDNA introgression are uncovered in this study, including an early invasion of the D. yakuba mitochondrial genome that fully replaced the D. santomea mtDNA native haplotypes and a more recent, ongoing event centred in the hybrid zone. Interestingly, this recent introgression event bears the signature of Darwinian natural selection, and the selective haplotype can be found at low frequency in Africa mainland populations of D. yakuba. We put forward the possibility that, because the effective population size of D. santomea is smaller than that of D. yakuba, the faster accumulation of mildly deleterious mutations associated with Muller's ratchet in the former species may have facilitated the replacement of the mutationally loaded mitochondrial genome of D.

On the Coyne and Orr-igin of species: effects of intrinsic postzygotic isolation, ecological differentiation, x chromosome size, and sympatry on Drosophila speciation

Coyne and Orr found that mating discrimination (premating isolation) evolves much faster between sympatric than allopatric Drosophila species pairs. Their meta-analyses established that this pattern, expected under reinforcement, is common and that Haldane's rule is ubiquitous in Drosophila species divergence. We examine three possible contributors to the reinforcement pattern: intrinsic postzygotic isolation, dichotomized as to whether hybrid males show complete inviability/sterility; host-plant divergence, as a surrogate for extrinsic postzygotic isolation; and X chromosome size, whether roughly 20% or 40% of the genome is X-linked. We focus on "young" species pairs with overlapping ranges, contrasted with allopatric pairs. Using alternative criteria for "sympatry" and tests that compare either level of prezygotic isolation in sympatry or frequency of sympatry, we find no statistically significant effects associated with X chromosome size or our coarse quantifications of intrinsic postzygotic isolation or ecological differentiation. Although sympatric speciation seems very rare in animals, the pervasiveness of the reinforcement pattern and the commonness of range overlap for close relatives indicate that speciation in Drosophila is often not purely allopatric. It remains to determine whether increased premating isolation with sympatry results from secondary contact versus parapatric speciation and what drives this pattern.