Phylogeographic inferences from chloroplast DNA: quantifying the effects of mutations in repetitive and non-repetitive sequences

Phylogeographic inference can be a powerful tool in reconstructing species’ evolutionary histories; however, although inferred phylogeographic patterns should depend in part on the underlying types and rates of mutations, the effects of different types of mutations have seldom been quantified. In this study we identified two chloroplast minisatellites in the common reed Phragmites australis, and showed that these are more variable than chloroplast microsatellites. We then recreated parsimony networks of the global phylogeography of P. australis based on data that either included or excluded repetitive sequences (minisatellites and microsatellites), thereby illustrating the influence that these repetitive sequences can have on large‐scale phylogeographic inference. The resulting networks differed in the numbers of mutational steps, degrees of uncertainty, and total numbers of haplotypes. In addition, the suggested ancestor‐descendant relationships among lineages changed substantially depending on whether repetitive sequences were included. We therefore caution against the inclusion of repetitive sequences in large‐scale networks because of their high potential for homoplasy. Nevertheless, we advocate the inclusion of repetitive sequences in other analyses: specifically, we show that the ratio of mutations in repetitive vs. non‐repetitive regions can provide insight into the relative ages of lineages.     

Patterns of gene flow and population genetic structure in the canyon treefrog, Hyla arenicolor (Cope)

Patterns of gene flow and genetic structuring were examined in the canyon treefrog, Hyla arenicolor (Cope). Hierarchical analysis of genetic variation was performed on mitochondrial cytochrome b haplotypes from 323 individuals, representing 32 populations from previously described phylogeographic regions. Results from AMOVA revealed that 60.4-78.9% of the recovered genetic variation was the result of differences in the appointment of genetic variation between subdivisions of the primary phylogeographic regions. In contrast, populations only contained between 13.9 and 30.1% of the observed haplotypic variation. Gene flow estimates based on calculations of phi ST revealed moderate levels of gene flow within phylogeographic regions, but there was no evidence of gene flow between these regions, suggesting that geographical boundaries were probably important in the formation of phylogeographic structure in H. arenicolor. Phylogeographic regions exhibited very different patterns of gene flow. One region showed evidence of recent colonization. Another region exhibited very limited gene flow. Moderate to high estimates of gene flow were obtained for populations from two distinct phylogeographic regions characterized by mesic and xeric environments. Isolation by distance was observed in both regions suggesting that these regions are in genetic equilibrium. Because gene flow is extremely unlikely between the populations in the xeric region, this result is interpreted as historical gene flow. These results indicate that isolation-by-distance effects may still be observed even when population genetic structure and gene flow are the result of historical association.     

Excessive variation in Y chromosomal DNA in Rumex acetosa (Polygonaceae)

Rumex acetosa is one of the few angiosperms that possesses sex chromosomes. The same types of abundant repetitive sequences cover both heterochromatic Y chromosomes present in males. The aim of this study was to investigate genetic variation in paternally inherited Y chromosomal DNA and in maternally inherited cpDNA, and to find out whether the examined genomic regions are suited to a phylogeographic study in R. acetosa. DNA sequence polymorphisms present in the 850-bp heterochromatic segment on the Y chromosomes were compared to variation in the 409-bp long chloroplast section (trnL- trnF spacer) in R. acetosa originating from several European locations and from the Altai mountains in Russia. A great amount of genetic variation was detected within the Y chromosomal region while only four chloroplast genotypes were detected. Although the chloroplast haplotypes possessed some geographic pattern, no clear phylogeographic pattern was detected based on the variable Y chromosomes. The mean Y chromosomal nucleotide diversity among all samples equaled 6.6 %, and the mean proportion of polymorphic sites per individual equaled 8.2 % among SNP sites and 1.7 % among all sites investigated. The high number of substitutions detected in the Y chromosomal DNA shows that this heterochromatic sequence has a high mutation rate. The diversity pattern indicates that gene flow via pollen is extensive and it blurs any geographical pattern in the Y chromosomal variation. The high number of repeats and uncertainty concerning the extent of recombination between the two Y chromosomes impair the usability of the Y chromosomal segment for phylogeographic or population genetic studies.     

