Temporal and spatial mosaics: deep host association and shallow geographic drivers shape genetic structure in a widespread pinworm,Rauschtineria eutamii (Nematoda: Oxyuridae)

Climate and host demographic cycling often shape both parasite genetic diversity and host distributions, processes that transcend a history of strict host–parasite association. We explored host associations and histories based on an evaluation of mitochondrial and nuclear sequences to reveal the underlying history and genetic structure of a pinworm, Rauschtineria eutamii, infecting ten species of western North American chipmunks (Rodentia:Tamias, subgenus Neotamias). Rauschtineria eutamii contains divergent lineages influenced by the diversity of hosts and variation across the complex topography of western North America. We recovered six reciprocally monophyletic R. eutamii mitochondrial clades, largely supported by a multilocus concordance tree, exhibiting divergence levels comparable with intraspecific variation reported for other nematodes. Phylogenetic relationships among pinworm clades suggest that R. eutamii colonized an ancestral lineage of western chipmunks and lineages persisted during historical isolation in diverging Neotamias species or species groups. Pinworm diversification, however, is incongruent and asynchronous relative to host diversification. Secondarily, patterns of shallow divergence were shaped by geography through events of episodic colonization reflecting an interaction of taxon pulses and ecological fitting among assemblages in recurrent sympatry. Pinworms occasionally infect geographically proximal host species; however, host switching may be unstable or ephemeral, as there is no signal of host switching in the deeper history of R. eutamii.     

First report on Wolbachia endosymbiosis in freshwater Aphelocheirus aestivalis (Heteroptera: Aphelocheiridae) and its potential impact on genetic diversity of host

Alphaproteobacteria Wolbachia have been described as endosymbionts of approximately half of all aquatic insect species. These bacteria might affect not only reproduction but also the genetic diversity of its hosts. In the present study we identified Wolbachia endosymbiosis in freshwater true bug Aphelocheirus aestivalis F., 1794 (Heteroptera: Aphelocheiridae). Despite the fact that A. aestivalis is widely distributed in Europe, it occurs rather locally, often in isolated populations. Taking into account that Wolbachia, close relationships and past demographic phenomena could affect the genetic diversity of its host, we analyzed mitochondrial (COI and 16S) and nuclear (internal transcribed spacer 2) markers determined for A. aestivalis individuals collected from five populations. Moreover, we compared obtained COI sequences with those deposited in GenBank. Analyses revealed low genetic differentiation among samples tested, as well as low variation among determined COI sequences and those downloaded from the database. Although Wolbachia infection could correlate with decreasing mitochondrial diversity of its host, we suggest that low genetic variation observed in tested A. aestivalis samples (at both mitochondrial and nuclear levels) is a result of populations’ close relationships, past demographic phenomena or is characteristic for this species. Detailed analysis of the wsp gene fragment revealed two distinct strains of Wolbachia infecting A. aestivalis. Both of them belong to supergroup A, also found in other arthropods.     

Population genetics of Wolbachia‐infected,parthenogenetic and uninfected,sexual populations of Tetrastichus coeruleus (Hymenoptera: Eulophidae)

Wolbachia are endosymbiotic bacteria known to manipulate the reproduction of their hosts. These manipulations are expected to have consequences on the population genetics of the host, such as heterozygosity levels, genetic diversity and gene flow. The parasitoid wasp Tetrastichus coeruleus has populations that are infected with parthenogenesis‐inducing Wolbachia and populations that are not infected. We studied the population genetics of T. coeruleus between and within Wolbachia‐infected and uninfected populations, using nuclear microsatellites and mitochondrial DNA. We expected reduced genetic diversity in both DNA types in infected populations. However, migration and gene flow could introduce new DNA variants into populations. We therefore paid special attention to individuals with unexpected (genetic) characteristics. Based on nuclear and mitochondrial DNA, two genetic clusters were evident: a thelytokous cluster containing all Wolbachia‐infected, parthenogenetic populations and an arrhenotokous cluster containing all uninfected, sexual populations. Nuclear and mitochondrial DNA did not exhibit concordant patterns of variation, although there was reduced genetic diversity in infected populations for both DNA types. Within the thelytokous cluster, there was nuclear DNA variation, but no mitochondrial DNA variation. This nuclear DNA variation may be explained by occasional sex between infected females and males, by horizontal transmission of Wolbachia, and/or by novel mutations. Several females from thelytokous populations were uninfected and/or heterozygous for microsatellite loci. These unexpected characteristics may be explained by migration, by inefficient transmission of Wolbachia, by horizontal transmission of Wolbachia, and/or by novel mutations. However, migration has not prevented the build‐up of considerable genetic differentiation between thelytokous and arrhenotokous populations.     

