Evolutionary and demographic processes shaping geographic patterns of genetic diversity in a keystone species,the African forest elephant (Loxodonta cyclotis)

The past processes that have shaped geographic patterns of genetic diversity may be difficult to infer from current patterns. However, in species with sex differences in dispersal, differing phylogeographic patterns between mitochondrial (mt) and nuclear (nu) DNA may provide contrasting insights into past events. Forest elephants (Loxodonta cyclotis) were impacted by climate and habitat change during the Pleistocene, which likely shaped phylogeographic patterns in mitochondrial (mt) DNA that have persisted due to limited female dispersal. By contrast, the nuclear (nu) DNA phylogeography of forest elephants in Central Africa has not been determined. We therefore examined the population structure of Central African forest elephants by genotyping 94 individuals from six localities at 21 microsatellite loci. Between forest elephants in western and eastern Congolian forests, there was only modest genetic differentiation, a pattern highly discordant with that of mtDNA. Nuclear genetic patterns are consistent with isolation by distance. Alternatively, male‐mediated gene flow may have reduced the previous regional differentiation in Central Africa suggested by mtDNA patterns, which likely reflect forest fragmentation during the Pleistocene. In species like elephants, male‐mediated gene flow erases the nuclear genetic signatures of past climate and habitat changes, but these continue to persist as patterns in mtDNA because females do not disperse. Conservation implications of these results are discussed.     

SEX-BIASED GENE FLOW IN SPECTACLED EIDERS (ANATIDAE): INFERENCES FROM MOLECULAR MARKERS WITH CONTRASTING MODES OF INHERITANCE

Abstract Genetic markers that differ in mode of inheritance and rate of evolution (a sex‐linked Z‐specific micro‐satellite locus, five biparentally inherited microsatellite loci, and maternally inherited mitochondrial [mtDNA] sequences) were used to evaluate the degree of spatial genetic structuring at macro‐ and microgeographic scales, among breeding regions and local nesting populations within each region, respectively, for a migratory sea duck species, the spectacled eider (Somateria fisheri). Disjunct and declining breeding populations coupled with sex‐specific differences in seasonal migratory patterns and life history provide a series of hypotheses regarding rates and directionality of gene flow among breeding populations from the Indigirka River Delta, Russia, and the North Slope and Yukon‐Kuskokwim Delta, Alaska. The degree of differentiation in mtDNA haplotype frequency among breeding regions and populations within regions was high (φ_CT= 0.189, P < 0.01; φ_SC= 0.059, P < 0.01, respectively). Eleven of 17 mtDNA haplotypes were restricted to a single breeding region. Genetic differences among regions were considerably lower for nuclear DNA loci (sex‐linked: φ_ST= 0.001, P > 0.05; biparentally inherited microsatellites: mean θ= 0.001, P > 0.05) than was observed for mtDNA. Using models explicitly designed for uniparental and biparentally inherited genes, estimates of spatial divergence based on nuclear and mtDNA data together with elements of the species' breeding ecology were used to estimate effective population size and degree of male and female gene flow. Differences in the magnitude and spatial patterns of gene correlations for maternally inherited and nuclear genes revealed that females exhibit greater natal philopatry than do males. Estimates of generational female and male rates of gene flow among breeding regions differed markedly (3.67 × 10^‐4 and 1.28 × 10^‐2, respectively). Effective population size for mtDNA was estimated to be at least three times lower than that for biparental genes (30,671 and 101,528, respectively). Large disparities in population sizes among breeding areas greatly reduces the proportion of total genetic variance captured by dispersal, which may accelerate rates of inbreeding (i.e., promote higher coancestries) within populations due to nonrandom pairing of males with females from the same breeding population.     

Spatial heterogeneity of mitochondrial DNA and allozymes among populations of white-tailed deer and mule deer.

A white-tailed deer (Odocoileus virginianus) population in northeastern Minnesota and a mule deer (O. hemionus) population in the Bridger Mountains Montana, have previously been shown to be spatially subdivided into contiguous subpopulations. We assessed the degree of genetic differentiation among subpopulations and tested the hypothesis that differentiation will be greater for mitochondrial DNA (mtDNA) than for nuclear-encoded allozymes. Differentiation of the white-tailed deer subpopulations was significant for two allozyme loci but not for mtDNA, and the overall degree of differentiation was low. Gene flow, recent founding of the subpopulations, and polygamous breeding structure may all have contributed to this pattern. Greater differentiation was evident among disjunct populations than between the contiguous subpopulations of white-tailed deer. The contiguous mule deer subpopulations were significantly differentiated for mtDNA and one allozyme locus. Differentiation was greater for mtDNA than for allozymes. These results are consistent with demographic data that indicate mule deer males disperse more than do females. Disjunct mule deer populations may be similar or dramatically different in mtDNA haplotype frequencies that do not necessarily vary with geographic distance. Current and historical gene flow and breeding structure will influence population genetic patterns.     

