iGLASS: an improvement to the GLASS method for estimating species trees from gene trees

Several methods have been designed to infer species trees from gene trees while taking into account gene tree/species tree discordance. Although some of these methods provide consistent species tree topology estimates under a standard model, most either do not estimate branch lengths or are computationally slow. An exception, the GLASS method of Mossel and Roch, is consistent for the species tree topology, estimates branch lengths, and is computationally fast. However, GLASS systematically overestimates divergence times, leading to biased estimates of species tree branch lengths. By assuming a multispecies coalescent model in which multiple lineages are sampled from each of two taxa at L independent loci, we derive the distribution of the waiting time until the first interspecific coalescence occurs between the two taxa, considering all loci and measuring from the divergence time. We then use the mean of this distribution to derive a correction to the GLASS estimator of pairwise divergence times. We show that our improved estimator, which we call iGLASS, consistently estimates the divergence time between a pair of taxa as the number of loci approaches infinity, and that it is an unbiased estimator of divergence times when one lineage is sampled per taxon. We also show that many commonly used clustering methods can be combined with the iGLASS estimator of pairwise divergence times to produce a consistent estimator of the species tree topology. Through simulations, we show that iGLASS can greatly reduce the bias and mean squared error in obtaining estimates of divergence times in a species tree.     

Phylogeny and evolutionary history of birch mice Sicista Griffith, 1827 (Sminthidae,Rodentia): Implications from a multigene study

Phylogeny of birch mice is estimated using sequences of ten nuclear genes and one mitochondrial gene. Based on the results of tree reconstructions and molecular dating, five major lineages are recognized: “tianschanica,” “concolor,” “caudata,” “betulina,” and “caucasica.” It is established that the three latter lineages constitute a clade and that the long‐tailed birch mouse Sicista caudata is the sister group of the “caucasica” lineage. The “tianschanica” lineage is placed as the sister branch to all other species, however, with insufficient support. The cytochrome b tree is generally concordant with the nuclear topology. The molecular clock results suggest that the radiation among the main lineages occurred in the Late Miocene–Early Pliocene (6.0–4.7 Mya). The correspondence between molecular dating and the fossil record is discussed. Based on nuclear data, a high level of divergence between cryptic species in the “tianschanica” lineage is confirmed. Mitochondrial and nuclear data suggest the existence of a potential cryptic species within Sicista strandi.     

Primate phylogenomics uncovers multiple rapid radiations and ancient interspecific introgression

Our understanding of the evolutionary history of primates is undergoing continual revision due to ongoing genome sequencing efforts. Bolstered by growing fossil evidence, these data have led to increased acceptance of once controversial hypotheses regarding phylogenetic relationships, hybridization and introgression, and the biogeographical history of primate groups. Among these findings is a pattern of recent introgression between species within all major primate groups examined to date, though little is known about introgression deeper in time. To address this and other phylogenetic questions, here, we present new reference genome assemblies for 3 Old World monkey (OWM) species: Colobus angolensis ssp. palliatus (the black and white colobus), Macaca nemestrina (southern pig-tailed macaque), and Mandrillus leucophaeus (the drill). We combine these data with 23 additional primate genomes to estimate both the species tree and individual gene trees using thousands of loci. While our species tree is largely consistent with previous phylogenetic hypotheses, the gene trees reveal high levels of genealogical discordance associated with multiple primate radiations. We use strongly asymmetric patterns of gene tree discordance around specific branches to identify multiple instances of introgression between ancestral primate lineages. In addition, we exploit recent fossil evidence to perform fossil-calibrated molecular dating analyses across the tree. Taken together, our genome-wide data help to resolve multiple contentious sets of relationships among primates, while also providing insight into the biological processes and technical artifacts that led to the disagreements in the first place.

Combining three newly sequenced primate genomes with other published genomes, this study adapts a little-known method for detecting ancient introgression to genome-scale data, revealing multiple previously unknown examples of hybridization between primate species.     

Coalescent histories for discordant gene trees and species trees

Given a gene tree and a species tree, a coalescent history is a list of the branches of the species tree on which coalescences in the gene tree take place. Each pair consisting of a gene tree topology and a species tree topology has some number of possible coalescent histories. Here we show that, for each n≥7, there exist a species tree topology S and a gene tree topology G≠S, both with n leaves, for which the number of coalescent histories exceeds the corresponding number of coalescent histories when the species tree topology is S and the gene tree topology is also S. This result has the interpretation that the gene tree topology G discordant with the species tree topology S can be produced by the evolutionary process in more ways than can the gene tree topology that matches the species tree topology, providing further insight into the surprising combinatorial properties of gene trees that arise from their joint consideration with species trees.     