Effects of vicariant barriers, habitat stability, population isolation and environmental features on species divergence in the south-western Australian coastal reptile community

Identifying explicit hypotheses regarding the factors determining genetic structuring within species can be difficult, especially in species distributed in historically dynamic regions. To contend with these challenges, we use a framework that combines species distribution models, environmental data and multi-locus genetic data to generate and explore phylogeographic hypotheses for reptile species occupying the coastal sand-dune and sand-plain habitats of the south-western Australian biodiversity hotspot, a community which has both a high diversity of endemics and has varied dramatically in spatial extent over time. We use hierarchical amova, summary statistic and distance-based analyses to explicitly test specific phylogeographic hypotheses. Namely, we test if biogeographic vicariance across barriers, habitat stability, population isolation along a linear habitat or fragmentation across different environments can explain genetic divergence within five co-distributed squamate reptile species. Our results show that patterns of genetic variation reflect complex and species-specific interactions related to the spatial distribution of habitats present currently and during repeated glacial minima, as opposed to being associated with historical factors such as habitat stability between glacial and inter-glacial periods or vicariant barriers. We suggest that the large impact of habitat characteristics over time (i.e. relative levels of habitat connectivity, climatic gradients and spatial heterogeneity of soil types) reflects the ecological restrictions of the sand-dune and sand-plain reptile communities and may explain the lack of concordance across taxa. The study demonstrates the general utility of the approach for assemblage-level, as well as single species, phylogeographic study, including its usefulness for exploring biologically informed hypotheses about what factors have influenced patterns of genetic variation.     

Comparison of whole chloroplast genome sequences to choose noncoding regions for phylogenetic studies in angiosperms: the tortoise and the hare III

Although the chloroplast genome contains many noncoding regions, relatively few have been exploited for interspecific phylogenetic and intraspecific phylogeographic studies. In our recent evaluation of the phylogenetic utility of 21 noncoding chloroplast regions, we found the most widely used noncoding regions are among the least variable, but the more variable regions have rarely been employed. That study led us to conclude that there may be unexplored regions of the chloroplast genome that have even higher relative levels of variability. To explore the potential variability of previously unexplored regions, we compared three pairs of single-copy chloroplast genome sequences in three disparate angiosperm lineages: Atropa vs. Nicotiana (asterids); Lotus vs. Medicago (rosids); and Saccharum vs. Oryza (monocots). These three separate sequence alignments highlighted 13 mutational hotspots that may be more variable than the best regions of our former study. These 13 regions were then selected for a more detailed analysis. Here we show that nine of these newly explored regions (rpl32-trnL((UAG)), trnQ((UUG))-5'rps16, 3'trnV((UAC))-ndhC, ndhF-rpl32, psbD-trnT((GGU)), psbJ-petA, 3'rps16-5'trnK((UUU)), atpI-atpH, and petL-psbE) offer levels of variation better than the best regions identified in our earlier study and are therefore likely to be the best choices for molecular studies at low taxonomic levels.     

Structure and Evolution of the Mitochondrial Control Region of Leaf Beetles (Coleoptera: Chrysomelidae): A Hierarchical Analysis of Nucleotide Sequence Variation

