Search

Interaction network structure reflects the ecological mechanisms acting within biological communities, which are affected by environmental conditions. In tropical forests, higher precipitation usually increases fruit production, which may lead frugivores to increase specialization, resulting in more modular and less nested animal-plant networks. In these ecosystems, El Niño is a major driver of precipitation, but we still lack knowledge of how species interactions change under this influence. To understand bat-plant network structure during an extreme El Niño-Southern Oscillation event, we determined the links between plantivorous bat species and the plants they consume by DNA barcoding seeds and pulp in bat faeces. These interactions were recorded in the dry forest and rainforest of Costa Rica, during the dry and the wet seasons of an extreme El Niño year. From these we constructed seasonal and whole-year bat-plant networks and analysed their structures and dissimilarities. In general, networks had low nestedness, had high modularity, and were dominated by one large compartment which included most species and interactions. Contrary to our expectations, networks were less nested and more modular in drier conditions, both in the comparison between forest types and between seasons. We suggest that increased competition, when resources are scarce during drier seasons and habitats, lead to higher resource partitioning among bats and thus higher modularity. Moreover, we have found similar network structures between dry and rainforests during El Niño and non-El Niño years. Finally, most interaction dissimilarity among networks occurred due to interaction rewiring among species, potentially driven by seasonal changes in resource availability., (© 2022 The Authors. Molecular Ecology published by John Wiley & Sons Ltd.)

Comprising more than 1,400 species, bats possess adaptations unique among mammals including powered flight, unexpected longevity, and extraordinary immunity. Some of the molecular mechanisms underlying these unique adaptations includes DNA repair, metabolism and immunity. However, analyses have been limited to a few divergent lineages, reducing the scope of inferences on gene family evolution across the Order Chiroptera. We conducted an exhaustive comparative genomic study of 37 bat species, one generated in this study, encompassing a large number of lineages, with a particular emphasis on multi-gene family evolution across immune and metabolic genes. In agreement with previous analyses, we found lineage-specific expansions of the APOBEC3 and MHC-I gene families, and loss of the proinflammatory PYHIN gene family. We inferred more than 1,000 gene losses unique to bats, including genes involved in the regulation of inflammasome pathways such as epithelial defence receptors, the natural killer gene complex and the interferon-gamma induced pathway. Gene set enrichment analyses revealed genes lost in bats are involved in defence response against pathogen-associated molecular patterns and damage-associated molecular patterns. Gene family evolution and selection analyses indicate bats have evolved fundamental functional differences compared to other mammals in both innate and adaptive immune system, with the potential to enhance antiviral immune response while dampening inflammatory signalling. In addition, metabolic genes have experienced repeated expansions related to convergent shifts to plant-based diets. Our analyses support the hypothesis that, in tandem with flight, ancestral bats had evolved a unique set of immune adaptations whose functional implications remain to be explored., (© 2021 John Wiley & Sons Ltd.)

Habitat degradation is pervasive across the tropics and is particularly acute in Southeast Asia, with major implications for biodiversity. Much research has addressed the impact of degradation on species diversity; however, little is known about how ecological interactions are altered, including those that constitute important ecosystem functions such as consumption of herbivores. To examine how rainforest degradation alters trophic interaction networks, we applied DNA metabarcoding to construct interaction networks linking forest-dwelling insectivorous bat species and their prey, comparing old-growth forest and forest degraded by logging in Sabah, Borneo. Individual bats in logged rainforest consumed a lower richness of prey than those in old-growth forest. As a result, interaction networks in logged forests had a less nested structure. These network structures were associated with reduced network redundancy and thus increased vulnerability to perturbations in logged forests. Our results show how ecological interactions change between old-growth and logged forests, with potentially negative implications for ecosystem function and network stability., (© 2021 The Authors. Molecular Ecology published by John Wiley & Sons Ltd.)