Understanding crown shyness from a 3-D perspective

Background and AimsCrown shyness describes the phenomenon whereby tree crowns avoid growing into each other, producing a puzzle-like pattern of complementary tree crowns in the canopy. Previous studies found that tree slenderness plays a role in the development of crown shyness. Attempts to quantify crown shyness have largely been confined to 2-D approaches. This study aimed to expand the current set of metrics for crown shyness by quantifying the characteristic of 3-D surface complementarity between trees displaying crown shyness, using LiDAR-derived tree point clouds. Subsequently, the relationship between crown surface complementarity and slenderness of trees was assessed.MethodsFourteen trees were scanned using a laser scanning device. Individual tree points clouds were extracted semi-automatically and manually corrected where needed. A metric that quantifies the surface complementarity (S_c) of a pair of protein molecules is applied to point clouds of pairs of adjacent trees. Then 3-D tree crown surfaces were generated from point clouds by computing their α shapes.Key ResultsTree pairs that were visually determined to have overlapping crowns scored significantly lower S_c values than pairs that did not overlap (n = 14, P < 0.01). Furthermore, average slenderness of pairs of trees correlated positively with their S_c score (R² = 0.484, P < 0.01), showing agreement with previous studies on crown shyness.ConclusionsThe characteristic of crown surface complementarity present in trees displaying crown shyness was succesfully quantified using a 3-D surface complementarity metric adopted from molecular biology. Crown surface complementarity showed a positive relationship to tree slenderness, similar to other metrics used for measuring crown shyness. The 3-D metric developed in this study revealed how trees adapt the shape of their crowns to those of adjacent trees and how this is linked to the slenderness of the trees.     

Evaluating evolutionary history in the face of high gene tree discordance in Australian Gehyra (Reptilia: Gekkonidae)

Species tree methods have provided improvements for estimating species relationships and the timing of diversification in recent radiations by allowing for gene tree discordance. Although gene tree discordance is often observed, most discordance is attributed to incomplete lineage sorting rather than other biological phenomena, and the causes of discordance are rarely investigated. We use species trees from multi-locus data to estimate the species relationships, evolutionary history and timing of diversification among Australian Gehyra—a group renowned for taxonomic uncertainty and showing a large degree of gene tree discordance. We find support for a recent Asian origin and two major clades: a tropically adapted clade and an arid adapted clade, with some exceptions, but no support for allopatric speciation driven by chromosomal rearrangement in the group. Bayesian concordance analysis revealed high gene tree discordance and comparisons of Robinson–Foulds distances showed that discordance between gene trees was significantly higher than that generated by topological uncertainty within each gene. Analysis of gene tree discordance and incomplete taxon sampling revealed that gene tree discordance was high whether terminal taxon or gene sampling was maximized, indicating discordance is due to biological processes, which may be important in contributing to gene tree discordance in many recently diversified organisms.     

Data Concatenation,Bayesian Concordance and Coalescent-Based Analyses of the Species Tree for the Rapid Radiation of Triturus Newts

The phylogenetic relationships for rapid species radiations are difficult to disentangle. Here we study one such case, namely the genus Triturus, which is composed of the marbled and crested newts. We analyze data for 38 genetic markers, positioned in 3-prime untranslated regions of protein-coding genes, obtained with 454 sequencing. Our dataset includes twenty Triturus newts and represents all nine species. Bayesian analysis of population structure allocates all individuals to their respective species. The branching patterns obtained by data concatenation, Bayesian concordance analysis and coalescent-based estimations of the species tree differ from one another. The data concatenation based species tree shows high branch support but branching order is considerably affected by allele choice in the case of heterozygotes in the concatenation process. Bayesian concordance analysis expresses the conflict between individual gene trees for part of the Triturus species tree as low concordance factors. The coalescent-based species tree is relatively similar to a previously published species tree based upon morphology and full mtDNA and any conflicting internal branches are not highly supported. Our findings reflect high gene tree discordance due to incomplete lineage sorting (possibly aggravated by hybridization) in combination with low information content of the markers employed (as can be expected for relatively recent species radiations). This case study highlights the complexity of resolving rapid radiations and we acknowledge that to convincingly resolve the Triturus species tree even more genes will have to be consulted.     

