The biogeography of mitochondrial and nuclear discordance in animals

Taking the discovery approach in integrative taxonomy: decrypting a complex of narrow-endemic Alpine harvestmen (Opiliones: Phalangiidae: Megabunus)

Species delimitation is fundamental for biological studies, yet precise delimitation is not an easy task, and every involved approach has an inherent failure rate. Integrative taxonomy, a method that merges multiple lines of evidence, can profoundly contribute to reliable alpha-taxonomy and shed light on the processes behind speciation. In this study, we explored and validated species limits in a group of closely related Megabunus harvestmen (Eupnoi, Phalangiidae) endemic to the European Alps. Without a priori species hypotheses, we used multiple sources of inference, including mitochondrial and multilocus nuclear DNA, morphometrics and chemistry. The results of these discovery approaches revealed morphological crypsis and multiple new species within two of the five hitherto known species. Based on our analyses, we discussed the most plausible evolutionary scenarios, invoked the most reasonable species hypotheses and validated the new species limits. Building upon the achieved rigour, three new species, Megabunus cryptobergomas Muster and Wachter sp. nov., Megabunus coelodonta Muster and Steiner sp. nov., and Megabunus lentipes Muster and Komposch sp. nov., are formally described. In addition, we provide a dichotomous morphological key to the Megabunus species of the Alps. Our work demonstrates the suitability of integrative, discovery-based approaches in combination with validation approaches to precisely characterize species and enabled us to implement nomenclatural consequences for this genus.

The phylogeny of pikas (Ochotona) inferred from a multilocus coalescent approach

The clarification of the systematics of pikas (genus Ochotona) has been hindered by largely overlapping morphological characters among species and the lack of a comprehensive molecular phylogeny. Here we estimate the first multilocus phylogeny of the genus to date, by analysing 12 nuclear DNA markers (total of 7.5Kb) in 11 species of pikas from the four classified subgenera (Pika, Ochotona, Lagotona and Conothoa) using a multispecies coalescent-based framework. The species-tree confirmed the subgeneric classification by retrieving as monophyletic the subgenera represented here by more than one species. Contrary to previous phylogenies based on mtDNA alone, Lagotona was found to be sister to Pika. Also, support for the monophyly of the alpina group was not strong, thus caution should be used in future analyses of this group. A relaxed molecular clock calibrated using the Ochotonidae-Leporidae divergence resulted in more recent estimates of divergence times relative to previous studies. Strong concordance with inferences based on fossil records was found, suggesting that the initial diversification of the genus took place by the end of late Miocene. Finally, this work sets up methodologies and gathers molecular markers that can be used to extend the understanding of the evolutionary history of the genus.

Evidence for high dispersal ability and mito-nuclear discordance in the small brown planthopper, Laodelphax striatellus

Understanding dispersal ability in pest species is critical for both theoretical aspects of evolutionary and population biology and from a practical standpoint, such as implementing effective forecasting systems. The small brown planthopper (SBPH), Laodelphax striatellus (Fallén), is an economically important pest, but few data exist on its dispersal ability. Here, we used mitochondrial and nuclear markers to elucidate the population genetic structure of SBPH and of the parasitic bacterium Wolbachia throughout temperate and subtropical China. Our results showed that the SBPH populations in China lack significant differences in genetic structure, suggesting extensive gene flow. Multilocus sequence typing revealed that Wolbachia infection was systematic and due to the same strain (wStri) within and across populations. However, the mtDNA haplogroups had a nonrandom distribution across the sampling localities, which correlated to latitudinal and climatic gradients. We explain this mito-nuclear discordance as a result of historical population recolonization or mitochondria adaptation to climate.

Differential introgression and effective size of marker type influence phylogenetic inference of a recently divergent avian group (Phasianidae: Tympanuchus)

Life history strategies can influence the effective population size (Ne) of loci differently based on their mode of inheritance. Recognizing how this may affect the rate of lineage sorting among marker types is important for studies focused on resolving phylogenetic relationships among recently divergent taxa. In this study, we use gene tree, coalescent-based species tree, and isolation-with-migration analyses to explore the differences between marker types (autosomal, Z-linked, and mitochondrial) in resolving phylogenetic relationships among North American prairie grouse (Tympanuchus). We found that Z-linked loci were more likely to identify monophyletic relationships among prairie grouse species compared to autosomal and mtDNA loci in both species and gene tree analyses, with species tree analyses outperforming gene trees. These results were further supported with isolation-with-migration analyses, where Z-linked loci largely followed a strict isolation model while autosomal loci were more likely to fit a model with gene flow between species following population divergence. While accounting for differences in inheritance pattern (or Ne) for marker type, results suggest that additional factors, such as strong sexual selection and sex-biased introgression (i.e., male-biased postzygotic hybrid behavioral isolation or "unsexy son"), may further explain the decreased diversity levels and increased rate of lineage sorting observed with the Z-linked loci relative to autosomal and mtDNA loci. In fact, to our knowledge no hybrid male prairie grouse have been observed breeding in the wild, yet hybrid females along with backcross females are known to produce viable offspring. Overall, this study highlights that more work is needed to determine how complex models of gene flow (i.e., sex biased introgression) and differences in the effective size among marker types based on differing life history strategies influence divergence date estimation and species delimitation.