Host species,and not environment,predicts variation in blood parasite prevalence,distribution, and diversity along a humidity gradient in northern South America

Environmental factors strongly influence the ecology and evolution of vector‐borne infectious diseases. However, our understanding of the influence of climatic variation on host–parasite interactions in tropical systems is rudimentary. We studied five species of birds and their haemosporidian parasites (Plasmodium and Haemoproteus) at 16 sampling sites to understand how environmental heterogeneity influences patterns of parasite prevalence, distribution, and diversity across a marked gradient in water availability in northern South America. We used molecular methods to screen for parasite infections and to identify parasite lineages. To characterize spatial heterogeneity in water availability, we used weather‐station and remotely sensed climate data. We estimated parasite prevalence while accounting for spatial autocorrelation, and used a model selection approach to determine the effect of variables related to water availability and host species on prevalence. The prevalence, distribution, and lineage diversity of haemosporidian parasites varied among localities and host species, but we found no support for the hypothesis that the prevalence and diversity of parasites increase with increasing water availability. Host species and host × climate interactions had stronger effects on infection prevalence, and parasite lineages were strongly associated with particular host species. Because climatic variables had little effect on the overall prevalence and lineage diversity of haemosporidian parasites across study sites, our results suggest that independent host–parasite dynamics may influence patterns in parasitism in environmentally heterogeneous landscapes.     

Out of the forest: past and present range expansion of a parthenogenetic weevil pest,or how to colonize the world successfully

Previous research revealed complex diversification patterns in the parthenogenetic weevil Naupactus cervinus. To understand the origin of clonal diversity and successful spreading of this weevil, we investigated its geographic origin and possible dispersal routes and whether parthenogens can persist in habitats under unsuitable environmental conditions. This study is based on samples taken throughout a broad area of the species’ range. We used both mitochondrial and nuclear markers and applied phylogenetic and network analyses to infer possible relationships between haplotypes. Bayesian phylogeographic analyses and ecological niche modeling were used to investigate the processes that shaped genetic diversity and enabled the colonization of new geographic areas. Southeastern Brazil emerges as the original distribution area of N. cervinus. We detected two range expansions, one along natural corridors during the Pleistocene and the other in countries outside South America during recent times. Isolation due to climate shifts during the early Pleistocene led to diversification in two divergent clades, which probably survived in different refugia of the Paranaense Forest and the Paraná River delta. The origin of the clonal diversity was probably a complex process including mutational diversification, hybridization, and secondary colonization. The establishment of N. cervinus in areas outside its native range may indicate adaptation to drier and cooler conditions. Parthenogenesis would be advantageous for the colonization of new environments by preventing the breakup of successful gene combinations. As in other insect pests, the present distribution of N. cervinus results from both its evolutionary history and its recent history related to human activities.     

First record of the bacterial endosymbiont Wolbachia for phytophagous hoverflies from genus Merodon (Diptera: Syrphidae)

