MITOCHONDRIAL INTROGRESSION AND INCOMPLETE LINEAGE SORTING THROUGH SPACE AND TIME: PHYLOGENETICS OF CROTAPHYTID LIZARDS

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.

Snake and bird predation drive the repeated convergent evolution of correlated life history traits and phenotype in the Izu Island Scincid lizard (Plestiodon latiscutatus)

Predation may create strong natural selection pressure on the phenotype and life history characteristics of prey species. The Izu scincid lizards (Plestiodon latiscutatus) that inhabit the four Japanese Izu Islands with only bird predators are drab brown, mature later, lay small clutches of large eggs, and hatch large neonates. In contrast, skinks on seven islands with both snake and bird predators are conspicuously colored, mature early, lay large clutches of small eggs, and hatch small neonates. We test the hypothesis that these suites of traits have evolved independently on each island via natural selection pressures from one of two predator regimes--birds-only and birds + snakes. Using two mtDNA genes and a nuclear locus, we infer a time-calibrated phylogeny of P. latiscutatus that reveals a basal split between Mikura and all islands south, and Miyake, all islands north, and the Izu Peninsula. Populations inhabiting Miyake, Niijima, Shikine, and Toshima are not monophyletic, suggesting either multiple colonizations or an artifact of incomplete lineage sorting (ILS). We therefore developed novel phylogenetic comparative analyses that assume either a multiple colonization or more restrictive single colonization ILS scenario and found 1) statistically significant support for the of different suites of phenotypic and life history characteristics with the presence of bird-only or bird + snake predator assemblages, and 2) strong phylogenetic support for at least two independent derivations of either the "bird-only" or "snakes + birds" phenotypes regardless of colonization scenario. Finally, our time-calibrated phylogeographic analysis supports the conclusion that the ancestor to modern Izu Island P. latiscutatus dispersed from the mainland to the Izu proto-islands between 3-7.6 million years ago (Ma). These lineages remained present in the area during successive formation of the islands, with one lineage re-colonizing the mainland 0.24-0.7 Ma.

Reproductive isolation between phylogeographic lineages scales with divergence

Phylogeographic studies frequently reveal multiple morphologically cryptic lineages within species. What is not yet clear is whether such lineages represent nascent species or evolutionary ephemera. To address this question, we compare five contact zones, each of which occurs between ecomorphologically cryptic lineages of skinks from the rainforests of the Australian Wet Tropics. Although the contacts probably formed concurrently in response to Holocene expansion from glacial refugia, we estimate that the divergence times (τ) of the lineage pairs range from 3.1 to 11.5 Ma. Multi-locus analyses of the contact zones yielded estimates of reproductive isolation that are tightly correlated with divergence time and, for lineages with older divergence times (τ > 5 Myr), substantial. These results show that phylogeographic splits of increasing depth represent stages along the speciation continuum, even in the absence of overt change in ecologically relevant morphology.

Molecular systematics and plumage coloration evolution of an enigmatic babbler (Pomatorhinus ruficollis) in East Asia

The streak-breasted scimitar babbler, Pomatorhinus ruficollis, is a polytypic and taxonomically enigmatic babbler common in southern, eastern, and southeastern Asia. To infer the phylogeny of the P. ruficollis, we examined the sequences of two complete mitochondrial genes (2184 bp in total) from fourteen of the fifteen known subspecies, and an additional five nuclear genes (2657 bp in total) from ten subspecies. The mitochondrial phylogeny indicates four major clades with large geographical identity in P. ruficollis and paraphyly of the P. ruficollis species complex, with the inclusion of the olivaceus group of congeneric P. schisticeps. Together with their interbreeding in northern Indochina, we propose to lump this group into P. ruficollis. Analysis of both multilocus networks and species-tree inference recovered poor phylogenetic structure among mainland/ Hainan subspecies and exclusive groupings of the Taiwanese subspecies, consistent with the recent taxonomic revision of its species status. Our analyses also suggest strong incongruence between the morphological-based classification and molecular systematics, implying the strength of multilocus data for taxonomy.