The application of metabarcoding to environmental and invertebrate-derived DNA (eDNA and iDNA) is a new and increasingly applied method for monitoring biodiversity across a diverse range of habitats. This approach is particularly promising for sampling in the biodiverse humid tropics, where rapid land-use change for agriculture means there is a growing need to understand the conservation value of the remaining mosaic and degraded landscapes. Here we use iDNA from blood-feeding leeches (Haemadipsa picta) to assess differences in mammalian diversity across a gradient of forest degradation in Sabah, Malaysian Borneo. We screened 557 individual leeches for mammal DNA by targeting fragments of the 16S rRNA gene and detected 14 mammalian genera. We recorded lower mammal diversity in the most heavily degraded forest compared to higher quality twice logged forest. Although the accumulation curves of diversity estimates were comparable across these habitat types, diversity was higher in twice logged forest, with more taxa of conservation concern. In addition, our analysis revealed differences between the community recorded in the heavily logged forest and that of the twice logged forest. By revealing differences in mammal diversity across a human-modified tropical landscape, our study demonstrates the value of iDNA as a noninvasive biomonitoring approach in conservation assessments., (© 2020 The Authors. Molecular Ecology published by John Wiley & Sons Ltd.)

Changes in behaviour may initiate shifts to new adaptive zones, with physical adaptations for novel environments evolving later. While new mutations are commonly considered engines of adaptive change, sensory evolution enabling access to new resources might also arise from standing genetic diversity, and even gene loss. We examine the relative contribution of molecular adaptations, measured by positive and relaxed selection, acting on eye-expressed genes associated with shifts to new adaptive zones in ecologically diverse bats from the superfamily Noctilionoidea. Collectively, noctilionoids display remarkable ecological breadth, from highly divergent echolocation to flight strategies linked to specialized insectivory, the parallel evolution of diverse plant-based diets (e.g., nectar, pollen and fruit) from ancestral insectivory, and-unusually for echolocating bats-often have large, well-developed eyes. We report contrasting levels of positive selection in genes associated with the development, maintenance and scope of visual function, tracing back to the origins of noctilionoids and Phyllostomidae (the bat family with most dietary diversity), instead of during shifts to novel diets. Generalized plant visiting was not associated with exceptional molecular adaptation, and exploration of these novel niches took place in an ancestral phyllostomid genetic background. In contrast, evidence for positive selection in vision genes was found at subsequent shifts to either nectarivory or frugivory. Thus, neotropical noctilionoids that use visual cues for identifying food and roosts, as well as for orientation, were effectively preadapted, with subsequent molecular adaptations in nectar-feeding lineages and the subfamily Stenodermatinae of fig-eating bats fine-tuning pre-existing visual adaptations for specialized purposes., (© 2020 The Authors. Molecular Ecology published by John Wiley & Sons Ltd.)

Cases of geographically restricted co-occurring sister taxa are rare and may point to potential divergence with gene flow. The two bat species Murina gracilis and Murina recondita are both endemic to Taiwan and are putative sister species. To test for nonallopatric divergence and gene flow in these taxa, we generated sequences using Sanger and next-generation sequencing, and combined these with microsatellite data for coalescent-based analyses. MtDNA phylogenies supported the reciprocally monophyletic sister relationship between M. gracilis and M. recondita; however, clustering of microsatellite genotypes revealed several cases of species admixture suggesting possible introgression. Sequencing of microsatellite flanking regions revealed that admixture signatures stemmed from microsatellite allele homoplasy rather than recent introgressive hybridization, and also uncovered an unexpected sister relationship between M. recondita and the continental species Murina eleryi, to the exclusion of M. gracilis. To dissect the basis of these conflicts between ncDNA and mtDNA, we analysed sequences from 10 anonymous ncDNA loci with *beast and isolation-with-migration and found two distinct clades of M. eleryi, one of which was sister to M. recondita. We conclude that Taiwan was colonized by the ancestor of M. gracilis first, followed by the ancestor of M. recondita after a period of allopatric divergence. After colonization, the mitochondrial genome of M. recondita was replaced by that of the resident M. gracilis. This study illustrates how apparent signatures of sympatric divergence can arise from complex histories of allopatric divergence, colonization and hybridization, thus highlighting the need for rigorous analyses to distinguish between such scenarios., (© 2015 John Wiley & Sons Ltd.)

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., (© 2014 John Wiley & Sons Ltd.)

