Deciphering patterns of transoceanic dispersal: the evolutionary origin and biogeography of coastal lizards (Cryptoblepharus) in the Western Indian Ocean region

Aim Cryptoblepharus is a genus of small arboreal or rock‐dwelling scincid lizards, widespread through the Indo‐Pacific and Australian regions, with a disjunct outlier in the Malagasy region. The taxonomy within this genus is controversial, with different authors ranking the different forms (now some 36) at various levels, from different species to subspecies of a single species, Cryptoblepharus boutonii. We investigated the biogeography and genetic differentiation of the Cryptoblepharus from the Western Indian Ocean region, in order to understand their origin and history. Location Western Indian Ocean region. Methods We analysed sequences of mitochondrial DNA (partial 12s and 16s rRNA genes, 766 bp) from 48 specimens collected in Madagascar, Mauritius, the four Comoros islands and East Africa, and also in New Caledonia, representing the Australo‐Pacific unit of the distribution. Results Pairwise sequence divergences of c. 3.1% were found between the New Caledonian forms and the ones from the Western Indian Ocean. Two clades were identified in Madagascar, probably corresponding to the recognized forms cognatus and voeltzkowi, and two clades were identified in the Comoro islands, where each island population formed a distinct haplotype clade. The East African samples form a monophyletic unit, with some variation existing between Pemba, Zanzibar and continental Tanzania populations. Individuals from Mauritius form a divergent group, more related to populations from Moheli and Grand Comore (Comoros islands) than to the others. Main conclusions The level of divergence between the populations from the Western Indian Ocean and Australian regions and the geographic coherence of the variation within the Western Indian Ocean group are concordant with the hypothesis of a colonization of this region by a natural transoceanic dispersal (from Australia or Indonesia). The group then may have diversified in Madagascar, from where it separately colonized the East African coast, the Comoros islands (twice), and Mauritius. The genetic divergence found is congruent with the known morphological variation, but its degree is much lower than typically seen between distinct species of reptiles.     

Northern Hemisphere origin,transoceanic dispersal,and diversification of Ranunculeae DC. (Ranunculaceae) in the Cenozoic

Aim The role of dispersal versus vicariance for plant distribution patterns has long been disputed. We study the temporal and spatial diversification of Ranunculeae, an almost cosmopolitan tribe comprising 19 genera, to understand the processes that have resulted in the present inter‐continental disjunctions. Location All continents (except Antarctica). Methods Based on phylogenetic analyses of nuclear and chloroplast DNA sequences for 18 genera and 89 species, we develop a temporal–spatial framework for the reconstruction of the biogeographical history of Ranunculeae. To estimate divergence dates, Bayesian uncorrelated rates analyses and four calibration points derived from geological, fossil and external molecular information were applied. Parsimony‐based methods for dispersal–vicariance analysis (diva and Mesquite ) and a maximum likelihood‐based method (Lagrange ) were used for reconstructing ancestral areas. Six areas corresponding to continents were delimited. Results The reconstruction of ancestral areas is congruent in the diva and maximum likelihood‐based analyses for most nodes, but Mesquite reveals equivocal results at deep nodes. Our study suggests a Northern Hemisphere origin for the Ranunculeae in the Eocene and a weakly supported vicariance event between North America and Eurasia. The Eurasian clade diversified between the early Oligocene and the late Miocene, with at least three independent migrations to the Southern Hemisphere. The North American clade diversified in the Miocene and dispersed later to Eurasia, South America and Africa. Main conclusions Ranunculeae diversified between the late Eocene and the late Miocene. During this time period, the main oceanic barriers already existed between continents and thus dispersal is the most likely explanation for the current distribution of the tribe. In the Southern Hemisphere, a vicariance model related to the break‐up of Gondwana is clearly rejected. Dispersals between continents could have occurred via migration over land bridges, such as the Bering Land Bridge, or via long‐distance dispersal.     

How marine currents influenced the widespread natural overseas dispersal of reptiles in the Western Indian Ocean region

In a recent contribution to this journal, Wilmé et al. (2016) proposed that the giant tortoises of the islands of the Western Indian Ocean (WIO: Aldabra, the Mascarenes, and the Granitic Seychelles) might have originated from translocation by early Austronesian sailors. Prompted by this paper we review recent literature and show that natural overseas dispersal was remarkably widespread in the colonization history of terrestrial reptiles in the WIO region. Almost 90% of the successful colonization events are supported by prevailing marine surface currents. However, these currents may change over geological (and evolutionary) time‐scales, and eddies and counter‐currents may facilitate transport against the main current direction. We review the cases of the extant and extinct WIO giant tortoises and suggest that the current distribution of all lineages can be convincingly explained by overseas dispersal.     

