Contrasted patterns of mitochondrial and nuclear structure among nursery colonies of the bat Myotis myotis

Thirteen nursery colonies of Myotis myotis were sampled in central Europe to investigate the dispersal behaviour of this bat species. Mitochondrial DNA sequences of 260 bats reveal the occurrence of three evolutionary lineages that have probably originated in distinct glacial refugia and meet in a contact zone near the Alps. Moreover, the strong haplotypic segregation (Φ_ST=0.540) suggests that breeding females are philopatric. Contrastingly, the low population structure at 15 microsatellite loci (F_ST=0.022), suggests the homogenizing effect of nuclear gene flow. The different perspectives given by these two markers are consistent with strong male‐biased dispersal. As a result of female philopatry, the local haplotypic variability seems to be largely influenced by historical processes of colonization. Conversely, the homogeneity of nuclear variability within roosts that are located north of the Alps seems to mainly reflect contemporary gene flow. Finally, despite the fact that females are faithful to their natal colony, movements of both males and females occur outside the breeding period. Mitochondrial survey of individuals sampled exclusively in nurseries may thus poorly reflect the metapopulation dynamics of this species.     

Italy as a major Ice Age refuge area for the bat Myotis myotis (Chiroptera: Vespertilionidae) in Europe

The distribution of biota from the temperate regions changed considerably during the climatic fluctuations of the Quaternary. This is especially true for many bat species that depend on warm roosts to install their nursery colonies. Surveys of genetic variation among European bats have shown that the southern peninsulas (Iberia and the Balkans) harbour endemic diversity, but to date, no such surveys have been conducted in the third potential glacial refuge area, the Apennine peninsula. We report here the phylogeographical analysis of 115 greater mouse-eared bats ( Myotis myotis ) sampled throughout Italy, and show that 15 of the 18 different haplotypes found in the mitochondrial control region of these bats were unique to the Apennine peninsula. Colonies within this region also showed substantial genetic structure at both mitochondrial ( Φ_ST  = 0.47, P  < 0.001) and nuclear markers ( F _ST = 0.039, P  < 0.001). Based on a comprehensive survey of 575 bats from Europe, these genetic markers further indicate that central Italian populations of M. myotis are more closely related to Greek samples from across the Adriatic Sea, than to other Italian bats. Mouse-eared bat populations from the Apennine peninsula thus represent a complex mixture of several endemic lineages, which evolved in situ , with others that colonized this region more recently along an Adriatic route. Our broad survey also confirms that the Alps represent a relatively impermeable barrier to gene flow for Apennine lineages, even for vagile animals such as bats. These results underline the conservation value of bats from this region and the need to include the Apennine peninsula in phylogeographical surveys in order to provide a more accurate view of the evolution of bats in Europe.     

Genetic isolation of insular populations of the Maghrebian bat,Myotis punicus,in the Mediterranean Basin

Aim We investigate the population genetic structure of the Maghrebian bat, Myotis punicus, between the mainland and islands to assess the island colonization pattern and current gene flow between nearby islands and within the mainland. Location North Africa and the Mediterranean islands of Corsica and Sardinia. Methods We sequenced part of the control region (HVII) of 79 bats across 11 colonies. The phylogeographical pattern was assessed by analysing molecular diversity indices, examining differentiation among populations and estimating divergence time. In addition, we genotyped 182 bats across 10 colonies at seven microsatellite loci. We used analysis of molecular variance and a Bayesian approach to infer nuclear population structure. Finally, we estimated sex‐specific dispersal between Corsica and Sardinia. Results Mitochondrial analyses indicated that colonies between Corsica, Sardinia and North Africa are highly differentiated. Within islands there was no difference between colonies, while at the continental level Moroccan and Tunisian populations were highly differentiated. Analyses with seven microsatellite loci showed a similar pattern. The sole difference was the lack of nuclear differentiation between populations in North Africa, suggesting a male‐biased dispersal over the continental area. The divergence time of Sardinian and Corsican populations was estimated to date back to the early and mid‐Pleistocene. Main conclusions Island colonization by the Maghrebian bats seems to have occurred in a stepping‐stone manner and certainly pre‐dated human colonization. Currently, open water seems to prevent exchange of bats between the two islands, despite their ability to fly and the narrowness of the strait of Bonifacio. Corsican and Sardinian populations are thus currently isolated from any continental gene pool and must therefore be considered as different evolutionarily significant units (ESU).     

Some genetic consequences of ice ages, and their role, in divergence and speciation