To assess the levels of variation at different evolutionary scales in the mitochondrial (mt) control region of leaf beetles, we sequenced and compared the full mt control region in two genera ( Chrysomela and Gonioctena), in two species within a genus ( Gonioctena olivacea and G. pallida), in individuals from distant populations of these species in Europe, and in individuals from populations separated by moderate (10- to 100-km) to short (<5-km) distances. In all individuals, a highly repetitive section consisting of the tandem repetition of 12 to 17 imperfect copies of a 107- to 159-bp-long core sequence was observed. This repetitive fragment accounts for roughly 50% of the full control-region length. The sequence variability among repeated elements within the control region of a given individual depends on the species considered: the variability within any G. olivacea individual is much higher than that within G. pallida individuals. Comparisons of the repeated elements, in a phylogenetic framework, within and among individuals of G. olivacea and G. pallida suggests that the repetitive section of the control region experienced recurrent duplications/deletions, leading to some degree of concerted evolution. Comparisons between Chrysomela and Gonioctena control regions revealed virtually no significant sequence similarity, except for two long stretches of A's and several [T(T)A(A)] repeats, all found in the control region of other insect orders. Our analyses allowed us to identify portions of the control region with enough variation for population genetic or phylogeographic studies.     

Comparative avian phylogeography of Cameroon and equatorial Guinea mountains: implications for conservation

We illustrate the use of Faith's 'Phylogenetic Diversity' measure to compare the phylogeographic structure of two bird species with patterns of avian endemism across six mountains in Cameroon and Equatorial Guinea. The Mountain Greenbul and Cameroon Blue-headed Sunbird showed phylogeographic patterns that together defined three biogeographic regions: Bioko, Mt. Cameroon, and the northern mountains of Cameroon. In contrast, the distributions of endemic species were largely a function of geographical distance, with close mountains sharing more endemic species than distant mountains. Moreover, for both species, populations on Mt. Cameroon were distinctive with respect to the ecologically relevant character bill size. Our results, while preliminary, illustrate the utility of a comparative approach for identifying geographical regions that harbour evolutionarily distinct populations and caution against using only the distributional patterns of endemics to prioritize regions for conservation. Results show that patterns of endemism may not be concordant with patterns of phylogenetic diversity nor morphological variation in a character important in fitness. While incorporation of additional species from unrelated taxa will be necessary to draw definitive conclusions about evolutionarily distinct regions, our preliminary results suggest a conservation approach for the Afromontane region of the Gulf of Guinea that would: (i) emphasize protection of both Bioko and Mt. Cameroon, thereby maximizing preservation of within-species phylogenetic and morphologic diversity; (ii) emphasize protection within the northern mountains to further conserve intraspecific phylogenetic diversity and maximize protection of endemic species.     

Fat frogs,mobile genes: unexpected phylogeographic patterns for the ornate chorus frog (Pseudacris ornata)

The southeastern coastal plain of the United States is a region marked by extraordinary phylogeographic congruence that is frequently attributed to the changing sea levels that occurred during the glacial‐interglacial cycles of the Pleistocene epoch. A phylogeographic break corresponding to the Apalachicola River has been suggested for many species studied to date that are endemic to this region. Here, we used this pattern of phylogeographic congruence to develop and test explicit hypotheses about the genetic structure in the ornate chorus frog (Pseudacris ornata). Using 1299 bp of mtDNA sequence and seven nuclear microsatellite markers in 13 natural populations of P. ornata, we found three clades corresponding to geographically distinct regions; one spans the Apalachicola River (Southern Clade), one encompasses Georgia and South Carolina (Central Clade) and a third comprises more northerly individuals (Northern Clade). However, it does not appear that typical phylogeographic barriers demarcate these clades. Instead, isolation by distance across the range of the entire species explained the pattern of genetic variation that we observed. We propose that P. ornata was historically widespread in the southeastern United States, and that a balance between genetic drift and migration was the root of the genetic divergence among populations. Additionally, we investigated fine‐scale patterns of genetic structure and found the spatial scale at which there was significant genetic structure varied among the regions studied. Furthermore, we discuss our results in light of other phylogeographic studies of southeastern coastal plain organisms and in relation to amphibian conservation and management.     

Intraspecific sequence variation of chloroplast DNA among the component species of evergreen broad-leaved forests in Japan

