The effects of cultural drift on geographic variation in echolocation calls of the Chinese rufous horseshoe bat (Rhinolophus sinicus)

Drift, selection, or their combined effects commonly drive geographic variation in traits. Clarifying the relative roles of each process is a long‐standing research goal in evolutionary biology. Acoustic signals of bats are a phenotypic characteristic that plays an important role in social organization and species recognition. We extensively sampled the Chinese rufous horseshoe bat (Rhinolophus sinicus) throughout China and Vietnam and reconstructed a species phylogeny to better understand the patterns and causes of the geographic variation of acoustic signals. Our results showed that the resting frequency (RF) of calls varied with latitude, sex, and distance among different colony locations. RF differences were not correlated with genetic distance (based on only one nuclear locus and the mitochondrial locus), climatic factors (mean annual temperature and mean annual precipitable water), or body size, although differences in calls increased with distance among various populations. This suggests that cultural drift may play more important roles than genetic drift and acoustic adaptation in shaping acoustic differences within regions in R. sinicus.     

Testing drivers of acoustic divergence in cicadas (Cicadidae: Tettigettalna)

Divergence in acoustic signals may have a crucial role in the speciation process of animals that rely on sound for intra-specific recognition and mate attraction. The acoustic adaptation hypothesis (AAH) postulates that signals should diverge according to the physical properties of the signalling environment. To be efficient, signals should maximize transmission and decrease degradation. To test which drivers of divergence exert the most influence in a speciose group of insects, we used a phylogenetic approach to the evolution of acoustic signals in the cicada genus Tettigettalna, investigating the relationship between acoustic traits (and their mode of evolution) and body size, climate and micro-/macro-habitat usage. Different traits showed different evolutionary paths. While acoustic divergence was generally independent of phylogenetic history, some temporal variables’ divergence was associated with genetic drift. We found support for ecological adaptation at the temporal but not the spectral level. Temporal patterns are correlated with micro- and macro-habitat usage and temperature stochasticity in ways that run against the AAH predictions, degrading signals more easily. These traits are likely to have evolved as an anti-predator strategy in conspicuous environments and low-density populations. Our results support a role of ecological selection, not excluding a likely role of sexual selection in the evolution of Tettigettalna calling songs, which should be further investigated in an integrative approach.     

Weather conditions determine attenuation and speed of sound: Environmental limitations for monitoring and analyzing bat echolocation

Echolocating bats are regularly studied to investigate auditory‐guided behaviors and as important bioindicators. Bioacoustic monitoring methods based on echolocation calls are increasingly used for risk assessment and to ultimately inform conservation strategies for bats. As echolocation calls transmit through the air at the speed of sound, they undergo changes due to atmospheric and geometric attenuation. Both the speed of sound and atmospheric attenuation, however, are variable and determined by weather conditions, particularly temperature and relative humidity. Changing weather conditions thus cause variation in analyzed call parameters, limiting our ability to detect, and correctly analyze bat calls. Here, I use real‐world weather data to exemplify the effect of varying weather conditions on the acoustic properties of air. I then present atmospheric attenuation and speed of sound for the global range of weather conditions and bat call frequencies to show their relative effects. Atmospheric attenuation is a nonlinear function of call frequency, temperature, relative humidity, and atmospheric pressure. While atmospheric attenuation is strongly positively correlated with call frequency, it is also significantly influenced by temperature and relative humidity in a complex nonlinear fashion. Variable weather conditions thus result in variable and unknown effects on the recorded call, affecting estimates of call frequency and intensity, particularly for high frequencies. Weather‐induced variation in speed of sound reaches up to about ±3%, but is generally much smaller and only relevant for acoustic localization methods of bats. The frequency‐ and weather‐dependent variation in atmospheric attenuation has a threefold effect on bioacoustic monitoring of bats: It limits our capability (1) to monitor bats equally across time, space, and species, (2) to correctly measure frequency parameters of bat echolocation calls, particularly for high frequencies, and (3) to correctly identify bat species in species‐rich assemblies or for sympatric species with similar call designs.     

Environmental implications of pollen spectra in bat droppings from southeastern Spain and potential for palaeoenvironmental reconstructions