The genetic effects of pleistocene ice ages are approached by deduction from paleoenvironmental information, by induction from the genetic structure of populations and species, and by their combination to infer likely consequences. (1) Recent palaeoclimatic information indicate rapid global reversals and changes in ranges of species which would involve elimination with spreading from the edge. Leading edge colonization during a rapid expansion would be leptokurtic and lead to homozygosity and spatial assortment of genomes. In Europe and North America, ice age contractions were into southern refugia, which would promote genome reorganization. (2) The present day genetic structure of species shows frequent geographic subdivision, with parapatric genomes, hybrid zones and suture zones. A survey of recent DNA phylogeographic information supports and extends earlier work. (3) The grasshopper Chorthippus parallelus is used to illustrate such data and processes. Its range in Europe is divided on DNA sequences into five parapatric races, with southern genomes showing greater haplotype diversity - probably due to southern mountain blocks acting as refugia and northern expansion reducing diversity. (4) Comparison with other recent studies shows a concordance of such phylogeographic data over pleistocene time scales. (5) The role that ice age range changes may have played in changing adaptations is explored, including the limits of range, rapid change in new invasions and refugial differentiation in a variety of organisms. (6) The effects of these events in causing divergence and speciation are explored using Chorthippus as a paradigm. Repeated contraction and expansion would accumulate genome differences and adaptations, protected from mixing by hybrid zones, and such a composite mode of speciation could apply to many organisms.     

The legacy of ice ages in mountain species: post‐glacial colonization of mountain tops rather than current range fragmentation determines mitochondrial genetic diversity in an endemic Pyrenean rock lizard

Aim The genetic impact of Quaternary climatic fluctuations on mountain endemic species has rarely been investigated. The Pyrenean rock lizard (Iberolacerta bonnali) is restricted to alpine habitats in the Pyrenees where it exhibits a highly fragmented distribution between massifs and between habitats within massifs. Using mitochondrial DNA markers, we set out: (1) to test whether several evolutionary units exist within the species; (2) to understand how the species persisted through the Last Glacial Maximum and whether the current range fragmentation originates from distribution shifts after the Last Glacial Maximum or from more ancient events; and (3) to investigate whether current mitochondrial diversity reflects past population history or current habitat fragmentation. Location The Pyrenees in south‐western France and northern Spain. Methods We used variation in the hypervariable left domain of the mitochondrial control region of 146 lizards collected in 15 localities, supplemented by cytochrome b sequences downloaded from GenBank to cover most of the species’ distribution range. Measures of population genetic diversity were contrasted with population isolation inferred from topography. Classical (F‐statistics) and coalescence‐based methods were used to assess the level of gene flow and estimate divergence time between populations. We used coalescence‐based simulations to test the congruence of our genetic data with a scenario of simultaneous divergence of current populations. Results Coalescence‐based analyses suggested that these peripheral populations diverged simultaneously at the end of the last glacial episode when their habitats became isolated on mountain summits. High mitochondrial diversity was found in peripheral, isolated populations, while the populations from the core of the species’ range were genetically impoverished. Where mitochondrial diversity has been retained, populations within the same massif exhibited high levels of genetic differentiation. Main conclusions As suggested for many other mountain species, the Pyrenean rock lizard survived glacial maxima through short‐distance range shifts instead of migration or contraction in distant southern refugia. Most of the main Pyrenean range has apparently been re‐colonized from a single or a few source populations, resulting in a loss of genetic diversity in re‐colonized areas. As a result, current levels of intra‐population mitochondrial diversity are better explained by post‐glacial population history than by current habitat fragmentation. Genetic population differentiation within massifs implies severe reduction in female‐mediated gene flow between patches of habitats.     

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.     

Paternity assessment and population subdivision in a natural population of the larger mouse-eared bat Myotis myotis

A hypervariable simple sequence locus and mitochondrial D-loop sequences were used to analyse genetically a natural population of the larger mouse-eared bat Myotis myotis in southern Bavaria. Tests for population subdivision and direct observations suggest that females return to their natal sites, while males disperse. The males present in female nursery colonies are not related to the females. Paternity assessment for 46 offspring from a particular nursery colony showed that there are no males that monopolize the reproduction, and that the resident males in the colony had only a small mating success. Instead, the results suggest that females actively seek matings outside their colony. Most interestingly, it appears that a group of males about 16 km away from the nursery colony had a relatively high mating success and that this group of males may be related to the females of the nursery colonies. If this finding can be confirmed in a larger study, it may have important consequences for future conservation strategies.     