Genetic population structure and introgression in Anopheles dirus mosquitoes in South-east Asia

Genetic structure and species relationships were studied in three closely related mosquito species, Anopheles dirus A, C and D in Thailand using 11 microsatellite loci and compared with previous mitochondrial DNA (mtDNA) data on the same populations. All three species were well differentiated from each other at the microsatellite loci. Given the almost complete absence of mtDNA differentiation between An. dirus A and D, this endorses the previous suggestion of mtDNA introgression between these species. The high degree of differentiation between the northern and southern population of An. dirus C (RST = 0.401), in agreement with mtDNA data, is suggestive of incipient species. The lack of genetic structure indicated by microsatellites in four populations of An. dirus A across northern Thailand also concurs with mtDNA data. However, in An. dirus D a limited but significant level of structure was detected by microsatellites over ~400 km in northern Thailand, whereas the mtDNA detected no population differentiation over a much larger area (>1200 km). There is prior evidence for population expansion in the mtDNA. If this is due to a selective sweep originating in An. dirus D, the microsatellite data may indicate greater barriers to gene flow within An. dirus D than in species A. Alternatively, there may have been historical introgression of mtDNA and subsequent demographic expansion which occurred first in An. dirus D so enabling it to accumulate some population differentiation. In the latter case the lack of migration-drift equilibrium precludes the inference of absolute or relative values of gene flow in An. dirus A and D.     

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

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-by-distance 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.     

Male-biased dispersal promotes large scale gene flow in a subterranean army ant, Dorylus (Typhlopone) fulvus

Sex-biased dispersal is a widespread phenomenon in the animal kingdom, which strongly influences gene flow and population structure. Particularly army ants, important key-stone predators in tropical ecosystems, are prone to population fragmentation and isolation due to their extraordinary mating system: queens are permanently wingless, propagate via colony fission, and only the males disperse in mating flights. Here we report on sex-biased dispersal and the genetic population structure of an African subterranean army ant, Dorylus (Typhlopone) fulvus. Using maternally inherited mtDNA markers and bi-parentally inherited nuclear microsatellites we found strong geographical structuring of mtDNA haplotypes, whereas the nuclear genetic population structure was less pronounced. Strong mtDNA (Φ _ST = 0.85), but significantly lower nuclear (F _ST = 0.23) genetic differentiation translated into a more than an order of magnitude larger male migration rate compared to that of queens, reflecting the low motility of queens and strong, promiscuous dispersal by males. Thus, the well flying D. fulvus males appear to be the sex to promote large scale gene flow, and D. fulvus is indeed a species in which sex specific dispersal patterns and the mating system profoundly affect the population structure and phylogeography.     

Male‐biased dispersal and the potential impact of human‐induced habitat modifications on the Neotropical bat Trachops cirrhosus

Gene flow, maintained through natal dispersal and subsequent mating events, is one of the most important processes in both ecology and population genetics. Among mammalian populations, gene flow is strongly affected by a variety of factors, including the species’ ability to disperse, and the composition of the environment which can limit dispersal. Information on dispersal patterns is thus crucial both for conservation management and for understanding the social system of a species. We used 16 polymorphic nuclear microsatellite loci in addition to mitochondrial DNA sequences (1.61 kbp) to analyse the population structure and the sex‐specific pattern of natal dispersal in the frog‐eating fringe‐lipped bat, Trachops cirrhosus, in Central Panama. Our study revealed that—unlike most of the few other investigated Neotropical bats—gene flow in this species is mostly male‐mediated. Nevertheless, distinct genetic clusters occur in both sexes. In particular, the presence of genetic differentiation in the dataset only consisting of the dispersing sex (males) indicates that gene flow is impeded within our study area. Our data are in line with the Panama Canal in connection with the widening of the Río Chagres during the canal construction acting as a recent barrier to gene flow. The sensitivity of T. cirrhosus to human‐induced habitat modifications is further indicated by an extremely low capture success in highly fragmented areas. Taken together, our genetic and capture data provide evidence for this species to be classified as less mobile and thus vulnerable to habitat change, information that is important for conservation management.     

