Sex-biased gene flow and colonization in the Formosan lesser horseshoe bat: inference from nuclear and mitochondrial markers

Estimating and mitigating amplification bias in qualitative and quantitative arthropod metabarcoding

Amplicon based metabarcoding promises rapid and cost-efficient analyses of species composition. However, it is disputed whether abundance estimates can be derived from metabarcoding due to taxon specific PCR amplification biases. PCR-free approaches have been suggested to mitigate this problem, but come with considerable increases in workload and cost. Here, we analyze multilocus datasets of diverse arthropod communities, to evaluate whether amplification bias can be countered by (1) targeting loci with highly degenerate primers or conserved priming sites, (2) increasing PCR template concentration, (3) reducing PCR cycle number or (4) avoiding locus specific amplification by directly sequencing genomic DNA. Amplification bias is reduced considerably by degenerate primers or targeting amplicons with conserved priming sites. Surprisingly, a reduction of PCR cycles did not have a strong effect on amplification bias. The association of taxon abundance and read count was actually less predictable with fewer cycles. Even a complete exclusion of locus specific amplification did not exclude bias. Copy number variation of the target loci may be another explanation for read abundance differences between taxa, which would affect amplicon based and PCR free methods alike. As read abundance biases are taxon specific and predictable, the application of correction factors allows abundance estimates.

Female philopatry and limited male-biased dispersal in the Ussuri tube-nosed bat, Murina ussuriensis

Male-biased dispersal and female philopatry are common traits among social mammals, often leading to elevated relatedness within social groups. However, exceptions do occur, with documented cases of female-biased dispersal, dispersal by both sexes, and philopatry of both sexes. In this study, we examined levels of dispersal and relatedness based on analyses of the mitochondrial D-loop and 12 nuclear microsatellite markers in the woodland specialist Ussuri tube-nosed bat ( Murina ussuriensis ), a relatively widespread, yet locally rare species that is thought to be threatened by loss of its forest habitats across its range. A total of 85 individuals were captured in lowland temperate forest on the island of Yakushima and 28 individuals in montane forest in Hokkaido, Japan. In the former, haplotypes showed extreme spatial clustering among females consistent with strong philopatry, but spatial mixing among males suggesting dispersal over short distances. These findings were broadly supported by microsatellite analyses, which indicated considerable genealogical structure within sampling locations but a lack of spatial structure, again indicating that some gene flow does occur in one or both sexes. The Hokkaido data, although limited, did not show these patterns and instead suggest that differences in environmental and behavioral variables may influence movement ecology of individuals at these 2 sites and consequently fine-scale genetic structure within this species.

Genetic connectivity among swarming sites in the wide ranging and recently declining little brown bat (Myotis lucifugus)

Characterizing movement dynamics and spatial aspects of gene flow within a species permits inference on population structuring. As patterns of structuring are products of historical and current demographics and gene flow, assessment of structure through time can yield an understanding of evolutionary dynamics acting on populations that are necessary to inform management. Recent dramatic population declines in hibernating bats in eastern North America from white-nose syndrome have prompted the need for information on movement dynamics for multiple bat species. We characterized population genetic structure of the little brown bat, Myotis lucifugus, at swarming sites in southeastern Canada using 9 nuclear microsatellites and a 292-bp region of the mitochondrial genome. Analyses of F ST, ΦST, and Bayesian clustering (STRUCTURE) found weak levels of genetic structure among swarming sites for the nuclear and mitochondrial genome (Global F ST = 0.001, P < 0.05, Global ΦST = 0.045, P < 0.01, STRUCTURE K = 1) suggesting high contemporary gene flow. Hierarchical AMOVA also suggests little structuring at a regional (provincial) level. Metrics of nuclear genetic structure were not found to differ between males and females suggesting weak asymmetries in gene flow between the sexes. However, a greater degree of mitochondrial structuring does support male-biased dispersal long term. Demographic analyses were consistent with past population growth and suggest a population expansion occurred from approximately 1250 to 12,500 BP, following Pleistocene deglaciation in the region. Our study suggests high gene flow and thus a high degree of connectivity among bats that visit swarming sites whereby mainland areas of the region may be best considered as one large gene pool for management and conservation.

Sensory trait variation in an echolocating bat suggests roles for both selection and plasticity

Background