Multilocus phylogeography of the common lizard Zootoca vivipara at the Ibero-Pyrenean suture zone reveals lowland barriers and high-elevation introgression

BACKGROUND

The geographic distribution of evolutionary lineages and the patterns of gene flow upon secondary contact provide insight into the process of divergence and speciation. We explore the evolutionary history of the common lizard Zootoca vivipara (= Lacerta vivipara) in the Iberian Peninsula and test the role of the Pyrenees and the Cantabrian Mountains in restricting gene flow and driving lineage isolation and divergence. We also assess patterns of introgression among lineages upon secondary contact, and test for the role of high-elevation trans-mountain colonisations in explaining spatial patterns of genetic diversity. We use mtDNA sequence data and genome-wide AFLP loci to reconstruct phylogenetic relationships among lineages, and measure genetic structure.

RESULTS

The main genetic split in mtDNA corresponds generally to the French and Spanish sides of the Pyrenees as previously reported, in contrast to genome-wide AFLP data, which show a major division between NW Spain and the rest. Both types of markers support the existence of four distinct and geographically congruent genetic groups, which are consistent with major topographic barriers. Both datasets reveal the presence of three independent contact zones between lineages in the Pyrenean region, one in the Basque lowlands, one in the low-elevation mountains of the western Pyrenees, and one in the French side of the central Pyrenees. The latter shows genetic evidence of a recent, high-altitude trans-Pyrenean incursion from Spain into France.

CONCLUSIONS

The distribution and age of major lineages is consistent with a Pleistocene origin and a role for both the Pyrenees and the Cantabrian Mountains in driving isolation and differentiation of Z. vivipara lineages at large geographic scales. However, mountain ranges are not always effective barriers to dispersal, and have not prevented a recent high-elevation trans-Pyrenean incursion that has led to asymmetrical introgression among divergent lineages. Cytonuclear discordance in patterns of genetic structure and introgression at contact zones suggests selection may be involved at various scales. Suture zones are important areas for the study of lineage formation and speciation, and our results show that biogeographic barriers can yield markedly different phylogeographic patterns in different vertebrate and invertebrate taxa.

Polyphyly and hidden species among Hawai'i's dominant mesophotic coral genera, Leptoseris and Pavona (Scleractinia: Agariciidae)

Widespread polyphyly in stony corals (order Scleractinia) has prompted efforts to revise their systematics through approaches that integrate molecular and micromorphological evidence. To date, these approaches have not been comprehensively applied to the dominant genera in mesophotic coral ecosystems (MCEs) because several species in these genera occur primarily at depths that are poorly explored and from which sample collections are limited. This study is the first integrated morphological and molecular systematic analysis of the genera Leptoseris and Pavona to examine material both from shallow-water reefs (<30 m) and from mid- to lower-MCEs (>60 m). Skeletal and tissue samples were collected throughout the Hawaiian Archipelago between 2–127 m. A novel mitochondrial marker (cox1-1-rRNA intron) was sequenced for 70 colonies, and the micromorphologies of 94 skeletons, plus selected type material, were analyzed. The cox1-1-rRNA intron resolved 8 clades, yet Leptoseris and Pavona were polyphyletic. Skeletal micromorphology, especially costal ornamentation, showed strong correspondence and discrete differences between mitochondrial groups. One putative new Leptoseris species was identified and the global depth range of the genus Pavona was extended to 89 m, suggesting that the diversity of mesophotic scleractinians has been underestimated. Examination of species’ depth distributions revealed a pattern of depth zonation: Species common in shallow-water were absent or rare >40 m, whereas others occurred only >60 m. These patterns emphasize the importance of integrated systematic analyses and more comprehensive sampling by depth in assessing the connectivity and diversity of MCEs.