Freshwater paths across the ocean: molecular phylogeny of the frog Ptychadena newtoni gives insights into amphibian colonization of oceanic islands

Aim Amphibians are a model group for studies of the biogeographical origins of salt‐intolerant taxa on oceanic islands. We used the Gulf of Guinea islands to explore the biogeographical origins of island endemism of one species of frog, and used this to gain insights into potential colonization mechanisms. Location São Tomé and Príncipe, two of the four major islands in the Gulf of Guinea, West Africa, are truly oceanic and have an exceptionally high biodiversity. Methods Mitochondrial DNA is used to test the endemic status of a frog from São Tomé and compare it with congeneric taxa from tropical Africa. Existing data on surface currents, surface salinity, atmospheric circulation and bird migration in the Gulf of Guinea are summarized to address hypotheses concerning colonization mechanisms. Results The endemic status of Ptychadena newtoni (Bocage) is supported here by mitochondrial DNA sequences, and analysis of this and other molecular data indicates that an East African species close to Ptychadena mascareniensis (Duméril and Bibron) is its nearest relative. We refute the possibility that this population was anthropogenically introduced, in favour of a natural dispersal mechanism. Main conclusions With six endemic frogs and one caecilian, the Gulf of Guinea islands harbour a diverse amphibian fauna. Five of these species appear to have their closest relatives in East Africa. Insufficient evidence exists for transportation by storms, birds or rafts alone. However, we propose a synergy of rafting, favourable surface currents and a reduction in salinity of surface waters. Catastrophic events, or wet periods in climatic history, could allow freshwater paths to open far enough to enable continental flora and fauna to reach these and other isolated oceanic islands.     

Biogeographical history of cuckoo‐shrikes (Aves: Passeriformes): transoceanic colonization of Africa from Australo‐Papua

Aim Cuckoo‐shrikes and allies (Campephagidae) form a radiation of birds widely distributed in the Indo‐Pacific and Africa. Recent studies on the group have been hampered by poor taxon sampling, causing inferences about systematics and biogeography to be rather speculative. With improved taxon sampling and analyses within an explicit spatiotemporal framework, we elucidate biogeographical patterns of dispersal and diversification within this diverse clade of passerine birds. Location Africa, Asia, Australo‐Papua, the Pacific, the Philippines and Wallacea. Methods We use model‐based phylogenetic methods (Mr Bayes and garli ) to construct a phylogenetic hypothesis of the core Campephagidae (Campephagidae with the exclusion of Pericrocotus). The phylogeny is used to assess the biogeographical history of the group with a newly developed Bayesian approach to dispersal–vicariance analysis (Bayes‐diva) . We also made use of a partitioned beast analysis, with several calibration points taken from island ages, passerine mitochondrial substitution rates and secondary calibration points for passerine birds, to assess the timing of diversification and dispersal. Results We present a robust molecular phylogeny that includes all genera and 84% of the species within the core Campephagidae. Furthermore, we estimate divergence dates and ancestral area relationships. We demonstrate that Campephagidae originated in Australo‐Papua with a single lineage (Pericrocotus) dispersing to Asia early. Later, there was further extensive transoceanic dispersal from Australo‐Papua to Africa involving lineages within the core Campephagidae radiation. Main conclusions The phylogenetic relationships, along with the results of the ancestral area analysis and the timing of dispersal events, support a transoceanic dispersal scenario from Australo‐Papua to Africa by the core Campephagidae. The sister group to core Campephagidae, Pericrocotus, dispersed to mainland Asia in the late Oligocene. Asia remained uncolonized by the core Campephagidae until the Pliocene. Transoceanic dispersal is by no means an unknown phenomenon, but our results represent a convincing case of colonization over a significant water gap of thousands of kilometres from Australo‐Papua to Africa.     

Multiple colonization and dispersal events hide the early origin and induce a lack of genetic structure of the moss Bryum argenteum in Antarctica

