Current geography masks dynamic history of gene flow during speciation in northern Australian birds

Genome divergence is greatly influenced by gene flow during early stages of speciation. As populations differentiate, geographic barriers can constrain gene flow and so affect the dynamics of divergence and speciation. Current geography, specifically disjunction and continuity of ranges, is often used to predict the historical gene flow during the divergence process. We test this prediction in eight meliphagoid bird species complexes codistributed in four regions. These regions are separated by known biogeographical barriers across northern Australia and Papua New Guinea. We find that bird populations currently separated by terrestrial habitat barriers within Australia and marine barriers between Australia and Papua New Guinea have a range of divergence levels and probability of gene flow not associated with current range connectivity. Instead, geographic distance and historical range connectivity better predict divergence and probability of gene flow. In this dynamic environmental context, we also find support for a nonlinear decrease of the probability of gene flow during the divergence process. The probability of gene flow initially decreases gradually after a certain level of divergence is reached. Its decrease then accelerates until the probability is close to zero. This implies that although geographic connectivity may have more of an effect early in speciation, other factors associated with higher divergence may play a more important role in influencing gene flow midway through and later in speciation. Current geographic connectivity may then mislead inferences regarding potential for gene flow during speciation under a complex and dynamic history of geographic and reproductive isolation.     

The relative influence of natural selection and geography on gene flow in guppies

Two general processes may influence gene flow among populations. One involves divergent selection, wherein the maladaptation of immigrants and hybrids impedes gene flow between ecological environments (i.e. ecological speciation). The other involves geographic features that limit dispersal. We determined the relative influence of these two processes in natural populations of Trinidadian guppies (Poecilia reticulata). If selection is important, gene flow should be reduced between different selective environments. If geography is important, gene flow should be impeded by geographic distance and physical barriers. We examined how genetic divergence, long-term gene flow, and contemporary dispersal within a watershed were influenced by waterfalls, geographic distance, predation, and habitat features. We found that waterfalls and geographic distance increased genetic divergence and reduced dispersal and long-term gene flow. Differences in predation or habitat features did not influence genetic divergence or gene flow. In contrast, differences in predation did appear to reduce contemporary dispersal. We suggest that the standard predictions of ecological speciation may be heavily nuanced by the mating behaviour and life history strategies of guppies.     

Refugial isolation versus ecological gradients

Hypotheses for divergence and speciation in rainforests generally fall into two categories: those emphasizing the role of geographic isolation and those emphasizing the role of divergent selection along gradients. While a majority of studies have attempted to infer mechanisms based on the pattern of species richness and congruence of geographic boundaries, relatively few have tried to simultaneously test alternative hypotheses for diversification. Here we discuss four examples, taken from our work on diversification of tropical rainforest vertebrates, in which we examine patterns of genetic and morphological variation within and between biogeographic regions to address two alternative hypotheses. By estimating morphological divergence between geographically contiguous and isolated populations under similar and different ecological conditions, we attempt to evaluate the relative roles of geographic isolation and natural selection in population divergence. Results suggest that natural selection, even in the presence of appreciable gene flow, can result in morphological divergence that is greater than that found between populations isolated for millions of years and, in some cases, even greater than that found between congeneric, but distinct, species. The relatively small phenotypic divergence that occurs among long-term geographic isolates in similar habitats suggests that morphological divergence via drift may be negligible and/or that selection is acting to produce similar phenotypes in populations occupying similar habitats. Our results demonstrate that significant phenotypic divergence: (1) is not necessarily coupled with divergence in neutral molecular markers; and (2) can occur without geographic isolation in the presence of gene flow.     

Divergence in coloration and ecological speciation in the Anolis marmoratus species complex

Adaptive divergence in coloration is expected to produce reproductive isolation in species that use colourful signals in mate choice and species recognition. Indeed, many adaptive radiations are characterized by differentiation in colourful signals, suggesting that divergent selection acting on coloration may be an important component of speciation. Populations in the Anolis marmoratus species complex from the Caribbean island of Guadeloupe display striking divergence in the colour and pattern of adult males that occurs over small geographic distances, suggesting strong divergent selection. Here we test the hypothesis that divergence in coloration results in reduced gene flow among populations. We quantify variation in adult male coloration across a habitat gradient between mesic and xeric habitats, use a multilocus coalescent approach to infer historical demographic parameters of divergence, and examine gene flow and population structure using microsatellite variation. We find that colour variation evolved without geographic isolation and in the face of gene flow, consistent with strong divergent selection and that both ecological and sexual selection are implicated. However, we find no significant differentiation at microsatellite loci across populations, suggesting little reproductive isolation and high levels of contemporary gene exchange. Strong divergent selection on loci affecting coloration probably maintains clinal phenotypic variation despite high gene flow at neutral loci, supporting the notion of a porous genome in which adaptive portions of the genome remain fixed whereas neutral portions are homogenized by gene flow and recombination. We discuss the impact of these findings for studies of colour evolution and ecological speciation.     

