The origin and maintenance of intraspecific variation in vocal signals is important for population divergence and speciation. Where vocalizations are transmitted by vertical cultural inheritance, similarity will reflect co-ancestry, and thus vocal divergence should reflect genetic structure. Horseshoe bats are characterized by echolocation calls dominated by a constant frequency component that is partly determined by maternal imprinting. Although previous studies showed that constant frequency calls are also influenced by some non-genetic factors, it is not known how frequency relates to genetic structure. To test this, we related constant frequency variation to genetic and non-genetic variables in the Formosan lesser horseshoe bat (Rhinolophus monoceros). Recordings of bats from across Taiwan revealed that females called at higher frequencies than males; however, we found no effect of environmental or morphological factors on call frequency. By comparison, variation showed clear population structure, with frequencies lower in the centre and east, and higher in the north and south. Within these regions, frequency divergence was directional and correlated with geographical distance, suggesting that call frequencies are subject to cultural drift. However, microsatellite clustering analysis showed that broad differences in constant frequency among populations corresponded to discontinuities in allele frequencies resulting from vicariant events. Our results provide evidence that the processes shaping genetic subdivision have concomitant consequences for divergence in echolocation call frequency.

Intraspecific phylogenies can provide useful insights into how populations have been shaped by historical and contemporary processes. Taiwan formed around 5 million years ago from tectonic uplift, and has been connected to mainland Asia several times since its emergence. A central mountain range runs north to south, bisecting the island, and potentially impedes gene flow along an east-west axis. The Formosan lesser horseshoe bat (Rhinolophus monoceros) is endemic to Taiwan, where it is found mainly at low altitude. To determine the population structure and the demographic and colonization history of this species, we examined variation in the mitochondrial DNA control region in 203 bats sampled at 26 sites. We found very high haplotype and nucleotide diversity, which decreased from the centre to the south and north. Population differentiation followed a pattern of isolation by distance, though most regional genetic variance was attributable to differences between the relatively isolated southern population and those from other regions. A haplotype network was consistent with these findings and also suggested a southward colonization, followed by subsequent secondary contact between the south and other regions. Mismatch distributions were used to infer a past population expansion predating the last glacial maximum, and a neighbour-joining tree showed that R. monoceros formed a monophyletic grouping with respect to its sister taxa. Taken together, our results suggest that this taxon arose from a single period of colonization, and that demographic growth followed in the late Pleistocene. Current genetic structure reflects limited gene flow, probably coupled with stepwise colonization in the past. We consider explanations for the persistence of the species through multiple glacial maxima.

No abstract

Sex-biased behaviours are expected to play an important role in partitioning genetic variance in animal populations. Comparing genetic structure at markers with different modes of inheritance provides a means of detecting these behaviours and their consequences for population genetic structure. In colonially breeding mammals, the common combination of female philopatry and male vagility can promote contrasting patterns of genetic differentiation between the sexes, both via their effects on recurrent gene flow and on colonization. We examined sex differences in gene flow and structure by comparing maternally inherited mitochondrial DNA (mtDNA) and biparentally inherited autosomal loci in the Formosan lesser horseshoe bat Rhinolophus monoceros. We found that genetic partitioning was higher at mtDNA than autosomal markers in both sexes, indicative of female-biased philopatry and male-biased dispersal. Across Taiwan, isolation-by-distance was detected for all sex/marker combinations but was steeper for mtDNA than for nuclear markers. We suggest that isolation-by-distance shown from mtDNA at large scales is likely to reflect the stepwise founding of new breeding colonies by females during colonization. In contrast, no isolation-bydistance was found at smaller distances of up to 100 km, indicating that gene flow and/or recent shared ancestry homogenises genetic structure among nearby sites. Our results highlight the value of an indirect genetic approach to understanding sex-biased behaviours and their consequences in a little-studied species.

Phylogeographic reconstructions of codistributed taxa can help reveal the interplay between abiotic factors, such as altitude and climate, and species-specific attributes, in shaping patterns of population genetic structure. Recent studies also demonstrate the value of both rangewide sampling and species distribution modelling (SDM) in comparative phylogeography. Here, we combine these approaches to study the population histories of four phylogenetically related forest-dependent bat species. All are endemic to the mountainous island of Taiwan but show differences in their tolerance to altitude, with Murina gracilis considered to be a high-altitude specialist, M. recondita and Kerivoula sp. low-altitude specialists, and M. puta an altitudinal generalist. We tested the prediction that contrasting habitat preferences would impact on patterns of past and contemporary gene flow and found broad concordance between the results of population genetic analyses and species distribution models based on the Model for Interdisciplinary Research on Climate. Both lowland species showed evidence of genetic divergence between the east and west of the island, consistent with SDMs that indicated the Central Mountain Range (CMR) has presented a long-term and continuous barrier to gene flow since before the Last Glacial Maximum. In contrast, Murina gracilis and M. puta showed lower degrees of historical isolation and genetic differentiation associated with the CMR, reflecting greater gene flow, possibly coupled with past population growth in M. puta. Together our results highlight the usefulness of combining distribution models with phylogeographic analyses to understand the drivers of genetic structure.