Over time, species evolve substantial phenotype differences. Yet, genetic analysis of these traits is limited by reproductive barriers to those phenotypes that distinguish closely related species. Here, we conduct a genome-wide non-complementation screen to identify genes that contribute to a major difference in thermal growth profile between two Saccharomyces species. S. cerevisiae is capable of growing at temperatures exceeding 40C, whereas S. uvarum cannot grow above 33C but outperforms S. cerevisiae at 4C. The screen revealed only a single nuclear-encoded gene with a modest contribution to heat tolerance, but a large effect of the species' mitochondrial DNA (mitotype). Furthermore, we found that, while the S. cerevisiae mitotype confers heat tolerance, the S. uvarum mitotype confers cold tolerance. Recombinant mitotypes indicate multiple genes contribute to thermal divergence. Mitochondrial allele replacements showed that divergence in the coding sequence of COX1 has a moderate effect on both heat and cold tolerance, but it does not explain the entire difference between the two mitochondrial genomes. Our results highlight a polygenic architecture for interspecific phenotypic divergence and point to the mitochondrial genome as an evolutionary hotspot for not only reproductive incompatibilities, but also thermal divergence in yeast.