Assessing phylogenetic relationships among galliformes: a multigene phylogeny with expanded taxon sampling in Phasianidae

Galliform birds (relatives of the chicken and turkey) have attracted substantial attention due to their importance to society and value as model systems. This makes understanding the evolutionary history of Galliformes, especially the species-rich family Phasianidae, particularly interesting and important for comparative studies in this group. Previous studies have differed in their conclusions regarding galliform phylogeny. Some of these studies have suggested that specific clades within this order underwent rapid radiations, potentially leading to the observed difficulty in resolving their phylogenetic relationships. Here we presented analyses of six nuclear intron sequences and two mitochondrial regions, an amount of sequence data larger than many previous studies, and expanded taxon sampling by collecting data from 88 galliform species and four anseriform outgroups. Our results corroborated recent studies describing relationships among the major families, and provided further evidence that the traditional division of the largest family, the Phasianidae into two major groups ("pheasants" and "partridges") is not valid. Within the Phasianidae, relationships among many genera have varied among studies and there has been little consensus for the placement of many taxa. Using this large dataset, with substantial sampling within the Phasianidae, we obtained strong bootstrap support to confirm some previously hypothesized relationships and we were able to exclude others. In addition, we added the first nuclear sequence data for the partridge and quail genera Ammoperdix, Caloperdix, Excalfactoria, and Margaroperdix, placing these taxa in the galliform tree of life with confidence. Despite the novel insights obtained by combining increased sampling of taxa and loci, our results suggest that additional data collection will be necessary to solve the remaining uncertainties.

Genetic drift or natural selection? Hybridization and asymmetric mitochondrial introgression in two Caribbean lizards (Anolis pulchellus and Anolis krugi)

Hybridization and gene introgression can occur frequently between closely related taxa, but appear to be rare phenomena among members of the species-rich West Indian radiation of Anolis lizards. We investigated the pattern and possible mechanism of introgression between two sister species from Puerto Rico, Anolis pulchellus and Anolis krugi, using mitochondrial (ND2) and nuclear (DNAH3, NKTR) DNA sequences. Our findings demonstrated extensive introgression of A. krugi mtDNA (k-mtDNA) into the genome of A. pulchellus in western Puerto Rico, to the extent that k-mtDNA has mostly or completely replaced the native mtDNA of A. pulchellus on this part of the island. We proposed two not mutually exclusive scenarios to account for the interspecific matings between A. pulchellus and A. krugi. We inferred that hybridization events occurred independently in several populations, and determined that k-mtDNA haplotypes harboured in individuals of A. pulchellus can be assigned to four of the five major mtDNA clades of A. krugi. Further, the spatial distribution of k-mtDNA clades in the two species is largely congruent. Based on this evidence, we concluded that natural selection was the probable driving mechanism for the extensive k-mtDNA introgression into A. pulchellus. Our two nuclear data sets yielded different results. DNAH3 showed reciprocal monophyly of A. pulchellus and A. krugi, indicating no effect of hybridization on this marker. In contrast, the two species shared nine NKTR alleles, probably due to incomplete lineage sorting. Our study system will provide an excellent opportunity to experimentally assess the behavioural and ecological mechanisms that can lead to hybridization in closely related taxa.

Tracing glacial refugia of Triturus newts based on mitochondrial DNA phylogeography and species distribution modeling

INTRODUCTION

The major climatic oscillations during the Quaternary Ice Age heavily influenced the distribution of species and left their mark on intraspecific genetic diversity. Past range shifts can be reconstructed with the aid of species distribution modeling and phylogeographical analyses. We test the responses of the different members of the genus Triturus (i.e. the marbled and crested newts) as the climate shifted from the previous glacial period (the Last Glacial Maximum, ~21 Ka) to the current interglacial.

RESULTS

