Microsatellites and mitochondrial DNA (mtDNA) have traditionally been used in population genetics because of their variability and presumed neutrality, whereas genes of the major histocompatibility complex (MHC) are increasingly of interest because strong selective pressures shape their standing variation. Despite the potential for MHC genes, microsatellites, and mtDNA sequences to complement one another in deciphering population history and demography, the three are rarely used in tandem. Here we report on MHC, microsatellite, and mtDNA variability in a single large population of the eastern tiger salamander (Ambystoma tigrinum tigrinum). We use the mtDNA mismatch distribution and, on microsatellite data, the imbalance index and bottleneck tests to infer aspects of population history and demography. Haplotype and allelic variation was high at all loci surveyed, and heterozygosity was high at the nuclear loci. We find concordance among neutral molecular markers that suggests our study population originated from post-Pleistocene expansions of multiple, fragmented sources that shared few migrants. Differences in N_e estimates derived from haploid and diploid genetic markers are potentially attributable to secondary contact among source populations that experienced rapid mtDNA divergence and comparatively low levels of nuclear DNA divergence. We find strong evidence of natural selection acting on MHC genes and estimate long-term effective population sizes (N_e) that are very large, making small selection intensities significant evolutionary forces in this population.