Genome‐wide differentiation in closely related populations: the roles of selection and geographic isolation

Molecular and acoustic evidence support the species status of Anthus rubescens rubescens and Anthus [rubescens] japonicus (Passeriformes: Motacillidae)

The Buff-bellied Pipit Anthus rubescens comprises two allopatric subspecies groups: A. r. rubescens and A. r. alticola in North America and A. [r.] japonicus in north-east Asia. Despite their great morphological resemblance in breeding plumage, most individuals can be assigned to one or the other subspecies group in non-breeding plumage. Allopatric distributions, morphological differentiation and previously reported molecular divergence suggested the need for additional taxonomic study to assess the rank of these two populations. To resolve the taxonomy of the Buff-bellied Pipit species complex we analysed i) two mitochondrial DNA (mtDNA) loci and ii) nine bioacoustic parameters across 69 sound recordings (338 flight calls) recovered from public databases using principal component analysis and Euclidean distance measures. By comparing our mtDNA and call divergence measures with similar values measured between long-recognised species pairs of the genus, we show that the level of mitochondrial and acoustic divergence between the two Buff-bellied Pipit subspecies groups is typical of species-level divergence in the genus Anthus. Therefore, we recommend splitting the Buff-bellied Pipit species complex into two species: Anthus rubescens (American Pipit) and Anthus japonicus (Siberian Pipit). Our results also suggest that the Water Pipit A. spinoletta deserves taxonomic reassessment as its lineages are highly divergent in acoustics and mtDNA, while mtDNA relationships suggest paraphyly relative to the Rock Pipit A. petrosus. Our work highlights the crucial importance of integrative approaches in taxonomy and the usefulness of bioacoustics in studying cryptic diversity.

A chromosome-level reference genome and pangenome for barn swallow population genomics

Insights into the evolution of non-model organisms are limited by the lack of reference genomes of high accuracy, completeness, and contiguity. Here, we present a chromosome-level, karyotype-validated reference genome and pangenome for the barn swallow (Hirundo rustica). We complement these resources with a reference-free multialignment of the reference genome with other bird genomes and with the most comprehensive catalog of genetic markers for the barn swallow. We identify potentially conserved and accelerated genes using the multialignment and estimate genome-wide linkage disequilibrium using the catalog. We use the pangenome to infer core and accessory genes and to detect variants using it as a reference. Overall, these resources will foster population genomics studies in the barn swallow, enable detection of candidate genes in comparative genomics studies, and help reduce bias toward a single reference genome.

The Mitogenome Relationships and Phylogeography of Barn Swallows (Hirundo rustica)

The barn swallow (Hirundo rustica) poses a number of fascinating scientific questions, including the taxonomic status of postulated subspecies. Here we obtained and assessed the sequence variation of 411 complete mitogenomes, mainly from the European H. r. rustica, but other subspecies as well. In almost every case, we observed subspecies-specific haplogroups, which we employed together with estimated radiation times to postulate a model for the geographical and temporal worldwide spread of the species. The female barn swallow carrying the Hirundo rustica ancestral mitogenome left Africa (or its vicinity) around 280 thousand years ago (kya), and her descendants expanded first into Eurasia and then, at least 51 kya, into the Americas, from where a relatively recent (< 20 kya) back migration to Asia took place. The exception to the haplogroup subspecies specificity is represented by the sedentary Levantine H. r. transitiva that extensively shares haplogroup A with the migratory European H. r. rustica and, to a lesser extent, haplogroup B with the Egyptian H. r. savignii. Our data indicate that rustica and transitiva most likely derive from a sedentary Levantine population source that split at the end of the Younger Dryas (11.7 kya). Since then, however, transitiva received genetic inputs from and admixed with both the closely related rustica and the adjacent savignii. Demographic analyses confirm this species' strong link with climate fluctuations and human activities making it an excellent indicator for monitoring and assessing the impact of current global changes on wildlife.

A migratory divide spanning two continents is associated with genomic and ecological divergence