The dispersal routes of taxa with transoceanic disjunctions remain poorly understood, with the potential roles of Antarctica not yet demonstrated. Mosses are suitable organisms to test direct intra‐Antarctic dispersal, as major component of the extant Antarctic flora, with the cosmopolitan moss Bryum argenteum as ideal target species. We analyzed the genetic structure of B. argenteum to provide an evolutionary time frame for its radiation and shed light into its historical biogeography in the Antarctic region. We tested two alternative scenarios: (a) intra‐Antarctic panmixia and (b) intra‐Antarctic genetic differentiation. Furthermore, we tested for evidence of the existence of specific intra‐Antarctic dispersal routes. Sixty‐seven new samples (40 collected in Antarctica) were sequenced for ITS nrDNA and rps4 cpDNA regions, and phylogenetic trees of B. argenteum were constructed, with a focus on its Southern Hemisphere. Combining our new nrDNA dataset with previously published datasets, we estimated time‐calibrated phylogenies based on two different substitution rates (derived from angiosperms and bryophytes) along with ancestral area estimations. Minimum spanning network and pairwise genetic distances were also calculated. B. argenteum was potentially distributed across Africa and Antarctica soon after its origin. Its earliest intra‐Antarctic dispersal and diversification occurred during a warming period in the Pliocene. On the same timescale, a radiation took place involving a dispersal event from Antarctica to the sub‐Antarctic islands. A more recent event of dispersal and diversification within Antarctica occurred during a warm period in the Pleistocene, creating favorable conditions also for its colonization outside the Antarctic continent worldwide. We provide evidence supporting the hypothesis that contemporary populations of B. argenteum in Antarctica integrate a history of both multiple long‐range dispersal events and local persistence combined with in situ diversification. Our data support the hypothesis that B. argenteum has been characterized by strong connectivity within Antarctica, suggesting the existence of intra‐Antarctic dispersal routes.     

Phylogeny of Gaertnera Lam. (Rubiaceae) based on multiple DNA markers: evidence of a rapid radiation in a widespread, morphologically diverse genus

Phylogenetic relationships among 28 of the 68 species of the paleotropical genus Gaertnera (Rubiaceae) and two related genera were inferred from nucleotide sequence variation in four nuclear DNA (nDNA) markers: the internal transcribed spacers of nuclear rDNA (ITS), the large and small copies of phosphoenolpyruvate carboxylase (PepC-large and PepC-small), and triose phosphate isomerase (Tpi). Phylogenetic analysis of the combined nDNA dataset suggested that Gaertnera is monophyletic, but genetic variation among species was insufficient to reconstruct well-supported relationships within the genus. This was counter to expectations based on the very distinct morphologies and widespread distribution of the genus (West Africa to Sulawesi). Molecular clock analyses suggested variable dates of origin for Gaertnera depending upon the calibration method used. The most plausible calibration implies that Gaertnera migrated to Africa during the early Tertiary, possibly via a boreotropical land bridge and suggests that Gaertnera started to radiate 5.21 +/- 0.14 million years ago. This implies that range expansion in the group has occurred via a number of long-distance dispersal events rather than vicariance. The molecular clock estimate in turn estimated an unusually rapid lineage diversification rate within the radiation of 0.717-0.832 species/million years, comparable to those estimated for radiations on oceanic islands. Although low interspecific competition levels may have contributed to the diversification of Gaertnera on Mauritius, the mechanisms driving the rapid radiation of the group in other parts of its range remain elusive.     

Comparative biology of pollination systems in the African-Malagasy genus Brownleea (Brownleeinae: Orchidaceae)

The role of pollinators in the evolution of the African-Malagasy orchid genus Brownleea (Brownleeinae) was investigated. Taxa show specialization for pollination by bees ( B. parviflora , B. recurvata ), short-proboscid flies ( B. galpinii ), and long-proboscid flies ( B. macroceras , B. coerulea ). All species in the genus produce nectar, but some ( B. coerulea , B. galpinii ) appear to mimic flowers of abundant sympatric species as an additional strategy to attract pollinators. Species investigated in terms of their breeding systems ( B. coerulea , B. parviflora , B. macroceras , B. galpinii ) are reliant on pollinator visits for seed production. Self-pollination results in strong inbreeding depression during embryo formation. A phylogeny of the genus, constructed using data from the ITS1-5.8S-ITS2 region and morphology, indicates that fly pollination is likely to be basal in the genus, and that there has been a single shift to bee pollination.  © 2008 The Linnean Society of London, Botanical Journal of the Linnean Society , 2008, 156 , 65–78.     

