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

Evolutionarily labile dispersal behavior and discontinuous habitats enhance population differentiation in island versus continentally distributed swallows

The causes of population divergence in vagile groups remain a paradox in evolutionary biology: dispersive species should be able to colonize new areas, a prerequisite for allopatric speciation, but dispersal also facilitates gene flow, which erodes population differentiation. Strong dispersal ability has been suggested to enhance divergence in patchy habitats and inhibit divergence in continuous landscapes, but empirical support for this hypothesis is lacking. Here we compared patterns of population divergence in a dispersive clade of swallows distributed across both patchy and continuous habitats. The Pacific Swallow (Hirundo tahitica) has an insular distribution throughout Southeast Asia and the Pacific, while its sister species, the Welcome Swallow (H. neoxena), has a continental distribution in Australia. We used whole-genome data to demonstrate strong genetic structure and limited introgression among insular populations, but not among continental populations. Demographic models show that historic changes in habitat connectivity have contributed to population structure within the clade. Swallows appear to exhibit evolutionarily labile dispersal behavior in which they reduce dispersal propensity after island colonization despite retaining strong flight ability. Our data support the hypothesis that fragmented habitats enhance population differentiation in vagile groups, and suggest that labile dispersal behavior is a key mechanism underlying this pattern.

Reduced representation approaches produce similar results to whole genome sequencing for some common phylogeographic analyses

When designing phylogeographic investigations researchers can choose to collect many different types of molecular markers, including mitochondrial genes or genomes, SNPs from reduced representation protocols, large sequence capture data sets, and even whole genomes. Given that the statistical power and accuracy of various analyses are expected to differ depending on both the type of marker and the amount of data collected, an exploration of the variance across methodological results as a function of marker type should provide valuable information to researchers. Here we collect mitochondrial Cytochrome b sequences, whole mitochondrial genomes, single nucleotide polymorphisms (SNP)s isolated using a genotype by sequencing (GBS) protocol, sequences from ultraconserved elements, and low-coverage nuclear genomes from the North American water vole (Microtus richardsoni). We estimate genetic distances, population genetic structure, and historical demography using data from each of these datasets and compare the results across markers. As anticipated, the results exhibit differences across marker types, particularly in terms of the resolution offered by different analyses. A cost-benefit analysis indicates that SNPs collected using a GBS protocol are the most cost-effective molecular marker, with inferences that mirror those collected from the whole genome data at a fraction of the cost per sample.

Temporal matches between monarch butterfly and milkweed population changes over the past 25,000 years

In intimate ecological interactions, the interdependency of species may result in correlated demographic histories. For species of conservation concern, understanding the long-term dynamics of such interactions may shed light on the drivers of population decline. Here, we address the demographic history of the monarch butterfly, Danaus plexippus, and its dominant host plant, the common milkweed Asclepias syriaca (A. syriaca), using broad-scale sampling and genomic inference. Because genetic resources for milkweed have lagged behind those for monarchs, we first release a chromosome-level genome assembly and annotation for common milkweed. Next, we show that despite its enormous geographic range across eastern North America, A. syriaca is best characterized as a single, roughly panmictic population. Using approximate Bayesian computation with random forests (ABC-RF), a machine learning method for reconstructing demographic histories, we show that both monarchs and milkweed experienced population expansion during the most recent recession of North American glaciers 10,000-20,000 years ago. Our data also identify concurrent population expansions in both species during the large-scale clearing of eastern forests (∼200 years ago). Finally, we find no evidence that either species experienced a reduction in effective population size over the past 75 years. Thus, the well-documented decline of monarch abundance over the past 40 years is not visible in our genomic dataset, reflecting a possible mismatch of the overwintering census population to effective population size in this species.

A multidimensional framework to quantify the effects of urbanization on avian breeding fitness

Urbanization has dramatically altered Earth's landscapes and changed a multitude of environmental factors. This has resulted in intense land-use change, and adverse consequences such as the urban heat island effect (UHI), noise pollution, and artificial light at night (ALAN). However, there is a lack of research on the combined effects of these environmental factors on life-history traits and fitness, and on how these interactions shape food resources and drive patterns of species persistence. Here, we systematically reviewed the literature and created a comprehensive framework of the mechanistic pathways by which urbanization affects fitness and thus favors certain species. We found that urbanization-induced changes in urban vegetation, habitat quality, spring temperature, resource availability, acoustic environment, nighttime light, and species behaviors (e.g., laying, foraging, and communicating) influence breeding choices, optimal time windows that reduce phenological mismatch, and breeding success. Insectivorous and omnivorous species that are especially sensitive to temperature often experience advanced laying behaviors and smaller clutch sizes in urban areas. By contrast, some granivorous and omnivorous species experience little difference in clutch size and number of fledglings because urban areas make it easier to access anthropogenic food resources and to avoid predation. Furthermore, the interactive effect of land-use change and UHI on species could be synergistic in locations where habitat loss and fragmentation are greatest and when extreme-hot weather events take place in urban areas. However, in some instances, UHI may mitigate the impact of land-use changes at local scales and provide suitable breeding conditions by shifting the environment to be more favorable for species' thermal limits and by extending the time window in which food resources are available in urban areas. As a result, we determined five broad directions for further research to highlight that urbanization provides a great opportunity to study environmental filtering processes and population dynamics.

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.