Migratory divides are contact zones between breeding populations with divergent migratory strategies during the nonbreeding season. These locations provide an opportunity to evaluate the role of seasonal migration in the maintenance of reproductive isolation, particularly the relationship between population structure and features associated with distinct migratory strategies. We combine light-level geolocators, genomic sequencing, and stable isotopes to investigate the timing of migration and migratory routes of individuals breeding on either side of a migratory divide coinciding with genomic differentiation across a hybrid zone between barn swallow (Hirundo rustica) subspecies in China. Individuals west of the hybrid zone, with H. r. rustica ancestry, had comparatively enriched stable-carbon and hydrogen isotope values and overwintered in eastern Africa, whereas birds east of the hybrid zone, with H. r. gutturalis ancestry, had depleted isotope values and migrated to southern India. The two subspecies took divergent migratory routes around the high-altitude Karakoram Range and arrived on the breeding grounds over 3 weeks apart. These results indicate that assortative mating by timing of arrival and/or selection against hybrids with intermediate migratory traits may maintain reproductive isolation between the subspecies, and that inhospitable geographic features may have contributed to the diversification of Asian avifauna by influencing migratory patterns.

Neo-sex chromosome evolution and phenotypic differentiation across an elevational gradient in horned larks (Eremophila Alpestris)

Genetic structure and phenotypic variation among populations is affected by both geographic distance and environmental variation across species' distributions. Understanding the relative contributions of isolation by distance (IBD) and isolation by environment (IBE) is important for elucidating population dynamics across habitats and ecological gradients. In this study, we compared phenotypic and genetic variation among Horned Lark (Eremophila alpestris) populations from 10 sites encompassing an elevational gradient from low-elevation desert scrub in Death Valley (285 a.s.l.) to high-elevation meadows in the White Mountains of the Sierra Nevada of California (greater than 3000 m a.s.l.). Using a ddRAD dataset of 28,474 SNPs aligned to a high-quality reference genome, we compared genetic structure with elevational, environmental, and spatial distance to quantify how different aspects of the landscape drive genomic and phenotypic differentiation in Horned Larks. We found larger-bodied birds were associated with sites that had less seasonality and higher annual precipitation, and longer spurs occurred in soils with more clay and silt content, less sand, and finer fragments. Larks have large neo-sex chromosomes, and we found that associations with elevation and environmental variation were much stronger among neo-sex chromosomes compared to autosomes. Furthermore, we found that putative chromosomal translocations, fusions, and inversions were associated with elevation and may underlie local adaptation across an elevational gradient in Horned Larks. Our results suggest that genetic variation in Horned Larks is affected more by IBD than IBE, but specific phenotypes and genomic regions-particually on neo-sex chromosomes-bear stronger associations with the environment.

Low genome-wide divergence between two lizard populations with high adaptive phenotypic differentiation

Usually, adaptive phenotypic differentiation is paralleled by genetic divergence between locally adapted populations. However, adaptation can also happen in a scenario of nonsignificant genetic divergence due to intense gene flow and/or recent differentiation. While this phenomenon is rarely published, findings on incipient ecologically driven divergence or isolation by adaptation are relatively common, which could confound our understanding about the frequency at which they actually occur in nature. Here, we explore genome-wide traces of divergence between two populations of the lacertid lizard Psammodromus algirus separated by a 600 m elevational gradient. These populations seem to be differentially adapted to their environments despite showing low levels of genetic differentiation (according to previously studies of mtDNA and microsatellite data). We performed a search for outliers (i.e., loci subject to selection) trying to identify specific loci with FST statistics significantly higher than those expected on the basis of overall, genome-wide estimates of genetic divergence. We find that local phenotypic adaptation (in terms of a wide diversity of characters) was not accompanied by genome-wide differentiation, even when we maximized the chances of unveiling such differentiation at particular loci with FST-based outlier detection tests. Instead, our analyses confirmed the lack of genome-wide differentiation on the basis of more than 70,000 SNPs, which is concordant with a scenario of local adaptation without isolation by environment. Our results add evidence to previous studies in which local adaptation does not lead to any kind of isolation (or early stages of ecological speciation), but maintains phenotypic divergence despite the lack of a differentiated genomic background.

Lack of subspecies-recognition in breeding Barn Swallows (Hirundo rustica transitiva)

Assortative social interactions based on (sub)species recognition can be a driving force in speciation processes. To determine whether breeding Barn Swallows Hirundo rustica transitiva in Israel behave differentially towards members of their own subspecies, relative to a different, transient subspecies H. r. rustica and two sympatrically breeding species (Sand Martin Riparia riparia and House Sparrow Passer domesticus), we conducted a territory intrusion experiment near active nests using taxidermy models. Females responded less to the models than males, and the patterns of the recorded behavioral response traits co-varied statistically with sub- or species identity of the models, but none showed patterns of response selectivity for con(sub)specific model types only. These results do not support a role for subspecies recognition in the territorial intrusion responses of H. r. transitiva.