We present the results of a dense mitochondrial DNA phylogeography (visualizing genetic diversity within and divergence among populations) and species distribution modeling (using two different climate simulations) for the nine Triturus species on composite maps.

CONCLUSIONS

The combined use of species distribution modeling and mitochondrial phylogeography provides insight in the glacial contraction and postglacial expansion of Triturus. The combined use of the two independent techniques yields a more complete understanding of the historical biogeography of Triturus than both approaches would on their own. Triturus newts generally conform to the 'southern richness and northern purity' paradigm, but we also find more intricate patterns, such as the absence of genetic variation and suitable area at the Last Glacial Maximum (T. dobrogicus), an 'extra-Mediterranean' refugium in the Carpathian Basin (T. cristatus), and areas where species displaced one another postglacially (e.g. T. macedonicus and western T. karelinii). We provide a biogeographical scenario for Triturus, showing the positions of glacial refugia, the regions that were postglacially colonized and the areas where species displaced one another as they shifted their ranges.

Testing hypotheses for genealogical discordance in a rainforest lizard

Genealogical discordance, or when different genes tell distinct stories although they evolved under a shared history, often emerges from either coalescent stochasticity or introgression. In this study, we present a strong case of mito-nuclear genealogical discordance in the Australian rainforest lizard species complex of Saproscincus basiliscus and S. lewisi. One of the lineages that comprises this complex, the Southern S. basiliscus lineage, is deeply divergent at the mitochondrial genome but shows markedly less divergence at the nuclear genome. By placing our results in a comparative context and reconstructing the lineages' demography via multilocus and coalescent-based approximate Bayesian computation methods, we test hypotheses for how coalescent variance and introgression contribute to this pattern. These analyses suggest that the observed genealogical discordance likely results from introgression. Further, to generate such strong discordance, introgression probably acted in concert with other factors promoting asymmetric gene flow between the mitochondrial and nuclear genomes, such as selection or sex-biased dispersal. This study offers a framework for testing sources of genealogical discordance and suggests that historical introgression can be an important force shaping the genetic diversity of species and their populations.

Life history and biogeographic diversification of an endemic western North American freshwater fish clade using a comparative species tree approach

The west coast of North America contains a number of biogeographic freshwater provinces which reflect an ever-changing aquatic landscape. Clues to understanding this complex structure are often encapsulated genetically in the ichthyofauna, though frequently as unresolved evolutionary relationships and putative cryptic species. Advances in molecular phylogenetics through species tree analyses now allow for improved exploration of these relationships. Using a comprehensive approach, we analyzed two mitochondrial and nine nuclear loci for a group of endemic freshwater fish (sculpin-Cottus) known for a wide ranging distribution and complex species structure in this region. Species delimitation techniques identified three novel cryptic lineages, all well supported by phylogenetic analyses. Comparative phylogenetic analyses consistently found five distinct clades reflecting a number of unique biogeographic provinces. Some internal node relationships varied by species tree reconstruction method, and were associated with either Bayesian or maximum likelihood statistical approaches or between mitochondrial, nuclear, and combined datasets. Limited cases of mitochondrial capture were also evident, suggestive of putative ancestral hybridization between species. Biogeographic diversification was associated with four major regions and revealed historical faunal exchanges across regions. Mapping of an important life-history character (amphidromy) revealed two separate instances of trait evolution, a transition that has occurred repeatedly in Cottus. This study demonstrates the power of current phylogenetic methods, the need for a comprehensive phylogenetic approach, and the potential for sculpin to serve as an indicator of biogeographic history for native ichthyofauna in the region.

Postglacial species displacement in Triturus newts deduced from asymmetrically introgressed mitochondrial DNA and ecological niche models

BACKGROUND

If the geographical displacement of one species by another is accompanied by hybridization, mitochondrial DNA can introgress asymmetrically, from the outcompeted species into the invading species, over a large area. We explore this phenomenon using the two parapatric crested newt species, Triturus macedonicus and T. karelinii, distributed on the Balkan Peninsula in south-eastern Europe, as a model.

