Phylogenetic structure of the Thomomys bottae–umbrinus complex in North America

The phylogeography of the Thomomys bottae–umbrinus complex in the United States and Mexico was assessed with sequences of the mitochondrial cytochrome b gene. These sequences were obtained from 225 individuals representing 108 locations over the range, including 56 sequences from GenBank. 110 (500 bp) sequences were used for Bayesian inference and neighbor-joining analyses, and 34 (1140 bp) specimens from the main clades obtained from the Bayesian inference were used in maximum-parsimony and maximum-likelihood analyses. The different analyses indicate significant variation within the species complex that averages 13% among major groups of genetic differences among Thomomys bottae–umbrinus. The overall pattern of geographic variation is not concordant with the current taxonomy. To the contrary, eight monophyletic groups are supported by all analyses and can be considered phylogenetic species. Overall divergence among these groups appears influenced by historical biogeographic events active during the Pliocene and Pleistocene.