Sex-linked genetic diversity and differentiation in a globally distributed avian species complex

Sex chromosomes often bear distinct patterns of genetic variation due to unique patterns of inheritance and demography. The processes of mutation, recombination, genetic drift and selection also influence rates of evolution on sex chromosomes differently than autosomes. Measuring such differences provides information about how these processes shape genomic variation and their roles in the origin of species. To test hypotheses and predictions about patterns of autosomal and sex-linked genomic diversity and differentiation, we measured population genetic statistics within and between populations and subspecies of the barn swallow (Hirundo rustica) and performed explicit comparisons between autosomal and Z-linked genomic regions. We first tested for evidence of low Z-linked genetic diversity and high Z-linked population differentiation relative to autosomes, then for evidence that the Z chromosome bears greater ancestry information due to faster lineage sorting. Finally, we investigated geographical clines across hybrid zones for evidence that the Z chromosome is resistant to introgression due to selection against hybrids. We found evidence that the barn swallow mating system, demographic history and linked selection each contribute to low Z-linked diversity and high Z-linked differentiation. While incomplete lineage sorting is rampant across the genome, our results indicate faster sorting of ancestral polymorphism on the Z. Finally, hybrid zone analyses indicate barriers to introgression on the Z chromosome, suggesting that sex-linked traits are important in reproductive isolation, especially in migratory divide regions. Our study highlights how selection, gene flow and demography shape sex-linked genetic diversity and underlines the relevance of the Z chromosome in speciation.

Genomic variation, population history and within-archipelago adaptation between island bird populations

Oceanic island archipelagos provide excellent models to understand evolutionary processes. Colonization events and gene flow can interact with selection to shape genetic variation at different spatial scales. Landscape-scale variation in biotic and abiotic factors may drive fine-scale selection within islands, while long-term evolutionary processes may drive divergence between distantly related populations. Here, we examine patterns of population history and selection between recently diverged populations of the Berthelot's pipit (Anthus berthelotii), a passerine endemic to three North Atlantic archipelagos. First, we use demographic trees and f3 statistics to show that genome-wide divergence across the species range is largely shaped by colonization and bottlenecks, with evidence of very weak gene flow between populations. Then, using a genome scan approach, we identify signatures of divergent selection within archipelagos at single nucleotide polymorphisms (SNPs) in genes potentially associated with craniofacial development and DNA repair. We did not detect within-archipelago selection at the same SNPs as were detected previously at broader spatial scales between archipelagos, but did identify signatures of selection at loci associated with similar biological functions. These findings suggest that similar ecological factors may repeatedly drive selection between recently separated populations, as well as at broad spatial scales across varied landscapes.

Phylogeography of the iconic Australian pink cockatoo, Lophochroa leadbeateri

The pink cockatoo (Lophochroa leadbeateri; or Major Mitchell’s cockatoo) is one of Australia’s most iconic bird species. Two subspecies based on morphology are separated by a biogeographical divide, the Eyrean Barrier. Testing the genetic basis for this subspecies delineation, clarifying barriers to gene flow and identifying any cryptic genetic diversity will likely have important implications for conservation and management. Here, we used genome-wide single nucleotide polymorphisms (SNPs) and mitochondrial DNA data to conduct the first range-wide genetic assessment of the species. The aims were to investigate the phylogeography of the pink cockatoo, to characterize conservation units and to reassess subspecies boundaries. We found consistent but weak genetic structure between the two subspecies based on nuclear SNPs. However, phylogenetic analysis of nuclear SNPs and mitochondrial DNA sequence data did not recover reciprocally monophyletic groups, indicating incomplete evolutionary separation between the subspecies. Consequently, we have proposed that the two currently recognized subspecies be treated as separate management units rather than evolutionarily significant units. Given that poaching is suspected to be a threat to this species, we assessed the utility of our data for wildlife forensic applications. We demonstrated that a subspecies identification test could be designed using as few as 20 SNPs.

Migration, pathogens and the avian microbiome: A comparative study in sympatric migrants and residents