Ecological reproductive isolation of coast and inland races of Mimulus guttatus

Adaptive divergence due to habitat differences is thought to play a major role in formation of new species. However it is rarely clear the extent to which individual reproductive isolating barriers related to habitat differentiation contribute to total isolation. Furthermore, it is often difficult to determine the specific environmental variables that drive the evolution of those ecological barriers, and the geographic scale at which habitat-mediated speciation occurs. Here, we address these questions through an analysis of the population structure and reproductive isolation between coastal perennial and inland annual forms of the yellow monkeyflower, Mimulus guttatus. We found substantial morphological and molecular genetic divergence among populations derived from coast and inland habitats. Reciprocal transplant experiments revealed nearly complete reproductive isolation between coast and inland populations mediated by selection against immigrants and flowering time differences, but not postzygotic isolation. Our results suggest that selection against immigrants is a function of adaptations to seasonal drought in inland habitat and to year round soil moisture and salt spray in coastal habitat. We conclude that the coast and inland populations collectively comprise distinct ecological races. Overall, this study suggests that adaptations to widespread habitats can lead to the formation of reproductively isolated species.     

Consequences of divergence and introgression for speciation in Andean cloud forest birds

Divergence with gene flow is well documented and reveals the influence of ecological adaptation on speciation. Yet, it remains intuitive that gene exchange inhibits speciation in many scenarios, particularly among ecologically similar populations. The influence of gene flow on the divergence of populations facing similar selection pressures has received less empirical attention than scenarios where differentiation is coupled with local environmental adaptation. I used a paired study design to test the influence of genomic divergence and introgression on plumage differentiation between ecologically similar allopatric replacements of Andean cloud forest birds. Through analyses of short‐read genome‐wide sequences from over 160 individuals in 16 codistributed lineages, I found that plumage divergence is associated with deep genetic divergence, implicating a prominent role of geographic isolation in speciation. By contrast, lineages that lack plumage divergence across the same geographic barrier are more recently isolated or exhibit a signature of secondary genetic introgression, indicating a negative relationship between gene flow and divergence in phenotypic traits important to speciation. My results suggest that the evolutionary outcomes of cycles of isolation and divergence in this important theatre of biotic diversification are sensitive to time spent in the absence of gene flow.     

Sexual selection drives speciation in an Amazonian frog

One proposed mechanism of speciation is divergent sexual selection, whereby divergence in female preferences and male signals results in behavioural isolation. Despite the appeal of this hypothesis, evidence for it remains inconclusive. Here, we present several lines of evidence that sexual selection is driving behavioural isolation and speciation among populations of an Amazonian frog (Physalaemus petersi). First, sexual selection has promoted divergence in male mating calls and female preferences for calls between neighbouring populations, resulting in strong behavioural isolation. Second, phylogenetic analysis indicates that populations have become fixed for alternative call types several times throughout the species' range, and coalescent analysis rejects genetic drift as a cause for this pattern, suggesting that this divergence is due to selection. Finally, gene flow estimated with microsatellite loci is an average of 30 times lower between populations with different call types than between populations separated by a similar geographical distance with the same call type, demonstrating genetic divergence and incipient speciation. Taken together, these data provide strong evidence that sexual selection is driving behavioural isolation and speciation, supporting sexual selection as a cause for speciation in the wild.     

Spatially explicit models of divergence and genome hitchhiking