Wolbachia is a widespread bacterial endosymbiont among arthropod species. It influences the reproduction of the host species and also mitochondrial DNA diversity. Until now there were only a few studies that detected Wolbachia infections in hoverflies (Diptera: Syrphidae), and this is the first broader study with the aim of examining the incidence of Wolbachia in the hoverfly genus Merodon. The obtained results indicate an infection rate of 96% and the presence of both Wolbachia supergroup A and B, which are characteristic for most of the infected arthropod species. Additionally, the presence of multiple Wolbachia strains in the Merodon aureus group species was detected and the mitochondrial DNA COI‐based relationships of the group are discussed in the light of infection. Finally, we discuss plant‐mediated horizontal transmission of Wolbachia strains among the studied hoverfly species.     

Monophyly of Wolbachia pipientis genomes within Drosophila melanogaster: geographic structuring,titre variation and host effects across five populations

Wolbachia pipientis is one of the most widely studied endosymbionts today, yet we know little about its short‐term adaptation and evolution. Here, using a set of 91 inbred Drosophila melanogaster lines from five populations, we explore patterns of diversity and recent evolution in the Wolbachia strain wMel. Within the D. melanogaster lines, we identify six major mitochondrial clades and four wMel clades. Concordant with past studies, the Wolbachia haplotypes contain an overall low level of nucleotide diversity, yet they still display geographic structuring. Using Bayesian analysis informed with demographic estimates of colonization times, we estimate that all extant D. melanogaster mitochondrial haplotypes coalesce to a Wolbachia‐infected ancestor approximately 2200 years ago. Finally, we measure wMel titre within the infected flies and find that titre varies across populations, an effect attributable to host genetic factors. This demonstration of local phenotypic divergence suggests that intraspecific host genetic variation plays a key role in shaping this model symbiotic system.     

Geographic and host‐mediated population genetic structure in a cestode parasite of the three‐spined stickleback

Comparative studies of genetic diversity and population structure can shed light on the ecological and evolutionary factors that influence host–parasite interactions. Here we examined whether geography, time and genetic variation in Alaskan three‐spined stickleback (Gasterosteus aculeatus Linneaus) hosts shape the population genetic structure of the diphyllobothridean cestode parasite Schistocephalus solidus (Müller, 1776). Host lineages and haplotypes were identified by sequencing the mitochondrial cytochrome b gene, and host population structure was assessed by Bayesian clustering analysis of allelic variation at 11 microsatellite loci. Parasite population structure was characterized according to allelic variation at eight microsatellite loci. Mantel tests and canonical redundancy analysis were conducted to evaluate the proportion of parasite genetic variation attributable to time and geography vs. host lineage, haplotype, and genotypic cluster. Host and parasite population structure were largely discordant across the study area, probably reflecting differences in gene flow, environmental influences external to the host, and genomic admixture among host lineages. We found that geography explained the greatest proportion of parasite genetic variation, but that variation also reflects time, host lineage, and host haplotype. Associations with host haplotypes suggest that one parasite genotypic cluster exhibits a narrower host range, predominantly infecting the most common host haplotypes, whereas the other parasite cluster infects all haplotypes equally, including rare haplotypes. Although experimental infection trials might prove otherwise, distributional differences in hosts preferentially infected by S. solidus could underlie the observed pattern of population structure.     

Gut microbial diversity across a contact zone for California voles: Implications for lineage divergence of hosts and mitonuclear mismatch in the assembly of the mammalian gut microbiome