RESULTS

We first delimit a ca. 54,000 km(2) area in which T. macedonicus contains T. karelinii mitochondrial DNA. This introgression zone bisects the range of T. karelinii, cutting off a T. karelinii enclave. The high similarity of introgressed mitochondrial DNA haplotypes with those found in T. karelinii suggests a recent transfer across the species boundary. We then use ecological niche modeling to explore habitat suitability of the location of the present day introgression zone under current, mid-Holocene and Last Glacial Maximum conditions. This area was inhospitable during the Last Glacial Maximum for both species, but would have been habitable at the mid-Holocene. Since the mid-Holocene, habitat suitability generally increased for T. macedonicus, whereas it decreased for T. karelinii.

CONCLUSION

The presence of a T. karelinii enclave suggests that T. karelinii was the first to colonize the area where the present day introgression zone is positioned after the Last Glacial Maximum. Subsequently, we propose T. karelinii was outcompeted by T. macedonicus, which captured T. karelinii mitochondrial DNA via introgressive hybridization in the process. Ecological niche modeling suggests that this replacement was likely facilitated by a shift in climate since the mid-Holocene. We suggest that the northwestern part of the current introgression zone was probably never inhabited by T. karelinii itself, and that T. karelinii mitochondrial DNA spread there through T. macedonicus exclusively. Considering the spatial distribution of the introgressed mitochondrial DNA and the signal derived from ecological niche modeling, we do not favor the hypothesis that foreign mitochondrial DNA was pulled into the T. macedonicus range by natural selection.

Multiple mitochondrial introgression events and heteroplasmy in trypanosoma cruzi revealed by maxicircle MLST and next generation sequencing

Background

Mitochondrial DNA is a valuable taxonomic marker due to its relatively fast rate of evolution. In Trypanosoma cruzi, the causative agent of Chagas disease, the mitochondrial genome has a unique structural organization consisting of 20–50 maxicircles (∼20 kb) and thousands of minicircles (0.5–10 kb). T. cruzi is an early diverging protist displaying remarkable genetic heterogeneity and is recognized as a complex of six discrete typing units (DTUs). The majority of infected humans are asymptomatic for life while 30–35% develop potentially fatal cardiac and/or digestive syndromes. However, the relationship between specific clinical outcomes and T. cruzi genotype remains elusive. The availability of whole genome sequences has driven advances in high resolution genotyping techniques and re-invigorated interest in exploring the diversity present within the various DTUs.

Methodology/Principal Findings

To describe intra-DTU diversity, we developed a highly resolutive maxicircle multilocus sequence typing (mtMLST) scheme based on ten gene fragments. A panel of 32 TcI isolates was genotyped using the mtMLST scheme, GPI, mini-exon and 25 microsatellite loci. Comparison of nuclear and mitochondrial data revealed clearly incongruent phylogenetic histories among different geographical populations as well as major DTUs. In parallel, we exploited read depth data, generated by Illumina sequencing of the maxicircle genome from the TcI reference strain Sylvio X10/1, to provide the first evidence of mitochondrial heteroplasmy (heterogeneous mitochondrial genomes in an individual cell) in T. cruzi.

Conclusions/Significance

mtMLST provides a powerful approach to genotyping at the sub-DTU level. This strategy will facilitate attempts to resolve phenotypic variation in T. cruzi and to address epidemiologically important hypotheses in conjunction with intensive spatio-temporal sampling. The observations of both general and specific incidences of nuclear-mitochondrial phylogenetic incongruence indicate that genetic recombination is geographically widespread and continues to influence the natural population structure of TcI, a conclusion which challenges the traditional paradigm of clonality in T. cruzi.