Strong barriers to genetic exchange can exist at divergently selected loci, whereas alleles at neutral loci flow more readily between populations, thus impeding divergence and speciation in the face of gene flow. However, ‘divergence hitchhiking’ theory posits that divergent selection can generate large regions of differentiation around selected loci. ‘Genome hitchhiking’ theory suggests that selection can also cause reductions in average genome‐wide rates of gene flow, resulting in widespread genomic divergence (rather than divergence only around specific selected loci). Spatial heterogeneity is ubiquitous in nature, yet previous models of genetic barriers to gene flow have explored limited combinations of spatial and selective scenarios. Using simulations of secondary contact of populations, we explore barriers to gene flow in various selective and spatial contexts in continuous, two‐dimensional, spatially explicit environments. In general, the effects of hitchhiking are strongest in environments with regular spatial patterning of starkly divergent habitat types. When divergent selection is very strong, the absence of intermediate habitat types increases the effects of hitchhiking. However, when selection is moderate or weak, regular (vs. random) spatial arrangement of habitat types becomes more important than the presence of intermediate habitats per se. We also document counterintuitive processes arising from the stochastic interplay between selection, gene flow and drift. Our results indicate that generalization of results from two‐deme models requires caution and increase understanding of the genomic and geographic basis of population divergence.     

Both Geography and Ecology Contribute to Mating Isolation in Guppies

Local adaptation to different environments can promote mating isolation – either as an incidental by-product of trait divergence, or as a result of selection to avoid maladaptive mating. Numerous recent empirical examples point to the common influence of divergent natural selection on speciation based largely on evidence of strong pre-mating isolation between populations from different habitat types. Accumulating evidence for natural selection''s influence on speciation is therefore no longer a challenge. The difficulty, rather, is in determining the mechanisms involved in the progress of adaptive divergence to speciation once barriers to gene flow are already present. Here, we present results of both laboratory and field experiments with Trinidadian guppies (Poecilia reticulata) from different environments, who do not show complete reproductive isolation despite adaptive divergence. We investigate patterns of mating isolation between populations that do and do not exchange migrants and show evidence for both by-product and reinforcement mechanisms depending on female ecology. Specifically, low-predation females discriminate against all high-predation males thus implying a by-product mechanism, whereas high-predation females only discriminate against low-predation males from further upstream in the same river, implying selection to avoid maladaptive mating. Our study thus confirms that mechanisms of adaptive speciation are not necessarily mutually exclusive and uncovers the complex ecology-geography interactions that underlie the evolution of mating isolation in nature.     

VARIABLE PROGRESS TOWARD ECOLOGICAL SPECIATION IN PARAPATRY: STICKLEBACK ACROSS EIGHT LAKE-STREAM TRANSITIONS

Divergent selection between contrasting habitats can sometimes drive adaptive divergence and the evolution of reproductive isolation in the face of initially high gene flow. "Progress" along this ecological speciation pathway can range from minimal divergence to full speciation. We examine this variation for threespine stickleback fish that evolved independently across eight lake-stream habitat transitions. By quantifying stickleback diets, we show that lake-stream transitions usually coincide with limnetic-benthic ecotones. By measuring genetically based phenotypes, we show that these ecotones often generate adaptive divergence in foraging morphology. By analyzing neutral genetic markers (microsatellites), we show that adaptive divergence is often associated with the presence of two populations maintaining at least partial reproductive isolation in parapatry. Coalescent-based simulations further suggest that these populations have diverged with gene flow within a few thousand generations, although we cannot rule out the possibility of phases of allopatric divergence. Finally, we find striking variation among the eight lake-stream transitions in progress toward ecological speciation. This variation allows us to hypothesize that progress is generally promoted by strong divergent selection and limited dispersal across the habitat transitions. Our study thus makes a case for ecological speciation in a parapatric context, while also highlighting variation in the outcome.     

Evolutionary radiation of "stone plants" in the genus Argyroderma (Aizoaceae): unraveling the effects of landscape, habitat, and flowering time

Recent phylogenetic evidence suggests that the extraordinary diversity of the Cape Floristic Kingdom in South Africa may be the result of widespread evolutionary radiation. Our understanding of the role of adaptive versus neutral processes in these radiations remains largely speculative. In this study we investigated factors involved in the diversification of Argyroderma, a genus within the most spectacular of the Cape radiations, that of the Ruschioid subfamily of the Aizoaceae. We used amplified fragment length polymorphisms and a suite of morphological traits to elucidate patterns of differentiation within and between species of Argyroderma across the range of the genus. We then used a matrix correlation approach to assess the influence of landscape structure, edaphic gradients, and flowering phenology on phenotypic and neutral genetic divergence in the system. We found evidence for strong spatial genetic isolation at all taxonomic levels. In addition, genetic differentiation occurs along a temporal axis, between sympatric species with divergent flowering times. Morphological differentiation, which previous studies suggest is adaptive, occurs along a habitat axis, between populations occupying different edaphic microenvironments. Morphological differentiation is in turn significantly associated with flowering time shifts. Thus we propose that diversification within Argyroderma has occurred through a process of adaptive speciation in allopatry. Spatially isolated populations diverge phenotypically in response to divergent habitat selection, which in turn leads to the evolution of reproductive isolation through divergence of flowering phenologies, perhaps as a correlated response to morphological divergence. Evidence suggests that diversification of the group has proceeded in two phases: the first involving divergence of allopatric taxa on varied microhabitats within a novel habitat type (the quartz gravel plains), and the second involving range expansion of an early flowering phenotype on the most extreme edaphic habitat and subsequent incomplete differentiation of allopatric populations of the early flowering group. These results point to adaptive speciation in allopatry as a likely model for the spectacular diversification of the ice-plant family in the dissected landscapes of the southern African winter rainfall deserts.     