Genetic and phenotypic characterization of a hybrid zone between polyandrous Northern and Wattled Jacanas in Western Panama

Background

Hybridization provides a unique perspective into the ecological, genetic and behavioral context of speciation. Hybridization is common in birds, but has not yet been reported among bird species with a simultaneously polyandrous mating system; a mating system where a single female defends a harem of males who provide nearly all parental care. Unlike simple polyandry, polyandrous mating is extremely rare in birds, with only 1% of bird species employing this mating system. Although it is classically held that females are “choosy” in avian hybrid systems, nearly-exclusive male parental care raises the possibility that female selection against heterospecific matings might be reduced compared to birds with other mating systems.

Results

We describe a narrow hybrid zone in southwestern Panama between two polyandrous freshwater waders: Northern Jacana, Jacana spinosa and Wattled Jacana, J. jacana. We document coincident cline centers for three phenotypic traits, mtDNA, and one of two autosomal introns. Cline widths for these six markers varied from seven to 142 km, with mtDNA being the narrowest, and five of the six markers having widths less than 100 km. Cline tails were asymmetrical, with greater introgression of J. jacana traits extending westward into the range of J. spinosa. Likewise, within the hybrid zone, the average hybrid index of phenotypic hybrids was significantly biased towards J. spinosa. Species distribution models indicate that the hybrid zone is located at the edge of a roughly 100 km wide overlap where habitat is predicted to be suitable for both species, with more westerly areas suitable only for spinosa and eastward habitats suitable only for J. jacana.

Conclusion

The two species of New World jacanas maintain a narrow, and persistent hybrid zone in western Panama. The hybrid zone may be maintained by the behavioral dominance of J. spinosa counterbalanced by unsuitable habitat for J. spinosa east of the contact zone. Although the two parental species are relatively young, mitochondrial cline width was extremely narrow. This result suggests strong selection against maternally-inherited markers, which may indicate either mitonuclear incompatibilities and/or female choice against heterospecific matings typical of avian hybrid systems, despite jacana sex role reversal.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-014-0227-7) contains supplementary material, which is available to authorized users.

Whole plastome sequencing reveals deep plastid divergence and cytonuclear discordance between closely related balsam poplars, Populus balsamifera and P. trichocarpa (Salicaceae)

As molecular phylogenetic analyses incorporate ever-greater numbers of loci, cases of cytonuclear discordance - the phenomenon in which nuclear gene trees deviate significantly from organellar gene trees - are being reported more frequently. Plant examples of topological discordance, caused by recent hybridization between extant species, are well known. However, examples of branch-length discordance are less reported in plants relative to animals. We use a combination of de novo assembly and reference-based mapping using short-read shotgun sequences to construct a robust phylogeny of the plastome for multiple individuals of all the common Populus species in North America. We demonstrate a case of strikingly high plastome divergence, in contrast to little nuclear genome divergence, in two closely related balsam poplars, Populus balsamifera and Populus trichocarpa (Populus balsamifera ssp. trichocarpa). Previous studies with nuclear loci indicate that the two species (or subspecies) diverged since the late Pleistocene, whereas their plastomes indicate deep divergence, dating to at least the Pliocene (6-7 Myr ago). Our finding is in marked contrast to the estimated Pleistocene divergence of the nuclear genomes, previously calculated at 75 000 yr ago, suggesting plastid capture from a 'ghost lineage' of a now-extinct North American poplar.

Hybridization and mitochondrial genome introgression between Rana chensinensis and R. kukunoris

Mitochondrial genome (mito-genome) introgression among metazoans is commonplace, and several biological processes may promote such introgression. We examined two proposed processes for the mito-genome introgression between Rana chensinensis and R. kukunoris: natural hybridization and sex-biased dispersal. We sampled 477 individuals from 28 sites in the potential hybrid zone in the western Tsinling Mountains. Mitochondrial gene (cyt-b) trees were used to examine the introgression events. Microsatellite DNA loci, cyt-b and morphological data were used to identify hybrids and to examine the extent of natural hybridization. We detected rampant bidirectional introgressions, both ancient and recent, between the two species. Furthermore, we found a wide hybrid zone, and frequent and asymmetric hybridization. The hybrid zone cline analysis revealed a clear mitochondrial-nuclear discordance; while most nuclear markers displayed similar and steep clines, cyt-b had a displaced cline centre and a more gradual and wider cline. We also detected strong and asymmetric historical maternal gene flow across the hybrid zone. This widespread hybridization and detected low mito-nuclear conflicts may, at least partially, explain the high frequency of introgression. Lastly, microsatellite data and population genetic methods were used to assess sex-biased dispersal. A weak pattern of female-biased dispersal was detected in both species, suggesting it may not play an important role in the observed introgression. Our data are consistent with the hybridization hypothesis, but support for the sex-biased dispersal hypothesis is weak. We further suggest that selective advantages of the R. kukunoris-type mito-genome in thermal adaptation may also contribute to the introgression between the two species.