Pollen was analysed from bat guano from nine caves in southeastern Spain and surface soils in their immediate surroundings. We compare the pollen spectra of 34 modern dung samples from the nine caves with one modern surface pollen sample from each cave. The contents suggest reasonable pollen diversity and richness, including anemophilous and zoophilous pollen types. Since the latter is usually under-represented in atmospheric pollen, the guano spectra therefore appear to reflect the vegetation more effectively than normal surface soil samples. Despite health hazards such as histoplasmosis, the difficulties of obtaining bat guano in deep caves and possible interpretational concerns relating to behaviour and feeding habits of different bat species, this material can be very useful in palaeoecological research provided that the dung was fossilized under favourable environmental conditions that allowed the preservation of pollen.     

Sex‐specific phenotypic selection and geographic variation in gender divergence in a gynodioecious shrub

In sexually polymorphic plant species the extent of gender divergence in floral morphology and phenology may be influenced by gender‐specific selection patterns imposed by pollinators, which may change geographically. Distribution margins are areas where changes in the pollinator fauna, and thus variation in gender divergence of floral traits, are expected. We tested for pollination‐driven geographic variation in the gender divergence in floral and phenological traits in the gynodioecious shrub Daphne laureola, in core and marginal areas differing in the identity of the main pollinator. Pollinators selected for longer corolla tubes in hermaphrodite individuals only in core populations, which in turn recorded higher fruit set. Consistent with these phenotypic selection patterns, gender divergence in flower corolla length was higher in core populations. Moreover, pollinators selected towards delayed flowering on hermaphrodite individuals only in marginal populations, where the two sexes differed more in flowering time. Our results support that a shift in main pollinators is able to contribute to geographic variation in the gender divergence of sexually polymorphic plant species.     

Ecological characterization of bat species distributions in Michoacán, México, using a geographic information system

Aim To investigate the relationship between chiropteran species distributions and four fundamental environmental factors — temperature, precipitation, soil and vegetation — and to construct a species richness prediction map based on the environmental factors. Location The state of Michoacán, México. Methods Bat specimens collected during a 2‐year project were combined with museum specimens to form a dataset of over 3200 entries pertaining to 71 species of bats. Coordinates of the collection localities were recorded with GPS receivers or determined from maps. ArcView GIS was used to characterize the distribution of the species relative to the four environmental factors by projecting coordinates of the collection sites onto digitized maps of those factors. Correspondence analysis (CA) was used to evaluate the relationship between species distributions and the environmental factors. Results The CA results indicated that the order of importance of these factors is (from highest to lowest): temperature, vegetation, precipitation and soil. A predicted distribution map was constructed for each species of bat, based on the result of the CA analysis, using correspondences of each species to climate, vegetation and precipitation. Soil types were excluded from the prediction model because soil type does not appear to carry high predictive value for bat species in Michoacán. Distribution maps of the 71 bat species were then overlaid to generate a map of bat species richness for the state of Michoacán. Main conclusions Neither family membership nor feeding guild affiliation appear to play important roles in chiropteran species distributions in Michoacán. The bat species richness prediction map will be a useful tool for conservation works in the region.     

Drivers of longitudinal telomere dynamics in a long‐lived bat species,Myotis myotis

Age‐related telomere shortening is considered a hallmark of the ageing process. However, a recent cross‐sectional ageing study of relative telomere length (rTL) in bats failed to detect a relationship between rTL and age in the long‐lived genus Myotis (M. myotis and M. bechsteinii), suggesting some other factors are responsible for driving telomere dynamics in these species. Here, we test if longitudinal rTL data show signatures of age‐associated telomere attrition in M. myotis and differentiate which intrinsic or extrinsic factors are likely to drive telomere length dynamics. Using quantitative polymerase chain reaction, rTL was measured in 504 samples from a marked population, from Brittany, France, captured between 2013 and 2016. These represent 174 individuals with an age range of 0 to 7+ years. We find no significant relationship between rTL and age (p = .762), but demonstrate that within‐individual rTL is highly variable from year to year. To investigate the heritability of rTL, a population pedigree (n = 1744) was constructed from genotype data generated from a 16‐microsatellite multiplex, designed from an initial, low‐coverage, Illumina genome for M. myotis. Heritability was estimated in a Bayesian, mixed model framework, and showed that little of the observed variance in rTL is heritable (h² = 0.01–0.06). Rather, correlations of first differences, correlating yearly changes in telomere length and weather variables, demonstrate that, during the spring transition, average temperature, minimum temperature, rainfall and windspeed correlate with changes in longitudinal telomere dynamics. As such, rTL may represent a useful biomarker to quantify the physiological impact of various environmental stressors in bats.     