High gene flow on a continental scale in the polyandrous Kentish plover Charadrius alexandrinus

Gene flow promotes genetic homogeneity of species in time and space. Gene flow can be modulated by sex‐biased dispersal that links population genetics to mating systems. We investigated the phylogeography of the widely distributed Kentish plover Charadrius alexandrinus. This small shorebird has a large breeding range spanning from Western Europe to Japan and exhibits an unusually flexible mating system with high female breeding dispersal. We analysed genetic structure and gene flow using a 427‐bp fragment of the mitochondrial (mtDNA) control region, 21 autosomal microsatellite markers and a Z microsatellite marker in 397 unrelated individuals from 21 locations. We found no structure or isolation‐by‐distance over the continental range. However, island populations had low genetic diversity and were moderately differentiated from mainland locations. Genetic differentiation based on autosomal markers was positively correlated with distance between mainland and each island. Comparisons of uniparentally and biparentally inherited markers were consistent with female‐biased gene flow. Maternally inherited mtDNA was less structured, whereas the Z‐chromosomal marker was more structured than autosomal microsatellites. Adult males were more related than females within genetic clusters. Taken together, our results suggest a prominent role for polyandrous females in maintaining genetic homogeneity across large geographic distances.     

Genetic population structure of the paper wasp Polistes olivaceus (Hymenoptera: Vespidae) in Bangladesh

Dispersal triggers gene flow, which in turn strongly affects the ensuing genetic population structure of a species. Using nuclear microsatellite loci and mitochondrial DNA (mtDNA), we estimated the genetic population structure of the wasp Polistes olivaceus throughout Bangladesh. The level of population differentiation using nuclear markers (F _ST) appeared to be much lower than that estimated using mtDNA haplotype sequences (Ф_ST), even after correcting for effective population size differences between the two markers. These results suggest a philopatric tendency, in which gynes disperse less than males. We observed no isolation by distance among the study populations at either the nuclear or mtDNA level, suggesting nonequilibrium between gene flow and drift as a result of very frequent interpopulation movement. For the nuclear markers, an individual assignment test showed no genetically and geographically distinct groups. Instead, phylogenetic analyses as well as a minimum spanning network using mtDNA haplotypes consistently revealed two distinct lineages. The distribution of haplotypes indicated western populations with a single lineage and offered clear evidence for restricted gene flow across the Jamuna–Padma–Upper Meghna river system. Mismatch distributions exhibited a unimodal distribution, which along with a starlike haplotype network, suggested a population expansion in lineage I but not in lineage II. Overall, these results suggest that gene flow among populations of P. olivaceus was affected by both female philopatry and a major river system across Bangladesh.     

Ecology,environment and evolutionary history influence genetic structure in five mammal species from the Italian Alps

The identification and evaluation of the ecological and environmental factors shaping patterns of natural genetic variation are fundamental goals of population and conservation genetics. Many studies focus on factors affecting single species, but it is also important to test whether some influential biotic and abiotic factors are common drivers of genetic diversity across species, or if species or species groups are each affected by different forces; a multi‐species analysis is necessary for this. Here we analysed the molecular variation from five mammal species (roe deer, red deer, chamois, mountain hare and European brown hare) at mtDNA and microsatellite loci from the eastern Italian Alps. We use phylogeographical and landscape‐level analyses to test the relative influence of large‐scale geographical history and contemporary environmental characteristics of the landscape on genetic diversity and differentiation. We found: (1) all study species except brown hare are strongly differentiated into two main groups, located west and east of a major river valley; (2) significant correlations between levels of within‐population diversity at both mtDNA and microsatellite loci, and several landscape features such as alpine grassland, water courses and anthropized areas. We conclude that heterogeneous landscape has some influence on within‐population diversity, but biogeographical history has probably had the stronger influence on current genetic patterns, despite an apparently large dispersal potential of certain species. However, our results for brown hare show that management actions such as stocking may alter these large‐scale patterns.     

Species‐level view of population structure and gene flow for a critically endangered primate (Varecia variegata)