Animals generally benefit from their gastrointestinal microbiome, but the factors that influence the composition and dynamics of their microbiota remain poorly understood. Studies of nonmodel host species can illuminate how microbiota and their hosts interact in natural environments. We investigated the role of migratory behaviour in shaping the gut microbiota of free-ranging barn swallows (Hirundo rustica) by studying co-occurring migrant and resident subspecies sampled during the autumn migration at a migratory bottleneck. We found that within-host microbial richness (α-diversity) was similar between migrant and resident microbial communities. In contrast, we found that microbial communities (β-diversity) were significantly different between groups regarding both microbes present and their relative abundances. Compositional differences were found for 36 bacterial genera, with 27 exhibiting greater abundance in migrants and nine exhibiting greater abundance in residents. There was heightened abundance of Mycoplasma spp. and Corynebacterium spp. in migrants, a pattern shared by other studies of migratory species. Screens for key regional pathogens revealed that neither residents nor migrants carried avian influenza viruses and Newcastle disease virus, suggesting that the status of these diseases did not underlie observed differences in microbiome composition. Furthermore, the prevalence and abundance of Salmonella spp., as determined from microbiome data and cultural assays, were both low and similar across the groups. Overall, our results indicate that microbial composition differs between migratory and resident barn swallows, even when they are conspecific and sympatrically occurring. Differences in host origins (breeding sites) may result in microbial community divergence, and varied behaviours throughout the annual cycle (e.g., migration) could further differentiate compositional structure as it relates to functional needs.

ELIXIR-IT HPC@CINECA: high performance computing resources for the bioinformatics community

BACKGROUND

The advent of Next Generation Sequencing (NGS) technologies and the concomitant reduction in sequencing costs allows unprecedented high throughput profiling of biological systems in a cost-efficient manner. Modern biological experiments are increasingly becoming both data and computationally intensive and the wealth of publicly available biological data is introducing bioinformatics into the "Big Data" era. For these reasons, the effective application of High Performance Computing (HPC) architectures is becoming progressively more recognized also by bioinformaticians. Here we describe HPC resources provisioning pilot programs dedicated to bioinformaticians, run by the Italian Node of ELIXIR (ELIXIR-IT) in collaboration with CINECA, the main Italian supercomputing center.

RESULTS

Starting from April 2016, CINECA and ELIXIR-IT launched the pilot Call "ELIXIR-IT HPC@CINECA", offering streamlined access to HPC resources for bioinformatics. Resources are made available either through web front-ends to dedicated workflows developed at CINECA or by providing direct access to the High Performance Computing systems through a standard command-line interface tailored for bioinformatics data analysis. This allows to offer to the biomedical research community a production scale environment, continuously updated with the latest available versions of publicly available reference datasets and bioinformatic tools. Currently, 63 research projects have gained access to the HPC@CINECA program, for a total handout of ~ 8 Millions of CPU/hours and, for data storage, ~ 100 TB of permanent and ~ 300 TB of temporary space.

CONCLUSIONS

Three years after the beginning of the ELIXIR-IT HPC@CINECA program, we can appreciate its impact over the Italian bioinformatics community and draw some considerations. Several Italian researchers who applied to the program have gained access to one of the top-ranking public scientific supercomputing facilities in Europe. Those investigators had the opportunity to sensibly reduce computational turnaround times in their research projects and to process massive amounts of data, pursuing research approaches that would have been otherwise difficult or impossible to undertake. Moreover, by taking advantage of the wealth of documentation and training material provided by CINECA, participants had the opportunity to improve their skills in the usage of HPC systems and be better positioned to apply to similar EU programs of greater scale, such as PRACE. To illustrate the effective usage and impact of the resources awarded by the program - in different research applications - we report five successful use cases, which have already published their findings in peer-reviewed journals.

Change in sexual signalling traits outruns morphological divergence across an ecological gradient in the post-glacial radiation of the songbird genus Junco