Divergent echolocation call frequencies in insular rhinolophids (Chiroptera): a case of character displacement?

Aims Because rhinolophids have been hypothesized to use echolocation call frequency to recognize conspecifics, sympatric species calling at similar frequencies should be subject to acoustic character displacement, i.e. a drift in frequency values to minimize the risk of misidentification of conspecifics. However, it has been proposed that insufficient geographical separation between populations in sympatry and allopatry may counter the establishment of frequency differences by character displacement. Here we tested the hypothesis that insular populations should exhibit acoustic divergence, and this should be revealed by comparing call frequencies with those observed in mainland, allopatric populations of conspecifics. We also tested whether the evolutionary pressure towards acoustic divergence should be especially strong at sites where rhinolophid species emitting similar call frequencies roost together in order to minimize interspecific frequency overlap. Location Sardinia and southern Italy (Campania, Lazio, Abruzzo). Methods Time‐expanded echolocation calls and body size were recorded from Sardinian populations of Rhinolophus mehelyi Matschie, Rhinolophus hipposideros (Bechstein) and Rhinolophus euryale Blasius. Both call frequencies and forearm length of insular R. hipposideros and R. euryale were compared with those of populations from mainland areas of Italy where R. mehelyi is absent, to explore the hypothesis that the presence of the latter species (which calls at frequency values intermediate between the other two) may determine acoustic divergence in the other species. For Sardinian R. mehelyi and R. euryale, we also carried out intraspecific comparisons of call frequencies between bats from monospecific colonies and those from mixed colonies. Results As hypothesized, Sardinian R. hipposideros and R. euryale called at frequencies higher and lower, respectively, than in the peninsula. In this way, overlap with R. mehelyi is avoided. Body size showed no difference between insular and peninsular populations, i.e. frequency differences are not a by‐product of difference in body size determined by insularity. Frequency values in Sardinian R. euryale from monospecific colonies did not differ from those of bats roosting together with R. mehelyi. However, R. mehelyi showed frequency values significantly higher when associated with R. euryale, possibly to minimize the risk of species misrecognition. Main conclusions At least under geographical isolation, character displacement may be a causal mechanism for shifts in call frequency of sympatric rhinolophids. Species recognition and facilitation of intraspecific communication (with possible implications for mate recognition) constitute the best candidate factors for the phenomenon we observed.     

The maintenance of phenotypic divergence through sexual selection: An experimental study in barn swallows Hirundo rustica

Previous studies have shown that sexual signals can rapidly diverge among closely related species. However, we lack experimental studies to demonstrate that differences in trait‐associated reproductive performance maintain sexual trait differences between closely related populations, in support for a role of sexual selection in speciation. Populations of Northern Hemisphere distributed barn swallows Hirundo rustica are closely related, yet differ in two plumage‐based traits: ventral color and length of the outermost tail feathers (streamers). Here we provide experimental evidence that manipulations of these traits result in different reproductive consequences in two subspecies of barn swallow: (H. r. erythrogaster in North America and H. r. transitiva in the East Mediterranean). Experimental results in Colorado, USA, demonstrate that males with (1) darkened ventral coloration and (2) shortened streamers gained paternity between two successive reproductive bouts. In contrast, exaggeration of both traits improved reproductive performance within H. r. transitiva in Israel: males with a combination treatment of darkened ventral coloration and elongated streamers gained paternity between two successive reproductive bouts. Collectively, these experimental results fill an important gap in our understanding for how divergent sexual selection maintains phenotype differentiation in closely related populations, an important aspect of the speciation process.     

The relative contribution of drift and selection to phenotypic divergence: A test case using the horseshoe bats Rhinolophus simulator and Rhinolophus swinnyi

Natural selection and drift can act on populations individually, simultaneously or in tandem and our understanding of phenotypic divergence depends on our ability to recognize the contribution of each. According to the quantitative theory of evolution, if an organism has diversified through neutral evolutionary processes (mutation and drift), variation of phenotypic characteristics between different geographic localities (B) should be directly proportional to the variation within localities (W), that is, B ∝ W. Significant deviations from this null model imply that non‐neutral forces such as natural selection are acting on a phenotype. We investigated the relative contributions of drift and selection to intraspecific diversity using southern African horseshoe bats as a test case. We characterized phenotypic diversity across the distributional range of Rhinolophus simulator (n = 101) and Rhinolophus swinnyi (n = 125) using several traits associated with flight and echolocation. Our results suggest that geographic variation in both species was predominantly caused by disruptive natural selection (B was not directly proportional to W). Evidence for correlated selection (co‐selection) among traits further confirmed that our results were not compatible with drift. Selection rather than drift is likely the predominant evolutionary process shaping intraspecific variation in traits that strongly impact fitness.     