The demographic history of Rhinolophus hipposideros (lesser horseshoe bat) was reconstructed across its European, North African and Middle-Eastern distribution prior to, during and following the most recent glaciations by generating and analysing a multimarker data set. This data set consisted of an X-linked nuclear intron (Bgn; 543 bp), mitochondrial DNA (cytb-tRNA-control region; 1630 bp) and eight variable microsatellite loci for up to 373 individuals from 86 localities. Using this data set of diverse markers, it was possible to determine the species' demography at three temporal stages. Nuclear intron data revealed early colonization into Europe from the east, which pre-dates the Quaternary glaciations. The mtDNA data supported multiple glacial refugia across the Mediterranean, the largest of which were found in the Ibero-Maghreb region and an eastern location (Anatolia/Middle East)-that were used by R. hipposideros during the most recent glacial cycles. Finally, microsatellites provided the most recent information on these species' movements since the Last Glacial Maximum and suggested that lineages that had diverged into glacial refugia, such as in the Ibero-Maghreb region, have remained isolated. These findings should be used to inform future conservation management strategies for R. hipposideros and show the power of using a multimarker data set for phylogeographic studies., (© 2013 John Wiley & Sons Ltd.)

Phylogenetic conflicts between genetic markers can help to disentangle complex histories of phylogeography and introgression among taxa. We previously proposed that the Chinese mainland subspecies of the intermediate horseshoe bat Rhinolophus affinis himalayanus colonized Hainan Island to form the subspecies R. a. hainanus. Subsequent recolonization of the mainland formed a third taxon, R. a macrurus, and a secondary contact zone with the ancestral himalayanus. To test for historical and recurrent genetic exchange between these mainland subspecies, we sampled populations of each from two parapatric zones and undertook analyses using one mtDNA marker, three nuclear genes and 14 microsatellites. Nuclear DNA, echolocation call frequencies and morphological data all recovered two taxa; however, a mtDNA phylogeny revealed two himalayanus clades, of which one clustered with macrurus, as well as some shared or related mtDNA haplotypes in eastern populations. Isolation-with-migration (IM) models suggested some mtDNA gene flow from macrurus to himalayanus. However, strong population structure in himalayanus raises the possibility that macrurus captured mtDNA from a coastal population of himalayanus that has since become rare or extinct. To reconcile these two sets of results, we suggest that the IM estimates might reflect historical mtDNA gene flow among populations of himalayanus, before mtDNA was subsequently captured by macrurus. Finally, microsatellite-based ABC analyses supported the island origin of macrurus; however, mtDNA-based ABC analyses suggest this taxon might have evolved on the mainland. Our findings highlight the importance of understanding population history and structure for interpreting hybridization and introgression events., (© 2012 Blackwell Publishing Ltd.)

The impact of ecology and social organization on genetic structure at landscape spatial scales, where gene dynamics shape evolution as well as determine susceptibility to habitat fragmentation, is poorly understood. Attempts to assess these effects must take into account the potentially confounding effects of history. We used microsatellites to compare genetic structure in seven bat species with contrasting patterns of roosting ecology and social organization, all of which are codistributed in an ancient forest habitat that has been exceptionally buffered from radical habitat shifts. Over one thousand individuals were captured at foraging sites and genotyped at polymorphic microsatellite loci. Analyses of spatially explicit genotype data revealed interspecies differences in the extent of movement and gene flow and genetic structure across continuous intact forest. Highest positive genetic structure was observed in tree-roosting taxa that roost either alone or in small groups. By comparison, a complete absence of genetic autocorrelation was noted in the cave-roosting colonial species across the study area. Our results thus reveal measurable interspecies differences in the natural limits of gene flow in an unmodified habitat, which we attribute to contrasting roosting ecology and social organization. The consequences of ecology and behaviour for gene flow have important implications for conservation. In particular, tree-roosting species characterized by lower vagility and thus gene flow will be disproportionally impacted by landscape-scale forest clearance and habitat fragmentation, which are prevalent in the study region. Our method also highlights the usefulness of rapid sampling of foraging bats for assaying genetic structure, particularly where roosting sites are not always known., (© 2011 Blackwell Publishing Ltd.)

