Genome-wide sequence data show no evidence of hybridization and introgression among pollinator wasps associated with a community of Panamanian strangler figs.

The specificity of pollinator host choice influences opportunities for reproductive isolation in their host plants. Similarly, host plants can influence opportunities for reproductive isolation in their pollinators. For example, in the fig and fig wasp mutualism, offspring of fig pollinator wasps mate inside the inflorescence that the mothers pollinate. Although often host specific, multiple fig pollinator species are sometimes associated with the same fig species, potentially enabling hybridization between wasp species. Here, we study the 19 pollinator species (Pegoscapus spp.) associated with an entire community of 16 Panamanian strangler fig species (Ficus subgenus Urostigma, section Americanae) to determine whether the previously documented history of pollinator host switching and current host sharing predicts genetic admixture among the pollinator species, as has been observed in their host figs. Specifically, we use genome-wide ultraconserved element (UCE) loci to estimate phylogenetic relationships and test for hybridization and introgression among the pollinator species. In all cases, we recover well-delimited pollinator species that contain high interspecific divergence. Even among pairs of pollinator species that currently reproduce within syconia of shared host fig species, we found no evidence of hybridization or introgression. This is in contrast to their host figs, where hybridization and introgression have been detected within this community, and more generally, within figs worldwide. Consistent with general patterns recovered among other obligate pollination mutualisms (e.g. yucca moths and yuccas), our results suggest that while hybridization and introgression are processes operating within the host plants, these processes are relatively unimportant within their associated insect pollinators.

Evolution of GSTD1 in Cactophilic Drosophila.

GSTD1 is an insect glutathione S-transferase that has received considerable attention because of its role in detoxification of xenobiotic compounds, specifically pesticides and plant allelochemicals involved in detoxification, or in the use of some substrates as a nutritional source. GSTD1 has been implicated in the adaptation to a new cactus host in Drosophila mojavensis and thus constitutes an interesting candidate to study ecological genetics of adaptation in Drosophila. We conducted population genetic and molecular evolution analyses of the GstD1 gene in the context of association with different cactus hosts (Opuntia sp. vs. Columnar) in nine Drosophila species from the repleta group. We observed strong evidence of selection in GstD1 from D. hexastigma. This species is associated with a diverse set of columnar cacti with very complex chemistries. GstD1 sequences from D. hexastigma show evidence of a recent selective sweep, and positive selection at one residue just outside of the active site of the enzyme. The substitution (Q116T) at the site under selection leads to a conformational change in the enzyme that could have important consequences for substrate binding efficiency. Our results suggest that GSTD1 from D. hexastigma may have evolved improved substrate binding in order to adapt to the diverse chemical environments that this species encounters in the wild.

Lack of genetic isolation by distance, similar genetic structuring but different demographic histories in a fig-pollinating wasp mutualism.

Historical abiotic factors such as climatic oscillations and extreme climatic events as well as biotic factors have shaped the structuring of species' genetic diversity. In obligate species-specific mutualisms, the biogeographic histories of the interacting species are tightly linked. This could be particularly true for nuclear genes in the Ficus-pollinating wasp mutualistic association as the insects disperse pollen from their natal tree. In this study, we compare spatial genetic structure of plant and pollinator for the Ficus hirta-Valisia javana association throughout southeast China including Hainan Island, for both nuclear and cytoplasmic markers. We show that dispersal of the insect leads to plant and insect presenting similar signatures of lack of genetic isolation by distance for nuclear genes on the continent over a distance of 1000 km. But we also show that the demographic histories of plant and insect are strikingly different. This is in agreement with extreme climatic events leading to transient regional extinctions of the insects, associated with local survival of the plants. We also observe evidence of genetic differentiation for both wasps and fig-tree between the continent and Hainan Island, although the Qiongzhou Strait is only on average 30 km wide, suggesting that geographic isolation by itself has not been sufficient to generate this differentiation. Hence, our results suggest that in highly dispersive mutualistic systems, isolation-by-dispersal limitation across a geographic barrier could be supplemented by isolation by adaptation, and maybe by coevolution, allowing further genetic divergence. In such systems, species may frequently be composed of a single population.