Importance of high quality early-successional habitats in managed forest landscapes to rare beetle species

Species adapted to early-successional forest habitats are in managed landscapes largely confined to clearcuts. To improve habitat quality on clearcuts, green tree and dead wood retention is widely applied in forestry; however, its effects on rare early-successional species have rarely been shown. We repeatedly surveyed two red-listed beetle species (Upis ceramboides and Platysoma minus) on clearcuts in a managed boreal forest landscape. We found that U. ceramboides decreased its occupancy over time while P. minus increased, indicating that red-listed species vary in their ability to successfully utilise managed habitats. We found no effect of connectivity on probability of occurrence, colonisation or extinction per clearcut. Trees retained alive improved habitat quality of clearcuts, since both species were more frequent in dead wood of such trees, in comparison to logging residues. We suggest that retention can be improved by protecting and creating dead wood as intact trees during harvesting. Rare specialist species require habitat of high quality, and consequently it is impossible to meet the requirements of these species on every clearcut. To preserve all early-successional species at a regional scale, we recommend focusing retention of green trees and dead wood to one or a few trees species on each clearcut and in each landscape.     

Incomplete lineage sorting impacts the inference of macroevolutionary regimes from molecular phylogenies when concatenation is employed: An analysis based on Cetacea

Interest in methods that estimate speciation and extinction rates from molecular phylogenies has increased over the last decade. The application of such methods requires reliable estimates of tree topology and node ages, which are frequently obtained using standard phylogenetic inference combining concatenated loci and molecular dating. However, this practice disregards population‐level processes that generate gene tree/species tree discordance. We evaluated the impact of employing concatenation and coalescent‐based phylogeny inference in recovering the correct macroevolutionary regime using simulated data based on the well‐established diversification rate shift of delphinids in Cetacea. We found that under scenarios of strong incomplete lineage sorting, macroevolutionary analysis of phylogenies inferred by concatenating loci failed to recover the delphinid diversification shift, while the coalescent‐based tree consistently retrieved the correct rate regime. We suggest that ignoring microevolutionary processes reduces the power of methods that estimate macroevolutionary regimes from molecular data.     

Ageing culturally significant relic trees in southeast Queensland to support bushfire management strategies

Appropriate fire management strategies are needed to protect forests and large old ecologically and culturally significant trees in natural landscapes. The aim of this study was to determine the age of large old and relic trees of cultural significance that included Cypress Pine (Callitris columellaris F. Muell.), a species that is sensitive to crown scorching fires in a fire‐prone landscape, and to calibrate a tree‐growth‐rate method for estimating tree age. Twelve large trees were dated using radiocarbon (¹⁴C) dating. The trees are located on North Stradbroke Island (Indigenous name: Minjerribah), southeast Queensland (Australia) in a fire‐prone landscape where recent wildfires have destroyed many large trees. The median tree ages ranged from 155 to 369 years. These results suggest an important role of past Indigenous land management practices in protecting Cypress Pine from crown scorching fires. The tree‐growth‐rate‐based method for estimating tree age generally overestimated the age derived from radiocarbon dating. Bias correction factors were developed for correcting various measures of periodic growth rates. This study provides evidence that appropriate low‐intensity fire strategies have the potential to contribute to the survival of forests and conserve large old trees.     

Wood density and its radial variation in six canopy tree species differing in shade-tolerance in western Thailand

Background and Aims

Wood density is a key variable for understanding life history strategies in tropical trees. Differences in wood density and its radial variation were related to the shade-tolerance of six canopy tree species in seasonally dry tropical forest in Thailand. In addition, using tree ring measurements, the influence of tree size, age and annual increment on radial density gradients was analysed.

Methods

Wood density was determined from tree cores using X-ray densitometry. X-ray films were digitized and images were measured, resulting in a continuous density profile for each sample. Mixed models were then developed to analyse differences in average wood density and in radial gradients in density among the six tree species, as well as the effects of tree age, size and annual increment on radial increases in Melia azedarach.

Key Results