Gut microbial diversity is thought to reflect the co‐evolution of microbes and their hosts as well as current host‐specific attributes such as genetic background and environmental setting. To explore interactions among these parameters, we characterized variation in gut microbiome composition of California voles (Microtus californicus) across a contact zone between two recently diverged lineages of this species. Because this contact zone contains individuals with mismatched mitochondrial‐nuclear genomes (cybrids), it provides an important opportunity to explore how different components of the genotype contribute to gut microbial diversity. Analyses of bacterial 16S rRNA sequences and joint species distribution modelling revealed that host genotypes and genetic differentiation among host populations together explained more than 50% of microbial community variation across our sampling transect. The ranked importance (most to least) of factors contributing to gut microbial diversity in our study populations were: genome‐wide population differentiation, local environmental conditions, and host genotypes. However, differences in microbial communities among vole populations (β‐diversity) did not follow patterns of lineage divergence (i.e., phylosymbiosis). Instead, among‐population variation was best explained by the spatial distribution of hosts, as expected if the environment is a primary source of gut microbial diversity (i.e., dispersal limitation hypothesis). Across the contact zone, several bacterial taxa differed in relative abundance between the two parental lineages as well as among individuals with mismatched mitochondrial and nuclear genomes. Thus, genetic divergence among host lineages and mitonuclear genomic mismatches may also contribute to microbial diversity by altering interactions between host genomes and gut microbiota (i.e., hologenome speciation hypothesis).     

Collembola,the biological species concept and the underestimation of global species richness

Despite its ancient origin, global distribution and abundance in nearly all habitats, the class Collembola is comprised of only 8000 described species and is estimated to number no more than 50 000. Many morphologically defined species have broad geographical ranges that span continents, and recent molecular work has revealed high genetic diversity within species. However, the evolutionary significance of this genetic diversity is unknown. In this study, we sample five morphological species of the globally distributed genus Lepidocyrtus from 14 Panamanian sampling sites to characterize genetic diversity and test morphospecies against the biological species concept. Mitochondrial and nuclear DNA sequence data were analysed and a total of 58 molecular lineages revealed. Deep lineage diversification was recovered, with 30 molecular lineages estimated to have established more than 10 million years ago, and the origin almost all contemporary lineages preceding the onset of the Pleistocene (~2 Mya). Thirty‐four lineages were sampled in sympatry revealing unambiguous cosegregation of mitochondrial and nuclear DNA sequence variation, consistent with biological species. Species richness within the class Collembola and the geographical structure of this diversity are substantially misrepresented components of terrestrial animal biodiversity. We speculate that global species richness of Collembola could be at least an order of magnitude greater than a previous estimate of 50 000 species.     

Gene flow and Andean uplift shape the diversification of Gasteracantha cancriformis (Araneae: Araneidae) in Northern South America

The Andean uplift has played a major role in shaping the current Neotropical biodiversity. However, in arthropods other than butterflies, little is known about how this geographic barrier has impacted species historical diversification. Here, we examined the phylogeography of the widespread color polymorphic spider Gasteracantha cancriformis to evaluate the effect of the northern Andean uplift on its divergence and assess whether its diversification occurred in the presence of gene flow. We inferred phylogenetic relationships and divergence times in G. cancriformis using mitochondrial and nuclear data from 105 individuals in northern South America. Genetic diversity, divergence, and population structure were quantified. We also compared multiple demographic scenarios for this species using a model‐based approach (Phrapl ) to determine divergence with or without gene flow. At last, we evaluated the association between genetic variation and color polymorphism. Both nuclear and mitochondrial data supported two well‐differentiated clades, which correspond to populations occurring on opposite sides of the Eastern cordillera of the Colombian Andes. The final uplift of this cordillera was identified as the most likely force that shaped the diversification of G. cancriformis in northern South America, resulting in a cis‐ and trans‐Andean phylogeographic structure for the species. We also found shared genetic variation between the cis‐ and trans‐Andean clades, which is better explained by a scenario of historical divergence in the face of gene flow. This has been likely facilitated by the presence of low‐elevation passes across the Eastern Colombian cordillera. Our work constitutes the first example in which the Andean uplift coupled with gene flow influenced the evolutionary history of an arachnid lineage.     

Tracing horizontal Wolbachia movements among bees (Anthophila): a combined approach using multilocus sequence typing data and host phylogeny