Genetic evidence for widespread population size expansion in North American boreal birds prior to the Last Glacial Maximum

Pleistocene climate cycles are well documented to have shaped contemporary species distributions and genetic diversity. Northward range expansions in response to deglaciation following the Last Glacial Maximum (LGM; approximately 21 000 years ago) are surmised to have led to population size expansions in terrestrial taxa and changes in seasonal migratory behaviour. Recent findings, however, suggest that some northern temperate populations may have been more stable than expected through the LGM. We modelled the demographic history of 19 co-distributed boreal-breeding North American bird species from full mitochondrial gene sets and species-specific molecular rates. We used these demographic reconstructions to test how species with different migratory strategies were affected by glacial cycles. Our results suggest that effective population sizes increased in response to Pleistocene deglaciation earlier than the LGM, whereas genetic diversity was maintained throughout the LGM despite shifts in geographical range. We conclude that glacial cycles prior to the LGM have most strongly shaped contemporary genetic diversity in these species. We did not find a relationship between historic population dynamics and migratory strategy, contributing to growing evidence that major switches in migratory strategy during the LGM are unnecessary to explain contemporary migratory patterns.

A field theoretic approach to non-equilibrium population genetics in the strong selection regime

Natural populations are virtually never observed in equilibrium, yet equilibrium approximations comprise the majority of our understanding of population genetics. Using standard tools from statistical physics, a formalism is presented that re-expresses the stochastic equations describing allelic evolution as a partition functional over all possible allelic trajectories ('paths') governed by selection, mutation, and drift. A perturbative field theory is developed for strong additive selection, relevant to disease variation, that facilitates the straightforward computation of closed-form approximations for time-dependent moments of the allele frequency distribution across a wide range of non-equilibrium scenarios; examples are presented for constant population size, exponential growth, bottlenecks, and oscillatory size, all of which align well to simulations and break down just above the drift barrier. Equilibration times are computed and, even for static population size, generically extend beyond the order 1/s timescale associated with exponential frequency decay. Though the mutation load is largely robust to variable population size, perturbative drift-based corrections to the deterministic trajectory are readily computed. Under strong selection, the variance of a new mutation's frequency (related to homozygosity) is dominated by drift-driven dynamics and a transient increase in variance often occurs prior to equilibrating. The excess kurtosis over skew squared is roughly constant (i.e., independent of selection, provided 2Ns ≳ 5) for static population size, and thus potentially sensitive to deviation from equilibrium. These insights highlight the value of such closed-form approximations, naturally generated from Feynman diagrams in a perturbative field theory, which can simply and accurately capture the parameter dependences describing a variety of non-equilibrium population genetic phenomena of interest.

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.

Population structure, genomic diversity and demographic history of Komodo dragons inferred from whole-genome sequencing

Population and conservation genetics studies have greatly benefited from the development of new techniques and bioinformatic tools associated with next-generation sequencing. Analysis of extensive data sets from whole-genome sequencing of even a few individuals allows the detection of patterns of fine-scale population structure and detailed reconstruction of demographic dynamics through time. In this study, we investigated the population structure, genomic diversity and demographic history of the Komodo dragon (Varanus komodoensis), the world's largest lizard, by sequencing the whole genomes of 24 individuals from the five main Indonesian islands comprising the entire range of the species. Three main genomic groups were observed. The populations of the Island of Komodo and the northern coast of Flores, in particular, were identified as two distinct conservation units. Degrees of genomic divergence among island populations were interpreted as a result of changes in sea level affecting connectivity across islands. Demographic inference suggested that Komodo dragons probably experienced a relatively steep population decline over the last million years, reaching a relatively stable N_e during the Saalian glacial cycle (400-150 thousand years ago) followed by a rapid N_e decrease. Genomic diversity of Komodo dragons was similar to that found in endangered or already extinct reptile species. Overall, this study provides an example of how whole-genome analysis of a few individuals per population can help define population structure and intraspecific demographic dynamics. This is particularly important when applying population genomics data to conservation of rare or elusive endangered species.

Effects of artificial light at night on the nest-site selection, reproductive success and behavior of a synanthropic bird