Mitochondrial sequences reveal a clear separation between Angolan and South African giraffe along a cryptic rift valley

BACKGROUND

The current taxonomy of the African giraffe (Giraffa camelopardalis) is primarily based on pelage pattern and geographic distribution, and nine subspecies are currently recognized. Although genetic studies have been conducted, their resolution is low, mainly due to limited sampling. Detailed knowledge about the genetic variation and phylogeography of the South African giraffe (G. c. giraffa) and the Angolan giraffe (G. c. angolensis) is lacking. We investigate genetic variation among giraffe matrilines by increased sampling, with a focus on giraffe key areas in southern Africa.

RESULTS

The 1,562 nucleotides long mitochondrial DNA dataset (cytochrome b and partial control region) comprises 138 parsimony informative sites among 161 giraffe individuals from eight populations. We additionally included two okapis as an outgroup. The analyses of the maternally inherited sequences reveal a deep divergence between northern and southern giraffe populations in Africa, and a general pattern of distinct matrilineal clades corresponding to their geographic distribution. Divergence time estimates among giraffe populations place the deepest splits at several hundred thousand years ago.

CONCLUSIONS

Our increased sampling in southern Africa suggests that the distribution ranges of the Angolan and South African giraffe need to be redefined. Knowledge about the phylogeography and genetic variation of these two maternal lineages is crucial for the development of appropriate management strategies.

On the adaptive value of cytoplasmic genomes in plants

Is DNA variation maintained in organelle genomes selectively neutral? The answer to this question has important implications for many aspects of ecology and evolution. While traditionally the answer has been 'yes', recent studies in animals have shown that, on the contrary, mitochondrial DNA polymorphism is frequently adaptive. In plants, however, the neutrality assumption has not been strongly challenged. Here, we begin with a critical evaluation of arguments in favour of this long-held view. We then discuss the latest empirical evidence for the opposing prediction that sequence variation in plant cytoplasmic genomes is frequently adaptive. While outstanding research progress is being made towards understanding this fundamental topic, we highlight the need for studies that combine information ranging from field experiments to physiology to molecular evolutionary biology. Such an interdisciplinary approach provides a means for determining the frequency, drivers and evolutionary significance of adaptive organelle DNA variation.

Temporal relationship between genetic and warning signal variation in the aposematic wood tiger moth (Parasemia plantaginis)

Many plants and animals advertise unpalatability through warning signals in the form of colour and shape. Variation in warning signals within local populations is not expected because they are subject to directional selection. However, mounting evidence of warning signal variation within local populations suggests that other selective forces may be acting. Moreover, different selective pressures may act on the individual components of a warning signal. At present, we have a limited understanding about how multiple selection processes operate simultaneously on warning signal components, and even less about their temporal and spatial dynamics. Here, we examined temporal variation of several wing warning signal components (colour, UV-reflectance, signal size and pattern) of two co-occurring colour morphs of the aposematic wood tiger moth (Parasemia plantaginis). Sampling was carried out in four geographical regions over three consecutive years. We also evaluated each morph's temporal genetic structure by analysing mitochondrial sequence data and nuclear microsatellite markers. Our results revealed temporal differences between the morphs for most signal components measured. Moreover, variation occurred differently in the fore- and hindwings. We found no differences in the genetic structure between the morphs within years and regions, suggesting single local populations. However, local genetic structure fluctuated temporally. Negative correlations were found between variation produced by neutrally evolving genetic markers and those of the different signal components, indicating a non-neutral evolution for most warning signal components. Taken together, our results suggest that differential selection on warning signal components and fluctuating population structure can be one explanation for the maintenance of warning signal variation in this aposematic species.

Extensive unidirectional introgression between two salamander lineages of ancient divergence and its evolutionary implications

Hybridization and introgression, contrary to previous beliefs, are now considered to be widespread processes even among animal species. Nonetheless, the range of their possible outcomes and roles in moulding biodiversity patterns are still far from being fully appraised. Here we investigated the pattern of hybridization and introgression between Salamandrina perspicillata and S. terdigitata, two salamanders endemic to the Italian peninsula. Using a set of diagnostic or differentiated genetic markers (9 nuclear and 1 mitochondrial), we documented extensive unidirectional introgression of S. terdigitata alleles into the S. perspicillata gene pool in central Italy, indicating that barriers against hybridization were permeable when they came into secondary contact, and despite their ancient divergence. Nonetheless, purebred S. terdigitata, as well as F1, F2, and backcrosses were not found within the hybrid zone. Moreover, Bayesian analyses of population structure identified admixed populations belonging to a differentiated gene pool with respect to both parental populations. Overall, the observed genetic structure, together with their geographic pattern of distribution, suggests that Salamandrina populations in central Italy could have entered a distinct evolutionary pathway. How far they have gone along this pathway will deserve future investigation.