Genetic structure of brown trout (Salmo trutta, L.) populations from south-western France: data from mitochondrial control region variability

The genetic relationships between several French Atlantic populations of brown trout were analysed using polymerase chain reaction-single-strand conformation polymorphism and sequencing of two parts of the mitochondrial control region. Ten mitochondrial haplotypes were identified, separated by a small number of mutations. The distribution of these haplotypes confirmed our hypothesis of the existence of two genetically well-differentiated groups of populations in this area, which has already been suggested using nuclear markers. Nuclear and mitochondrial markers, however, show different patterns of genetic differentiation: some populations which belong to one group according to allozyme results appear to be part of the other when using mitochondrial DNA. The origin of these different forms is discussed; they seem to belong to the Atlantic grouping. Nevertheless, the hypothesis of an ancestral and a modern form suggested for northern populations does not appear valid in this particular case.     

A new perspective on the zoogeography of the sibling mouse-eared bat species Myotis myotis and Myotis blythii: morphological, genetical and ecological evidence

The actual geographic distribution of the two sibling mouse-eared bat species Myotis myotis and Myotis blythii , which occur widely sympatrically in the western Palaearctic region, remains largely controversial. This concerns particularly the specific attribution of marginal populations from the Mediterranean islands and from adjacent areas of North Africa and Asia, which are morphologically intermediate between continental M. myotis and M. blythii from Europe. This study attempts to clarify this question by using four different approaches: cranial morphology, external morphology, genetics and trophic ecology. The three latter methods show unambiguously that North Africa, Malta, Sardinia and Corsica are presently inhabited by monospecific populations of M. myotis . In contrast, cranial morphometrics do not yield conclusive results. These results contradict all recent studies, which attribute North African and Maltese mouse-eared bats to M. blythii and consider that Sardinia and Corsica harbour sympatric populations of the two species. As concerns south-eastem populations, doubts are also expressed about the attribution of the subspecific taxon omari which may actually refer to M. myotis instead of M. blythii. Protein electrophoresis is presently the only absolute method available for determining M. myotis and M. blythii throughout their distribution ranges. However, species identification may be approached by relying on less sophisticated morphometrical methods as presented in this study. Species-specific habitat specializations are probably responsible for the differences observed between the geographic distributions of M. myotis and M. blythii , as they provide a logical groundwork for a coherent model of speciation for these two bat species.     

Riverine effects on mitochondrial structure of Bornean orangutans (Pongo pygmaeus) at two spatial scales

We examined mitochondrial DNA control region sequences of 73 Kinabatangan orang-utans to test the hypothesis that the phylogeographical structure of the Bornean orang-utan is influenced by riverine barriers. The Lower Kinabatangan Wildlife Sanctuary contains one of the most northern populations of orang-utans ( Pongo pygmaeus ) on Borneo and is bisected by the Kinabatangan River, the longest river in Sabah. Orang-utan samples on either side of the river were strongly differentiated with a high Φ_ST value of 0.404 ( P <  0.001). Results also suggest an east–west gradient of genetic diversity and evidence for population expansion along the river, possibly reflecting a postglacial colonization of the Kinabatangan floodplain. We compared our data with previously published sequences of Bornean orang-utans in the context of river catchment structure on the island and evaluated the general relevance of rivers as barriers to gene flow in this long-lived, solitary arboreal ape.     

Genetic consequences of polygyny and social structure in an Indian fruit bat, Cynopterus sphinx. I. Inbreeding, outbreeding, and population subdivision

Population subdivision into behaviorally cohesive kin groups influences rates of inbreeding and genetic drift and has important implications for the evolution of social behavior. Here we report the results of a study designed to test the hypothesis that harem social structure promotes inbreeding and genetic subdivision in a population with overlapping generations. Genetic consequences of harem social structure were investigated in a natural population of a highly polygynous fruit bat, Cynopterus sphinx (Chiroptera: Pteropodidae), in western India. The partitioning of genetic variance within and among breeding groups was assessed using 10-locus microsatellite genotypes for 431 individually marked bats. Genetic analysis of the C. sphinx study population was integrated with field data on demography and social structure to determine the specific ways in which mating, dispersal, and new social group formation influenced population genetic structure. Microsatellite data revealed striking contrasts in genetic structure between consecutive offspring cohorts and between generations. Relative to the 1998 (dry-season) offspring cohort, the 1997 (wet-season) cohort was characterized by a more extensive degree of within-group heterozygote excess (F(IS) = -0.164 vs. -0.050), a greater degree of among-group subdivision (F(ST) = 0.123 vs. 0.008), and higher average within-group relatedness (r = 0.251 vs. 0.017). Differences in genetic structure between the two offspring cohorts were attributable to seasonal differences in the number and proportional representation of male parents. Relative to adult age-classes, offspring cohorts were characterized by more extensive departures from allelic and genotypic equilibria and a greater degree of genetic subdivision. Generational differences in F-statistics indicated that genetic structuring of offspring cohorts was randomized by natal dispersal prior to recruitment into the breeding population. Low relatedness among harem females (r = 0.002-0.005) was primarily attributable to high rates of natal dispersal and low rates of juvenile survivorship. Kin selection is therefore an unlikely explanation for the formation and maintenance of behaviorally cohesive breeding groups in this highly social mammal.     