Chagas disease, caused by the protozoan parasite Trypanosoma cruzi, is an important public health problem in Latin America. While molecular techniques can differentiate the major T. cruzi genetic lineages, few have sufficient resolution to describe diversity among closely related strains. The online availability of three mitochondrial genomes allowed us to design a multilocus sequence typing (mtMLST) scheme to exploit these rapidly evolving markers. We compared mtMLST with current nuclear typing tools using isolates belonging to the oldest and most widely occurring lineage TcI. T. cruzi is generally believed to reproduce clonally. However, in this study, distinct branching patterns between mitochondrial and nuclear phylogenetic trees revealed multiple incidences of genetic exchange within different geographical populations and major lineages. We also examined Illumina sequencing data from the TcI genome strain which revealed multiple different mitochondrial genomes within an individual parasite (heteroplasmy) that were, however, not sufficiently divergent to represent a major source of typing error. We strongly recommend this combined nuclear and mitochondrial genotyping methodology to reveal cryptic diversity and genetic exchange in T. cruzi. The level of resolution that this mtMLST provides should greatly assist attempts to elucidate the complex interactions between parasite genotype, clinical outcome and disease distribution.

Direction and extent of organelle DNA introgression between two spruce species in the Qinghai-Tibetan Plateau

A recent model has shown that, during range expansion of one species in a territory already occupied by a related species, introgression should take place preferentially from the resident species towards the invading species and genome components experiencing low rates of gene flow should introgress more readily than those experiencing high rates of gene flow. Here, we use molecular markers from two organelle genomes with contrasted rates of gene flow to test these predictions by examining genetic exchanges between two morphologically distinct spruce Picea species growing in the Qinghai-Tibetan Plateau. The haplotypes from both mitochondrial (mt) DNA and chloroplast (cp) DNA cluster into two distinct lineages that differentiate allopatric populations of the two species. By contrast, in sympatry, the species share the same haplotypes, suggesting interspecific genetic exchanges. As predicted by the neutral model, all sympatric populations of the expanding species had received their maternally inherited mtDNA from the resident species, whereas for paternally inherited cpDNA introgression is more limited and not strictly unidirectional. Our results underscore cryptic introgressions of organelle DNAs in plants and the importance of considering rates of gene flow and range shifts to predict direction and extent of interspecific genetic exchanges.

Bonobos fall within the genomic variation of chimpanzees

To gain insight into the patterns of genetic variation and evolutionary relationships within and between bonobos and chimpanzees, we sequenced 150,000 base pairs of nuclear DNA divided among 15 autosomal regions as well as the complete mitochondrial genomes from 20 bonobos and 58 chimpanzees. Except for western chimpanzees, we found poor genetic separation of chimpanzees based on sample locality. In contrast, bonobos consistently cluster together but fall as a group within the variation of chimpanzees for many of the regions. Thus, while chimpanzees retain genomic variation that predates bonobo-chimpanzee speciation, extensive lineage sorting has occurred within bonobos such that much of their genome traces its ancestry back to a single common ancestor that postdates their origin as a group separate from chimpanzees.

Phylogeography and demographic history of Lacerta lepida in the Iberian Peninsula: multiple refugia, range expansions and secondary contact zones

BACKGROUND

The Iberian Peninsula is recognized as an important refugial area for species survival and diversification during the climatic cycles of the Quaternary. Recent phylogeographic studies have revealed Iberia as a complex of multiple refugia. However, most of these studies have focused either on species with narrow distributions within the region or species groups that, although widely distributed, generally have a genetic structure that relates to pre-Quaternary cladogenetic events. In this study we undertake a detailed phylogeographic analysis of the lizard species, Lacerta lepida, whose distribution encompasses the entire Iberian Peninsula. We attempt to identify refugial areas, recolonization routes, zones of secondary contact and date demographic events within this species.

RESULTS