Weak population structure in European roe deer (Capreolus capreolus) and evidence of introgressive hybridization with Siberian roe deer (C. pygargus) in northeastern Poland

We investigated contemporary and historical influences on the pattern of genetic diversity of European roe deer (Capreolus capreolus). The study was conducted in northeastern Poland, a zone where vast areas of primeval forests are conserved and where the European roe deer was never driven to extinction. A total of 319 unique samples collected in three sampling areas were genotyped at 16 microsatellites and one fragment (610 bp) of mitochondrial DNA (mtDNA) control region. Genetic diversity was high, and a low degree of genetic differentiation among sampling areas was observed with both microsatellites and mtDNA. No evidence of genetic differentiation between roe deer inhabiting open fields and forested areas was found, indicating that the ability of the species to exploit these contrasting environments might be the result of its phenotypic plasticity. Half of the studied individuals carried an mtDNA haplotype that did not belong to C. capreolus, but to a related species that does not occur naturally in the area, the Siberian roe deer (C. pygargus). No differentiation between individuals with Siberian and European mtDNA haplotypes was detected at microsatellite loci. Introgression of mtDNA of Siberian roe deer into the genome of European roe deer has recently been detected in eastern Europe. Such introgression might be caused by human-mediated translocations of Siberian roe deer within the range of European roe deer or by natural hybridization between these species in the past.

Are genes faster than crabs? Mitochondrial introgression exceeds larval dispersal during population expansion of the invasive crab Carcinus maenas

Biological invasions offer unique opportunities to investigate evolutionary dynamics at the peripheries of expanding populations. Here, we examine genetic patterns associated with admixture between two distinct invasive lineages of the European green crab, Carcinus maenas L., independently introduced to the northwest Atlantic. Previous investigations based on mitochondrial DNA sequences demonstrated that larval dispersal driven by advective currents could explain observed southward displacement of an admixture zone between the two invasions. Comparison of published mitochondrial results with new nuclear data from nine microsatellite loci, however, reveals striking discordance in their introgression patterns. Specifically, introgression of mitochondrial genomes relative to nuclear background suggests that demographic processes such as sex-biased reproductive dynamics and population size imbalances-and not solely larval dispersal-play an important role in driving the evolution of the genetic cline. In particular, the unpredicted introgression of mitochondrial alleles against the direction of mean larval dispersal in the region is consistent with recent models invoking similar demographic processes to explain movements of genes into invading populations. These observations have important implications for understanding historical shifts in C. maenas range limits, and more generally for inferences of larval dispersal based on genetic data.

Evolutionary neutrality of mtDNA introgression: evidence from complete mitogenome analysis in roe deer

Introgressive hybridization offers a unique platform for studying the molecular basis of natural selection acting on mitogenomes. Most of the mtDNA protein-coding genes are extremely conserved; however, some of the observed variations have potentially adaptive significance. Here, we evaluated whether the evolution of mtDNA in closely related roe deer species affected by widespread mtDNA introgression is neutral or adaptive. We characterized and compared 16 complete mitogenomes of European (Capreolus capreolus) and Siberian (C. pygargus) roe deer, including four of Siberian origin introgressed into European species. The average sequence divergence of species-specific lineages was estimated at 2.8% and varied across gene classes. Only 21 of 315 fixed differences identified in protein-coding genes represented nonsynonymous changes. Only three of them were determined to have arisen in the C. pygargus lineage since the time to the most recent common ancestor (TMRCA) of both Capreolus species, reflecting a decelerated evolutionary ratio. The almost four-fold higher dN /dS ratio described for the European roe deer lineage is constrained by overall purifying selection, especially pronounced in the ND4 and ND5 genes. We suggest that the highly divergent C. capreolus lineage could have maintained a capability for genomic incorporation of the well-preserved and almost ancestral type of mtDNA present in C. pygargus. Our analyses did not indicate any signs of positive selection for Siberian roe deer mtDNA, suggesting that the present widespread introgression is evolutionarily neutral.

A mitogenomic phylogeny and genetic history of sable (Martes zibellina)