Average wood density generally reflected differences in shade-tolerance, varying by nearly a factor of two. Radial gradients occurred in all species, ranging from an increase of (approx. 70%) in the shade-intolerant Melia azedarach to a decrease of approx. 13% in the shade-tolerant Neolitsea obtusifolia, but the slopes of radial gradients were generally unrelated to shade-tolerance. For Melia azedarach, radial increases were most-parsimoniously explained by log-transformed tree age and annual increment rather than by tree size.

Conclusions

The results indicate that average wood density generally reflects differences in shade-tolerance in seasonally dry tropical forests; however, inferences based on wood density alone are potentially misleading for species with complex life histories. In addition, the findings suggest that a ‘whole-tree’ view of life history and biomechanics is important for understanding patterns of radial variation in wood density. Finally, accounting for wood density gradients is likely to improve the accuracy of estimates of stem biomass and carbon in tropical trees.Key words: Radial gradients, shade-tolerance, tree biomass estimates, tree rings, tropical trees, wood density     

TESTING MITOCHONDRIAL CAPTURE AND DEEP COALESCENCE IN AMAZONIAN CICHLID FISHES (CICHLIDAE: CICHLA)

Hybridization and introgression have important consequences in evolution, such as increasing the genetic diversity and adaptive potential of a species. One of their most conspicuous footprints is discordance among gene trees or between genes and phenotypes. However, most studies that report introgression fail to disprove the null hypothesis that genetic incongruence may result from stochastic sorting of ancestral allelic polymorphisms. In the case of ancient introgression, these two processes may be especially difficult to distinguish topologically, but they make different predictions about the patterns of coalescence among loci. Here we apply three methods, molecular dating, multispecies coalescent models, and gene tree simulation under coalescence, to compare these two hypotheses that explain the polyphyletic mtDNA of the butterfly peacock bass, Cichla orinocensis. In comparison with a species tree based on 20 unlinked nuclear loci, we determined that mtDNA divergences were too recent to be explained by ancestral polymorphism. Similarly, coalescent species tree branches were significantly shorter when putative introgressed mtDNA was incorporated, and simulations showed the mtDNA topology to be unlikely under lineage sorting only. We conclude that introgression approximately 1.5 million years ago resulted in capture by C. orinocensis of an mtDNA lineage ancestral to the modern subspecies C. oc. monoculus.     

CALIBRATING DIVERGENCE TIMES ON SPECIES TREES VERSUS GENE TREES: IMPLICATIONS FOR SPECIATION HISTORY OF APHELOCOMA JAYS

Estimates of the timing of divergence are central to testing the underlying causes of speciation. Relaxed molecular clocks and fossil calibration have improved these estimates; however, these advances are implemented in the context of gene trees, which can overestimate divergence times. Here we couple recent innovations for dating speciation events with the analytical power of species trees, where multilocus data are considered in a coalescent context. Divergence times are estimated in the bird genus Aphelocoma to test whether speciation in these jays coincided with mountain uplift or glacial cycles. Gene trees and species trees show general agreement that diversification began in the Miocene amid mountain uplift. However, dates from the multilocus species tree are more recent, occurring predominately in the Pleistocene, consistent with theory that divergence times can be significantly overestimated with gene‐tree based approaches that do not correct for genetic divergence that predates speciation. In addition to coalescent stochasticity, Haldane's rule could account for some differences in timing estimates between mitochondrial DNA and nuclear genes. By incorporating a fossil calibration applied to the species tree, in addition to the process of gene lineage coalescence, the present approach provides a more biologically realistic framework for dating speciation events, and hence for testing the links between diversification and specific biogeographic and geologic events.     

On the influence of tree size on the climate–growth relationship of New Zealand kauri (Agathis australis): insights from annual, monthly and daily growth patterns

Many tree-ring-based climate reconstructions are based on the assumption that the climate reaction of trees is independent of their size. Here, we test this assumption for New Zealand kauri (Agathis australis), one of the longest tree ring-based proxies for the El Niño-Southern Oscillation (ENSO). The most recent kauri chronology contains a large amount of archaeological material, e.g. timber for which the original tree size is often unknown. We analyzed the climate–growth relationship of different-sized kauri in a pristine forest using different temporal scales, i.e. annually, monthly and daily data on tree growth and climate conditions. Trees of different life stages exhibited approximately the same seasonal growth peaks during austral spring (October and November). The dormancy period overlaps with the period where weekly air temperature maxima are below ca. 17–18 °C, and where the corresponding daily minima are below ca. 8 °C. However, both correlation functions between annual growth and seasonal climate as well as Kalman filter regressions between daily growth and climate conditions suggest an influence of tree size on the climate–growth relationship for kauri. Smaller trees (DBH < 40 cm) contain weaker climate signals than larger trees. Therefore, the precautionary stripping of near-pith material (first 20 cm) from kauri chronologies may result in more uniform responses to climate forcing and thus enhance the reliability of long-term climate reconstructions.     