The endosymbiotic bacterium Wolbachia enhances its spread via vertical transmission by generating reproductive effects in its hosts, most notably cytoplasmic incompatibility (CI). Additionally, frequent interspecific horizontal transfer is evident from a lack of phylogenetic congruence between Wolbachia and its hosts. The mechanisms of this lateral transfer are largely unclear. To identify potential pathways of Wolbachia movements, we performed multilocus sequence typing of Wolbachia strains from bees (Anthophila). Using a host phylogeny and ecological data, we tested various models of horizontal endosymbiont transmission. In general, Wolbachia strains seem to be randomly distributed among bee hosts. Kleptoparasite‐host associations among bees as well as other ecological links could not be supported as sole basis for the spread of Wolbachia. However, cophylogenetic analyses and divergence time estimations suggest that Wolbachia may persist within a host lineage over considerable timescales and that strictly vertical transmission and subsequent random loss of infections across lineages may have had a greater impact on Wolbachia strain distribution than previously estimated. Although general conclusions about Wolbachia movements among arthropod hosts cannot be made, we present a framework by which precise assumptions about shared evolutionary histories of Wolbachia and a host taxon can be modelled and tested.     

Reproductive isolation due to prezygotic isolation and postzygotic cytoplasmic incompatibility in parasitoid wasps

The reproductive barriers that prevent gene flow between closely related species are a major topic in evolutionary research. Insect clades with parasitoid lifestyle are among the most species‐rich insects and new species are constantly described, indicating that speciation occurs frequently in this group. However, there are only very few studies on speciation in parasitoids. We studied reproductive barriers in two lineages of Lariophagus distinguendus (Chalcidoidea: Hymenoptera), a parasitoid wasp of pest beetle larvae that occur in human environments. One of the two lineages occurs in households preferably attacking larvae of the drugstore beetle Stegobium paniceum (“DB‐lineage”), the other in grain stores with larvae of the granary weevil Sitophilus granarius as main host (“GW‐lineage”). Between two populations of the DB‐lineage, we identified slight sexual isolation as intraspecific barrier. Between populations from both lineages, we found almost complete sexual isolation caused by female mate choice, and postzygotic isolation, which is partially caused by cytoplasmic incompatibility induced by so far undescribed endosymbionts which are not Wolbachia or Cardinium. Because separation between the two lineages is almost complete, they should be considered as separate species according to the biological species concept. This demonstrates that cryptic species within parasitoid Hymenoptera also occur in Central Europe in close contact to humans.     

Speciation with gene flow in whiptail lizards from a Neotropical xeric biome

Two main hypotheses have been proposed to explain the diversification of the Caatinga biota. The riverine barrier hypothesis (RBH) claims that the São Francisco River (SFR) is a major biogeographic barrier to gene flow. The Pleistocene climatic fluctuation hypothesis (PCH) states that gene flow, geographic genetic structure and demographic signatures on endemic Caatinga taxa were influenced by Quaternary climate fluctuation cycles. Herein, we analyse genetic diversity and structure, phylogeographic history, and diversification of a widespread Caatinga lizard (Cnemidophorus ocellifer) based on large geographical sampling for multiple loci to test the predictions derived from the RBH and PCH. We inferred two well‐delimited lineages (Northeast and Southwest) that have diverged along the Cerrado–Caatinga border during the Mid‐Late Miocene (6–14 Ma) despite the presence of gene flow. We reject both major hypotheses proposed to explain diversification in the Caatinga. Surprisingly, our results revealed a striking complex diversification pattern where the Northeast lineage originated as a founder effect from a few individuals located along the edge of the Southwest lineage that eventually expanded throughout the Caatinga. The Southwest lineage is more diverse, older and associated with the Cerrado–Caatinga boundaries. Finally, we suggest that C. ocellifer from the Caatinga is composed of two distinct species. Our data support speciation in the presence of gene flow and highlight the role of environmental gradients in the diversification process.     

Phylogeographic structure and deep lineage diversification of the red alga Chondrus ocellatus Holmes in the Northwest Pacific