Biogeography of the Indonesian archipelago: mitochondrial DNA variation in the fruit bat, Eonycteris spelaea

The fruit bat, Eonycteris spelaea , occurs from India through the Philippines to the southeast limit of its distribution in the Lesser Sunda islands of Indonesia. Mitochondrial DNA (mtDNA) variation was examined in Indonesian E. spelaea island populations by amplification of the D-loop and digestion with restriction endonucleases. In addition, microgeographic variation was assessed by investigation of three cave populations within one island. A total of 24 genotypes, comprising two broad clades, was detected. The pattern of mtDNA variation reflects the colonization history of E. spelaea with estimates of haplotype and sequence diversity highest in the older western populations and lowest at the eastern periphery of the species' distribution. These findings may also be associated with an environmental cline from west to east. There is also evidence that genetic distance between populations reflects geographic relationships, especially historical connectedness, as measured by Pleistocene sea-crossing distances. At the microgeographic level, cave populations were heterogeneous and composed of diverse lineages suggesting restricted local interchange.     

Genetic structure among continental and island populations of gyrfalcons

Little is known about the possible influence that past glacial events have had on the phylogeography and population structure of avian predators in the Arctic and sub-Arctic. In this study, we use microsatellite and mitochondrial control region DNA variation to investigate the population genetic structure of gyrfalcons (Falco rusticolus) throughout a large portion of their circumpolar distribution. In most locations sampled, the mtDNA data revealed little geographic structure; however, five out of eight mtDNA haplotypes were unique to a particular geographic area (Greenland, Iceland, or Alaska) and the Iceland population differed from others based on haplotype frequency differences (F(ST)). With the microsatellite results, significant population structure (F(ST), principal components analysis, and cluster analysis) was observed identifying Greenland and Iceland as separate populations, while Norway, Alaska and Canada were identified as a single population consistent with contemporary gene flow across Russia. Within Greenland, differing levels of gene flow between western and eastern sampling locations was indicated with apparent asymmetric dispersal in western Greenland from north to south. This dispersal bias is in agreement with the distribution of plumage colour variants with white gyrfalcons in much higher proportion in northern Greenland. Lastly, because the mtDNA control region sequence differed by only one to four nucleotides from a common haplotype among all gyrfalcons, we infer that the observed microsatellite population genetic structure has developed since the last glacial maximum. This conclusion is further supported by our finding that a closely related species, the saker falcon (Falco cherrug), has greater genetic heterogeneity, including mtDNA haplotypes differing by 1-16 nucleotide substitutions from a common gyrfalcon haplotype. This is consistent with gyrfalcons having expanded rapidly from a single glacial-age refugium to their current circumpolar distribution. Additional sampling of gyrfalcons from Fennoscandia and Russia throughout Siberia is necessary to test putative gene flow between Norway and Alaska and Canada as suggested by this study.     

Genetic variation across the historical range of the wild turkey (Meleagris gallopavo)

Genetic differences within and among naturally occurring populations of wild turkeys (Meleagris gallopavo) were characterized across five subspecies' historical ranges using amplified fragment length polymorphism (AFLP) analysis, microsatellite loci and mitochondrial control region sequencing. Current subspecific designations based on morphological traits were generally supported by these analyses, with the exception of the eastern (M. g. silvestris) and Florida (M. g. osceola) subspecies, which consistently formed a single unit. The Gould's subspecies was both the most genetically divergent and the least genetically diverse of the subspecies. These genetic patterns were consistent with current and historical patterns of habitat continuity. Merriam's populations showed a positive association between genetic and geographical distance, Rio Grande populations showed a weaker association and the eastern populations showed none, suggesting differing demographic forces at work in these subspecies. We recommend managing turkeys to maintain subspecies integrity, while recognizing the importance of maintaining regional population structure that may reflect important adaptive variation.     