The relative roles of natural and sexual selection in promoting evolutionary lineage divergence remains controversial and difficult to assess in natural systems. Local adaptation through natural selection is known to play a central role in promoting evolutionary divergence, yet secondary sexual traits can vary widely among species in recent radiations, suggesting that sexual selection may also be important in the early stages of speciation. Here, we compare rates of divergence in ecologically relevant traits (morphology) and sexually selected signalling traits (coloration) relative to neutral structure in genome-wide molecular markers and examine patterns of variation in sexual dichromatism to explore the roles of natural and sexual selection in the diversification of the songbird genus Junco (Aves: Passerellidae). Juncos include divergent lineages in Central America and several dark-eyed junco (J. hyemalis) lineages that diversified recently as the group recolonized North America following the last glacial maximum (ca. 18,000 years ago). We found an accelerated rate of divergence in sexually selected characters relative to ecologically relevant traits. Moreover, sexual dichromatism measurements suggested a positive relationship between the degree of colour divergence and the strength of sexual selection when controlling for neutral genetic distance. We also found a positive correlation between dichromatism and latitude, which coincides with the geographic axis of decreasing lineage age in juncos but also with a steep ecological gradient. Finally, we found significant associations between genome-wide variants linked to functional genes and proxies of both sexual and natural selection. These results suggest that the joint effects of sexual and ecological selection have played a prominent role in the junco radiation.

Divergent sexual signals reflect costs of local parasites

Many closely related populations are distinguished by variation in sexual signals and this variation is hypothesized to play an important role in reproductive isolation and speciation. Within populations, there is considerable evidence that sexual signals provide information about the incidence and severity of parasite infections, but it remains unclear if variation in parasite communities across space could play a role in initiating or maintaining sexual trait divergence. To test for variation in parasite-associated selection, we compared three barn swallow subspecies with divergent sexual signals. We found that parasite community structure and host tolerance to ecologically similar parasites varied between subspecies. Across subspecies we also found that different parasites were costly in terms of male survival and reproductive success. For each subspecies, the preferred sexual signal(s) were associated with the most costly local parasite(s), indicating that divergent signals are providing relevant information to females about local parasite communities. Across subspecies, the same traits were often associated with different parasites, indicating that parasite-sexual signal links are quite flexible and may evolve relatively quickly. This study provides evidence for (1) variation in parasite communities and (2) different parasite-sexual signal links among three closely related subspecies with divergent sexual signal traits, suggesting that parasites may play an important role in initiating and/or maintaining the divergence of sexual signals among these closely related, yet geographically isolated populations.

Senecio as a model system for integrating studies of genotype, phenotype and fitness

Two major developments have made it possible to use examples of ecological radiations as model systems to understand evolution and ecology. First, the integration of quantitative genetics with ecological experiments allows detailed connections to be made between genotype, phenotype, and fitness in the field. Second, dramatic advances in molecular genetics have created new possibilities for integrating field and laboratory experiments with detailed genetic sequencing. Combining these approaches allows evolutionary biologists to better study the interplay between genotype, phenotype, and fitness to explore a wide range of evolutionary processes. Here, we present the genus Senecio (Asteraceae) as an excellent system to integrate these developments, and to address fundamental questions in ecology and evolution. Senecio is one of the largest and most phenotypically diverse genera of flowering plants, containing species ranging from woody perennials to herbaceous annuals. These Senecio species exhibit many growth habits, life histories, and morphologies, and they occupy a multitude of environments. Common within the genus are species that have hybridized naturally, undergone polyploidization, and colonized diverse environments, often through rapid phenotypic divergence and adaptive radiation. These diverse experimental attributes make Senecio an attractive model system in which to address a broad range of questions in evolution and ecology.

SMRT long reads and Direct Label and Stain optical maps allow the generation of a high-quality genome assembly for the European barn swallow (Hirundo rustica rustica)

Background

The barn swallow (Hirundo rustica) is a migratory bird that has been the focus of a large number of ecological, behavioral, and genetic studies. To facilitate further population genetics and genomic studies, we present a reference genome assembly for the European subspecies (H. r. rustica).

Findings

As part of the Genome10K effort on generating high-quality vertebrate genomes (Vertebrate Genomes Project), we have assembled a highly contiguous genome assembly using single molecule real-time (SMRT) DNA sequencing and several Bionano optical map technologies. We compared and integrated optical maps derived from both the Nick, Label, Repair, and Stain technology and from the Direct Label and Stain (DLS) technology. As proposed by Bionano, DLS more than doubled the scaffold N50 with respect to the nickase. The dual enzyme hybrid scaffold led to a further marginal increase in scaffold N50 and an overall increase of confidence in the scaffolds. After removal of haplotigs, the final assembly is approximately 1.21 Gbp in size, with a scaffold N50 value of more than 25.95 Mbp.

Conclusions

This high-quality genome assembly represents a valuable resource for future studies of population genetics and genomics in the barn swallow and for studies concerning the evolution of avian genomes. It also represents one of the very first genomes assembled by combining SMRT long-read sequencing with the new Bionano DLS technology for scaffolding. The quality of this assembly demonstrates the potential of this methodology to substantially increase the contiguity of genome assemblies.