We assessed phylogeny of sable (Martes zibellina, Linnaeus, 1758) by sequence analysis of nearly complete, new mitochondrial genomes in 36 specimens from different localities in northern Eurasia (Primorye, Khabarovsk and Krasnoyarsk regions, the Kamchatka Peninsula, the Kuril Islands and the Urals). Phylogenetic analysis of mtDNA sequences demonstrates that two clades, A and BC, radiated about 200-300 thousandyears ago (kya) according to results of Bayesian molecular clock and RelTime analyses of different mitogenome alignments (nearly complete mtDNA sequences, protein-coding region, and synonymous sites), while the age estimates of clades A, B and C fall within the Late Pleistocene (~50-140 kya). Bayesian skyline plots (BSPs) of sable population size change based on analysis of nearly complete mtDNAs show an expansion around 40 kya in the warm Karganian time, without a decline of population size around the Last Glacial Maximum (21 kya). The BSPs based on synonymous clock rate indicate that M. zibellina experienced demographic expansions later, approximately 22 kya. The A2a clade that colonized Kamchatka ~23-50 kya (depending on the mutation rate used) survived the last glaciation there as demonstrated by the BSP analysis. In addition, we have found evidence of positive selection acting at ND4 and cytochrome b genes, thereby suggesting adaptive evolution of the A2a clade in Kamchatka.

Distinguishing between incomplete lineage sorting and genomic introgressions: complete fixation of allospecific mitochondrial DNA in a sexually reproducing fish (Cobitis; Teleostei), despite clonal reproduction of hybrids

Distinguishing between hybrid introgression and incomplete lineage sorting causing incongruence among gene trees in that they exhibit topological differences requires application of statistical approaches that are based on biologically relevant models. Such study is especially challenging in hybrid systems, where usual vectors mediating interspecific gene transfers--hybrids with Mendelian heredity--are absent or unknown. Here we study a complex of hybridizing species, which are known to produce clonal hybrids, to discover how one of the species, Cobitis tanaitica, has achieved a pattern of mito-nuclear mosaic genome over the whole geographic range. We appplied three distinct methods, including the method using solely the information on gene tree topologies, and found that the contrasting mito-nuclear signal might not have resulted from the retention of ancestral polymorphism. Instead, we found two signs of hybridization events related to C. tanaitica; one concerning nuclear gene flow and the other suggested mitochondrial capture. Interestingly, clonal inheritance (gynogenesis) of contemporary hybrids prevents genomic introgressions and non-clonal hybrids are either absent or too rare to be detected among European Cobitis. Our analyses therefore suggest that introgressive hybridizations are rather old episodes, mediated by previously existing hybrids whose inheritance was not entirely clonal. Cobitis complex thus supports the view that the type of resulting hybrids depends on a level of genomic divergence between sexual species.

Speciation in Western Scrub-Jays, Haldane's rule, and genetic clines in secondary contact

Background

Haldane’s Rule, the tendency for the heterogametic sex to show reduced fertility in hybrid crosses, can obscure the signal of gene flow in mtDNA between species where females are heterogametic. Therefore, it is important when studying speciation and species limits in female-heterogametic species like birds to assess the signature of gene flow in the nuclear genome as well. We studied introgression of microsatellites and mtDNA across a secondary contact zone between coastal and interior lineages of Western Scrub-Jays (Aphelocoma californica) to test for a signature of Haldane’s Rule: a narrower cline of introgression in mtDNA compared to nuclear markers.

Results

Our initial phylogeographic analysis revealed that there is only one major area of contact between coastal and interior lineages and identified five genetic clusters with strong spatial structuring: Pacific Slope, Interior US, Edwards Plateau (Texas), Northern Mexico, and Southern Mexico. Consistent with predictions from Haldane’s Rule, mtDNA showed a narrower cline than nuclear markers across a transect through the hybrid zone. This result is not being driven by female-biased dispersal because neutral diffusion analysis, which included estimates of sex-specific dispersal rates, also showed less diffusion of mtDNA. Lineage-specific plumage traits were associated with nuclear genetic profiles for individuals in the hybrid zone, indicating that these differences are under genetic control.

Conclusions

This study adds to a growing list of studies that support predictions of Haldane’s Rule using cline analysis of multiple loci of differing inheritance modes, although alternate hypotheses like selection on different mtDNA types cannot be ruled out. That Haldane’s Rule appears to be operating in this system suggests a measure of reproductive isolation between the Pacific Slope and interior lineages. Based on a variety of evidence from the phenotype, ecology, and genetics, we recommend elevating three lineages to species level: A. californica (Pacific Slope); A. woodhouseii (Interior US plus Edwards Plateau plus Northern Mexico); A. sumichrasti (Southern Mexico). The distinctive Edwards Plateau population in Texas, which was monophyletic in mtDNA except for one individual, should be studied in greater detail given habitat threat.

Mito-nuclear discord in six congeneric lineages of Holarctic ducks (genus Anas)