How males gain access to mates and the potential for female choice will determine whether polygyny can operate at several levels, from within litters and groups to the wider population. Female lesser flat-headed bats (Tylonycteris pachypus) form maternity groups in bamboo stems. Unusually for bats, they are multiparous, providing the opportunity to test whether multi-level polygyny differs among males depending on whether they roost with females, with males or are solitary. We genotyped 662 individuals from 54 internodes and analysed parentage of 165 litters. Our results revealed 170 sets of paternal twins/triplets, of which 96 were full-sibs and 74 were half-sibs. We found that males captured roosting with females typically sired more offspring overall than did other males and also showed a greater tendency to monopolize paternity within both litters and roosting groups. In comparison, males that sired fewer full-sibs were assigned more maternal half-sibs. These latter individuals, which included solitary males and those from all-male groups, might gain copulations either via roaming with furtive mating or during visits by females. Indeed, female lesser flat-headed bats store sperm, so could benefit from multiple mating to reduce genetic incompatibilities. At the same time, however, we found no evidence of outbreeding. Finally, relatedness and mtDNA analyses revealed that polygyny also operated within matrilineal kin, suggesting a system that might promote social cohesiveness. Future studies of individual movements will help to determine the extent to which mixed paternities in litters, matrilines and groups are driven by male or female behaviour., (© 2011 Blackwell Publishing Ltd.)

The repeated formation and loss of land-bridges during the Pleistocene have had lasting impacts on population genetic structure. In the tropics, where island populations persisted through multiple glacial cycles, alternating periods of isolation and contact are expected to have driven population and taxonomic divergence. Here, we combine mitochondrial and nuclear sequence data with microsatellites to dissect the impact of Pleistocene climate change on intra-specific diversification in the horseshoe bat Rhinolophus affinis. This taxon shows considerable morphological and acoustic variation: two parapatric subspecies (himalayanus and macrurus) occur on mainland China and a third (hainanus) on Hainan Island. Our phylogeographic reconstruction and coalescent analyses suggest the island subspecies formed from an ancestral population of himalayanus via two colonization events c. 800,000 years before present. R. a. hainanus then recolonized the mainland, forming macrurus and thus a secondary contact zone with himalayanus. Finally, macrurus recolonized Hainan following the LGM. We found that all three biological events corresponded to known periods of land-bridge formation. Evidence of introgression was detected between macrurus and both its sister taxa, with geographical proximity rather than length of separation appearing to be the biggest determinant of subsequent genetic exchange. Our study highlights the important role of climate-mediated sea level changes have had in shaping current processes and patterns of population structure and taxonomic diversification.

Instances of hybridization between mammalian taxa in the wild are rarely documented. To test for introgression between sibling species of horseshoe bat (Rhinolophus yunanensis and R. pearsoni) and two subspecies of the latter (R. p. pearsoni and R. p. chinensis), we sequenced two mtDNA and two ncDNA markers in individuals sampled from multiple localities within their overlapping ranges. The interspecific mtDNA gene tree corresponded to the expected taxonomic divisions, and coalescent-based analyses suggested divergence occurred around 4 MYA. However, these relationships strongly conflicted with those recovered from two independent nuclear gene trees, in which R. yunanensis clustered with R. p. pearsoni to the exclusion of R. p. chinensis. This geographically widespread discordance is best explained by large-scale historical introgression of ncDNA from R. yunanensis to R. pearsoni by male-mediated exchange in mixed species colonies during Pleistocene glacial periods, when ranges may have contracted and overlapped more than at present. Further species tree-gene tree conflicts were detected between R. p. pearsoni and R. p. chinensis, also indicating past and/or current introgression in their overlapping regions. However, here the patterns point to asymmetric mtDNA introgression without ncDNA introgression. Analyses of coalescence times indicate this exchange has occurred subsequent to the divergence of these subspecies from their common ancestor. Our work highlights the importance of using multiple data sets for reconstructing phylogeographic histories and resolving taxonomic relationships.