Results support the existence of 6 evolutionary lineages (phylogroups) with a strong association between genetic variation and geography, suggesting a history of allopatric divergence in different refugia. Diversification within phylogroups is concordant with the onset of the Pleistocene climatic oscillations. The southern regions of several phylogroups show a high incidence of ancestral alleles in contrast with high incidence of recently derived alleles in northern regions. All phylogroups show signs of recent demographic and spatial expansions. We have further identified several zones of secondary contact, with divergent mitochondrial haplotypes occurring in narrow zones of sympatry.

CONCLUSIONS

The concordant patterns of spatial and demographic expansions detected within phylogroups, together with the high incidence of ancestral haplotypes in southern regions of several phylogroups, suggests a pattern of contraction of populations into southern refugia during adverse climatic conditions from which subsequent northern expansions occurred. This study supports the emergent pattern of multiple refugia within Iberia but adds to it by identifying a pattern of refugia coincident with the southern distribution limits of individual evolutionary lineages. These areas are important in terms of long-term species persistence and therefore important areas for conservation.

Birth of a hotspot of intraspecific genetic diversity: notes from the underground

Hotspots of intraspecific diversity have been observed in most species, often within areas of putative Pleistocene refugia. They have thus mostly been viewed as the outcome of prolonged stability of large populations within the refugia. However, recent evidence has suggested that several other microevolutionary processes could also be involved in their formation. Here, we investigate the contribution of these processes to current range-wide patterns of genetic diversity in the Italian endemic mole Talpa romana, using both nuclear (30 allozyme loci) and mitochondrial markers (cytochrome b sequences). Southern populations of this species showed an allozyme variation that is amongst the highest observed in small mammals (most populations had an expected heterozygosity of 0.10 or above), which was particularly unexpected for a subterranean species. Population genetic, phylogeographic and historical demographic analyses indicated that T. romana populations repeatedly underwent allopatric differentiations followed by secondary admixture within the refugial range in southern Italy. A prolonged demographic stability was reliably inferred from the mitochondrial DNA data only for a population group located north and east of the Calabrian peninsula, showing comparatively lower levels of allozyme variability, and lacking evidence of secondary admixture with other groups. Thus, our results point to the admixture between differentiated lineages as the main cause of the higher levels of diversity of refugial populations. When compared with the Pleistocene evolutionary history recently inferred for species from both the same and other geographic regions, these results suggest the need for a reappraisal of the role of gene exchange in the formation of intraspecific hotspots of genetic diversity.

A young clade repeating an old pattern: diversity in Nothonotus darters (Teleostei: Percidae) endemic to the Cumberland River

Hypotheses of diversification in eastern North American freshwater fishes have focused primarily on allopatric distributions of species between disjunct highland areas and major river systems. However, these hypotheses do not fully explain the rich diversity of species within highland regions and river systems. Relatively old diversification events at small geographic scales have been observed in the Barcheek Darter subclade that occurs in the Cumberland River drainage (CRD) in Kentucky and Tennessee, United States of America, but it is unknown if this pattern is consistent in other darter subclades. We explored phylogeographic diversity in two species of Nothonotus darters, N. microlepidus and N. sanguifluus, endemic to the CRD to compare phylogenetic patterns between Barcheek Darters and species of Nothonotus. We collected sequence data for a mitochondrial gene (cytb) and three nuclear genes (MLL, S7 and RAG1) from 19 N. microlepidus and 35 N. sanguifluus specimens. Gene trees were estimated using maximum likelihood and Bayesian methods, and a 'species tree' was inferred using a Bayesian method. These trees indicate that species diversity in Nothonotus is underestimated. Five distinct lineages were evident, despite retained ancestral polymorphism and unsampled extirpated populations. Comparison of chronograms for Barcheek Darters and Nothonotus revealed that microendemism resulting from species diversification at small geographic scales in the CRD is a consistent pattern in both old and young darter subclades. Our analyses reveal that geographic isolating mechanisms that result in similar phylogeographic patterns in the CRD are persistent through long expanses of evolutionary time.