Many species have Holarctic distributions that extend across Europe, Asia and North America. Most genetics research on these species has examined only mitochondrial (mt) DNA, which has revealed wide variance in divergence between Old World (OW) and New World (NW) populations, ranging from shallow, unstructured genealogies to deeply divergent lineages. In this study, we sequenced 20 nuclear introns to test for concordant patterns of OW-NW differentiation between mtDNA and nuclear (nu) DNA for six lineages of Holarctic ducks (genus Anas). Genetic differentiation for both marker types varied widely among these lineages (idiosyncratic population histories), but mtDNA and nuDNA divergence within lineages was not significantly correlated. Moreover, compared with the association between mtDNA and nuDNA divergence observed among different species, OW-NW nuDNA differentiation was generally lower than mtDNA divergence, at least for lineages with deeply divergent mtDNA. Furthermore, coalescent estimates indicated significantly higher rates of gene flow for nuDNA than mtDNA for four of the six lineages. Thus, Holarctic ducks show prominent mito-nuclear discord between OW and NW populations, and we reject differences in sorting rates as the sole cause of the within-species discord. Male-mediated intercontinental gene flow is likely a leading contributor to this discord, although selection could also cause increased mtDNA divergence relative to weak nuDNA differentiation. The population genetics of these ducks contribute to growing evidence that mtDNA can be an unreliable indicator of stage of speciation and that more holistic approaches are needed for species delimitation.

Speciation, population structure, and demographic history of the Mojave Fringe-toed Lizard (Uma scoparia), a species of conservation concern

The North American deserts were impacted by both Neogene plate tectonics and Quaternary climatic fluctuations, yet it remains unclear how these events influenced speciation in this region. We tested published hypotheses regarding the timing and mode of speciation, population structure, and demographic history of the Mojave Fringe-toed Lizard (Uma scoparia), a sand dune specialist endemic to the Mojave Desert of California and Arizona. We sampled 109 individual lizards representing 22 insular dune localities, obtained DNA sequences for 14 nuclear loci, and found that U. scoparia has low genetic diversity relative to the U. notata species complex, comparable to that of chimpanzees and southern elephant seals. Analyses of genotypes using Bayesian clustering algorithms did not identify discrete populations within U. scoparia. Using isolation-with-migration (IM) models and a novel coalescent-based hypothesis testing approach, we estimated that U. scoparia diverged from U. notata in the Pleistocene epoch. The likelihood ratio test and the Akaike Information Criterion consistently rejected nested speciation models that included parameters for migration and population growth of U. scoparia. We reject the Neogene vicariance hypothesis for the speciation of U. scoparia and define this species as a single evolutionarily significant unit for conservation purposes.

What triggers the rising of an intraspecific biodiversity hotspot? Hints from the agile frog

Hotspots of genetic diversity are regions of utmost importance for species survival and conservation, and their intimate link with the geographic location of glacial refugia has been well established. Nonetheless, the microevolutionary processes underlying the generation of hotspots in such regions have only recently become a fervent field of research. We investigated the phylogeographic and population genetic structure of the agile frog, Rana dalmatina, within its putative refugium in peninsular Italy. We found this region to harbour far more diversity, phylogeographic structure, and lineages of ancient origin than that by the rest of the species' range in Europe. This pattern appeared to be well explained by climate-driven microevolutionary processes that occurred during both glacial and interglacial epochs. Therefore, the inferred evolutionary history of R. dalmatina in Italy supports a view of glacial refugia as 'factories' rather than as repositories of genetic diversity, with significant implications for conservation strategies for hotspots.

Evolutionary basis of mitonuclear discordance between sister species of mole salamanders (Ambystoma sp.)

Distinct genetic markers should show similar patterns of differentiation between species reflecting their common evolutionary histories, yet there are increasing examples of differences in the biogeographic distribution of species-specific nuclear (nuDNA) and mitochondrial DNA (mtDNA) variants within and between species. Identifying the evolutionary processes that underlie these anomalous patterns of genetic differentiation is an important goal. Here, we analyse the putative mitonuclear discordance observed between sister species of mole salamanders (Ambystoma barbouri and A. texanum) in which A. barbouri-specific mtDNA is found in animals located within the range of A. texanum. We test three hypotheses for this discordance (undetected range expansion, mtDNA introgression, and hybridization) using nuDNA and mtDNA data analysed with methods that varied in the parameters estimated and the timescales measured. Results from a Bayesian clustering technique (structure), bidirectional estimates of gene flow (migrate-n and IMa2) and phylogeny-based methods (*beast, bucky) all support the conclusion that the discordance is due to geographically restricted mtDNA introgression from A. barbouri into A. texanum. Limited data on species-specific tooth morphology match this conclusion. Significant differences in environmental conditions exist between sites where A. texanum with and without A. barbouri-like mtDNA occur, suggesting a possible role for selection in the process of introgression. Overall, our study provides a general example of the value of using complimentary analyses to make inferences of the directionality, timescale, and source of mtDNA introgression in animals.

Extensive gene flow characterizes the phylogeography of a North American migrant bird: Black-headed Grosbeak (Pheucticus melanocephalus)