Anomalous ring identification in two Australian subtropical Araucariaceae species permits annual ring dating and growth-climate relationship development

Almost all Australian tropical and subtropical regions lack annually-resolved long-term (multi-decadal to centennial scale) instrumental climate records. Reconstructing climate in these regions requires the use of sparse climate proxy records such as tree rings. Tree rings often archive annually-resolved centennial-scale climate information. However, many tropical and subtropical species have short life-spans, the timbers are poorly preserved, and there is a belief that the proxy records of these species are often compromised by ring anomalies. Additionally, for many species the relationship between climate (e.g. temperature and/or rainfall) and tree growth has not been established. These factors have led to tree-ring data being underutilized in the Australian subtropics. Trees in the Araucariaceae family, a common family in northern and eastern Australia, are both longer lived than many species in the Australian subtropics, present growth rings that are annual in nature, and their growth is known to vary with climate. In this study we examine two subtropical Araucariaceae species, Araucaria cunninghamii and Araucaria bidwillii, and quantify the relationship between their radial growth and climate variability. Ring anomalies including false, faint, locally absent, and pinching rings, are found to be present in these species, however, bomb-pulse radiocarbon dating of A. cunninghamii samples together with a whole tree approach helped to identify annual growth patterns despite such anomalous ring boundaries. Additionally, to determine which climate variables most influence growth in these species, dendrometers were installed at two locations in subtropical Southeast Queensland, Australia. We found that rainfall variability drives annual ring growth, while temperature constrains the onset and conclusion of the growth season each year. Our results demonstrate that through the use of A. cunninghamii and A. bidwillii trees which demonstrate annual growth in relation to climate variables there is potential to develop centennial scale climate reconstructions from the Australian subtropics. We provide recommendations on how to best identify ring anomalies in these species to help in the future development of long-term chronologies and climate reconstructions.     

On the optimality of the neighbor-joining algorithm

The popular neighbor-joining (NJ) algorithm used in phylogenetics is a greedy algorithm for finding the balanced minimum evolution (BME) tree associated to a dissimilarity map. From this point of view, NJ is "optimal" when the algorithm outputs the tree which minimizes the balanced minimum evolution criterion. We use the fact that the NJ tree topology and the BME tree topology are determined by polyhedral subdivisions of the spaces of dissimilarity maps to study the optimality of the neighbor-joining algorithm. In particular, we investigate and compare the polyhedral subdivisions for n ≤ 8. This requires the measurement of volumes of spherical polytopes in high dimension, which we obtain using a combination of Monte Carlo methods and polyhedral algorithms. Our results include a demonstration that highly unrelated trees can be co-optimal in BME reconstruction, and that NJ regions are not convex. We obtain the l ₂ radius for neighbor-joining for n = 5 and we conjecture that the ability of the neighbor-joining algorithm to recover the BME tree depends on the diameter of the BME tree.     

Sancho,the oldest known Iberian shrub

Cliffs are refuges for old trees and shrubs. In the Mediterranean Basin most dendroclimatic reconstructions have focused on high-elevation forests where tree radial growth is constrained by low temperatures in addition to drought stress. Old shrubs may provide longer ring-width series of hydroclimate proxies in low-elevation, drought-prone Mediterranean ecosystems where old trees are rare. To fill this research gap we investigated the maximum age and climate sensitivity of young, old, and recently dead Phoenician junipers (Juniperus phoenicea L.), growing on calcareous cliffs and nearby plains, in the Guara Natural Park (northeast Spain). The oldest living juniper was ¹⁴C-dated to be 927 years old, and it was named “Sancho” after Don Quixote’s squire. Based on ring counts, the maximum age was 655 years. The difference in age estimates between the ¹⁴C-dates and ring counts was 39 years indicating that ring counts underestimate age. This was due to missing and wedging rings making the cross-dating of old junipers unfeasible. Cool and wet conditions from May to July enhanced radial growth of young junipers. Old shrubs have a high dendroecological potential in Mediterranean sites where their growth is constrained by warm-dry conditions during the growing-season. Further techniques combining dendrochronological and wiggle-match ¹⁴C dating may allow reconstructing long-term hydroclimate in low-elevation Mediterranean areas.     