Genome-wide signals of drift and local adaptation during rapid lineage divergence in a songbird

The formation of independent evolutionary lineages involves neutral and selective factors, and understanding their relative roles in population divergence is a fundamental goal of speciation research. Correlations between allele frequencies and environmental variability can reveal the role of selection, yet the relative contribution of drift can be difficult to establish. Recently diversified taxa like the Oregon junco (Aves, Passerellidae, Junco hyemalis oreganus) of western North America provide ideal scenarios to apply genetic-environment association analyses (GEA) while controlling for population structure. Analysis of genome-wide SNP loci revealed marked genetic structure consisting of differentiated populations in isolated, dry southern mountain ranges, and less divergent, recently expanded populations in humid northern latitudes. We used correlations between genomic and environmental variance to test for three specific modes of evolutionary divergence: (a) drift in geographic isolation, (b) differentiation along continuous selective gradients and (c) isolation-by-adaptation. We found evidence of strong drift in southern mountains, but also signals of local adaptation driven by temperature, precipitation, elevation and vegetation, especially when controlling for population history. We identified numerous variants under selection scattered across the genome, suggesting that local adaptation can promote rapid differentiation when acting over multiple independent loci.

Demographic inference in barn swallows using whole-genome data shows signal for bottleneck and subspecies differentiation during the Holocene

Accounting for historical demographic features is vital for many types of evolutionary inferences, including the estimation of divergence times between closely related populations. In barn swallow, Hirundo rustica, inferring historical population sizes and subspecies divergence times can shed light on the recent co-evolution of this species with humans. Pairwise sequentially Markovian coalescent uncovered population growth beginning on the order of one million years ago-which may reflect the radiation of the broader Hirundo genus-and a more recent population decline. Additionally, we used approximate Bayesian computation to evaluate hypotheses about recent timescale barn swallow demography, including population growth due to human commensalism, and a potential founder event associated with the onset of nesting on human structures. We found signal for a bottleneck event approximately 7,700 years ago, near the time that humans began building substantial structures, although there was considerable uncertainty associated with this estimate. Subspecies differentiation and subsequent growth occurred after the bottleneck in the best-supported model, an order of magnitude more recently than previous estimates in this system. We also compared results obtained from whole-genome sequencing versus reduced representation sequencing, finding many similar results despite substantial allelic dropout in the reduced representation data, which may have affected estimates of some parameters. This study presents the first genetic evidence of a potential barn swallow founder effect and subspecies divergence coinciding with the Holocene, which is an important step in analysing the biogeographical history of a well-known human commensal species.

Environmental correlates of geographic divergence in a phenotypic trait: A case study using bat echolocation

Divergence in phenotypic traits may arise from the interaction of different evolutionary forces, including different kinds of selection (e.g., ecological), genetic drift, and phenotypic plasticity. Sensory systems play an important role in survival and reproduction, and divergent selection on such systems may result in lineage diversification. Such diversification could be largely influenced by selection in different environments as a result of isolation by environment (IbE). We investigated this process using geographic variation in the resting echolocation frequency of the horseshoe bat species, Rhinolophus damarensis, as a test case. Bats were sampled along a latitudinal gradient ranging from 16°S to 32°S in the arid western half of southern Africa. We measured body size and peak resting frequencies (RF) from handheld individual bats. Three hypotheses for the divergence in RF were tested: (1) James' Rule, (2) IbE, and (3) genetic drift through isolation by distance (IbD) to isolate the effects of body size, local climatic conditions, and geographic distance, respectively, on the resting frequency of R. damarensis. Our results did not support genetic drift because there was no correlation between RF variation and geographic distance. Our results also did not support James' Rule because there was no significant relationship between (1) geographic distances and RF, (2) body size and RF, or (3) body size and climatic variables. Instead, we found support for IbE in the form of a correlation between RF and both region and annual mean temperature, suggesting that RF variation may be the result of environmental discontinuities. The environmental discontinuities coincided with previously reported genetic divergence. Climatic gradients in conjunction with environmental discontinuities could lead to local adaptation in sensory signals and directed dispersal such that gene flow is restricted, allowing lineages to diverge. However, our study cannot exclude the role of processes like phenotypic plasticity in phenotypic variation.