Phylogeography and postglacial dispersal of arctic charr Salvelinus alpinus in North America

Restriction fragment length polymorphism analysis of mitochondrial DNA (mtDNA) was used to reconstruct postglacial dispersal routes of arctic charr Salvelinus alpinus in North America. Twelve of 35 restriction enzymes detected polymorphisms among representative populations, revealing two distinct lineages with an estimated nucleotide divergence of 1.32%. Subsequent screening of 869 fish from 54 populations with four diagnostic restriction enzymes showed that these lineages have largely allopatric distributions, suggesting their dispersal from separate northern and eastern glacial refugia. In addition, geographical and genetic structure among eastern populations suggested the existence of a second eastern refuge. Among the three lineages, the most divergent (Arctic) lineage occurred from Alaska east to northern Labrador. Quebec, New Brunswick, and New England were colonized by a second (Laurentian) lineage, and Labrador by a third group. Contact between refugial groups was only detected in two Labrador populations. The Arctic lineage was highly differentiated from eastern North American and European haplotypes, and probably diverged during the early Pleistocene. By contrast, the Laurentian and Labrador groups were similar to Old World charr, suggesting a shared ancestry during the mid-Pleistocene. In addition, the close relationship between Labrador and Laurentian charr indicates their probable divergence during the Wisconsinan glaciation.     

Vicariance and dispersal effects on phylogeographic structure and speciation in a widespread estuarine invertebrate

Vicariance and dispersal can strongly influence population genetic structure and allopatric speciation, but their importance in the origin of marine biodiversity is unresolved. In transitional estuarine environments, habitat discreteness and dispersal barriers could enhance divergence and provide insight to evolutionary mechanisms underlying marine and freshwater biodiversity. We examined this by assessing phylogeographic structure in the widespread amphipod Gammarus tigrinus across 13 estuaries spanning its northwest Atlantic range from Quebec to Florida. Mitochondrial cytochrome c oxidase I and nuclear internal transcribed spacer 1 phylogenies supported deep genetic structure consistent with Pliocene separation and cryptic northern and southern species. This break occurred across the Virginian-Carolinian coastal biogeographic zone, where an oceanographic discontinuity may restrict gene flow. Ten estuarine populations of the northern species occurred in four distinct clades, supportive of Pleistocene separation. Glaciation effects on genetic structure of estuarine populations are largely unknown, but analysis of molecular variance (AMOVA) supported a phylogeographic break among clades in formerly glaciated versus nonglaciated areas across Cape Cod, Massachusetts. This finding was concordant with patterns in other coastal species, though there was no significant relationship between latitude and genetic diversity. This supports Pleistocene vicariance events and divergence of clades in different northern glacial refugia. AMOVA results and private haplotypes in most populations support an allopatric distribution across estuaries. Clade mixture zones are consistent with historical colonization and human-mediated transfer. An isolation-by-distance model of divergence was detected after we excluded a suspected invasive haplotype in the St. Lawrence estuary. The occurrence of cryptic species and divergent population structure support limited dispersal, dispersed habitat distribution, and historical factors as important determinants of estuarine speciation and diversification.     

Intraspecific variation of seed floating ability in Sparganium emersum suggests a bimodal dispersal strategy

Water-mediated spread of seeds (hydrochory) plays an important role in the dispersal of aquatic plants. In this study we investigate intraspecific variation in floating ability and germination capacity of Sparganium emersum seeds in relation to seed mass, within three natural populations along the Rur River (the Netherlands–Germany). Our results suggest that S. emersum produces two types of seeds: (i) short-floating seeds (SFS) that sink within 4 weeks (approximately 71% of all seeds), and (ii) long-floating seeds (LFS) that float at least for 6 months (approximately 28% of all seeds). Our study further shows that short-floating seeds display a significantly higher germination (%) (SFS = 89.9% vs LFS = 32.6%), a faster germination rate (SFS = 8.71 ± 3.3 vs LFS = 9.32 ± 3.1 days to germination) and a higher mean seed mass (SFS = 15.17 ± 4.5 vs LFS = 11.25 ± 3.8 mg), compared to long-floating seeds. It is argued that the production of these two types of seeds by S. emersum plants, each type with a different potential for water-mediated dispersal, represents a bimodal hydrochoric dispersal strategy.     

Linking seed dispersal and genetic structure of trees: a biogeographical approach