Genetic divergences and intraspecific variation in corvids of the genus Corvus (Aves: Passeriformes: Corvidae) – a first survey based on museum specimens

The first comprehensive overview of intra‐ and interspecific variation within the genus Corvus as well as first insights into the phylogenetic relationships of its species is presented. DNA sequences of the mitochondrial control region were obtained from 34 of the 40 described species (including subspecies: 56 taxa). As the study was based mainly on museum material, several specimens did not yield the full length marker sequence. In these cases, only a short section of the control region could be analysed. Nevertheless, even these individuals could be assigned tentatively to clades established on the full length marker sequence. Inclusion of sequences of other corvid genera as available in GenBank clearly confirmed the monophyly of the genus Corvus. Within the Corvus clade several distinct subclades can be distinguished. Some represent lineages of single species or species pairs while other clades are composed of many species. In general, the composition of the clades reflects geographical contiguousness and confirms earlier assumptions of a Palearctic origin of the genus Corvus with several independent colonizations of the Nearctic and the Aethiopis. The Australasian radiation seems to be derived from a single lineage. The distribution of plumage colour in the phylogenetic tree indicates that the pale markings evolved several times independently. The white/grey plumage colour pattern – which is found also in other genera of the family Corvidae, for example, in Pica– occurs already in the species pair representing the first split within the genus Corvus (Corvus monedula, Corvus dauuricus). Thus, reversal to full black colour seems to have occurred as well. The use of colour traits as a phylogenetic marker within Corvus should be considered with severe caution.     

Population genetics of fruit bat reservoir informs the dynamics,distribution and diversity of Nipah virus

The structure and connectivity of wildlife host populations may influence zoonotic disease dynamics, evolution and therefore spillover risk to people. Fruit bats in the genus Pteropus, or flying foxes, are the primary natural reservoir for henipaviruses—a group of emerging paramyxoviruses that threaten livestock and public health. In Bangladesh, Pteropus medius is the reservoir for Nipah virus—and viral spillover has led to human fatalities nearly every year since 2001. Here, we use mitochondrial DNA and nuclear microsatellite markers to measure the population structure, demographic history and phylogeography of P. medius in Bangladesh. We combine this with a phylogeographic analysis of all known Nipah virus sequences and strains currently available to better inform the dynamics, distribution and evolutionary history of Nipah virus. We show that P. medius is primarily panmictic, but combined analysis of microsatellite and morphological data shows evidence for differentiation of two populations in eastern Bangladesh, corresponding to a divergent strain of Nipah virus also found in bats from eastern Bangladesh. Our demographic analyses indicate that a large, expanding population of flying foxes has existed in Bangladesh since the Late Pleistocene, coinciding with human population expansion in South Asia, suggesting repeated historical spillover of Nipah virus likely occurred. We present the first evidence of mitochondrial introgression, or hybridization, between P. medius and flying fox species found in South‐East Asia (P. vampyrus and P. hypomelanus), which may help to explain the distribution of Nipah virus strains across the region.     

Microsatellite and mitochondrial DNA analyses of the genetic structure of blacktip shark (Carcharhinus limbatus) nurseries in the northwestern Atlantic, Gulf of Mexico, and Caribbean Sea

Abstract We investigated the genetic structure of blacktip shark (Carcharhinus limbatus) continental nurseries in the northwestern Atlantic Ocean, Gulf of Mexico, and Caribbean Sea using mitochondrial DNA control region sequences and eight nuclear microsatellite loci scored in neonate and young-of-the-year sharks. Significant structure was detected with both markers among nine nurseries (mitochondrial PhiST = 0.350, P < 0.001; nuclear PhiST = 0.007, P < 0.001) and sharks from the northwestern Atlantic, eastern Gulf of Mexico, western Gulf of Mexico, northern Yucatan, and Belize possessed significantly different mitochondrial DNA haplotype frequencies. Microsatellite differentiation was limited to comparisons involving northern Yucatan and Belize sharks with nuclear genetic homogeneity throughout the eastern Gulf of Mexico, western Gulf of Mexico, and northwestern Atlantic. Differences in the magnitude of maternal vs. biparental genetic differentiation support female philopatry to northwestern Atlantic, Gulf of Mexico, and Caribbean Sea natal nursery regions with higher levels of male-mediated gene flow. Philopatry has produced multiple reproductive stocks of this commercially important shark species throughout the range of this study.