For the purpose of phylogeographic study of lucidophyllous (evergreen broad-leaved) forests in Japan, we surveyed intraspecific chloroplast DNA (cpDNA) variation in 41 component species of such forests. Intraspecific cpDNA variations were detected in 14 species. In 15 species and one species group, 16 non-coding cpDNA regions were examined to find intraspecific sequence variation. The extent of variation in these regions was compared. The largest amount of intraspecific variation was detected in the rps16 region. A relatively large amount of intraspecific variation was detected in the petD-rpoA, rpl16, and trnL-F regions. It is suggested that these regions of cpDNA would be useful for detecting intraspecific variation in plant species, and could provide valuable information for various research purposes.     

Phylogeography and spatial genetic structure of the Southern torrent salamander: implications for conservation and management

The Southern torrent salamander (Rhyacotriton variegatus) was recently found not warranted for listing under the US Endangered Species Act due to lack of information regarding population fragmentation and gene flow. Found in small-order streams associated with late-successional coniferous forests of the US Pacific Northwest, threats to their persistence include disturbance related to timber harvest activities. We conducted a study of genetic diversity throughout this species' range to 1) identify major phylogenetic lineages and phylogeographic barriers and 2) elucidate regional patterns of population genetic and spatial phylogeographic structure. Cytochrome b sequence variation was examined for 189 individuals from 72 localities. We identified 3 major lineages corresponding to nonoverlapping geographic regions: a northern California clade, a central Oregon clade, and a northern Oregon clade. The Yaquina River may be a phylogeographic barrier between the northern Oregon and central Oregon clades, whereas the Smith River in northern California appears to correspond to the discontinuity between the central Oregon and northern California clades. Spatial analyses of genetic variation within regions encompassing major clades indicated that the extent of genetic structure is comparable among regions. We discuss our results in the context of conservation efforts for Southern torrent salamanders.     

Phylogeographic signal variation in mitochondrial DNA among geographically isolated grassland butterflies

Aim We surveyed mitochondrial DNA (mtDNA) sequence variation among regionally isolated populations of 10 grassland‐associated butterfly species to determine: (1) the utility of phylogeographic comparisons among multiple species for assessing recent evolutionary patterns, and (2) the respective roles of isolation attributable to range disjunction versus isolation attributable solely to geographic distance in establishing divergence patterns. Location The Peace River grasslands of northern Alberta and British Columbia, Canada, which are isolated by 300+ km from similar communities to the south. Methods We sequenced mtDNA (1420 bp of cytochrome c oxidase subunit I) from five grassland‐restricted butterfly species that have geographically disjunct populations and from five ecologically broader species that have more continuous distributions across the same regions. Using analysis of molecular variance (AMOVA), Mantel and partial Mantel tests, and haplotype networks, we compared population structure within and between species in order to assess the validity of single‐species phylogeographic characterizations. We then contrasted variance components between disjunct and continuously distributed species to assess whether divergences were correlated more with disjunction or with geographic distance. Results Single‐species analyses varied substantially within both the disjunct and the continuous groups. One species in each of these groups had mtDNA with unusually deep intraspecific mitochondrial lineage divergences. On the whole, however, the five species with disjunct ranges exhibited greater divergence between geographically distant populations than did the five species with continuous distributions. Comparison of variance components between disjunct and continuous species indicated that isolation attributable only to geographic distance was responsible for up to half of the total sequence variation between disjunct populations of grassland butterflies. Main conclusions Our findings show that single‐species phylogeographic analyses of post‐Pleistocene butterfly distributions are inadequate for characterizing regional biogeographic divergence histories. However, comparison of mtDNA sequence divergences between groups of disjunct and continuously distributed species can allow isolation attributable to range interruption to be quantitatively distinguished from isolation attributable solely to gene flow attenuation over the same geographic area.     