We describe range-wide phylogeographic variation in the Black-headed Grosbeak (Pheucticus melanocephalus), a songbird that is widely distributed across North American scrublands and forests. Phylogenetic analysis of mitochondrial DNA (mtDNA, n=424) revealed three geographically structured clades. One widespread clade occurs throughout the Rocky Mountains, Great Basin, and Mexican Plateau, a second clade is found on the Pacific coast and in coastal ranges; and, a third in the Sierra Madre del Sur of Oaxaca and Guerrero. Some geographical structuring occurs in Mexican Plateau and Sierra Madre Oriental mtDNA clade, presumably because these populations have been more stable over time than northern populations. Multiple mitochondrial groups are found sympatrically in the Okanogan River Valley in Washington, the eastern Sierra Nevada, and the Transvolcanic Belt across central Mexico, indicating that there is a potential for introgression. Analyses of 12 nuclear loci did not recover the same geographically structured clades. Population analyses show high levels of gene flow in nucDNA from the Interior into the Sierra Madre del Sur and Pacific population groups, possibly indicating expansion of the Interior population at the expense of peripheral populations.

Past hybridization between two East Asian long-tailed tits (Aegithalos bonvaloti and A. fuliginosus)

Introduction

Incomplete lineage sorting and hybridization are two major nonexclusive causes of haplotype sharing between species. Distinguishing between these two processes is notoriously difficult as they can generate similar genetic signatures. Previous studies revealed that the mitochondrial DNA (mtDNA) differentiation between two East Asian long-tailed tits (Aegithalos bonvaloti and A. fuliginosus) was extremely low, even lower than intraspecific differentiation in some other long-tailed tits. Using a combination of multilocus and coalescent analyses, we explored the causes of the anomalous lack of mtDNA differentiation between the two species.

Results

The mtDNA divergence between the two species was shallow, while the nuclear DNA (nuDNA) divergence was considerably deeper. The IMa analyses based on the mtDNA dataset suggested relatively high gene flow from A. fuliginosus to A. bonvaloti, while negligible gene flow in the opposite direction. In contrast to mtDNA, the migration rates at autosomal and Z-linked nuDNA loci were negligible or much lower. The NEWHYBRIDS analysis assigned all individuals except one to pure parental species with high posterior probability. The Bayesian skyline plot showed that both species underwent population expansions during the Last Glacial Maximum (LGM), and the ecological niche modelling suggested that their ranges overlapped more during the LGM than at present.

Conclusions

We suggest that historical hybridization, in combination with selective sweep and/or genetic drift might be the main causes of the extremely low mtDNA differentiation between the two species. The hybridization probably occurred mainly between A. fuliginosus females and A. bonvaloti males. The LGM distribution expansion might have facilitated hybridization, while the post-LGM distribution contraction could have facilitated some mtDNA sorting. Ongoing hybridization between the two species might be very limited, but further studies with more samples from the contact zone are needed to test this conclusion.

Analyzing mosquito (Diptera: culicidae) diversity in Pakistan by DNA barcoding

BACKGROUND

Although they are important disease vectors mosquito biodiversity in Pakistan is poorly known. Recent epidemics of dengue fever have revealed the need for more detailed understanding of the diversity and distributions of mosquito species in this region. DNA barcoding improves the accuracy of mosquito inventories because morphological differences between many species are subtle, leading to misidentifications.

METHODOLOGY/PRINCIPAL FINDINGS

Sequence variation in the barcode region of the mitochondrial COI gene was used to identify mosquito species, reveal genetic diversity, and map the distribution of the dengue-vector species in Pakistan. Analysis of 1684 mosquitoes from 491 sites in Punjab and Khyber Pakhtunkhwa during 2010-2013 revealed 32 species with the assemblage dominated by Culex quinquefasciatus (61% of the collection). The genus Aedes (Stegomyia) comprised 15% of the specimens, and was represented by six taxa with the two dengue vector species, Ae. albopictus and Ae. aegypti, dominant and broadly distributed. Anopheles made up another 6% of the catch with An. subpictus dominating. Barcode sequence divergence in conspecific specimens ranged from 0-2.4%, while congeneric species showed from 2.3-17.8% divergence. A global haplotype analysis of disease-vectors showed the presence of multiple haplotypes, although a single haplotype of each dengue-vector species was dominant in most countries. Geographic distribution of Ae. aegypti and Ae. albopictus showed the later species was dominant and found in both rural and urban environments.

CONCLUSIONS

As the first DNA-based analysis of mosquitoes in Pakistan, this study has begun the construction of a barcode reference library for the mosquitoes of this region. Levels of genetic diversity varied among species. Because of its capacity to differentiate species, even those with subtle morphological differences, DNA barcoding aids accurate tracking of vector populations.

Divergence with gene flow within the recent chipmunk radiation (Tamias)

