Molecular systematics, hybridization, and phylogeography of the Bufo americanus complex in Eastern North America

We reconstruct phylogenetic relationships among a well-studied group of toads and find relationships that differ greatly from the current taxonomic understanding. We use mitochondrial sequences encoding ND1, tRNA^Leu(UUR), and part of 16S to infer relationships among members of the Bufo americanus complex. Focusing on the four taxa that historically have been most problematic due to morphological similarity and hybridization in sympatry, we sample 150 individuals from multiple populations across each species' geographic range. Our evidence conflicts with previous taxonomic hypotheses that were based on ability to hybridize, geographic distribution, and call variation. First, sequences from B. fowleri do not comprise the sister clade to sequences of B. woodhousii; therefore the previous classifications of B. fowleri as sister species to, or eastern subspecies of, B. woodhousii are both called into question. Second, sequences from B. americanus are more closely related to those of B. woodhousii than to those of B. terrestris, indicating that similar advertisement call characteristics evolved independently. Third, sequences of B. fowleri are paraphyletic, with sequences of B. terrestris embedded within. Lastly, sequences from B. fowleri cluster into three distinct mitochondrial clades, with some divergences corresponding to greater than 2 mya. These clades are somewhat geographically structured, suggesting divergence in allopatry during the Pleistocene. These mitochondrial divergences are not accompanied by known phenotypic differences, however, suggesting either evolutionary stasis in morphology and behavior, cryptic phenotypic evolution, or that hybridization in secondary contact has homogenized phenotypic differences that may have arisen in allopatry.