Aim Natural and human‐induced differences in frugivore assemblages can influence the seed dispersal distances of trees. An important issue in seed dispersal systems is to understand whether differences in seed dispersal distances also affect the genetic structure of mature trees. One possible approach to test for a relationship between seed dispersal and the genetic structure of mature trees is to compare the genetic structure of two closely related tree species between two biogeographical regions that differ in frugivore assemblages and seed dispersal distances. Previous studies on two Commiphora species revealed that Commiphora guillauminii in Madagascar has a much lower seed dispersal distance than Commiphora harveyi in South Africa. We tested whether the lower seed dispersal distance might have caused decreased gene flow, resulting in a stronger genetic structure in Madagascar than in South Africa. Location Madagascar and South Africa. Methods Using amplified fragment length polymorphism markers we investigated the genetic structure of 134 trees in Madagascar and 158 trees in South Africa at a local and a regional spatial scale. Results In concordance with our hypothesis, kinship analysis suggests that gene flow was restricted mostly to 3 km in Madagascar and to 30 km in South Africa. At the local spatial scale, the genetic differentiation among groups of trees within sample sites was marginally significantly higher in Madagascar (F_ST = 0.069) than in South Africa (F_ST = 0.021). However, at a regional spatial scale genetic differentiation was lower in Madagascar (F_ST = 0.053) than in South Africa (F_ST = 0.163). Main conclusions Our results show that lower seed dispersal distances of trees were linked to higher genetic differentiation of trees only at a local spatial scale. This suggests that seed dispersal affects the genetic population structure of trees at a local, but not at a regional, spatial scale.     

Historical transoceanic dispersal of a freshwater shrimp: the colonization of the South Pacific by the genus Paratya (Atyidae)

Aim To infer the phylogenetic relationships within the freshwater shrimp genus Paratya Miers, 1882 (Atyidae) and to use these data to answer biogeographical questions about the location, timing and form of evolution of this genus in the South Pacific. Location Paratya are spread throughout various freshwater habitats in the western Pacific, with a disjunct northern range in the North Pacific (Japan, Korea, Ryukyu Islands, Siberia) and South Pacific (Australia, New Zealand, New Caledonia, Lord Howe, Norfolk Island). Methods Specimens were obtained from throughout its range. Mitochondrial sequences of cytochrome oxidase subunit I and 16S ribosomal DNA were analysed using phylogenetic techniques to identify whether landmasses are monophyletic and what the relationships are between landmasses. Molecular clock dating methods were used to date divergences between taxa. Results Each landmass was recovered as monophyletic. Japan/Ryukyu Islands is the most basal group, followed by New Zealand. Australian specimens form a sister group to a clade made up of two groups (New Caledonia and Lord Howe/Norfolk Island). The oldest divergence within the genus (between North and South Pacific) took place 12–19 Ma. Main conclusions The geographical origin of the genus (either Gondwana or Laurasia) is unclear. Dispersal occurred between the North and South Pacific long after the split up of Gondwana. Dispersal likely explains the presence of Paratya on each landmass in the South Pacific, from continent to isolated oceanic island. This dispersal is conjectured to have taken place through oceanic currents because of the amphidromous life cycle of some taxa of Paratya, given that amphyidromy is plesiomorphic in atyid shrimp.     

Phylogeography of the widespread African puff adder (Bitis arietans) reveals multiple Pleistocene refugia in southern Africa

Evidence from numerous Pan‐African savannah mammals indicates that open‐habitat refugia existed in Africa during the Pleistocene, isolated by expanding tropical forests during warm and humid interglacial periods. However, comparative data from other taxonomic groups are currently lacking. We present a phylogeographic investigation of the African puff adder (Bitis arietans), a snake that occurs in open‐habitat formations throughout sub‐Saharan Africa. Multiple parapatric mitochondrial clades occur across the current distribution of B. arietans, including a widespread southern African clade that is subdivided into four separate clades. We investigated the historical processes responsible for generating these phylogeographic patterns in southern Africa using species distribution modelling and genetic approaches. Our results show that interior regions of South Africa became largely inhospitable for B. arietans during glacial maxima, whereas coastal and more northerly areas remained habitable. This corresponds well with the locations of refugia inferred from mitochondrial data using a continuous phylogeographic diffusion model. Analysis of data from five anonymous nuclear loci revealed broadly similar patterns to mtDNA. Secondary admixture was detected between previously isolated refugial populations. In some cases, this is limited to individuals occurring near mitochondrial clade contact zones, but in other cases, more extensive admixture is evident. Overall, our study reveals a complex history of refugial isolation and secondary expansion for puff adders and a mosaic of isolated refugia in southern Africa. We also identify key differences between the processes that drove isolation in B. arietans and those hypothesized for sympatric savannah mammals.     