Interspecies physiological variation as a tool for cross-species assessments of global warming-induced endangerment: validation of an intrinsic determinant of macroecological and phylogeographic structure

Global warming is now recognized as the dominant threat to biodiversity because even protected populations and habitats are susceptible. Nonetheless, current criteria for evaluating species' relative endangerment remain purely ecological, and the accepted conservation strategies of habitat preservation and population management assume that species can mount ecological responses if afforded protection. The insidious threat from climate change is that it will attenuate or preclude ecological responses by species that are physiologically constrained; yet, quantitative, objective criteria for assessing relative susceptibility of diverse taxa to warming-induced stress are wanting. We explored the utility of using interspecies physiological variation for this purpose by relating species' physiological phenotypes to landscape patterns of ecological and genetic exchange. Using a salamander model system in which ecological, genetic and physiological diversity are well characterized, we found strong quantitative relationships of basal metabolic rates (BMRs) to both macroecological and phylogeographic patterns, with decreasing BMR leading to dispersal limitation (small contemporary ranges with marked phylogeographic structure). Measures of intrinsic physiological tolerance, which vary systematically with macroecological and phylogeographic patterns, afford objective criteria for assessing endangerment across a wide range of species and should be incorporated into conservation assessment criteria that currently rely exclusively upon ecological predictors.     

基于叶绿体基因序列rbcL-accD研究莛子藨的分子遗传多样性和谱系地理结构

莛子藨是莛子藨属的一种多年生草本植物,主要分布在青藏高原东北部及邻近高山地区和日本,为东亚特有种。本文基于叶绿体DNA基因间隔区rbcL-accD序列,对莛子藨物种水平的序列特征、分子遗传多样性和谱系地理结构进行研究。分析结果发现该区域莛子藨种群的遗传变异主要来自于居群间（72.28%）,优势单倍型（H2）广泛分布,特有单倍型分散在多个居群。居群间遗传分化系数N_STST,没有检测到显著谱系地理结构,而且居群间基因流较低（N_m=0.096）,青藏高原及其邻近地区特殊地形和生态环境造成的地理隔离或者生境异质化可能是现有单倍型分布格局形成的原因。种群动态历史分析显示莛子藨经历过不显著的种群选择或者扩张事件,综合分析推测莛子藨居群应该是在间冰期或者冰期后从不同避难所扩散形成现在的分布格局,具体避难所还需进一步深入研究才能确认。     

A Passerine Bird's evolution corroborates the geologic history of the island of New Guinea

New Guinea is a biologically diverse island, with a unique geologic history and topography that has likely played a role in the evolution of species. Few island-wide studies, however, have examined the phylogeographic history of lowland species. The objective of this study was to examine patterns of phylogeographic variation of a common and widespread New Guinean bird species (Colluricincla megarhyncha). Specifically, we test the mechanisms hypothesized to cause geographic and genetic variation (e.g., vicariance, isolation by distance and founder-effect with dispersal). To accomplish this, we surveyed three regions of the mitochondrial genome and a nuclear intron and assessed differences among 23 of the 30 described subspecies from throughout their range. We found support for eight highly divergent lineages within C. megarhyncha. Genetic lineages were found within continuous lowland habitat or on smaller islands, but all individuals within clades were not necessarily structured by predicted biogeographic barriers. There was some evidence of isolation by distance and potential founder-effects. Mitochondrial DNA sequence divergence among lineages was at a level often observed among different species or even genera of birds (5-11%), suggesting lineages within regions have been isolated for long periods of time. When topographical barriers were associated with divergence patterns, the estimated divergence date for the clade coincided with the estimated time of barrier formation. We also found that dispersal distance and range size are positively correlated across lineages. Evidence from this research suggests that different phylogeographic mechanisms concurrently structure lineages of C. megarhyncha and are not mutually exclusive. These lineages are a result of evolutionary forces acting at different temporal and spatial scales concordant with New Guinea's geological history.     