Genomic divergence across ecological gradients in the Central African rainforest songbird (Andropadus virens)

The little greenbul, a common rainforest passerine from sub‐Saharan Africa, has been the subject of long‐term evolutionary studies to understand the mechanisms leading to rainforest speciation. Previous research found morphological and behavioural divergence across rainforest–savannah transition zones (ecotones), and a pattern of divergence with gene flow suggesting divergent natural selection has contributed to adaptive divergence and ecotones could be important areas for rainforests speciation. Recent advances in genomics and environmental modelling make it possible to examine patterns of genetic divergence in a more comprehensive fashion. To assess the extent to which natural selection may drive patterns of differentiation, here we investigate patterns of genomic differentiation among populations across environmental gradients and regions. We find compelling evidence that individuals form discrete genetic clusters corresponding to distinctive environmental characteristics and habitat types. Pairwise F_ST between populations in different habitats is significantly higher than within habitats, and this differentiation is greater than what is expected from geographic distance alone. Moreover, we identified 140 SNPs that showed extreme differentiation among populations through a genomewide selection scan. These outliers were significantly enriched in exonic and coding regions, suggesting their functional importance. Environmental association analysis of SNP variation indicates that several environmental variables, including temperature and elevation, play important roles in driving the pattern of genomic diversification. Results lend important new genomic evidence for environmental gradients being important in population differentiation.     

Rapid ecological speciation in three‐spined stickleback Gasterosteus aculeatus from Middleton Island,Alaska: the roles of selection and geographic isolation

Morphological analysis of three‐spined stickleback Gasterosteus aculeatus collected in Middleton Island, Alaska, was conducted in order to study how gene flow and selection interact during divergence. Middleton Island was uplifted by 3·4 m during the Great Alaska Earthquake in 1964; this event formed a series of new freshwater sites, triggering rapid evolution, and probably rapid speciation, in G. aculeatus populations that colonized them. The level of hybridization between the anadromous and the resident freshwater populations is reflected by the level of morphological variance of the resident freshwater G. aculeatus. Therefore, geographic isolation of the sites from the sea (approximating gene flow) and ionic concentration of the water (reflecting selection pressures) were correlated with morphological variance of the resident freshwater populations. Geographic isolation was negatively correlated with morphological variance in a majority of the analysed traits. Both selection and gene flow surrogates were found to be important influences on variance in morphology, though selection had a larger effect, especially on armour traits. It was concluded that gene flow appeared to constrain ecological speciation, but even in the presence of gene flow the strong selection in the freshwater environment was apparently leading to rapid divergence.     

DOES GENE FLOW CONSTRAIN ADAPTIVE DIVERGENCE OR VICE VERSA? A TEST USING ECOMORPHOLOGY AND SEXUAL ISOLATION IN TIMEMA CRISTINAE WALKING‐STICKS

Population differentiation often reflects a balance between divergent natural selection and the opportunity for homogenizing gene flow to erode the effects of selection. However, during ecological speciation, trait divergence results in reproductive isolation and becomes a cause, rather than a consequence, of reductions in gene flow. To assess both the causes and the reproductive consequences of morphological differentiation, we examined morphological divergence and sexual isolation among 17 populations of Timema cristinae walking-sticks. Individuals from populations adapted to using Adenostoma as a host plant tended to exhibit smaller overall body size, wide heads, and short legs relative to individuals using Ceonothus as a host. However, there was also significant variation in morphology among populations within host-plant species. Mean trait values for each single population could be reliably predicted based upon host-plant used and the potential for homogenizing gene flow, inferred from the size of the neighboring population using the alternate host and mitochondrial DNA estimates of gene flow. Morphology did not influence the probability of copulation in between-population mating trials. Thus, morphological divergence is facilitated by reductions in gene flow, but does not cause reductions in gene flow via the evolution of sexual isolation. Combined with rearing data indicating that size and shape have a partial genetic basis, evidence for parallel origins of the host-associated forms, and inferences from functional morphology, these results indicate that morphological divergence in T. cristinae reflects a balance between the effects of host-specific natural selection and gene flow. Our findings illustrate how data on mating preferences can help determine the causal associations between trait divergence and levels of gene flow.     