Measuring genome divergence in bacteria: a case study using chlamydian data

We have studied the relative contribution of inversions, transpositions, deletions, and nucleotide substitutions to the evolution of Chlamydia trachomatis and Chlamydia pneumoniae. The minimal number of rearrangement events required for converting the gene order structure of one genome into that of the other was estimated to 59 +/- 6 events, including 13% inversions, 38% short inversions, and 49% transpositions. In contrast to previous findings, no examples of horizontal gene transfer subsequent to species divergence were identified, nor any evidence for an excessive number of tandem gene duplications. A statistical model was used to compare nucleotide frequencies for a set of genes uniquely present in one species to a set of orthologous genes present in both species. The two data sets were not significantly different, which is indicative of a low frequency of horizontal gene transfer events. This is based on the assumption that a foreign gene of different nucleotide content will not have become completely ameliorated, as verified by simulations of the amelioration rate at twofold and fourfold degenerate codon sites. The frequencies of nucleotide substitutions at twofold and fourfold degenerate sites, deletions, inversions, and translocations were estimated to 1.42, 0.62, 0.18, 0.01, and 0.01 per site, respectively.     

Genome-wide patterns of divergence during speciation: the lake whitefish case study

The nature, size and distribution of the genomic regions underlying divergence and promoting reproductive isolation remain largely unknown. Here, we summarize ongoing efforts using young (12 000 yr BP) species pairs of lake whitefish (Coregonus clupeaformis) to expand our understanding of the initial genomic patterns of divergence observed during speciation. Our results confirmed the predictions that: (i) on average, phenotypic quantitative trait loci (pQTL) show higher F(ST) values and are more likely to be outliers (and therefore candidates for being targets of divergent selection) than non-pQTL markers; (ii) large islands of divergence rather than small independent regions under selection characterize the early stages of adaptive divergence of lake whitefish; and (iii) there is a general trend towards an increase in terms of numbers and size of genomic regions of divergence from the least (East L.) to the most differentiated species pair (Cliff L.). This is consistent with previous estimates of reproductive isolation between these species pairs being driven by the same selective forces responsible for environment specialization. Altogether, dwarf and normal whitefish species pairs represent a continuum of both morphological and genomic differentiation contributing to ecological speciation. Admittedly, much progress is still required to more finely map and circumscribe genomic islands of speciation. This will be achieved through the use of next generation sequencing data but also through a better quantification of phenotypic traits moulded by selection as organisms adapt to new environmental conditions.     

河南济源人工引水渠隧道蝙蝠冬眠生态学特征

2017-2020年期间,每年1月份对河南省济源市邵原镇布袋沟水库人工引水渠隧道内蝙蝠进行冬眠生态学特征调查,共发现2科5属7种蝙蝠在此冬眠,包括马铁菊头蝠（Rhinolophus ferrumequinum）、小菊头蝠（R.pusillus）、华南水鼠耳蝠（Myotis laniger）、白腹管鼻蝠（Murina leucogaster）、金管鼻蝠（Mu.aurata）、奥氏长耳蝠（Plecotus ognevi）和亚洲宽耳蝠（Barbastella leucomelas）。马铁菊头蝠是优势种（约52%-73%的冬眠个体）,其次是小菊头蝠（约19%-37%）、华南水鼠耳蝠（约5%-8%）,其余蝙蝠物种数量不足3%。2017-2020年冬眠蝙蝠个体总数呈增长趋势,但仍少于早期报道的数量。有42个隧道每年均有蝙蝠冬眠,而且不同年度冬眠数量也不尽相同。通过多元线性回归分析发现,隧道长度可能是影响蝙蝠冬眠栖息场所选择的主要影响因子（Adjusted R²=0.208,P=0.001）。每个隧道内,蝙蝠具有不同的冬眠栖点位置,约4/5的蝙蝠选择温暖且环境相对稳定的隧道深处（> 30 m）作为冬眠栖点,超过95%的个体选择长度> 60 m的隧道冬眠。蝙蝠具有不同的冬眠方式,绝大多数个体采用独栖方式进行冬眠（> 90%）,少数采用聚集方式。不同的冬眠栖点和冬眠方式可能有利于冬眠成本优化。此外,栖点温度与蝙蝠体温之间呈显著正相关（R²=0.98,P < 0.001）,而且蝙蝠冬眠期间的栖点温度具有种内和种间差异。研究结果为我国蝙蝠种群保护和冬眠场所管理提供科学依据。     