Testing the role of ecology and life history in structuring genetic variation across a landscape: a trait‐based phylogeographic approach

Hypotheses to explain phylogeographic structure traditionally invoke geographic features, but often fail to provide a general explanation for spatial patterns of genetic variation. Organisms' intrinsic characteristics might play more important roles than landscape features in determining phylogeographic structure. We developed a novel comparative approach to explore the role of ecological and life‐history variables in determining spatial genetic variation and tested it on frog communities in Panama. We quantified spatial genetic variation within 31 anuran species based on mitochondrial DNA sequences, for which hierarchical approximate Bayesian computation analyses rejected simultaneous divergence over a common landscape. Regressing ecological variables, on genetic divergence allowed us to test the importance of individual variables revealing that body size, current landscape resistance, geographic range, biogeographic origin and reproductive mode were significant predictors of spatial genetic variation. Our results support the idea that phylogeographic structure represents the outcome of an interaction between organisms and their environment, and suggest a conceptual integration we refer to as trait‐based phylogeography.     

Cranial morphological variation of Dromiciops gliroides (Microbiotheria) along its geographical distribution in south-central Chile: A three-dimensional analysis

We analyzed the variation in cranial morphology of the marsupial Dromiciops gliroides along its distribution in south-central Chile. We evaluated whether the cranial morphological variation is congruent with the phylogeographic structure previously observed in this species. We built three-dimensional models of 69 crania on which we digitized 30 landmarks. We used standard geometric morphometric methods to extract and analyze the shape and size components of the crania. Our data showed a subtle but consistent cranial size and shape variation along the studied distributional range, suggesting a geographic variation pattern rather than a phylogeographic structuring. Indeed, our multivariate analyses recovered a subtle morphological differentiation between island and mainland populations, contrary to what is suggested by a former phylogeographic study. We detected that either the cranial size variation, as well as the insularity and the latitude could be important factors underlying the cranial shape changes. We suggest that an interplay of historical and contemporary processes could be shaping the morphological pattern observed in this marsupial.     

Sequence capture and next‐generation sequencing of ultraconserved elements in a large‐genome salamander

Amidst the rapid advancement in next‐generation sequencing (NGS) technology over the last few years, salamanders have been left behind. Salamanders have enormous genomes—up to 40 times the size of the human genome—and this poses challenges to generating NGS data sets of quality and quantity similar to those of other vertebrates. However, optimization of laboratory protocols is time‐consuming and often cost prohibitive, and continued omission of salamanders from novel phylogeographic research is detrimental to species facing decline. Here, we use a salamander endemic to the southeastern United States, Plethodon serratus, to test the utility of an established protocol for sequence capture of ultraconserved elements (UCEs) in resolving intraspecific phylogeographic relationships and delimiting cryptic species. Without modifying the standard laboratory protocol, we generated a data set consisting of over 600 million reads for 85 P. serratus samples. Species delimitation analyses support recognition of seven species within P. serratus sensu lato, and all phylogenetic relationships among the seven species are fully resolved under a coalescent model. Results also corroborate previous data suggesting nonmonophyly of the Ouachita and Louisiana regions. Our results demonstrate that established UCE protocols can successfully be used in phylogeographic studies of salamander species, providing a powerful tool for future research on evolutionary history of amphibians and other organisms with large genomes.     

Structure and Evolution of the Mitochondrial Control Region of Leaf Beetles (Coleoptera: Chrysomelidae): A Hierarchical Analysis of Nucleotide Sequence Variation