The Complex Landscape of Structural Divergence Between the Drosophila pseudoobscura and D. persimilis Genomes.

Structural genomic variants are key drivers of phenotypic evolution. They can span hundreds to millions of base pairs and can thus affect large numbers of genetic elements. Although structural variation is quite common within and between species, its characterization depends upon the quality of genome assemblies and the proportion of repetitive elements. Using new high-quality genome assemblies, we report a complex and previously hidden landscape of structural divergence between the genomes of Drosophila persimilis and D. pseudoobscura, two classic species in speciation research, and study the relationships among structural variants, transposable elements, and gene expression divergence. The new assemblies confirm the already known fixed inversion differences between these species. Consistent with previous studies showing higher levels of nucleotide divergence between fixed inversions relative to collinear regions of the genome, we also find a significant overrepresentation of INDELs inside the inversions. We find that transposable elements accumulate in regions with low levels of recombination, and spatial correlation analyses reveal a strong association between transposable elements and structural variants. We also report a strong association between differentially expressed (DE) genes and structural variants and an overrepresentation of DE genes inside the fixed chromosomal inversions that separate this species pair. Interestingly, species-specific structural variants are overrepresented in DE genes involved in neural development, spermatogenesis, and oocyte-to-embryo transition. Overall, our results highlight the association of transposable elements with structural variants and their importance in driving evolutionary divergence.

Field sampling of fig pollinator wasps across host species and host developmental phase: Implications for host recognition and specificity.

Previous genetic studies of pollinator wasps associated with a community of strangler figs (Ficus subgenus Urostigma, section Americana) in Central Panama suggest that the wasp species exhibit a range in host specificity across their host figs. To better understand factors that might contribute to this observed range of specificity, we used sticky traps to capture fig-pollinating wasp individuals at 13 Ficus species, sampling at different phases of the reproductive cycle of the host figs (e.g., trees with receptive inflorescences, or vegetative trees, bearing only leaves). We also sampled at other tree species, using them as non-Ficus controls. DNA barcoding allowed us to identify the wasps to species and therefore assign their presence and abundance to host fig species and the developmental phase of that individual tree. We found: (1) wasps were only very rarely captured at non-Ficus trees; (2) nonetheless, pollinators were captured often at vegetative individuals of some host species; (3) overwhelmingly, wasp individuals were captured at receptive host fig trees representing the fig species from which they usually emerge. Our results indicate that wasp occurrence is not random either spatially or temporally within the forest and across these hosts, and that wasp specificity is generally high, both at receptive and vegetative host trees. Therefore, in addition to studies that show chemicals produced by receptive fig inflorescences attract pollinator wasps, we suggest that other cues (e.g., chemicals produced by the leaves) can also play a role in host recognition. We discuss our results in the context of recent findings on the role of host shifts in diversification processes in the Ficus genus.

Pollinator and host sharing lead to hybridization and introgression in Panamanian free-standing figs, but not in their pollinator wasps.

Obligate pollination mutualisms, in which plant and pollinator lineages depend on each other for reproduction, often exhibit high levels of species specificity. However, cases in which two or more pollinator species share a single host species (host sharing), or two or more host species share a single pollinator species (pollinator sharing), are known to occur in current ecological time. Further, evidence for host switching in evolutionary time is increasingly being recognized in these systems. The degree to which departures from strict specificity differentially affect the potential for hybridization and introgression in the associated host or pollinator is unclear. We addressed this question using genome-wide sequence data from five sympatric Panamanian free-standing fig species (Ficus subgenus Pharmacosycea, section Pharmacosycea) and their six associated fig-pollinator wasp species (Tetrapus). Two of the five fig species, F. glabrata and F. maxima, were found to regularly share pollinators. In these species, ongoing hybridization was demonstrated by the detection of several first-generation (F1) hybrid individuals, and historical introgression was indicated by phylogenetic network analysis. By contrast, although two of the pollinator species regularly share hosts, all six species were genetically distinct and deeply divergent, with no evidence for either hybridization or introgression. This pattern is consistent with results from other obligate pollination mutualisms, suggesting that, in contrast to their host plants, pollinators appear to be reproductively isolated, even when different species of pollinators mate in shared hosts.