A major goal of phylogeographic analysis using molecular markers is to understand the ecological and historical variables that influence genetic diversity within a species. Here, we used sequences of the mitochondrial Cox1 gene and nuclear internal transcribed spacer to reconstruct its phylogeography and demographic history of the intertidal red seaweed Chondrus ocellatus over most of its geographical range in the Northwest Pacific. We found three deeply separated lineages A, B and C, which diverged from one another in the early Pliocene–late Miocene (c. 4.5–7.7 Ma). The remarkably deep divergences, both within and between lineages, appear to have resulted from ancient isolations, accelerated by random drift and limited genetic exchange between regions. The disjunct distributions of lineages A and C along the coasts of Japan may reflect divergence during isolation in scattered refugia. The distribution of lineage B, from the South China Sea to the Korean Peninsula, appears to reflect postglacial recolonizations of coastal habitats. These three lineages do not coincide with the three documented morphological formae in C. ocellatus, suggesting that additional cryptic species may exist in this taxon. Our study illustrates the interaction of environmental variability and demographic processes in producing lineage diversification in an intertidal seaweed and highlights the importance of phylogeographic approaches for discovering cryptic marine biodiversity.     

The global impact of Wolbachia on mitochondrial diversity and evolution

The spread of maternally inherited microorganisms, such as Wolbachia bacteria, can induce indirect selective sweeps on host mitochondria, to which they are linked within the cytoplasm. The resulting reduction in effective population size might lead to smaller mitochondrial diversity and reduced efficiency of natural selection. While documented in several host species, it is currently unclear if such a scenario is common enough to globally impact the diversity and evolution of mitochondria in Wolbachia‐infected lineages. Here, we address this question using a mapping of Wolbachia acquisition/extinction events on a large mitochondrial DNA tree, including over 1000 species. Our analyses indicate that on a large phylogenetic scale, other sources of variation, such as mutation rates, tend to hide the effects of Wolbachia. However, paired comparisons between closely related infected and uninfected taxa reveal that Wolbachia is associated with a twofold reduction in silent mitochondrial polymorphism, and a 13% increase in nonsynonymous substitution rates. These findings validate the conjecture that the widespread distribution of Wolbachia infections throughout arthropods impacts the effective population size of mitochondria. These effects might in part explain the disconnection between genetic diversity and demographic population size in mitochondria, and also fuel red‐queen‐like cytonuclear co‐evolution through the fixation of deleterious mitochondrial alleles.     

Limited female dispersal predicts the incidence of Wolbachia across ants (Hymenoptera: Formicidae)

The endosymbiotic bacterium Wolbachia is perhaps the greatest panzootic in the history of life on Earth, yet remarkably little is known regarding the factors that determine its incidence across species. One possibility is that Wolbachia more easily invades species with structured populations, due to the increased strength of genetic drift and higher initial frequency of infection. This should enable strains that induce mating incompatibilities to more easily cross the threshold prevalence above which they spread to either fixation or a stable equilibrium infection prevalence. Here, we provide empirical support for this hypothesis by analysing the relationship between female dispersal (as a proxy for population structure) and the incidence of Wolbachia across 250 species of ants. We show that species in which the dispersal of reproductive females is limited are significantly more likely to be infected with Wolbachia than species whose reproductive ecology is consistent with significant dispersal of females, and that this relationship remains after controlling for host phylogeny. We suggest that structured host populations, in this case resulting from limited female dispersal, may be an important feature determining how easily Wolbachia becomes successfully established in a novel host, and thus its occurrence across a wide diversity of invertebrate hosts.     

Widespread mito‐nuclear discordance with evidence for introgressive hybridization and selective sweeps in Lycaeides