Abstract To assess the levels of variation at different evolutionary scales in the mitochondrial (mt) control region of leaf beetles, we sequenced and compared the full mt control region in two genera (Chrysomela and Gonioctena), in two species within a genus (Gonioctena olivacea and G. pallida), in individuals from distant populations of these species in Europe, and in individuals from populations separated by moderate (10- to 100-km) to short (<5-km) distances. In all individuals, a highly repetitive section consisting of the tandem repetition of 12 to 17 imperfect copies of a 107- to 159-bp-long core sequence was observed. This repetitive fragment accounts for roughly 50% of the full control-region length. The sequence variability among repeated elements within the control region of a given individual depends on the species considered: the variability within any G. olivacea individual is much higher than that within G. pallida individuals. Comparisons of the repeated elements, in a phylogenetic framework, within and among individuals of G. olivacea and G. pallida suggests that the repetitive section of the control region experienced recurrent duplications/deletions, leading to some degree of concerted evolution. Comparisons between Chrysomela and Gonioctena control regions revealed virtually no significant sequence similarity, except for two long stretches of A's and several [T(T)A(A)] repeats, all found in the control region of other insect orders. Our analyses allowed us to identify portions of the control region with enough variation for population genetic or phylogeographic studies.     

Genetic variation in the complete MgPa operon and its repetitive chromosomal elements in clinical strains of Mycoplasma genitalium

Mycoplasma genitalium has been increasingly recognized as an important microbe not only because of its significant association with human genital tract diseases but also because of its utility as a model for studying the minimum set of genes necessary to sustain life. Despite its small genome, 4.7% of the total genome sequence is devoted to making the MgPa adhesin operon and its nine chromosomal repetitive elements (termed MgPars). The MgPa operon, along with 9 MgPars, is believed to play an important role in pathogenesis of M. genitalium infection and has also served as the main target for development of diagnostic tools. However, genetic variation in the complete MgPa operon and MgPars among clinical strains of M. genitalium has not been addressed. In this study we examined the genetic variation in the complete MgPa operon (approximately 8.5 kb) and full or partial MgPar sequences (0.4-2.6 kb) in 15 geographically diverse strains of M. genitalium. Extensive variation was present in four repeat regions of the MgPa operon (with homology to MgPars) among and within strains while the non-repeat regions (without homology to MgPars) showed low-level variation among strains and no variation within strains. MgPars showed significant variation among strains but were highly homogeneous within strains, supporting gene conversion as the likely recombination mechanism. When applying our sequence data to evaluate published MgPa operon-based diagnostic PCR assays and genotyping systems, we found that 11 of 19 primers contain up to 19 variable nucleotides and that the target for one of two typing systems is located in a hypervariable repeat region, suggesting the likelihood of false results with some of these assays. This study not only provides new insights into the role of the MgPa operon in the pathogenesis of M. genitalium infection but has important implications for the development of diagnostic tools.     

Phylogeography and introgressive hybridization: chipmunks (genus Tamias) in the northern Rocky Mountains

Abstract. If phylogeographic studies are to be broadly used for assessing population-level processes relevant to speciation and systematics, the ability to identify and incorporate instances of hybridization into the analytical framework is essential. Here, we examine the evolutionary history of two chipmunk species, Tamias ruficaudus and Tamias amoenus , in the northern Rocky Mountains by integrating multivariate morphometrics of bacular (os penis) variation, phylogenetic estimation, and nested clade analysis with regional biogeography. Our results indicate multiple examples of mitochondrial DNA introgression layered within the evolutionary history of these nonsister species. Three of these events are most consistent with recent and/or ongoing asymmetric introgression of mitochondrial DNA across morphologically defined secondary contact zones. In addition, we find preliminary evidence where a fourth instance of nonconcordant characters may represent complete fixation of introgressed mitochondrial DNA via a more ancient hybridization event, although alternative explanations of convergence or incomplete sorting of ancestral polymorphisms cannot be dismissed with these data. The demonstration of hybridization among chipmunks with strongly differentiated bacular morphology contradicts long-standing assumptions that variation within this character is diagnostic of complete reproductive isolation within Tamias . Our results illustrate the utility of phylogeographic analyses for detecting instances of reticulate evolution and for incorporating this and other information in the inference of the evolutionary history of species.