Lemurs are among the world's most threatened mammals. The critically endangered black‐and‐white ruffed lemur (Varecia variegata), in particular, has recently experienced rapid population declines due to habitat loss, ecological sensitivities to habitat degradation, and extensive human hunting pressure. Despite this, a recent study indicates that ruffed lemurs retain among the highest levels of genetic diversity for primates. Identifying how this diversity is apportioned and whether gene flow is maintained among remnant populations will help to diagnose and target conservation priorities. We sampled 209 individuals from 19 sites throughout the remaining V. variegata range. We used 10 polymorphic microsatellite loci and ~550 bp of mtDNA sequence data to evaluate genetic structure and population dynamics, including dispersal patterns and recent population declines. Bayesian cluster analyses identified two distinct genetic clusters, which optimally partitioned data into populations occurring on either side of the Mangoro River. Localities north of the Mangoro were characterized by greater genetic diversity, greater gene flow (lower genetic differentiation) and higher mtDNA haplotype and nucleotide diversity than those in the south. Despite this, genetic differentiation across all sites was high, as indicated by high average F_ST (0.247) and Φ_ST (0.544), and followed a pattern of isolation‐by‐distance. We use these results to suggest future conservation strategies that include an effort to maintain genetic diversity in the north and restore connectivity in the south. We also note the discordance between patterns of genetic differentiation and current subspecies taxonomy, and encourage a re‐evaluation of conservation management units moving forward.     

Genetic structuring and recent demographic history of red pandas (Ailurus fulgens) inferred from microsatellite and mitochondrial DNA

Clarification of the genetic structure and population history of a species can shed light on the impacts of landscapes, historical climate change and contemporary human activities and thus enables evidence‐based conservation decisions for endangered organisms. The red panda (Ailurus fulgens) is an endangered species distributing at the edge of the Qinghai‐Tibetan Plateau and is currently subject to habitat loss, fragmentation and population decline, thus representing a good model to test the influences of the above‐mentioned factors on a plateau edge species. We combined nine microsatellite loci and 551 bp of mitochondrial control region (mtDNA CR) to explore the genetic structure and demographic history of this species. A total of 123 individuals were sampled from 23 locations across five populations. High levels of genetic variation were identified for both mtDNA and microsatellites. Phylogeographic analyses indicated little geographic structure, suggesting historically wide gene flow. However, microsatellite‐based Bayesian clustering clearly identified three groups (Qionglai‐Liangshan, Xiaoxiangling and Gaoligong‐Tibet). A significant isolation‐by‐distance pattern was detected only after removing Xiaoxiangling. For mtDNA data, there was no statistical support for a historical population expansion or contraction for the whole sample or any population except Xiaoxiangling where a signal of contraction was detected. However, Bayesian simulations of population history using microsatellite data did pinpoint population declines for Qionglai, Xiaoxiangling and Gaoligong, demonstrating significant influences of human activity on demography. The unique history of the Xiaoxiangling population plays a critical role in shaping the genetic structure of this species, and large‐scale habitat loss and fragmentation is hampering gene flow among populations. The implications of our findings for the biogeography of the Qinghai‐Tibetan Plateau, subspecies classification and conservation of red pandas are discussed.     

PREDICTING NUCLEAR GENE COALESCENCE FROM MITOCHONDRIAL DATA: THE THREE-TIMES RULE

Abstract.— Coalescence theory predicts when genetic drift at nuclear loci will result in fixation of sequence differences to produce monophyletic gene trees. However, the theory is difficult to apply to particular taxa because it hinges on genetically effective population size, which is generally unknown. Neutral theory also predicts that evolution of monophyly will be four times slower in nuclear than in mitochondrial genes primarily because genetic drift is slower at nuclear loci. Variation in mitochondrial DNA (mtDNA) within and between species has been studied extensively, but can these mtDNA data be used to predict coalescence in nuclear loci? Comparison of neutral theories of coalescence of mitochondrial and nuclear loci suggests a simple rule of thumb. The “three‐times rule” states that, on average, most nuclear loci will be monophyletic when the branch length leading to the mtDNA sequences of a species is three times longer than the average mtDNA sequence diversity observed within that species. A test using mitochondrial and nuclear intron data from seven species of whales and dolphins suggests general agreement with predictions of the three‐times rule. We define the coalescence ratio as the mitochondrial branch length for a species divided by intraspecific mtDNA diversity. We show that species with high coalescence ratios show nuclear monophyly, whereas species with low ratios have polyphyletic nuclear gene trees. As expected, species with intermediate coalescence ratios show a variety of patterns. Especially at very high or low coalescence ratios, the three‐times rule predicts nuclear gene patterns that can help detect the action of selection. The three‐times rule may be useful as an empirical benchmark for evaluating evolutionary processes occurring at multiple loci.     