We investigated the extent and potential cause(s) of mitochondrial introgression within the polytypic North American Lycaeides species complex (Lepidoptera). By comparing population genetic structure based on mitochondrial DNA (COI and COII) and nuclear DNA (251 polymorphic amplified fragment length polymorphism markers), we detected substantial mito‐nuclear discordance, primarily involving a single mitochondrial haplotype (h01), which is likely due to mitochondrial introgression between differentiated Lycaeides populations and/or species. We detected reduced mitochondrial genetic diversity relative to nuclear genetic diversity in populations where mitochondrial haplotype h01 occurs, suggesting that the spread of this haplotype was facilitated by selection. We found no evidence that haplotype h01 is associated with increased fitness (in terms of survival to eclosion, fresh adult weight, and adult longevity) in a polymorphic Lycaeides melissa population. However, we did find a positive association between mitochondrial haplotype h01 and infection by the endoparasitic bacterium Wolbachia in one out of three lineages tested. Linkage disequilibrium between mitochondrial haplotype h01 and Wolbachia infection status may have resulted in indirect selection favouring the spread of haplotype h01 in at least one lineage of North American Lycaeides. These results illustrate the potential for introgressive hybridization to produce substantial mito‐nuclear discordance and demonstrate that an individual's mitochondrial and nuclear genome may have strikingly different evolutionary histories resulting from non‐neutral processes and intrinsic differences in the inheritance and biology of these genomes.     

Radiation of the polymorphic Little Devil poison frog (Oophaga sylvatica) in Ecuador

Some South American poison frogs (Dendrobatidae) are chemically defended and use bright aposematic colors to warn potential predators of their unpalatability. Aposematic signals are often frequency‐dependent where individuals deviating from a local model are at a higher risk of predation. However, extreme diversity in the aposematic signal has been documented in poison frogs, especially in Oophaga. Here, we explore the phylogeographic pattern among color‐divergent populations of the Little Devil poison frog Oophaga sylvatica by analyzing population structure and genetic differentiation to evaluate which processes could account for color diversity within and among populations. With a combination of PCR amplicons (three mitochondrial and three nuclear markers) and genome‐wide markers from a double‐digested RAD (ddRAD) approach, we characterized the phylogenetic and genetic structure of 199 individuals from 13 populations (12 monomorphic and 1 polymorphic) across the O. sylvatica distribution. Individuals segregated into two main lineages by their northern or southern latitudinal distribution. A high level of genetic and phenotypic polymorphism within the northern lineage suggests ongoing gene flow. In contrast, low levels of genetic differentiation were detected among the southern lineage populations and support recent range expansions from populations in the northern lineage. We propose that a combination of climatic gradients and structured landscapes might be promoting gene flow and phylogenetic diversification. Alternatively, we cannot rule out that the observed phenotypic and genomic variations are the result of genetic drift on near or neutral alleles in a small number of genes.     

Geographical features are the predominant driver of molecular diversification in widely distributed North American whipsnakes

Allopatric divergence following the formation of geographical features has been implicated as a major driver of evolutionary diversification. Widespread species complexes provide opportunities to examine allopatric divergence across varying degrees of isolation in both time and space. In North America, several geographical features may play such a role in diversification, including the Mississippi River, Pecos River, Rocky Mountains, Cochise Filter Barrier, Gulf of California and Isthmus of Tehuantepec. We used thousands of nuclear single nucleotide polymorphisms (SNPs) and mitochondrial DNA from several species of whipsnakes (genera Masticophis and Coluber) distributed across North and Central America to investigate the role that these geographical features have played on lineage divergence. We hypothesize that these features restrict gene flow and separate whipsnakes into diagnosable genomic clusters. We performed genomic clustering and phylogenetic reconstructions at the species and population levels using Bayesian and likelihood analyses and quantified migration levels across geographical features to assess the degree of genetic isolation due to allopatry. Our analyses suggest that (i) major genetic divisions are often consistent with isolation by geographical features, (ii) migration rates between clusters are asymmetrical across major geographical features, and (iii) areas that receive proportionally more migrants possess higher levels of genetic diversity. Collectively, our findings suggest that multiple features of the North American landscape contributed to allopatric divergence in this widely distributed snake group.