Across heterogeneous environments selection and gene flow interact to influence the rate and extent of adaptive trait evolution. This complex relationship is further influenced by the rarely considered role of phenotypic plasticity in the evolution of adaptive population variation. Plasticity can be adaptive if it promotes colonization and survival in novel environments and in doing so may increase the potential for future population differentiation via selection. Gene flow between selectively divergent environments may favour the evolution of phenotypic plasticity or conversely, plasticity itself may promote gene flow, leading to a pattern of trait differentiation in the presence of gene flow. Variation in sensory traits is particularly informative in testing the role of environment in trait and population differentiation. Here we test the hypothesis of ‘adaptive differentiation with minimal gene flow’ in resting echolocation frequencies (RF) of Cape horseshoe bats (Rhinolophus capensis) across a gradient of increasingly cluttered habitats.

Results

Our analysis reveals a geographically structured pattern of increasing RF from open to highly cluttered habitats in R. capensis; however genetic drift appears to be a minor player in the processes influencing this pattern. Although Bayesian analysis of population structure uncovered a number of spatially defined mitochondrial groups and coalescent methods revealed regional-scale gene flow, phylogenetic analysis of mitochondrial sequences did not correlate with RF differentiation. Instead, habitat discontinuities between biomes, and not genetic and geographic distances, best explained echolocation variation in this species. We argue that both selection for increased detection distance in relatively less cluttered habitats and adaptive phenotypic plasticity may have influenced the evolution of matched echolocation frequencies and habitats across different populations.

Conclusions

Our study reveals significant sensory trait differentiation in the presence of historical gene flow and suggests roles for both selection and plasticity in the evolution of echolocation variation in R. capensis. These results highlight the importance of population level analyses to i) illuminate the subtle interplay between selection, plasticity and gene flow in the evolution of adaptive traits and ii) demonstrate that evolutionary processes may act simultaneously and that their relative influence may vary across different environments.

Phylogeography of the Rickett's big-footed bat, Myotis pilosus (Chiroptera: Vespertilionidae): a novel pattern of genetic structure of bats in China

Background

China is characterized by complex topographic structure and dramatic palaeoclimatic changes, making species biogeography studies particularly interesting. Previous researchers have also demonstrated multiple species experienced complex population histories, meanwhile multiple shelters existed in Chinese mainland. Despite this, species phylogeography is still largely unexplored. In the present study, we used a combination of microsatellites and mitochondrial DNA (mtDNA) to investigate the phylogeography of the east Asian fish-eating bat (Myotis pilosus).

Results

Phylogenetic analyses showed that M. pilosus comprised three main lineages: A, B and C, which corresponded to distinct geographic populations of the Yangtze Plain (YTP), Sichuan Basin (SCB) and North and South of China (NSC), respectively. The most recent common ancestor of M. pilosus was dated as 0.25 million years before present (BP). Population expansion events were inferred for populations of Clade C, North China Plain region, Clade B and YunGui Plateau region at 38,700, 15,900, 4,520 and 4,520 years BP, respectively. Conflicting results were obtained from mtDNA and microsatellite analyses; strong population genetic structure was obtained from mtDNA data but not microsatellite data. The microsatellite data indicated that genetic subdivision fits an isolation-by-distance (IBD) model, but the mtDNA data failed to support this model.

Conclusions

Our results suggested that Pleistocene climatic oscillations might have had a profound influence on the demographic history of M. pilosus. Spatial genetic structures of maternal lineages that are different from those observed in other sympatric bats species may be as a result of interactions among special population history and local environmental factors. There are at least three possible refugia for M. pilosus during glacial episodes. Apparently contradictory genetic structure patterns of mtDNA and microsatellite could be explained by male-mediated gene flow among populations. This study also provides insights on the necessity of conservation of M. pilosus populations to conserve this genetic biodiversity, especially in the areas of YTP, SCB and NSC regions.

Demographic loss, genetic structure and the conservation implications for Indian tigers

India is home to approximately 60 per cent of the world's remaining wild tigers, a species that has declined in the last few centuries to occupy less than 7 per cent of its former geographical range. While Indian tiger numbers have somewhat stabilized in recent years, they remain low and populations are highly fragmented. Therefore, the application of evidence-based demographic and genetic management to enhance the remaining populations is a priority. In this context, and using genetic data from historical and modern tigers, we investigated anthropogenic impacts on genetic variation in Indian tigers using mitochondrial and nuclear genetic markers. We found a very high number of historical mitochondrial DNA variants, 93 per cent of which are not detected in modern populations. Population differentiation was higher in modern tigers. Simulations incorporating historical data support population decline, and suggest high population structure in extant populations. Decreased connectivity and habitat loss as a result of ongoing fragmentation in the Indian subcontinent has therefore resulted in a loss of genetic variants and increased genetic differentiation among tiger populations. These results highlight that anthropogenic fragmentation and species-specific demographic processes can interact to alter the partitioning of genetic variation over very short time scales. We conclude that ongoing strategies to maximize the size of some tiger populations, at the expense of losing others, is an inadequate conservation strategy, as it could result in a loss of genetic diversity that may be of adaptive significance for this emblematic species.