Evaluating summary statistics used to test for incomplete lineage sorting: mito‐nuclear discordance in the reef sponge Callyspongia vaginalis

Conflicting patterns of population differentiation between the mitochondrial and nuclear genomes (mito‐nuclear discordance) have become increasingly evident as multilocus data sets have become easier to generate. Incomplete lineage sorting (ILS) of nucDNA is often implicated as the cause of such discordance, stemming from the large effective population size of nucDNA relative to mtDNA. However, selection, sex‐biased dispersal and historical demography can also lead to mito‐nuclear discordance. Here, we compare patterns of genetic diversity and subdivision for six nuclear protein‐coding gene regions to those for mtDNA in a common Caribbean coral reef sponge, Callyspongia vaginalis, along the Florida reef tract. We also evaluated a suite of summary statistics to determine which are effective metrics for comparing empirical and simulated data when testing drivers of mito‐nuclear discordance in a statistical framework. While earlier work revealed three divergent and geographically subdivided mtDNACOI haplotypes separated by 2.4% sequence divergence, nuclear alleles were admixed with respect to mitochondrial clade and geography. Bayesian analysis showed that substitution rates for the nuclear loci were up to 7 times faster than for mitochondrial COI. Coalescent simulations and neutrality tests suggested that mito‐nuclear discordance in C. vaginalis is not the result of ILS in the nucDNA or selection on the mtDNA but is more likely caused by changes in population size. Sperm‐mediated gene flow may also influence patterns of population subdivision in the nucDNA.     

Spatial structure of lemming populations (Dicrostonyx groenlandicus) fluctuating in density

The pattern and scale of the genetic structure of populations provides valuable information for the understanding of the spatial ecology of populations, including the spatial aspects of density fluctuations. In the present paper, the genetic structure of periodically fluctuating lemmings (Dicrostonyx groenlandicus) in the Canadian Arctic was analysed using mitochondrial DNA (mtDNA) control region sequences and four nuclear microsatellite loci. Low genetic variability was found in mtDNA, while microsatellite loci were highly variable in all localities, including localities on isolated small islands. For both genetic markers the genetic differentiation was clear among geographical regions but weaker among localities within regions. Such a pattern implies gene flow within regions. Based on theoretical calculations and population census data from a snap-trapping survey, we argue that the observed genetic variability on small islands and the low level of differentiation among these islands cannot be explained without invoking long distance dispersal of lemmings over the sea ice. Such dispersal is unlikely to occur only during population density peaks.     

Population genetic structure and male-biased dispersal in the queenless ant Diacamma cyaneiventre

In this study we investigated the population genetic structure of the queenless ant Diacamma cyaneiventre. This species, lacking winged queens, is likely to have a restricted female dispersal. We used both mitochondrial and microsatellite markers to assess the consequence of such restricted female dispersal at three geographical scales: within a given locality (< 1 km), between localities within a given region (< 10 km) and between regions (> 36 km). Within a locality, a strong population structure was observed for mitochondrial DNA (mtDNA) whereas weak or nonexistent population genetic structure was observed for the microsatellites (around 5% of the value for mtDNA). Male gene flow was estimated to be about 20-30 times higher than female gene flow at this scale. At a larger spatial scale, very strong genetic differentiation for both markers was observed between localities - even within a single region. Female dispersal is nonexistent at these scales and male dispersal is very restricted, especially between regions. The phylogeographical structure of the mtDNA haplotypes as well as the very low genetic diversity of mtDNA within localities indicate that new sites are colonized by a single migration event from adjacent localities, followed by successive colony fissions. These patterns of genetic variability and differentiation agree with what is theoretically expected when colonization events are kin-structured and when, following colonization, dispersion is mainly performed by males.     

Marked genetic structuring and extreme dispersal limitation in the Pyrenean brook newt Calotriton asper (Amphibia: Salamandridae) revealed by genome-wide AFLP but not mtDNA