Transcriptome of the adult female malaria mosquito vector Anopheles albimanus.

BACKGROUND: Human Malaria is transmitted by mosquitoes of the genus Anopheles. Transmission is a complex phenomenon involving biological and environmental factors of humans, parasites and mosquitoes. Among more than 500 anopheline species, only a few species from different branches of the mosquito evolutionary tree transmit malaria, suggesting that their vectorial capacity has evolved independently. Anopheles albimanus (subgenus Nyssorhynchus) is an important malaria vector in the Americas. The divergence time between Anopheles gambiae, the main malaria vector in Africa, and the Neotropical vectors has been estimated to be 100 My. To better understand the biological basis of malaria transmission and to develop novel and effective means of vector control, there is a need to explore the mosquito biology beyond the An. gambiae complex. RESULTS: We sequenced the transcriptome of the An. albimanus adult female. By combining Sanger, 454 and Illumina sequences from cDNA libraries derived from the midgut, cuticular fat body, dorsal vessel, salivary gland and whole body, we generated a single, high-quality assembly containing 16,669 transcripts, 92% of which mapped to the An. darlingi genome and covered 90% of the core eukaryotic genome. Bidirectional comparisons between the An. gambiae, An. darlingi and An. albimanus predicted proteomes allowed the identification of 3,772 putative orthologs. More than half of the transcripts had a match to proteins in other insect vectors and had an InterPro annotation. We identified several protein families that may be relevant to the study of Plasmodium-mosquito interaction. An open source transcript annotation browser called GDAV (Genome-Delinked Annotation Viewer) was developed to facilitate public access to the data generated by this and future transcriptome projects. CONCLUSIONS: We have explored the adult female transcriptome of one important New World malaria vector, An. albimanus. We identified protein-coding transcripts involved in biological processes that may be relevant to the Plasmodium lifecycle and can serve as the starting point for searching targets for novel control strategies. Our data increase the available genomic information regarding An. albimanus several hundred-fold, and will facilitate molecular research in medical entomology, evolutionary biology, genomics and proteomics of anopheline mosquito vectors. The data reported in this manuscript is accessible to the community via the VectorBase website (http://www.vectorbase.org/Other/AdditionalOrganisms/).

Enrichment of mRNA-like noncoding RNAs in the divergence of Drosophila males.

With the advent of transcriptome data, it has become clear that mRNA-like noncoding RNAs (mlncRNAs) are widespread in eukaryotes. Although their functions are poorly understood, these transcripts may play an important role in development and could thus be involved in determining developmental complexity and phenotypic diversification. However, few studies have assessed their potential roles in the divergence of closely related species. Here, we identify and study patterns of sequence and expression divergence in ten novel candidate mlncRNAs from Drosophila pseudoobscura and its close relative D. persimilis. The candidate mlncRNAs were identified by randomly sequencing a group of 734 cDNA clones from a microarray that showed either no difference in expression (187 clones) or differential expression (547 clones) in comparisons between D. pseudoobscura and D. persimilis and between these two species and their F(1) hybrids. Candidate mlncRNAs are overrepresented among differentially expressed transcripts between males of D. pseudoobscura and D. persimilis, and although they have high sequence conservation between these two species, seven of them have no putative homologs in any of the other ten Drosophila species whose genomes have been sequenced. Expression of eight of the ten candidate mlncRNAs was detected either in whole bodies (adults) or testes using a custom-designed oligonucleotide microarray. Three of the ten candidate mlncRNAs are highly expressed (in the top 4% of the male transcriptome), differentially expressed between species, and show extreme levels of sex-bias, with one transcript having the highest level of male bias in the whole transcriptome. Proteomic data from testes show no traces of any predicted peptides from the candidate mlncRNAs. Our results suggest that these mlncRNAs may be important in male-specific processes related to sexual dimorphism and species divergence in this species group.

Anonymous and EST-based microsatellite DNA markers that transfer broadly across the fig tree genus (Ficus, Moraceae).

PREMISE OF THE STUDY: We developed a set of microsatellite markers for broad utility across the species-rich pantropical tree genus Ficus (fig trees). The markers were developed to study population structure, hybridization, and gene flow in neotropical species. METHODS AND RESULTS: We developed seven novel primer sets from expressed sequence tag (EST) libraries of F. citrifolia and F. popenoei (subgen. Urostigma sect. Americana) and optimized five previously developed anonymous loci for cross-species amplification. The markers were successfully tested on four species from the basal subgenus Pharmacosycea sect. Pharmacosycea (F. insipida, F. maxima, F. tonduzii, and F. yoponensis) and seven species of the derived subgenus Urostigma (F. citrifolia, F. colubrinae, F. costaricana, F. nymphaeifolia, F. obtusifolia, F. pertusa, and F. popenoei). The 12 markers amplified consistently and displayed polymorphism in all the species. CONCLUSIONS: This set of microsatellite markers is transferable across the phylogenetic breadth of Ficus, and should therefore be useful for studies of population structure and gene flow in approximately 750 fig species worldwide.

Culture-free survey reveals diverse and distinctive fungal communities associated with developing figs (Ficus spp.) in Panama.

The ancient association of figs (Ficus spp.) and their pollinating wasps (fig wasps; Chalcidoidea, Hymenoptera) is one of the most interdependent plant-insect mutualisms known. In addition to pollinating wasps, a diverse community of organisms develops within the microcosm of the fig inflorescence and fruit. To better understand the multipartite context of the fig-fig wasp association, we used a culture-free approach to examine fungal communities associated with syconia of six species of Ficus and their pollinating wasps in lowland Panama. Diverse fungi were recovered from surface-sterilized flowers of all Ficus species, including gall- and seed flowers at four developmental stages. Fungal communities in syconia and on pollinating wasps were similar, dominated by diverse and previously unknown Saccharomycotina, and distinct from leaf- and stem endophyte communities in the same region. Before pollination, fungal communities were similar between gall- and seed flowers and among Ficus species. However, fungal communities differed significantly in flowers after pollination vs. before pollination, and between anciently diverged lineages of Ficus with active vs. passive pollination syndromes. Within groups of relatively closely related figs, there was little evidence for strict-sense host specificity between figs and particular fungal species. Instead, mixing of fungal communities among related figs, coupled with evidence for possible transfer by pollinating wasps, is consistent with recent suggestions of pollinator mixing within syconia. In turn, changes in fungal communities during fig development and ripening suggest an unexplored role of yeasts in the context of the fig-pollinator wasp mutualism.

Phylogenetic diversity of two common Trypanosoma cruzi lineages in the Southwestern United States.

Trypanosoma cruzi is the causative agent of Chagas disease, a devastating parasitic disease endemic to Central and South America, Mexico, and the USA. We characterized the genetic diversity of Trypanosoma cruzi circulating in five triatomine species (Triatoma gerstaeckeri, T. lecticularia, T.indictiva, T. sanguisuga and T. recurva) collected in Texas and Southern Arizona using multilocus sequence typing (MLST) with four single-copy loci (cytochrome oxidase subunit II- NADH dehydrogensase subunit 1 region (COII-ND1), mismatch-repair class 2 (MSH2), dihydrofolate reductase-thymidylate synthase (DHFR-TS) and a nuclear gene with ID TcCLB.506529.310). All T. cruzi variants fall in two main genetic lineages: 75% of the samples corresponded to T. cruzi Discrete Typing Unit (DTU) I (TcI), and 25% to a North American specific lineage previously labelled TcIV-USA. Phylogenetic and sequence divergence analyses of our new data plus all previously published sequence data from those four loci collected in the USA, show that TcIV-USA is significantly different from any other previously defined T. cruzi DTUs. The significant level of genetic divergence between TcIV-USA and other T. cruzi DTUs should lead to an increased focus on understanding the epidemiological importance of this DTU, as well as its geographical range and pathogenicity in humans and domestic animals. Our findings further corroborate the fact that there is a high genetic diversity of the parasite in North America and emphasize the need for appropriate surveillance and vector control programs for Chagas disease in southern USA and Mexico.

Inversions shape the divergence of Drosophila pseudoobscura and Drosophila persimilis on multiple timescales.

By shaping meiotic recombination, chromosomal inversions can influence genetic exchange between hybridizing species. Despite the recognized importance of inversions in evolutionary processes such as divergence and speciation, teasing apart the effects of inversions over time remains challenging. For example, are their effects on sequence divergence primarily generated through creating blocks of linkage disequilibrium prespeciation or through preventing gene flux after speciation? We provide a comprehensive look into the influence of inversions on gene flow throughout the evolutionary history of a classic system: Drosophila pseudoobscura and Drosophila persimilis. We use extensive whole-genome sequence data to report patterns of introgression and divergence with respect to chromosomal arrangements. Overall, we find evidence that inversions have contributed to divergence patterns between D. pseudoobscura and D. persimilis over three distinct timescales: (1) segregation of ancestral polymorphism early in the speciation process, (2) gene flow after the split of D. pseudoobscura and D. persimilis, but prior to the split of D. pseudoobscura subspecies, and (3) recent gene flow between sympatric D. pseudoobscura and D. persimilis, after the split of D. pseudoobscura subspecies. We discuss these results in terms of our understanding of evolution in this classic system and provide cautions for interpreting divergence measures in other systems.

Genomic evidence of prevalent hybridization throughout the evolutionary history of the fig-wasp pollination mutualism.

Ficus (figs) and their agaonid wasp pollinators present an ecologically important mutualism that also provides a rich comparative system for studying functional co-diversification throughout its coevolutionary history (~75 million years). We obtained entire nuclear, mitochondrial, and chloroplast genomes for 15 species representing all major clades of Ficus. Multiple analyses of these genomic data suggest that hybridization events have occurred throughout Ficus evolutionary history. Furthermore, cophylogenetic reconciliation analyses detect significant incongruence among all nuclear, chloroplast, and mitochondrial-based phylogenies, none of which correspond with any published phylogenies of the associated pollinator wasps. These findings are most consistent with frequent host-switching by the pollinators, leading to fig hybridization, even between distantly related clades. Here, we suggest that these pollinator host-switches and fig hybridization events are a dominant feature of fig/wasp coevolutionary history, and by generating novel genomic combinations in the figs have likely contributed to the remarkable diversity exhibited by this mutualism.

Molecular mechanisms of mutualistic and antagonistic interactions in a plant-pollinator association.

Many insects metamorphose from antagonistic larvae into mutualistic adult pollinators, with reciprocal adaptation leading to specialized insect-plant associations. It remains unknown how such interactions are established at molecular level. Here we assemble high-quality genomes of a fig species, Ficus pumila var. pumila, and its specific pollinating wasp, Wiebesia pumilae. We combine multi-omics with validation experiments to reveal molecular mechanisms underlying this specialized interaction. In the plant, we identify the specific compound attracting pollinators and validate the function of several key genes regulating its biosynthesis. In the pollinator, we find a highly reduced number of odorant-binding protein genes and an odorant-binding protein mainly binding the attractant. During antagonistic interaction, we find similar chemical profiles and turnovers throughout the development of galled ovules and seeds, and a significant contraction of detoxification-related gene families in the pollinator. Our study identifies some key genes bridging coevolved mutualists, establishing expectations for more diffuse insect-pollinator systems.

Polytene chromosomal maps of 11 Drosophila species: the order of genomic scaffolds inferred from genetic and physical maps.

The sequencing of the 12 genomes of members of the genus Drosophila was taken as an opportunity to reevaluate the genetic and physical maps for 11 of the species, in part to aid in the mapping of assembled scaffolds. Here, we present an overview of the importance of cytogenetic maps to Drosophila biology and to the concepts of chromosomal evolution. Physical and genetic markers were used to anchor the genome assembly scaffolds to the polytene chromosomal maps for each species. In addition, a computational approach was used to anchor smaller scaffolds on the basis of the analysis of syntenic blocks. We present the chromosomal map data from each of the 11 sequenced non-Drosophila melanogaster species as a series of sections. Each section reviews the history of the polytene chromosome maps for each species, presents the new polytene chromosome maps, and anchors the genomic scaffolds to the cytological maps using genetic and physical markers. The mapping data agree with Muller's idea that the majority of Drosophila genes are syntenic. Despite the conservation of genes within homologous chromosome arms across species, the karyotypes of these species have changed through the fusion of chromosomal arms followed by subsequent rearrangement events.

Host specificity, phenotype matching and the evolution of reproductive isolation in a coevolved plant-pollinator mutualism.

Coevolutionary interactions between plants and their associated pollinators and seed dispersers are thought to have promoted the diversification of flowering plants (Raven 1977; Regal 1977; Stebbins 1981). The actual mechanisms by which pollinators could drive species diversification in plants are not fully understood. However, it is thought that pollinator host specialization can influence the evolution of reproductive isolation among plant populations because the pollinator's choice of host is what determines patterns of gene flow in its host plant, and host choice may also have important consequences on pollinator and host fitness (Grant 1949; Bawa 1992). In this issue of Molecular Ecology, Smith et al. (2009) present a very interesting study that addresses how host specialization affects pollinator fitness and patterns of gene flow in a plant host. Several aspects of this study match elements of a seminal mathematical model of plant-pollinator codivergence (Kiester et al. 1984) suggesting that reciprocal selection for matched plant and pollinator reproductive traits may lead to speciation in the host and its pollinator when there is strong host specialization and a pattern of geographic subdivision. Smith et al.'s study represents an important step to fill the gap in our understanding of how reciprocal selection may lead to speciation in coevolved plant-pollinator mutualisms.

Molecular dating and biogeography of fig-pollinating wasps.

Figs and fig-pollinating wasps are obligate mutualists that have coevolved for over 60 million years. But when and where did pollinating fig wasps (Agaonidae) originate? Some studies suggest that agaonids arose in the Late Cretaceous and the current distribution of fig-wasp faunas can be explained by the break-up of the Gondwanan landmass. However, recent molecular-dating studies suggest divergence time estimates that are inconsistent with the Gondwanan vicariance hypothesis and imply that long distance oceanic dispersal could have been an important process for explaining the current distribution of both figs and fig wasps. Here, we use a combination of phylogenetic and biogeographical data to infer the age, the major period of diversification, and the geographic origin of pollinating fig wasps. Age estimates ranged widely depending on the molecular-dating method used and even when using the same method but with slightly different constraints, making it difficult to assess with certainty a Gondwanan origin of agaonids. The reconstruction of ancestral areas suggests that the most recent common ancestor of all extant fig-pollinating wasps was most likely Asian, although a southern Gondwana origin cannot be rejected. Our analysis also suggests that dispersal has played a more important role in the development of the fig-wasp biota than previously assumed.

Inferring processes of coevolutionary diversification in a community of Panamanian strangler figs and associated pollinating wasps.

The fig and pollinator wasp obligate mutualism is diverse (∼750 described species), ecologically important, and ancient (∼80 Ma). Once thought to be an example of strict one-to-one cospeciation, current thinking suggests genera of pollinator wasps codiversify with corresponding sections of figs, but the degree to which cospeciation or other processes contribute to the association at finer scales is unclear. Here, we use genome-wide sequence data from a community of Panamanian strangler figs and associated wasp pollinators to estimate the relative contributions of four evolutionary processes generating cophylogenetic patterns in this mutualism: cospeciation, host switching, pollinator speciation, and pollinator extinction. Using a model-based approach adapted from the study of gene family evolution, our results demonstrate the importance of host switching of pollinator wasps at this fine phylogenetic and regional scale. Although we estimate a modest amount of cospeciation, simulations reveal the number of putative cospeciation events to be consistent with what would be expected by chance. Additionally, model selection tests identify host switching as a critical parameter for explaining cophylogenetic patterns in this system. Our study demonstrates a promising approach through which the history of evolutionary association between interacting lineages can be rigorously modeled and tested in a probabilistic phylogenetic framework.

Evaluation of the genomic extent of effects of fixed inversion differences on intraspecific variation and interspecific gene flow in Drosophila pseudoobscura and D. persimilis.

There is increasing evidence that chromosomal inversions may facilitate the formation or persistence of new species by allowing genetic factors conferring species-specific adaptations or reproductive isolation to be inherited together and by reducing or eliminating introgression. However, the genomic domain of influence of the inverted regions on introgression has not been carefully studied. Here, we present a detailed study on the consequences that distance from inversion breakpoints has had on the inferred level of gene flow and divergence between Drosophila pseudoobscura and D. persimilis. We identified the locations of the inversion breakpoints distinguishing D. pseudoobscura and D. persimilis in chromosomes 2, XR, and XL. Population genetic data were collected at specific distances from the inversion breakpoints of the second chromosome and at two loci inside the XR and XL inverted regions. For loci outside the inverted regions, we found that distance from the nearest inversion breakpoint had a significant effect on several measures of divergence and gene flow between D. pseudoobscura and D. persimilis. The data fitted a logarithmic relationship, showing that the suppression of crossovers in inversion heterozygotes also extends to loci located outside the inversion but close to it (within 1-2 Mb). Further, we detected a significant reduction in nucleotide variation inside the inverted second chromosome region of D. persimilis and near one breakpoint, consistent with a scenario in which this inversion arose and was fixed in this species by natural selection.

Divergence between the Drosophila pseudoobscura and D. persimilis genome sequences in relation to chromosomal inversions.

As whole-genome sequence assemblies accumulate, a challenge is to determine how these can be used to address fundamental evolutionary questions, such as inferring the process of speciation. Here, we use the sequence assemblies of Drosophila pseudoobscura and D. persimilis to test hypotheses regarding divergence with gene flow. We observe low differentiation between the two genome sequences in pericentromeric and peritelomeric regions. We interpret this result as primarily a remnant of the correlation between levels of variation and local recombination rate observed within populations. However, we also observe lower differentiation far from the fixed chromosomal inversions distinguishing these species and greater differentiation within and near these inversions. This finding is consistent with models suggesting that chromosomal inversions facilitate species divergence despite interspecies gene flow. We also document heterogeneity among the inverted regions in their degree of differentiation, suggesting temporal differences in the origin of each inverted region consistent with the inversions arising during a process of divergence with gene flow. While this study provides insights into the speciation process using two single-genome sequences, it was informed by lower throughput but more rigorous examinations of polymorphism and divergence. This reliance highlights the need for complementary genomic and population genetic approaches for tackling fundamental evolutionary questions such as speciation.