Ernst Mayr and the integration of geographic and ecological factors in speciation

Mayr's best recognized scientific contributions include the biological species concept and the theory of geographic speciation. In the latter, reproductive isolation evolves as an incidental by‐product of genetic divergence between allopatric populations. Mayr noted that divergent natural selection could accelerate speciation, but also argued that gene flow so strongly retards divergence that, even with selection, non‐allopatric speciation is unlikely. However, current theory and data demonstrate that substantial divergence, and even speciation, in the face of gene flow is possible. Here, I attempt to connect some opposing views about speciation by integrating Mayr's ideas about the roles of ecology and geography in speciation with current data and theory. My central premise is that the speciation process (i.e. divergence) is often continuous, and that the opposing processes of selection and gene flow interact to determine the degree of divergence (i.e. the degree of progress towards the completion of speciation). I first establish that, in the absence of gene flow, divergent selection often promotes speciation. I then discuss how population differentiation in the face of gene flow is common when divergent selection occurs. However, such population differentiation does not always lead to the evolution of discontinuities, strong reproductive isolation, and thus speciation per se. I therefore explore the genetic and ecological circumstances that facilitate speciation in the face of gene flow. For example, particular genetic architectures or ecological niches may tip the balance between selection and gene flow strongly in favour of selection. The circumstances allowing selection to overcome gene flow to the extent that a discontinuity develops, and how often these circumstances occur, are major remaining questions in speciation research. © 2008 The Linnean Society of London, Biological Journal of the Linnean Society, 2008, 95 , 26–46.     

Elevational speciation in action? Restricted gene flow associated with adaptive divergence across an altitudinal gradient

Evolutionary theory predicts that divergent selection pressures across elevational gradients could cause adaptive divergence and reproductive isolation in the process of ecological speciation. Although there is substantial evidence for adaptive divergence across elevation, there is less evidence that this restricts gene flow. Previous work in the boreal chorus frog (Pseudacris maculata) has demonstrated adaptive divergence in morphological, life history and physiological traits across an elevational gradient from approximately 1500–3000 m in the Colorado Front Range, USA. We tested whether this adaptive divergence is associated with restricted gene flow across elevation – as would be expected if incipient speciation were occurring – and, if so, whether behavioural isolation contributes to reproductive isolation. Our analysis of 12 microsatellite loci in 797 frogs from 53 populations revealed restricted gene flow across elevation, even after controlling for geographic distance and topography. Calls also varied significantly across elevation in dominant frequency, pulse number and pulse duration, which was partly, but not entirely, due to variation in body size and temperature across elevation. However, call variation did not result in strong behavioural isolation: in phonotaxis experiments, low‐elevation females tended to prefer an average low‐elevation call over a high‐elevation call, and vice versa for high‐elevation females, but this trend was not statistically significant. In summary, our results show that adaptive divergence across elevation restricts gene flow in P. maculata, but the mechanisms for this potential incipient speciation remain open.     

Gene flow in Dubautia arborea and D. ciliolata: the roles of ecology and isolation by distance in maintaining species boundaries despite ongoing hybridization

The relative roles of gene flow and natural selection in maintaining species differentiation have been a subject of debate for some time. The traditional view is that gene flow constrains adaptive divergence and maintains species cohesiveness. Alternatively, ecological speciation posits that the reverse is true: that adaptive ecological differentiation constrains gene flow. In this study, we examine gene flow and population differentiation among populations of two species of the Hawaiian silversword alliance, Dubautia arborea and D. ciliolata. We compare divergence in putatively neutral microsatellite markers with divergence in leaf morphometric traits, which may be selectively important or physiologically linked to selectively important traits. Gene flow between populations was found to be significant in only one of the two species, D. arborea. Leaf morphometric differentiation between species was significant, though not among populations within species. No evidence of effective genetic introgression was observed between apparently 'pure' populations of these species. Gene flow as measured by microsatellites was not correlated with geographic distance between populations, but was correlated with the linear placement of the widest part of the leaf. Because these two species are interfertile, as demonstrated by the presence of active hybrid zone, the lack of genetic introgression and the maintenance of species boundaries may be associated with natural selection on differential habitat.     