Phylogeographical studies are typically based on haplotype data, occasionally on nuclear markers such as microsatellites, but rarely combine both. This is unfortunate because the use of markers with contrasting modes of inheritance and rates of evolution might provide a more accurate and comprehensive understanding of a species' history. Here we present a detailed study of the phylogeography of the greater horseshoe bat, Rhinolophus ferrumequinum, using 1098 bp of the mitochondrial ND2 gene from 45 localities from across its Palaearctic range to infer population history. In addition, we re-analysed a large microsatellite data set available for this species and compared the results of both markers to infer population relationships and the historical processes influencing them. We show that mtDNA, the most popular marker in phylogeography studies, yielded a misleading result, and would have led us to conclude erroneously that a single expansion had taken place in Europe. Only by combining the mitochondrial and microsatellite data sets are we able to reconstruct the species' history and show two colonization events in Europe, one before the Last Glacial Maximum (LGM) and one after it. Combining markers also revealed the importance of Asia Minor as an ancient refugium for this species and a source population for the expansion of the greater horseshoe bat into Europe before the LGM.

The distribution of genetic variability across a species' range can provide valuable insights into colonization history. To assess the relative importance of European and Asian refugia in shaping current levels of genetic variation in the greater horseshoe bats, we applied a microsatellite-based approach to data collected from 56 localities ranging from the UK to Japan. A decline in allelic richness from west Asia to the UK and analyses of F(ST) both imply a northwestward colonization across Europe. However, sharp discontinuities in gene frequencies within Europe and between the Balkans and west Asia (Syria/Russia) are consistent with suture zones following expansion from multiple refugia, and a lack of recent gene flow from Asia Minor. Together, these results suggest European populations originated from west Asia in the ancient past, and experienced a more recent range expansion since the Last Glacial Maximum. Current populations in central Europe appear to originate from the Balkans and those from west Europe from either Iberia and/or Italy. Comparisons of R(ST )and F(ST) suggest that stepwise mutation has contributed to differentiation between island and continental populations (France/UK and China/Japan) and also among distant samples. However, pairwise R(ST) values between distant populations appear to be unreliable, probably due to size homoplasy. Our findings also highlight two priorities for conservation. First, stronger genetic subdivision within the UK than across 4000 km of continental Eurasia is most likely the result of population fragmentation and highlights the need to maintain gene flow in this species. Second, deep splits within China and between Europe and China are indicative of cryptic taxonomic divisions which need further investigation.

Most mammalian groups are characterized by male-biased sexual size dimorphism, in which size-dependent male-male competition and reproductive skew are tightly linked. By comparison, little is known about the opportunity for sexual selection in mammalian systems without male-biased dimorphism, where the traits under sexual selection might be less obvious. We examined 10 years of parentage data in a colony of greater horseshoe bats (Rhinolophus ferrumequinum) to determine the magnitude of male reproductive skew and the opportunity for sexual selection in a mammal in which females are the larger sex. Annual paternity success was weakly skewed but consistent patterns led to strong longitudinal paternity skew among breeders. Just three males accounted for a third of all paternity assignments, representing at least a fifth of all colony offspring born in a decade. Paternity success was in part determined by age but was not influenced by dispersal status. Our results show that paternity skew and the opportunity for sexual selection in a species with reversed sexual size dimorphism can approach levels reported for classical examples of species with polygyny and male-biased dimorphism, even where the traits under sexual selection are not known.

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.

Following a dramatic decline last century, the British population of the endangered greater horseshoe bat Rhinolophus ferrumequinum is highly fragmented. To examine the consequences of fragmentation and limited dispersal on patterns of genetic structure and variation, we used microsatellite markers to screen bats from around 50% of the known maternity colonies in Britain, and two areas from continental Europe. Analyses revealed that Welsh and English colonies were genetically isolated. This, and lower variability in Britain than north France, may result from either genetic drift, or the species' colonization history. Gene flow among most neighbouring colonies was not generally restricted, with one exception. These findings have important implications for the ongoing conservation management of this species.