Metastatic breast cancer is an extremely complex disease with limited treatment options due to the lack of information about the major characteristics of metastatic disease. There is an urgent need, therefore, to understand the changes in cellular complexity and dynamics that occur during metastatic progression. In the current study, we analyzed the cellular and molecular differences between primary tumors and paired lung metastases using a syngeneic p53-null mammary tumor model of basal-like breast cancer. Distinct subpopulations driven by the Wnt- and/or STAT3 signaling pathways were detected in vivo using a lentiviral Wnt- and STAT3 signaling reporter system. A significant increase in the overlapping populations driven by both the Wnt- and STAT3 signaling pathways was observed in the lung metastases as compared to the primary tumors. Furthermore, the overlapping populations showed a higher metastatic potential relative to the other populations and pharmacological inhibition of both signaling pathways was shown to markedly reduce the metastatic lesions in established lung metastases. An analysis of the unique molecular features of the lung metastases revealed a significant association with immune response signatures. Specifically, Foxp3 gene expression was markedly increased and elevated levels of Foxp3 + Treg cells were detected in close proximity to lung metastases. Collectively, these studies illustrate the importance of analyzing intratumoral heterogeneity, changes in population dynamics, and the immune microenvironment during metastatic progression.