Direct estimation of dispersal rates at large geographic scales can be technically and logistically challenging, especially in small animals of low vagility like amphibians. The use of molecular markers to reveal patterns of genetic structure provides an indirect way to infer dispersal rates and patterns of recent and historical gene flow among populations. Here, we use mitochondrial DNA (mtDNA) sequence data and genome-wide amplified fragment length polymorphism markers to examine population structure in the Pyrenean brook newt ( Calotriton asper ) across four main drainages in the French Pyrenees. mtDNA sequence data (2040 bp) revealed three phylogroups shallowly differentiated and with low genetic diversity. In sharp contrast, variation in 382 amplified fragment length polymorphism loci was high and revealed a clear pattern of isolation by distance consistent with long-term restriction of gene flow at three spatial scales: (i) among all four main drainages, (ii) between sites within drainages, and (iii) even between adjacent populations separated by less than 4 km. The high pairwise F _ST values between localities across numerous loci, together with the high frequency of fixed alleles in several populations, suggests a combination of marked geographic isolation, small population sizes and very limited dispersal in C. asper . The contrasting lack of variation detected in mtDNA sequence data is intriguing and underscores the importance of multilocus approaches to detect true patterns of gene flow in natural populations of amphibians.     

Selection and sex‐biased dispersal in a coastal shark: the influence of philopatry on adaptive variation

Sex‐biased dispersal is expected to homogenize nuclear genetic variation relative to variation in genetic material inherited through the philopatric sex. When site fidelity occurs across a heterogeneous environment, local selective regimes may alter this pattern. We assessed spatial patterns of variation in nuclear‐encoded, single nucleotide polymorphisms (SNPs) and sequences of the mitochondrial control region in bonnethead sharks (Sphyrna tiburo), a species thought to exhibit female philopatry, collected from summer habitats used for gestation. Geographic patterns of mtDNA haplotypes and putatively neutral SNPs confirmed female philopatry and male‐mediated gene flow along the northeastern coast of the Gulf of Mexico. A total of 30 outlier SNP loci were identified; alleles at over half of these loci exhibited signatures of latitude‐associated selection. Our results indicate that in species with sex‐biased dispersal, philopatry can facilitate sorting of locally adaptive variation, with the dispersing sex facilitating movement of potentially adaptive variation among locations and environments.     

Extensive phenotypic diversification coexists with little genetic divergence and a lack of population structure in the White Wagtail subspecies complex (Motacilla alba)

Geographically clustered phenotypes often demonstrate consistent patterns in molecular markers, particularly mitochondrial DNA (mtDNA) traditionally used in phylogeographic studies. However, distinct evolutionary trajectories among traits and markers can lead to their discordance. First, geographic structure in phenotypic traits and nuclear molecular markers can be co‐aligned but inconsistent with mtDNA (mito‐nuclear discordance). Alternatively, phenotypic variation can have little to do with patterns in neither mtDNA nor nuclear markers. Disentangling between these distinct patterns can provide insight into the role of selection, demography and gene flow in population divergence. Here, we examined a previously reported case of strong inconsistency between geographic structure in mtDNA and plumage traits in a widespread polytypic bird species, the White Wagtail (Motacilla alba). We tested whether this pattern is due to mito‐nuclear discordance or discrepancy between morphological evolution and both nuclear and mtDNA markers. We analysed population differentiation and structure across six out of nine commonly recognized subspecies using 17 microsatellite loci and a combination of microsatellites and plumage indices in a comprehensively sampled region of a contact between two subspecies. We did not find support for the mito‐nuclear discordance hypothesis: nuclear markers indicated a subtle signal of genetic clustering only partially consistent with plumage groups, similar to previous findings that relied on mtDNA. We discuss evolutionary factors that could have shaped the intricate patterns of phenotypic diversification in the White wagtail and the role that repeated selection on plumage ‘hotspots’ and hybridization may have played.     

Phylogeography and ecological niche modeling unravel the evolutionary history of the African green toad,Bufotes boulengeri boulengeri (Amphibia: Bufonidae), through the Quaternary

Recent integration of ecological niche models in phylogeographic studies is improving our understanding of the processes structuring genetic variation across landscapes. Previous studies on the amphibian Bufotes boulengeri boulengeri uncovered a surprisingly weak intraspecific differentiation across the Maghreb region. We widely sampled this species from Morocco to Egypt and sequenced one nuclear and three mitochondrial (mtDNA) genes to determine the level of genetic variability across its geographic range. We evaluated these data with ecological niche modeling to reveal its evolutionary history in response to climate change during the Quaternary. Our results highlight some mtDNA phylogeographic structure within this species, with one haplogroup endemic to coastal Morocco, and one haplogroup widely distributed throughout North Africa. No or little genetic differentiation is observed between isolated populations from the Hoggar Mountains, the Sabha district and the islands of Kerkennah and Lampedusa, compared to others populations. This can be explained by the expansion of the distribution range of B. b. boulengeri during glacial periods. This might have facilitated the species’ dispersal and subsequent gene flow between most North African localities.