Dispersal, mating events and fine-scale genetic structure in the lesser flat-headed bats

Population genetic structure has important consequences in evolutionary processes and conservation genetics in animals. Fine-scale population genetic structure depends on the pattern of landscape, the permanent movement of individuals, and the dispersal of their genes during temporary mating events. The lesser flat-headed bat (Tylonycteris pachypus) is a nonmigratory Asian bat species that roosts in small groups within the internodes of bamboo stems and the habitats are fragmented. Our previous parentage analyses revealed considerable extra-group mating in this species. To assess the spatial limits and sex-biased nature of gene flow in the same population, we used 20 microsatellite loci and mtDNA sequencing of the ND2 gene to quantify genetic structure among 54 groups of adult flat-headed bats, at nine localities in South China. AMOVA and F(ST) estimates revealed significant genetic differentiation among localities. Alternatively, the pairwise F(ST) values among roosting groups appeared to be related to the incidence of associated extra-group breeding, suggesting the impact of mating events on fine-scale genetic structure. Global spatial autocorrelation analyses showed positive genetic correlation for up to 3 km, indicating the role of fragmented habitat and the specialized social organization as a barrier in the movement of individuals among bamboo forests. The male-biased dispersal pattern resulted in weaker spatial genetic structure between localities among males than among females, and fine-scale analyses supported that relatedness levels within internodes were higher among females than among males. Finally, only females were more related to their same sex roost mates than to individuals from neighbouring roosts, suggestive of natal philopatry in females.

Mitochondrial DNA paradox: sex-specific genetic structure in a marine mussel--despite maternal inheritance and passive dispersal

Background

When genetic structure is identified using mitochondrial DNA (mtDNA), but no structure is identified using biparentally-inherited nuclear DNA, the discordance is often attributed to differences in dispersal potential between the sexes.

Results

We sampled the intertidal rocky shore mussel Perna perna in a South African bay and along the nearby open coast, and sequenced maternally-inherited mtDNA (there is no evidence for paternally-inherited mtDNA in this species) and a biparentally-inherited marker. By treating males and females as different populations, we identified significant genetic structure on the basis of mtDNA data in the females only.

Conclusions

This is the first study to report sex-specific differences in genetic structure based on matrilineally-inherited mtDNA in a passively dispersing species that lacks social structure or sexual dimorphism. The observed pattern most likely stems from females being more vulnerable to selection in habitats from which they did not originate, which also manifests itself in a male-biased sex ratio. Our results have three important implications for the interpretation of population genetic data. First, even when mtDNA is inherited exclusively in the female line, it also contains information about males. For that reason, using it to identify sex-specific differences in genetic structure by contrasting it with biparentally-inherited markers is problematic. Second, the fact that sex-specific differences were found in a passively dispersing species in which sex-biased dispersal is unlikely highlights the fact that significant genetic structure is not necessarily a function of low dispersal potential or physical barriers. Third, even though mtDNA is typically used to study historical demographic processes, it also contains information about contemporary processes. Higher survival rates of males in non-native habitats can erase the genetic structure present in their mothers within a single generation.

A tale of two genomes: contrasting patterns of phylogeographic structure in a widely distributed bat

One of the most widely distributed bats in the New World, the big brown bat (Eptesicus fuscus) exhibits well-documented geographic variation in morphology and life history traits, suggesting the potential for significant phylogeographic structure as well as adaptive differentiation among populations. In a pattern broadly consistent with morphologically defined subspecies, we found deeply divergent mitochondrial lineages restricted to different geographic regions. In contrast, sequence data from two nuclear loci suggest a general lack of regional genetic structure except for peripheral populations in the Caribbean and Mexico/South America. Coalescent analyses suggest that the striking difference in population structure between genomes cannot be attributed solely to different rates of lineage sorting, but is likely due to male-mediated gene flow homogenizing nuclear genetic diversity across most of the continental range. Despite this ongoing gene flow, selection has apparently been effective in producing and maintaining adaptive differentiation among populations, while strong female site fidelity, maintained over the course of millions of years, has produced remarkably deep divergence among geographically isolated matrilines. Our results highlight the importance of evaluating multiple genetic markers for a more complete understanding of population structure and history.

Ecology driving genetic variation: a comparative phylogeography of jungle cat (Felis chaus) and leopard cat (Prionailurus bengalensis) in India

BACKGROUND