Environmental correlates of body size influence range size and extinction risk: A global study in rodents

Aim

Species-level traits, such as body and range sizes, are important correlates of extinction risk. However, both are often related and are driven by environmental factors. Here, we elucidated links between environmental factors, body size, range size and susceptibility to extinction, across the whole order of rodents.

Location

Global.

Time period

Current.

Major taxa studied

Rodents (order Rodentia).

Methods

We compiled an unprecedentedly large database of rodent morphology, phylogeny, range size, conservation status, global climate and normalized difference vegetation index (NDVI), comprising >86% of all described species. Using phylogenetic regressions, we initially explored the environmental factors driving body size. Next, we modelled the relationship between body size and range size. From this relationship, we computed and mapped (at the assemblage level) an index of relative range size, corresponding to the deviation from the expected range size of each species, given its body size. Finally, we tested whether relative range was correlated with the risk of extinction of the species derived from an assessment by the International Union for Conservation of Nature.

Results

We found that, contrary to the expectations of Bergmann's rule, the body size of rodents was mostly influenced by variation in NDVI (rather than latitude/temperature). Body size, in turn, imposed a constraint on species range size, as evidenced by a triangular relationship that was segmented at the lower bound. The relative species range size derived from this relationship highlighted four geographical regions where rodents with small relative range were concentrated globally. We demonstrated that lower relative range size was associated with increased risk of extinction.

Main conclusions

Species that, given their body size, are distributed across ranges that are smaller than expected have elevated extinction risk. Therefore, investigating the relationships between environmental drivers, body size and range size might help to detect species that could become threatened in the near future.     

Heard but not seen: Comparing bat assemblages and study methods in a mosaic landscape in the Western Ghats of India

We used capture (mist‐netting) and acoustic methods to compare the species richness, abundance, and composition of a bat assemblage in different habitats in the Western Ghats of India. In the tropics, catching bats has been more commonly used as a survey method than acoustic recordings. In our study, acoustic methods based on recording echolocation calls detected greater bat activity and more species than mist‐netting. However, some species were detected more frequently or exclusively by capture. Ideally, the two methods should be used together to compensate for the biases in each. Using combined capture and acoustic data, we found that protected forests, forest fragments, and shade coffee plantations hosted similar and diverse species assemblages, although some species were recorded more frequently in protected forests. Tea plantations contained very few species from the overall bat assemblage. In riparian habitats, a strip of forested habitat on the river bank improved the habitat for bats compared to rivers with tea planted up to each bank. Our results show that shade coffee plantations are better bat habitat than tea plantations in biodiversity hotspots. However, if tea is to be the dominant land use, forest fragments and riparian corridors can improve the landscape considerably for bats. We encourage coffee growers to retain traditional plantations with mature native trees, rather than reverting to sun grown coffee or coffee shaded by a few species of timber trees.     

澳门翼手类物种多样性调查

2009 ～2012 年,对澳门翼手目(蝙蝠)物种多样性进行了调查。结果共捕捉到10 个物种,属5 科8 属,其中包括澳门原来记载的2 个物种,即蹄蝠科的大蹄蝠(Hipposideros armiger)和蝙蝠科的东亚伏翼(Pipistrellus abramus);本研究新增加8 个物种,即狐蝠科的犬蝠(Cynopterus sphinx) 和棕果蝠(Rousettus leschenaulti),鞘尾蝠科的黑髯墓蝠(Taphozous melanopogon),菊头蝠科的菲菊头蝠(Rhinolophus pusillus),以及蝙蝠科的大足鼠耳蝠(Myotis ricketti)、普通伏翼(P. pipistrellus)、普通长翼蝠(Miniopterus schreibersi) 和南长翼蝠(M. pusillus)。另外,通过野外录音和分析,并与已发表物种声音特征比较核对,发现菊头蝠科和蹄蝠科各一种,前者可能是泰国菊头蝠(R. siamensis)或者中菊头蝠(R. affinis),后者可能是果树蹄蝠(H. pomona) 或者三叶蹄蝠(Aselliscus stoliczkanus)。本文对已捕捉10 种蝙蝠的分布、形态特征和回声定位叫声特征进行报道,同时对其种群数量和保护现状进行了讨论。保护蝙蝠栖息生境(洞穴、古老建筑和蒲葵树等) 对保护澳门蝙蝠物种多样性至关重要。     