QUANTIFYING THE PLEISTOCENE HISTORY OF THE OAK GALL PARASITOID CECIDOSTIBA FUNGOSA USING TWENTY INTRON LOCI

The longitudinal spread of temperate organisms into refugial populations in Southern Europe is generally assumed to predate the last interglacial. However, few studies have attempted to quantify this process in nonmodel organisms using explicit models and multilocus data. We used sequence data for 20 intron‐spanning loci (12 kb per individual) to resolve the history of refugial populations of a widespread western Palaearctic oak gall parasitoid Cecidostiba fungosa (Pteromalidae). Using maximum likelihood and Bayesian methods we assess alternative population tree topologies and estimate divergence times and ancestral population sizes under a model of divergence between three refugia (Middle East, Balkans and Iberia). Both methods support an “Out of the East” history for C. fungosa, matching the pattern previously inferred for their gallwasp hosts. However, coalescent‐based estimates of the ages of population divides are much more recent (coinciding with the Eemian interglacial) than nodal ages of single gene trees for C. fungosa and other species. We also find that increasing the sample size from one haploid sequence per refugial population to three only marginally improves parameter estimates. Our results suggest that there is significant information in the minimal samples currently analyzable with maximum likelihood methods, and that similar methods could be applied to multiple species to test alternative models of assemblage evolution.     

Examining the Role of Effective Population Size on Mitochondrial and Multilocus Divergence Time Discordance in a Songbird

Estimates of speciation times are subject to a number of potential errors. One source of bias is that effective population size (N_e) has been shown to influence substitution rates. This issue is of particular interest for phylogeographic studies because population sizes can vary dramatically among genetically structured populations across species’ ranges. In this study, we used multilocus data to examine temporal phylogeographic patterns in a widespread North American songbird, the Northern Cardinal (Cardinalis cardinalis). Species tree estimation indicated that the phylogeographic structure of C. cardinalis was comprised of four well-supported mainland lineages with large population sizes (large N_e) and two island lineages comprised of much smaller populations (small N_e). We inferred speciation times from mtDNA and multilocus data and found there was discordance between events that represented island-mainland divergences, whereas both estimates were similar for divergences among mainland lineages. We performed coalescent simulations and found that the difference in speciation times could be attributed to stochasticity for a recently diverged island lineage. However, the magnitude of the change between speciation times estimated from mtDNA and multilocus data of an older island lineage was substantially greater than predicted by coalescent simulations. For this divergence, we found the discordance in time estimates was due to a substantial increase in the mtDNA substitution rate in the small island population. These findings indicate that in phylogeographic studies the relative tempo of evolution between mtDNA and nuclear DNA can become highly discordant in small populations.     

Fine-Scale Phylogenetic Discordance across the House Mouse Genome

Population genetic theory predicts discordance in the true phylogeny of different genomic regions when studying recently diverged species. Despite this expectation, genome-wide discordance in young species groups has rarely been statistically quantified. The house mouse subspecies group provides a model system for examining phylogenetic discordance. House mouse subspecies are recently derived, suggesting that even if there has been a simple tree-like population history, gene trees could disagree with the population history due to incomplete lineage sorting. Subspecies of house mice also hybridize in nature, raising the possibility that recent introgression might lead to additional phylogenetic discordance. Single-locus approaches have revealed support for conflicting topologies, resulting in a subspecies tree often summarized as a polytomy. To analyze phylogenetic histories on a genomic scale, we applied a recently developed method, Bayesian concordance analysis, to dense SNP data from three closely related subspecies of house mice: Mus musculus musculus, M. m. castaneus, and M. m. domesticus. We documented substantial variation in phylogenetic history across the genome. Although each of the three possible topologies was strongly supported by a large number of loci, there was statistical evidence for a primary phylogenetic history in which M. m. musculus and M. m. castaneus are sister subspecies. These results underscore the importance of measuring phylogenetic discordance in other recently diverged groups using methods such as Bayesian concordance analysis, which are designed for this purpose.