Genetic structure among Charadrius plovers on the African mainland and islands of Madagascar and St Helena

Colonization of islands by long-distance dispersers has great impact on genetic diversification among populations and may spearhead speciation events. We investigated intra- and interspecific divergence in Charadrius plovers with populations on mainland Africa, Madagascar and St Helena. We analysed microsatellite loci and sequence data from four nuclear and two mitochondrial gene regions. Charadrius plovers are shorebirds with high dispersal and mobility. Our results confirmed genetic differentiation between Madagascar and mainland populations of three plover species (White-fronted Plover Charadrius marginatus, Kittlitz's Plover Charadrius pecuarius and, based on sequence data only, Three-banded Plover Charadrius tricollaris) but highlight substantial variation in levels of intraspecific divergence among the three species. Namely, the Kittlitz's Plover, a dispersive habitat generalist with a polygamous mating system, exhibited lower island–mainland differentiation (0.05% COI sequence divergence) compared with the two monogamous species, the White-fronted Plover (0.6% COI divergence) and Three-banded Plover (1.6% COI divergence). In addition, past colonization of the islands of St Helena and Madagascar by ancestors of today's Kittlitz's Plover has led to the evolution of two endemic island species, the Madagascar Plover Charadrius thoracicus and the more closely-related St Helena Plover Charadrius sanctaehelenae. We discuss the factors driving species differences in island–mainland divergence and highlight the importance of conserving genetically unique island populations and island habitats to safeguard future evolutionary potential.     

Evolutionary history of Scinax treefrogs on land‐bridge islands in south‐eastern Brazil

Aim We investigated how Pleistocene refugia and recent (c. 12,000 years ago) sea level incursions shaped genetic differentiation in mainland and island populations of the Scinax perpusillus treefrog group. Location Brazilian Atlantic Forest, São Paulo state, south‐eastern Brazil. Methods Using mitochondrial and microsatellite loci, we examined population structure and genetic diversity in three species from the S. perpusillus group, sampled from three land‐bridge islands and five mainland populations, in order to understand the roles of Pleistocene forest fragmentation and sea level incursions on genetic differentiation. We calculated metrics of relatedness and genetic diversity to assess whether island populations exhibit signatures of genetic drift and isolation. Two of the three island populations in this study have previously been described as new species based on a combination of distinct morphological and behavioural characters, thus we used the molecular datasets to determine whether phenotypic change is consistent with genetic differentiation. Results Our analyses recovered three distinct lineages or demes composed of northern mainland São Paulo populations, southern mainland São Paulo populations, and one divergent island population. The two remaining island populations clustered with samples from adjacent mainland populations. Estimates of allelic richness were significantly lower, and estimates of relatedness were significantly higher, in island populations relative to their mainland counterparts. Main conclusions Fine‐scale genetic structure across mainland populations indicates the possible existence of local refugia within São Paulo state, underscoring the small geographic scale at which populations diverge in this species‐rich region of the Atlantic Coastal Forest. Variation in genetic signatures across the three islands indicates that the populations experienced different demographic processes after marine incursions fragmented the distribution of the S. perpusillus group. Genetic signatures of inbreeding and drift in some island populations indicate that small population sizes, coupled with strong ecological selection, may be important evolutionary forces driving speciation on land‐bridge islands.     

A question of time: the land snail Murella muralis (Gastropoda: Pulmonata) reveals constraints on past ecological speciation

The lively debate about speciation currently focuses on the relative importance of factors driving population differentiation. While many studies are increasingly producing results on the importance of selection, little is known about the interaction between drift and selection. Moreover, there is still little knowledge on the spatial‐temporal scales at which speciation occurs, that is, arrangement of habitat patches, abruptness of habitat transitions, climate and habitat changes interacting with selective forces. To investigate these questions, we quantified variation on a fine geographical scale analysing morphological (shell) and genetic data sets coupled with environmental data in the land snail Murella muralis, endemic to the Mediterranean island of Sicily. Analysis of a fragment of the mitochondrial DNA cytochrome oxidase I gene (COI) and eight nuclear microsatellite loci showed that genetic variation is highly structured at a very fine spatial scale by local palaeogeographical events and historical population dynamics. Molecular clock estimates, calibrated here specifically for Tyrrhenian land snails, provided a framework of palaeogeographical events responsible for the observed geographical variations and migration routes. Finally, we showed for the first time well‐documented lines of evidence of selection in the past, which explains divergence of land snail shell shapes. We suggest that time and palaeogeographical history acted as constraints in the progress along the ecological speciation continuum. Our study shows that testing for correlation among palaeogeography, morphology and genetic data on a fine geographical scale provides information fundamental for a detailed understanding of ecological speciation processes.     

Inter‐island divergence within Drosophila mauritiana,a species of the D. simulans complex: Past history and/or speciation in progress?

Speciation with gene flow may be more common than generally thought, which makes detailed understanding of the extent and pattern of genetic divergence between geographically isolated populations useful. Species of the Drosophila simulans complex provide a good model for speciation and evolutionary studies, and hence understanding their population genetic structure will increase our understanding of the context in which speciation has occurred. Here, we describe genetic diversity and genetic differentiation of two distant populations of D. mauritiana (Mauritius and Rodrigues Islands) at mitochondrial and nuclear loci. We surveyed the two populations for their mitochondrial haplotypes, eight nuclear genes and 18 microsatellite loci. A new mitochondrial type is fixed in the Rodrigues population of D. mauritiana. The two populations are highly differentiated, their divergence appears relatively ancient (100,000 years) compared to the origin of the species, around 0.25MYA, and they exhibit very limited gene flow. However, they have similar levels of divergence from their sibling, D. simulans. Both nuclear genes and microsatellites revealed contrasting demographic histories between the two populations, expansion for the Mauritius population and stable population size for the Rodrigues Island population. The discovery of pronounced geographic structure within D. mauritiana combined to genetic structuring and low gene flow between the two island populations illuminates the evolutionary history of the species and clearly merits further attention in the broad context of speciation.     

Rejecting strictly allopatric speciation on a continental island: prolonged postdivergence gene flow between Taiwan (Leucodioptron taewanus,Passeriformes Timaliidae) and Chinese (L. canorum canorum) hwameis

Allopatry is conventionally considered the geographical mode of speciation for continental island organisms. However, strictly allopatric speciation models that assume the lack of postdivergence gene flow seem oversimplified given the recurrence of land bridges during glacial periods since the late Pliocene. Here, to evaluate whether a continental island endemic, the Taiwan hwamei (Leucodioptron taewanus, Passeriformes Timaliidae) speciated in strict allopatry, we used weighted‐regression‐based approximate Bayesian computation (ABC) to analyse the genetic polymorphism of 18 neutral nuclear loci (total length: 8500 bp) in Taiwan hwamei and its continental sister species, the Chinese hwamei (L. canorum canorum). The nonallopatry model was found to fit better with observed genetic polymorphism of the two hwamei species (posterior possibility = 0.82). We also recovered unambiguous signals of nontrivial bidirectional postdivergence gene flow (N_em » 1) between Chinese hwamei and Taiwan hwamei until 0.5 Ma. Divergence time was estimated to be 3.5 to 2 million years earlier than that estimated from mitochondrial cytochrome b sequences. Finally, using the inferred nonallopatry model to simulate genetic variation at 24 nuclear genes examined showed that the adiponectin receptor 1 gene may be under divergent adaptation. Our findings imply that the role of geographical barrier may be less prominent for the speciation of continental island endemics, and suggest a shift in speciation studies from simply correlating geographical barrier and genetic divergence to examining factors that facilitate and maintain divergence, e.g. differential selection and sexual selection, especially in the face of interpopulation gene flow.     

Experimental evidence of genome‐wide impact of ecological selection during early stages of speciation‐with‐gene‐flow

Theory predicts that speciation‐with‐gene‐flow is more likely when the consequences of selection for population divergence transitions from mainly direct effects of selection acting on individual genes to a collective property of all selected genes in the genome. Thus, understanding the direct impacts of ecologically based selection, as well as the indirect effects due to correlations among loci, is critical to understanding speciation. Here, we measure the genome‐wide impacts of host‐associated selection between hawthorn and apple host races of Rhagoletis pomonella (Diptera: Tephritidae), a model for contemporary speciation‐with‐gene‐flow. Allele frequency shifts of 32 455 SNPs induced in a selection experiment based on host phenology were genome wide and highly concordant with genetic divergence between co‐occurring apple and hawthorn flies in nature. This striking genome‐wide similarity between experimental and natural populations of R. pomonella underscores the importance of ecological selection at early stages of divergence and calls for further integration of studies of eco‐evolutionary dynamics and genome divergence.     

Speciation of silverside Chirostoma attenuatum (Pisces: Atheriniformes) in Central Mexico

The high speciation rate of Chirostoma in Central Mexico has been associated with allopatric speciation events promoted by the emergence of vicariant barriers in freshwater habitats, as well as by sympatric ecological segregation, common in those species inhabiting lacustrine ecosystems. Through nuclear and mitochondrial markers, this study revealed a speciation process within Chirostoma attenuatum resulting in two evolutionary independent units that coincide with their morphological differentiation, indicating that Chirostoma attenuatum and Chirostoma zirahuen may be considered separate species. This process was the result of vicariance associated with geological dynamics of the region. Phylogeographic findings indicated two speciation stages: early allopatric isolation, during which the isolated populations accumulated unique adaptations, and secondary contact with low migration rate and the maintenance of the evolutionary trajectory. Historical demographic analysis indicated that the two well‐differentiated lineages underwent independent evolutionary histories in their respective lakes. Chirostoma zirahuen from Zirahuen and C. attenuatum from Patzcuaro represent unique and irreplaceable genetic diversity that must to be conserved.     

Genetic differentiation among populations of a Hispaniolan trunk anole that exhibit geographical variation in dewlap colour

Long neglected by classic island biogeographical theory, speciation within and among islands is increasingly recognized as a major contributor to insular diversity. Although the factors responsible for island speciation remain poorly understood, this process appears critically dependent on geographical variation and speciation in allopatry or parapatry. Here, we investigate geographical variation and speciation in a complex of Hispaniolan trunk anoles (Anolis distichus), where populations with strikingly distinct dewlap colours and patterns correspond with deeply divergent mtDNA structure. Using a multilocus, population‐level analysis, we investigate whether these phenotypically and mitochondrially distinct populations exhibit the type of nuclear differentiation expected among species or incipient species. Along a transect that extends across a recently recessed marine barrier, our results are consistent with the persistence of an abrupt phenotypic and mitochondrial transition following secondary contact, in spite of little or no evidence for a reduction in nuclear gene flow. Along a second transect extending across a steep environmental gradient, our phenotypic and microsatellite data suggest a sharp genetic break with little or no admixture, whereas mtDNA recovers a signature of extensive unidirectional introgression. Together, these results are consistent with previous studies of Lesser Antillean anoles, suggesting that allopatric divergence alone is insufficient for speciation, whereas reduced gene flow and partial reproductive isolation may accumulate in the presence of ecological gradients.     

Historical introgression and the persistence of ghost alleles in the intermediate horseshoe bat (Rhinolophus affinis)

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

Morphotypes or distinct species? A multilocus assessment of two East Asian scimitar babblers (Aves,Timaliidae)

In birds, widespread species complexes often exhibit dramatic plumage differences across their distribution, which can give rise to discordance between morphotypes and phylogroups. Accurate identification of species diversity may require an integrated approach in which multilocus genetic data are used for inference and further corroborated by ecological, morphometric or behavioural data. Pomatorhinus gravivox and P. swinhoei, which were formerly considered component races of P. erythrogenys, were recently split using a delimitation system that mainly relied on the quantification of differences in phenotypic traits. We therefore carried out a reassessment of this taxonomic recommendation by conducting phylogenetic and population genetic analyses with multilocus genetic data. A deep mitochondrial split with a Kimura 2-parameter distance of 0.061 was observed, mainly corresponding to the two morphologically defined species. However, one individual from P. swinhoei harboured a haplotype of P. gravivox. Individual-based analyses of nuclear loci identified two distinct clusters that were exactly congruent with the two species. BPP delimitation also provided support for the separation. These congruencies support the notion that these two taxa are best regarded as two separate species. The presence of the P. gravivox mitochondrial haplotype in P. swinhoei was most likely a result of hybridization, due to the clear separation of nuclear loci. The speciation might be attributed to paleoclimatic changes but requires further evaluation due to the likelihood of ecological speciation. This study was in accordance with the results inferred from the quantitative system but highlighted the importance of sampling various data, especially in contact zones, in the study of taxonomy and speciation history.     

ALLOPATRIC DIVERGENCE AND SPECIATION IN CORAL REEF FISH: THE THREE‐SPOT DASCYLLUS,DASCYLLUS TRIMACULATUS,SPECIES COMPLEX

Long pelagic larval phases and the absence of physical barriers impede rapid speciation and contrast the high diversity observed in marine ecosystems such as coral reefs. In this study, we used the three‐spot dascyllus (Dascyllus trimaculatus) species complex to evaluate speciation modes at the spatial scale of the Indo‐Pacific. The complex includes four recognized species and four main color morphs that differ in distribution. Previous studies of the group using mitochondrial DNA revealed a noncongruence between color morphs and genetic groupings; with two of the color morphs grouped together and one color morph separated into three clades. Using extensive geographic sampling of 563 individuals and a combination of mitochondrial DNA sequences and 13 nuclear microsatellites, we defined population/species boundaries and inferred different speciation modes. The complex is composed of seven genetically distinct entities, some of which are distinct morphologically. Despite extensive dispersal abilities and an apparent lack of barriers, observed genetic partitions are consistent with allopatric speciation. However, ecological pressure, assortative mating, and sexual selection, were likely important during periods of geographical isolation. This study therefore suggests that primarily historical factors later followed by ecological factors caused divergence and speciation in this group of coral reef fish.     

A comparative analysis of the mechanisms underlying speciation on Lord Howe Island

On Lord Howe Island, speciation is thought to have taken place in situ in a diverse array of distantly related plant taxa (Metrosideros, Howea and Coprosma; Proc. Natl Acad. Sci. USA 108 , 2011, 13188). We now investigate whether the speciation processes were driven by divergent natural selection in each genus by examining the extent of ecological and genetic divergence. We present new and extensive, ecological and genetic data for all three genera. Consistent with ecologically driven speciation, outlier loci were detected using genome scan methods. This mechanism is supported by individual‐based analyses of genotype–environment correlations within species, demonstrating that local adaptation is currently widespread on the island. Genetic analyses show that prezygotic isolating barriers within species are currently insufficiently strong to allow further population differentiation. Interspecific hybridization was found in both Howea and Coprosma, and species distribution modelling indicates that competitive exclusion may result in selection against admixed individuals. Colonization of new niches, partly fuelled by the rapid generation of new adaptive genotypes via hybridization, appears to have resulted in the adaptive radiation in Coprosma – supporting the ‘Syngameon hypothesis’.     

Evidence for ecological speciation via a host shift in the holly leaf miner,Phytomyza glabricola (Diptera: Agromyzidae)

Evolutionary radiations have been well documented in plants and insects, and natural selection may often underly these radiations. If radiations are adaptive, the diversity of species could be due to ecological speciation in these lineages. Agromyzid flies exhibit patterns of repeated host‐associated radiations. We investigated whether host‐associated population divergence and evidence of divergent selection exist in the leaf miner Phytomyza glabricola on its sympatric host plants, the holly species, Ilex coriacea and I. glabra. Using AFLPs and nuclear sequence data, we found substantial genetic divergence between host‐associated populations of these flies throughout their geographic range. Genome scans using the AFLP data identified 13 loci under divergent selection, consistent with processes of ecological speciation. EF‐1α data suggest that I. glabra is the original host of P. glabricola and that I. coriacea is the novel host, but the AFLP data are ambiguous with regard to directionality of the host shift.     

Late Pleistocene speciation of three closely related tree peonies endemic to the Qinling–Daba Mountains,a major glacial refugium in Central China

Determining the factors promoting speciation is a major task in ecological and evolutionary research and can be aided by phylogeographic analysis. The Qinling–Daba Mountains (QDM) located in central China form an important geographic barrier between southern subtropical and northern temperate regions, and exhibit complex topography, climatic, and ecological diversity. Surprisingly, few phylogeographic analyses and studies of plant speciation in this region have been conducted. To address this issue, we investigated the genetic divergence and evolutionary histories of three closely related tree peony species (Paeonia qiui, P. jishanensis, and P. rockii) endemic to the QDM. Forty populations of the three tree peony species were genotyped using 22 nuclear simple sequence repeat markers (nSSRs) and three chloroplast DNA sequences to assess genetic structure and phylogenetic relationships, supplemented by morphological characterization and ecological niche modeling (ENM). Morphological and molecular genetic analyses showed the three species to be clearly differentiated from each other. In addition, coalescent analyses using DIYABC conducted on nSSR variation indicated that the species diverged from each other in the late Pleistocene, while ecological niche modeling (ENM) suggested they occupied a larger area during the Last Glacial Maximum (LGM) than at present. The combined genetic evidence from nuclear and chloroplast DNA and the results of ENM indicate that each species persisted through the late Pleistocene in multiple refugia in the Qinling, Daba, and Taihang Mountains with divergence favored by restricted gene flow caused by geographic isolation, ecological divergence, and limited pollen and seed dispersal. Our study contributes to a growing understanding of the origin and population structure of tree peonies and provides insights into the high level of plant endemism present in the Qinling–Daba Mountains of Central China.     

Evolutionary insights into Rhinolophus episcopus (Chiroptera,Rhinolophidae) in China: Isolation by distance,environment, or sensory system?

Evolutionary processes can be influenced by several factors, such as geographic isolation, environmental selection, and sensory variation. For most nocturnal bats, echolocation is the primary sensory system used to prey and communicate, and plays important roles in chiropteran diversification and evolution. Understanding the relative contribution of geography, the environment, and this sensory system to population genetic divergence can elucidate the processes involved in bat incipient speciation and evolution. In this study, we collected spatial and environmental information, echolocation calls, as well as the previously published genetic data (six microsatellite loci and the mitochondrial cytochrome b gene) of widely distributed Rhinolophus episcopus populations to test three hypotheses for nuclear and mitochondrial divergence (isolation by distance, isolation by environment, and isolation by sensory variation) and unveil the factors that drive intraspecific genetic differentiation. The moderate level of nuclear differentiation was correlated with geographic/spatial distance and acoustic variation, whereas the relatively high level of mitochondrial differentiation was mainly associated with acoustic divergence. No significant correlation was observed between genetic divergence and environmental variables. Among the three factors, acoustic divergence explained the highest percentage of both nuclear and mitochondrial divergence. Thus, our results indicate that sensory variation may have played important roles in driving population isolation early in bat speciation, which is consistent with the hypothesis of isolation by sensory variation. Our study emphasizes the need to consider more factors, especially sensory traits, and combine multiple statistical methods in landscape genetic studies to test their potential contributions to driving population divergence.     

Nuclear locus divergence at the early stages of speciation in the Orchard Oriole complex

As two lineages diverge from one another, mitochondrial DNA should evolve fixed differences more rapidly than nuclear DNA due to its smaller effective population size and faster mutation rate. As a consequence, molecular systematists have focused on the criteria of reciprocal monophyly in mitochondrial DNA for delimiting species boundaries. However, mitochondrial gene trees do not necessarily reflect the evolutionary history of the taxa in question, and even mitochondrial loci are not expected to be reciprocally monophyletic when the speciation event happened very recently. The goal of this study was to examine mitochondrial paraphyly within the Orchard Oriole complex, which is composed of Icterus spurius (Orchard Oriole) and Icterus fuertesi (Fuertes' Oriole). We increased the geographic sampling, added four nuclear loci, and used a range of population genetic and coalescent methods to examine the divergence between the taxa. With increased taxon sampling, we found evidence of clear structure between the taxa for mitochondrial DNA. However, nuclear loci showed little evidence of population structure, indicating a very recent divergence between I. spurius and I. fuertesi. Another goal was to examine the genetic variation within each taxon to look for evidence of a past founder event within the I. fuertesi lineage. Based on the high amounts of genetic variation for all nuclear loci, we found no evidence of such an event – thus, we found no support for the possible founding of I. fuertesi through a change in migratory behavior, followed by peripheral isolates speciation. Our results demonstrate that these two taxa are in the earliest stages of speciation, at a point when they have fixed differences in plumage color that are not reflected in monophyly of the mitochondrial or nuclear DNA markers in this study. This very recent divergence makes them ideal for continued studies of species boundaries and the earliest stages of speciation.     

DIVERSIFICATION AND GENE FLOW IN NASCENT LINEAGES OF ISLAND AND MAINLAND NORTH AMERICAN TREE SQUIRRELS (TAMIASCIURUS)

Pleistocene climate cycles and glaciations had profound impacts on taxon diversification in the Boreal Forest Biome. Using population genetic analyses with multilocus data, we examined diversification, isolation, and hybridization in two sibling species of tree squirrels (Tamiasciurus douglasii and Tamiasciurus hudsonicus) with special attention to the geographically and genetically enigmatic population of T. hudsonicus on Vancouver Island, Canada. The two species differentiated only about 500,000 years ago, in the Late Pleistocene. The island population is phylogenetically nested within T. hudsonicus according to our nuclear analysis but within T. douglasii according to mitochondrial DNA. This conflict is more likely due to historical hybridization than to incomplete lineage sorting, and it appears that bidirectional gene flow occurred between the island population and both species on the mainland. This interpretation of our genetic analyses is consistent with our bioclimatic modeling, which demonstrates that both species were able to occupy this region throughout the Late Pleistocene. The divergence of the island population 40,000 years ago suggests that tree squirrels persisted in a refugium on Vancouver Island at the last glacial maximum, 20,000 years ago. Our observations demonstrate how Pleistocene climate change and habitat shifts have created incipient divergence in the presence of gene flow.     

Sky island diversification meets the multispecies coalescent – divergence in the spruce‐fir moss spider (Microhexura montivaga,Araneae, Mygalomorphae) on the highest peaks of southern Appalachia

Microhexura montivaga is a miniature tarantula‐like spider endemic to the highest peaks of the southern Appalachian mountains and is known only from six allopatric, highly disjunct montane populations. Because of severe declines in spruce‐fir forest in the late 20th century, M. montivaga was formally listed as a US federally endangered species in 1995. Using DNA sequence data from one mitochondrial and seven nuclear genes, patterns of multigenic genetic divergence were assessed for six montane populations. Independent mitochondrial and nuclear discovery analyses reveal obvious genetic fragmentation both within and among montane populations, with five to seven primary genetic lineages recovered. Multispecies coalescent validation analyses [guide tree and unguided Bayesian Phylogenetics and Phylogeography (BPP), Bayes factor delimitation (BFD)] using nuclear‐only data congruently recover six or seven distinct lineages; BFD analyses using combined nuclear plus mitochondrial data favour seven or eight lineages. In stark contrast to this clear genetic fragmentation, a survey of secondary sexual features for available males indicates morphological conservatism across montane populations. While it is certainly possible that morphologically cryptic speciation has occurred in this taxon, this system may alternatively represent a case where extreme population genetic structuring (but not speciation) leads to an oversplitting of lineage diversity by multispecies coalescent methods. Our results have clear conservation implications for this federally endangered taxon and illustrate a methodological issue expected to become more common as genomic‐scale data sets are gathered for taxa found in naturally fragmented habitats.     

Pleistocene origin and population history of a neoendemic alpine butterfly

Alpine environments underwent dramatic transformation during glacial–interglacial cycles, with the consequence that geographical, ecological and demographic changes of alpine populations provided the opportunity for formation of neoendemic species. Several biogeographical models have been proposed to account for the unique history of alpine populations, with different expectations of genetic divergence and speciation. The expanding alpine archipelago model proposes that alpine populations expand spatially and demographically during glacial events, dispersing between mountain ranges. Under this model, alpine populations are unlikely to diverge in isolation due to substantial interpopulation gene flow. In contrast, the alpine archipelago refuge model proposes that gene flow during glacial phases is limited and populations expand demographically during interglacial phases, increasing genetic isolation and the likelihood of speciation. We assess these models by reconstructing the evolutionary history of Colias behrii, a morphologically and ecologically distinct alpine butterfly restricted to the California Sierra Nevada. C. behrii exhibits very low genetic diversity at mitochondrial and nuclear loci, limited population structure and evidence of population expansion. C. behrii and Rocky Mountain C. meadii share identical mitochondrial haplotypes, while in contrast, nuclear data indicate common ancestry between C. behrii and Cascades Range Colias pelidne. The conflict in gene genealogies may be a result of recent expansion in North American Colias, but an isolation with migration analysis indicates that genetic patterns in C. behrii might result from differential introgression following hybridization. Based on the timing of population expansion and gene flow between mountain ranges, the expanding alpine archipelago model is supported in C. behrii.     

Sequential adaptive introgression of the mitochondrial genome in Drosophila yakuba and Drosophila santomea

Interspecific hybridization provides the unique opportunity for species to tap into genetic variation present in a closely related species and potentially take advantage of beneficial alleles. It has become increasingly clear that when hybridization occurs, mitochondrial DNA (mtDNA) often crosses species boundaries, raising the possibility that it could serve as a recurrent target of natural selection and source of species' adaptations. Here we report the sequences of 46 complete mitochondrial genomes of Drosophila yakuba and Drosophila santomea, two sister species known to produce hybrids in nature (~3%). At least two independent events of mtDNA introgression are uncovered in this study, including an early invasion of the D. yakuba mitochondrial genome that fully replaced the D. santomea mtDNA native haplotypes and a more recent, ongoing event centred in the hybrid zone. Interestingly, this recent introgression event bears the signature of Darwinian natural selection, and the selective haplotype can be found at low frequency in Africa mainland populations of D. yakuba. We put forward the possibility that, because the effective population size of D. santomea is smaller than that of D. yakuba, the faster accumulation of mildly deleterious mutations associated with Muller's ratchet in the former species may have facilitated the replacement of the mutationally loaded mitochondrial genome of D. santomea by that of D. yakuba.