Understanding how artificial light at night (ALAN) impacts wildlife is increasingly important because more and more species are colonizing urban areas. As most of the bird studies on ALAN use controlled light set inside or around nest-boxes, the ecological effect of ALAN resulting from in situ streetlight on birds remains contentious. The barn swallow (Hirundo rustica) often builds open nests on buildings, which are directly exposed to varying intensity of ALAN, and thus provides a good system to examine the effect of in situ ALAN on birds. By examining the nest-site selection, reproductive success and behavior of barn swallows under various ALAN intensity in Taipei City, we found a positive effect of ALAN on their fledging success; nonetheless, such effect was only found in the swallows' first brood, but not second one. We also found that parent birds in the nests with higher ALAN intensity had higher feeding rates and more extended feeding time past sunset, which were likely stimulated by the increased begging behavior of their chicks. The night-feeding behavior might contribute to the increased fledging success, especially at the early breeding season. Interestingly, despite of the reproductive benefits obtained from ALAN, we found that the barn swallows did not select nest sites regarding ALAN intensity. The weak nest-site selection perhaps result from the complex life history interactions involving ALAN and/or confounding factors associated with ALAN in cities. This study improves our understanding of how urban birds, especially open-nesting ones, respond to in situ ALAN and provides useful information for developing urban conservation strategies.

Nesting strategy shapes territorial aggression but not testosterone: A comparative approach in female and male birds

Our understanding of the proximate and ultimate mechanisms shaping competitive reproductive phenotypes primarily stems from research on male-male competition for mates, even though competition is widespread in both sexes. We evaluate the hypothesis that the restricted nature of a resource required for reproduction, i.e. nest site, is a key variable driving territorial competition and testosterone secretion in female and male birds. Obligate secondary cavity-nesting has evolved repeatedly across avian lineages, providing a useful comparative context to explore how competition over limited nest cavities shapes aggression and its underlying mechanisms across species. Although evidence from one or another cavity-nesting species suggests that territorial aggression is adaptive in both females and males, this has not yet been tested in a comparative framework. We predicted that cavity-nesting generates more robust territorial aggression, in comparison to close relatives with less restrictive nesting strategies. Our focal species were two obligate secondary cavity-nesting species and two related species with more flexible nesting strategies in the same avian family: tree swallow (Tachycineta bicolor) vs. barn swallow (Hirundo rustica); Eastern bluebird (Sialia sialis) vs. American robin (Turdus migratorius). We assayed conspecific aggression using simulated territorial intrusion and found that cavity-nesting species displayed greater territorial aggression than their close relatives. This pattern held for both females and males. Because territorial aggression is often associated with elevated testosterone, we also hypothesized that cavity-nesting species would exhibit higher testosterone levels in circulation. However, cavity-nesting species did not have higher testosterone in circulation for either sex, despite some correlative evidence that testosterone is associated with higher rates of physical attack in female tree swallows. Our focus on a context that is relevant to both sexes - competition over essential breeding resources - provides a useful framework for co-consideration of proximate and ultimate drivers of reproductive competition in females and males.

Whole-genome re-sequencing data to infer historical demography and speciation processes in land snails: the study of two Candidula sister species

Despite the global biodiversity of terrestrial gastropods and their ecological and economic importance, the genomic basis of ecological adaptation and speciation in land snail taxa is still largely unknown. Here, we combined whole-genome re-sequencing with population genomics to evaluate the historical demography and the speciation process of two closely related species of land snails from western Europe, Candidula unifasciata and C. rugosiuscula. Historical demographic analysis indicated fluctuations in the size of ancestral populations, probably driven by Pleistocene climatic fluctuations. Although the current population distributions of both species do not overlap, our approximate Bayesian computation model selection approach on several speciation scenarios suggested that gene flow has occurred throughout the divergence process until recently. Positively selected genes diverging early in the process were associated with intragenomic and cyto-nuclear incompatibilities, respectively, potentially fostering reproductive isolation as well as ecological divergence. Our results suggested that the speciation between species entails complex processes involving both gene flow and ecological speciation, and that further research based on whole-genome data can provide valuable understanding on species divergence. This article is part of the Theo Murphy meeting issue 'Molluscan genomics: broad insights and future directions for a neglected phylum'.

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.

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.

Analysing phenotypic variation in barn swallows (Hirundo rustica) across China to assess subspecies status

Phenotypic variation is often used to delineate species and subspecies boundaries. Statistical analysis of phenotypic variation within a species is helpful both for understanding biodiversity and for its conservation. At least three named subspecies of barn swallows (Hirundo rustica) have distributions that span China, yet, to date, no systematic study of phenotypic differentiation has been applied to understand the delineation of these subspecies. In this study, we collected 510 samples of barn swallow from 23 populations in China, including two recognized subspecies, H. r. rustica and H. r. gutturalis, and one potential subspecies, H. r. mandschurica. With these samples, we examined and found morphometric and colour differences among different populations. Western Chinese barn swallows (H. r. rustica) have larger body size and could be clearly differentiated from eastern H. r. gutturalis and H. r. mandschurica, while north-eastern populations (named H. r. mandschurica) have darker, redder ventral plumage than H. r. rustica and H. r. gutturalis. However, we inferred that although there were phenotypic differences between H. r. mandschurica and H. r. gutturalis, they were not sufficiently distinct to assign them to separate subspecies based on the 75% rule for defining subspecies.