Cryptic speciation in two widespread subterranean amphipod genera reflects historical drainage patterns in an ancient landscape

The landscape of the Pilbara region of Western Australia has been relatively unchanged for 100 million years. The ancient river systems of this region might be expected to be sources of isolation and divergence for aquatic species. Hence, the occurrence of widespread groundwater taxa in this landscape offers the opportunity to examine associations between genetic diversity and drainage patterns. Pilbarus and Chydaekata are two widespread genera of subterranean amphipods endemic to the Pilbara, each occupying multiple tributaries. We used molecular data to examine the roles of drainage patterns in structuring genetic diversity. Gene flow within a tributary may be facilitated by the occasional occurrence of these amphipods in springs, which results in their downstream dispersal during episodic flooding. However, tributary boundaries may form hydrological barriers to gene flow, resulting in localised isolation of populations and divergence. Samples of both genera, collected throughout three river basins, were examined for sequence divergence in the cytochrome c oxidase I mitochondrial gene. There was no evidence of contemporary gene flow among populations of either genus, and each tributary contained highly divergent lineages, which were not associated with similar morphological differentiation. This suggests cryptic speciation has occurred, and similar phylogenetic signals in both taxa imply similar evolutionary histories. Surface populations may have been driven into subterranean refugia by the cessation of flow in the rivers, associated with Tertiary climate change, while morphological evolution may have been constrained by stabilising selection. The lack of congruence between molecular diversity and morphology raises important practical issues for conservation.     

Evidence of recent and continuous speciation in a biodiversity hotspot: a population genetic approach in southern African gladioli (Gladiolus; Iridaceae)

There has been much debate over the origin of species diversity in biodiversity hotspots, particularly the rate of speciation over extinction and the geographic mode of speciation. Here, we looked at speciation with varying degrees of sympatry in a biodiversity hotspot, focusing on a distinct morphological clade in the Cape Floristic Region in southern Africa, the Gladiolus carinatus species complex (Iridaceae). We investigate the mechanisms involved in population and species differentiation through a combination of ecological and genomic approaches. We estimated spatial and phenological overlap, differences in floral morphology, genetic isolation and genomic selection. A genetic coalescent analysis estimated that the time of divergence between lineages followed the establishment of available habitat in the Cape littoral plain where these species currently overlap geographically. Marked shifts in flowering time and morphology, which act as barriers to gene flow, have developed to varying degrees over the last 0.3–1.4 million years. An amplified fragment length polymorphism genome scan revealed signatures of divergent and balancing selection, although half of the loci consistently behaved neutrally. Divergent species outliers (1%) and floral morph outliers (3%) represent a small proportion of the genome, but these loci produced clear genetic clusters of species and significant associations with floral traits. These results indicate that the G. carinatus complex represents a continuum of recent speciation. We provide further evidence for ecological adaptation in the face of gene flow.     

Factors Affecting Geographic Variation in Echolocation Calls of the Endemic Myotis davidii in China

The sensory drive hypothesis of speciation predicts that divergence in communication systems will occur when environments differ and that this sensory divergence can ultimately promote speciation. The factors affecting geographic evolution in acoustic signals remain poorly understood, especially in the contexts of high gene flow. This study investigated variation patterns in peak frequency emitted by the Chinese endemic Myotis davidii on a broad geographic scale by evaluating the relative importance of morphological, environmental, geographic, and genetic variables. Significant variation in peak frequency was observed among regions, but peak frequencies among populations within region had some percentage of similarity. Differences in peak frequency were not associated with morphological difference, genetic structure, and geographic distance among regions, which suggested that peak frequency divergences in M. davidii were not the primary driver of regions' isolation in a context of weak gene flow. Within the Middle East Plain (MEP), one of the regions delineated in this study, peak frequency differences of M. davidii were not significantly correlated with genetic distance and geographic distance among populations, suggesting that peak frequency was not be subject to cultural drift within MEP. Our results provide evidence that geographic variation in echolocation call design may evolve as a consequence of local adaptation to climate conditions.