Restricted dispersal reduces the strength of spatial density dependence in a tropical bird population

Spatial processes could play an important role in density-dependent population regulation because the disproportionate use of poor quality habitats as population size increases is widespread in animal populations-the so-called buffer effect. While the buffer effect patterns and their demographic consequences have been described in a number of wild populations, much less is known about how dispersal affects distribution patterns and ultimately density dependence. Here, we investigated the role of dispersal in spatial density dependence using an extraordinarily detailed dataset from a reintroduced Mauritius kestrel (Falco punctatus) population with a territorial (despotic) breeding system. We show that recruitment rates varied significantly between territories, and that territory occupancy was related to its recruitment rate, both of which are consistent with the buffer effect theory. However, we also show that restricted dispersal affects the patterns of territory occupancy with the territories close to release sites being occupied sooner and for longer as the population has grown than the territories further away. As a result of these dispersal patterns, the strength of spatial density dependence is significantly reduced. We conclude that restricted dispersal can modify spatial density dependence in the wild, which has implications for the way population dynamics are likely to be impacted by environmental change.     

Adding energy gradients and long‐distance dispersal to a neutral model improves predictions of Madagascan bird diversity

Macroecological patterns are likely the result of both stochastically neutral mechanisms and deterministic differences between species. In Madagascar, the simplest stochastically neutral hypothesis – the mid‐domain effects (MDE) hypothesis – has already been rejected. However, rejecting the MDE hypothesis does not necessarily refute the existence of all other neutral mechanisms. Here, we test whether adding complexity to a basic neutral model improves predictions of biodiversity patterns. The simplest MDE model assumes that: (1) species' ranges are continuous and unfragmented, (2) are randomly located throughout the landscape, and (3) can be stacked independently and indefinitely. We designed a simulation based on neutral theory that allowed us to weaken each of these assumptions incrementally by adjusting the habitat capacity as well as the likelihood of short‐ and long‐distance dispersal. Simulated outputs were compared to four empirical patterns of bird diversity: the frequency distributions of species richness and range size, the within‐island latitudinal diversity gradient, and the distance‐decay of species compositional similarity. Neutral models emulated empirical diversity patterns for Madagascan birds accurately. The frequency distribution of range size, latitudinal diversity gradient, and the distance‐decay of species compositional similarity could be attributed to stochastic long‐distance migration events and zero‐sum population dynamics. However, heterogenous environmental gradients improved predictions of the frequency distribution of species richness. Patterns of bird diversity in Madagascar can broadly be attributed to stochastic long‐distance migration events and zero‐sum population dynamics. This implies that rejecting simple hypotheses, such as MDE, does not serve as evidence against stochastic processes in general. However, environmental gradients were necessary to explain patterns of species richness and deterministic differences between species are probably important for explaining the distributions of narrow‐range and endemic species.     

Effects of habitat fragmentation and changes of dispersal behaviour after a recent population decline on the genetic variability of noncoding and coding DNA of a monogamous Malagasy rodent

While interactions among demography, behaviour and genetic structure are well-documented for neutral genetic markers, the role of these parameters and the effects of genetic drift and selection are considerably less well understood in functional genes, such as the major histocompatibility complex (MHC). In this study, the consequences of habitat fragmentation and the effects of a current population decline on noncoding (mitochondrial DNA) and two coding MHC loci (DQA, DRB) with different functional importance were investigated in the small remnant subdivided population of the endangered Malagasy giant jumping rat (Hypogeomys antimena). Both neutral and selective markers revealed a significant genetic differentiation between the two remnant subpopulations. The FST values were much lower in the MHC DQA and DRB genes than in the mitochondrial data. The MHC DRB loci display the effects of both balancing selection (high sequence diversity, four times higher nonsynonymous than synonymous substitutions in the functionally important antigen-binding site positions, twice the average heterozygosity of individual amino acids at the positions identified as part of the antigen-binding site (ABS) than those outside the ABS and nonselective forces including genetic drift. Simultaneously with a current population decline offspring reduced their dispersal distances. No substantial effects were detected within the first 6 years of reduced gene flow in either mitochondrial or MHC markers.