Qualitative population divergence in proximate determination of a sexually selected trait in the collared flycatcher

We examined proximate determination of sexually selected forehead patch size in a Central‐European population of Ficedula albicollis, the collared flycatcher, using a 9‐year database, and compared our results with those obtained in other populations of the same and the sister species. Between‐individual variation of forehead patch size was large, its repeatability larger than, and heritability similar to the Swedish population. Unlike in the other populations, the trait proved unaffected by body condition, and only very slightly influenced by age. There was no relationship between forehead patch size and breeding lifespan, and a marginal negative association with survivorship in adult males. Our results suggest that additive genetic variance of the trait in this population is large, but genes act independently of body condition, and there is no viability indicator value of the trait. This is the first report of a qualitative intraspecific difference in proximate determination of a sexually selected trait.     

Leader preference in Neoconocephalus ensiger katydids: a female preference for a nonheritable male trait

Female preferences for males producing their calls just ahead of their neighbours, leader preferences, are common in acoustically communicating insects and anurans. While these preferences have been well studied, their evolutionary origins remain unclear. We tested whether females gain a fitness benefit by mating with leading males in Neoconocephalus ensiger katydids. We mated leading and following males with random females and measured the number and quality of F₁, the number of F₂ and the heritability of the preferred male trait. We found that females mating with leaders and followers did not differ in the number of F₁ or F₂ offspring. Females mating with leading males had offspring that were in better condition than those mating with following males suggesting a benefit in the form of higher quality offspring. We found no evidence that the male trait, the production of leading calls, was heritable. This suggests that there is no genetic correlate for the production of leading calls and that the fitness benefit gained by females must be a direct benefit, potentially mediated by seminal proteins. The presence of benefits indicates that leader preference is adaptive in N. ensiger, which may explain the evolutionary origin of leader preference; further tests are required to determine whether fitness benefits can explain the phylogenetic distribution of leader preference in Neoconocephalus. The absence of heritability will prevent leader preference from becoming coupled with or exaggerating the male trait and prevent females from gaining a ‘sexy‐sons’ benefit, weakening the overall selection for leader preference.     

Biogeographical basis of recent phenotypic divergence among birds: a global study of subspecies richness

Theory predicts that biogeographic factors should play a central role in promoting population divergence and speciation. Previous empirical studies into biogeography and diversification have been relatively restricted in terms of the geographical area, phylogenetic scope, and the range of biogeographic factors considered. Here we present a global analysis of allopatric phenotypic divergence (measured as subspecies richness) across more than 9600 bird species. The main aim of this study was to examine the extent to which biogeographical factors can explain patterns of phenotypic divergence. Analysis of the taxonomic distribution of subspecies among species suggests that subspecies formation and extinction have occurred at a considerably faster rate than has species formation. However, the observed distribution departs from the expectation under a random birth-death model of diversification. Across 19 phylogenetic trees, we find no significant linear relationship between species age and subspecies richness, implying that species age is a poor predictor of subspecies richness. Both subspecies richness and subspecies diversification rate are found to exhibit low phylogenetic signal, meaning that closely related species do not tend to possess similar numbers of subspecies. As predicted by theory, high subspecies richness was associated with large breeding range size, island dwelling, inhabitation of montane regions, habitat heterogeneity, and low latitude. Of these factors, breeding range size was the variable that explained the most variation. Unravelling whether species that have invaded previously glacial areas have more or fewer subspecies than expected proves to be complicated due to a covariation between the postglacial colonization, latitude, geographic range size, and subspecies richness. However, the effect of postglacial colonization on subspecies richness appears to be small. Mapping the distribution of species' subspecies richness globally reveals geographical patterns that correspond to many of the predictions of the statistical models, but may also reflect geographical variation in taxonomic practice. Overall, we demonstrate that biogeographic models can explain about 30% of the global variation in subspecies richness in birds.