Increasing data have supported the importance of divergence with gene flow (DGF) in the generation of biological diversity. In such cases, lineage divergence occurs on a shorter timescale than does the completion of reproductive isolation. Although it is critical to explore the mechanisms driving divergence and preventing homogenization by hybridization, it is equally important to document cases of DGF in nature. Here we synthesize data that have accumulated over the last dozen or so years on DGF in the chipmunk (Tamias) radiation with new data that quantify very high rates of mitochondrial DNA (mtDNA) introgression among para- and sympatric species in the T. quadrivittatus group in the central and southern Rocky Mountains. These new data (188 cytochrome b sequences) bring the total number of sequences up to 1871; roughly 16% (298) of the chipmunks we have sequenced exhibit introgressed mtDNA. This includes ongoing introgression between subspecies and between both closely related and distantly related taxa. In addition, we have identified several taxa that are apparently fixed for ancient introgressions and in which there is no evidence of ongoing introgression. A recurrent observation is that these introgressions occur between ecologically and morphologically diverged, sometimes non-sister taxa that engage in well-documented niche partitioning. Thus, the chipmunk radiation in western North America represents an excellent mammalian example of speciation in the face of recurrent gene flow among lineages and where biogeography, habitat differentiation and mating systems suggest important roles for both ecological and sexual selection.

Unexpected population genetic structure of European roe deer in Poland: an invasion of the mtDNA genome from Siberian roe deer

Introgressive hybridization is a widespread evolutionary phenomenon which may lead to increased allelic variation at selective neutral loci and to transfer of fitness-related traits to introgressed lineages. We inferred the population genetic structure of the European roe deer (Capreolus capreolus) in Poland from mitochondrial (CR and cyt b) and sex-linked markers (ZFX, SRY, DBY4 and DBY8). Analyses of CR mtDNA sequences from 452 individuals indicated widespread introgression of Siberian roe deer (C. pygargus) mtDNA in the European roe deer genome, 2000 km from the current distribution range of C. pygargus. Introgressed individuals constituted 16.6% of the deer studied. Nearly 75% of them possessed haplotypes belonging to the group which arose 23 kyr ago and have not been detected within the natural range of Siberian roe deer, indicating that majority of present introgression has ancient origin. Unlike the mtDNA results, sex-specific markers did not show signs of introgression. Species distribution modelling analyses suggested that C. pygargus could have extended its range as far west as Central Europe after last glacial maximum. The main hybridization event was probably associated with range expansion of the most abundant European roe deer lineage from western refugia and took place in Central Europe after the Younger Dryas (10.8-10.0 ka BP). Initially, introgressed mtDNA variants could have spread out on the wave of expansion through the mechanism of gene surfing, reaching high frequencies in European roe deer populations and leading to observed asymmetrical gene flow. Human-mediated introductions of C. pygargus had minimal effect on the extent of mtDNA introgression.

Exploring the effect of asymmetric mitochondrial DNA introgression on estimating niche divergence in morphologically cryptic species

If potential morphologically cryptic species, identified based on differentiated mitochondrial DNA, express ecological divergence, this increases support for their treatment as distinct species. However, mitochondrial DNA introgression hampers the correct estimation of ecological divergence. We test the hypothesis that estimated niche divergence differs when considering nuclear DNA composition or mitochondrial DNA type as representing the true species range. We use empirical data of two crested newt species (Amphibia: Triturus) which possess introgressed mitochondrial DNA from a third species in part of their ranges. We analyze the data in environmental space by determining Fisher distances in a principal component analysis and in geographical space by determining geographical overlap of species distribution models. We find that under mtDNA guidance in one of the two study cases niche divergence is overestimated, whereas in the other it is underestimated. In the light of our results we discuss the role of estimated niche divergence in species delineation.

Molecular and ecological signs of mitochondrial adaptation: consequences for introgression?

The evolution of the mitochondrial genome and its potential adaptive impact still generates vital debates. Even if mitochondria have a crucial functional role, as they are the main cellular energy suppliers, mitochondrial DNA (mtDNA) introgression is common in nature, introducing variation in populations upon which selection may act. Here we evaluated whether the evolution of mtDNA in a rodent species affected by mtDNA introgression is explained by neutral expectations alone. Variation in one mitochondrial and six nuclear markers in Myodes glareolus voles was examined, including populations that show mtDNA introgression from its close relative, Myodes rutilus. In addition, we modelled protein structures of the mtDNA marker (cytochrome b) and estimated the environmental envelopes of mitotypes. We found that massive mtDNA introgression occurred without any trace of introgression in the analysed nuclear genes. The results show that the native glareolus mtDNA evolved under past positive selection, suggesting that mtDNA in this system has selective relevance. The environmental models indicate that the rutilus mitotype inhabits colder and drier habitats than the glareolus one that can result from local adaptation or from the geographic context of introgression. Finally, homology models of the cytochrome b protein revealed a substitution in rutilus mtDNA in the vicinity of the catalytic fraction, suggesting that differences between mitotypes may result in functional changes. These results suggest that the evolution of mtDNA in Myodes may have functional, ecological and adaptive significance. This work opens perspective onto future experimental tests of the role of natural selection in mtDNA introgression in this system.