Comparative phylogeography links historical population processes to current/ecological processes through congruent/incongruent patterns of genetic variation among species/lineages. Despite high biodiversity, India lacks a phylogeographic paradigm due to limited comparative studies. We compared the phylogenetic patterns of Indian populations of jungle cat (Felis chaus) and leopard cat (Prionailurus bengalensis). Given similarities in their distribution within India, evolutionary histories, body size and habits, congruent patterns of genetic variation were expected.

METHODOLOGY/PRINCIPAL FINDINGS

We collected scats from various biogeographic zones in India and analyzed mtDNA from 55 jungle cats (460 bp NADH5, 141 bp cytochrome b) and 40 leopard cats (362 bp NADH5, 202 bp cytochrome b). Jungle cats revealed high genetic variation, relatively low population structure and demographic expansion around the mid-Pleistocene. In contrast, leopard cats revealed lower genetic variation and high population structure with a F(ST) of 0.86 between North and South Indian populations. Niche-model analyses using two approaches (BIOCLIM and MaxEnt) support absence of leopard cats from Central India, indicating a climate associated barrier. We hypothesize that high summer temperatures limit leopard cat distribution and that a rise in temperature in the peninsular region of India during the LGM caused the split in leopard cat population in India.

CONCLUSIONS/SIGNIFICANCE

Our results indicate that ecological variables describing a species range can predict genetic patterns. Our study has also resolved the confusion over the distribution of the leopard cat in India. The reciprocally monophyletic island population in the South mandates conservation attention.

Determinants of echolocation call frequency variation in the Formosan lesser horseshoe bat (Rhinolophus monoceros)

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.

Why the Indian subcontinent holds the key to global tiger recovery

With only ∼3,000 wild individuals surviving restricted to just 7% of their historical range, tigers are now a globally threatened species. Therefore, conservation efforts must prioritize regions that harbor more tigers, as well try to capture most of the remaining genetic variation and habitat diversity. Only such prioritization based on demographic, genetic, and ecological considerations can ensure species recovery and retention of evolutionary flexibility in the face of ongoing global changes. Although scientific understanding of ecological and demographic aspects of extant wild tiger populations has improved recently, little is known about their genetic composition and variability. We sampled 73 individual tigers from 28 reserves spread across a diversity of habitats in the Indian subcontinent to obtain 1,263 bp of mitochondrial DNA and 10 microsatellite loci. Our analyses reveals that Indian tigers retain more than half of the extant genetic diversity in the species. Coalescent simulations attribute this high genetic diversity to a historically large population size of about 58,200 tigers for peninsular India south of the Gangetic plains. Furthermore, our analyses indicate a precipitous, possibly human-induced population crash ∼200 years ago in India, which is in concordance with historical records. Our results suggest that only 1.7% (with an upper limit of 13% and a lower limit of 0.2%) of tiger numbers in historical times remain now. In the global conservation context our results suggest that, based on genetic, demographic, and ecological considerations, the Indian subcontinent holds the key to global survival and recovery of wild tigers.

Tiger range and numbers have collapsed globally despite substantial conservation efforts. Genetic data quantifying variation from 73 wild tigers in 28 reserves in the Indian subcontinent suggests historically high numbers for tigers, and simulations reveal a signature of a 200-year-old, possibly human-induced decline. Simulations suggest that only 1.7% of historical tiger numbers now persist in peninsular Indian. Our data also reveal that tigers of the Indian subcontinent retain most of the species' genetic diversity, besides this region harbouring maximum diversity of tiger habitats. Overall, the Indian subcontinent appears to be a global hotspot holding the key to any future recovery of wild tigers from both an ecological and genetic perspective.

Asymmetrical male-mediated gene flow between harbor seal (Phoca vitulina) populations in Alaska

Harbor seals ( Phoca vitulina richardii (Gray, 1864)) in Alaska are currently treated as three distinct management stocks. Previous genetic analyses using mitochondrial DNA suggested that these stocks are differentiated genetically. We studied populations in Glacier Bay (GB; Southeast Alaska Stock), where harbor seals are declining, and Prince William Sound (PWS; Gulf of Alaska Stock), where the population has recently stabilized. Using six pairs of hypervariable microsatellite primers, we determined that these populations are a single panmictic unit with estimated migration rates of 22 animals/generation (PWS to GB) and 63 animals/generation (GB to PWS). The asymmetrical gene flow between GB and PWS is likely driven in part by a recent increase in competitors and predators of seals in GB. In contrast with males, emigration of females from PWS to GB (8.3 seals/generation) is higher than emigration of females from GB to PWS (3.3 seals/generation), likely because females use glacial ice as pupping habitat. Despite the high gene flow, the number of migrants per year (0.02% of the Gulf of Alaska population) is likely too low to influence the demographics of harbor seals in PWS, and the two populations may best be managed as separate stocks.