Phylogenomic Data Reveal Widespread Introgression Across the Range of an Alpine and Arctic Specialist

Understanding how gene flow affects population divergence and speciation remains challenging. Differentiating one evolutionary process from another can be difficult because multiple processes can produce similar patterns, and more than one process can occur simultaneously. Although simple population models produce predictable results, how these processes balance in taxa with patchy distributions and complicated natural histories is less certain. These types of populations might be highly connected through migration (gene flow), but can experience stronger effects of genetic drift and inbreeding, or localized selection. Although different signals can be difficult to separate, the application of high-throughput sequence data can provide the resolution necessary to distinguish many of these processes. We present whole-genome sequence data for an avian species group with an alpine and arctic tundra distribution to examine the role that different population genetic processes have played in their evolutionary history. Rosy-finches inhabit high elevation mountaintop sky islands and high-latitude island and continental tundra. They exhibit extensive plumage variation coupled with low levels of genetic variation. Additionally, the number of species within the complex is debated, making them excellent for studying the forces involved in the process of diversification, as well as an important species group in which to investigate species boundaries. Total genomic variation suggests a broadly continuous pattern of allele frequency changes across the mainland taxa of this group in North America. However, phylogenomic analyses recover multiple distinct, well supported, groups that coincide with previously described morphological variation and current species-level taxonomy. Tests of introgression using D-statistics and approximate Bayesian computation reveal significant levels of introgression between multiple North American taxa. These results provide insight into the balance between divergent and homogenizing population genetic processes and highlight remaining challenges in interpreting conflict between different types of analytical approaches with whole-genome sequence data. [ABBA-BABA; approximate Bayesian computation; gene flow; phylogenomics; speciation; whole-genome sequencing.].

Population history of the golden eagle inferred from whole-genome sequencing of three of its subspecies

The application of evolutionary genetic research to investigate the potential for endangered species to adapt to changing environments is important for conservation biology. Effective population size (Ne) is informative for understanding adaptive potential as it refers to the genetic variation in breeding individuals who have contributed to contemporary and historic population diversity. We reconstruct fluctuations in Ne in three golden eagle subspecies (Japanese, Scottish, North American) using the pairwise sequential Markovian coalescent (PSMC) model based on whole-genome sequence data. Our results indicate the timing of subspeciation events and suggest significant ongoing demographic reductions since the start of the Last Glacial Period. Importantly, we find evidence for gene flow from continental populations into the ancestral Japanese population resulting in a short, sharp recovery in genetic diversity. Timing agrees with the palaeogeographic estimates of land bridge connections between the Japanese archipelago and Asian continent and matches a similar Ne spike in the Scottish population, but not in the North American population. Given contemporary declines in isolated Japanese and UK island populations, our study highlights a concerning loss of local genetic diversity, but also indicates the likely response of populations to genetic reinforcement from neighbouring subspecies, increasing management options and encouraging a range-wide species conservation approach.

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.

Lack of fine-tuned egg rejection adjustment in barn swallows with variable local abundance of common cuckoos

Tracking variation in hosts' responses to parasitism across space in a metapopulation is critical to assess the current status of parasitism/anti-parasitism in a host-parasite system, which is also helpful to infer its coevolutionary history. The barn swallow Hirundo rustica is a tractable bird species to understand potential fine-tuned adaptations to avian brood parasitism across small geographic scales, both in the context of variation in its foreign-egg rejection pattern, and its widespread distribution in cuckoo-free urban areas and in cuckoo-present rural habitats, including reedbeds. Here we tested whether variation in co-occurrence between the host and its rare brood parasite, the common cuckoo Cuculus canorus (i.e., high level of co-occurrence in reed habitats, low in town habitats) at the metapopulation level predicts patterns of antiparasitic egg rejection behaviors in barn swallows in response to different types of model and natural eggs. Contrary to our predictions, higher parasite detectability in the reed habitat did not translate into higher parasitism rate and, in turn, we also found similarly low egg rejection rates across both sampled habitat types. These patterns implied a lack of fine-tuned increase of egg rejection rate in the reed-breeding population of barn swallows as a response to the increased perceived cuckoo encounter rate, perhaps because higher potential parasitism threat did not transfer into greater actual parasitism rate. It remains to be assessed whether the lack of small-scale geographic variation in barn swallows' egg rejection rates persists because this species responds to selection by parasitism as a spatially contiguous evolutionary unit.