Species distribution models identify regions with ideal environmental characteristics for the establishment and proliferation of species. The chayote is a crop that originated and domesticated in Mexico; however, it is cultivated in different parts of the world due to its nutritional and pharmaceutical importance. The objective of this research was to locate the potential distribution of S. edule in Japan supported on seven machine learning models, to also determine which bioclimatic variables influence its distribution, and which are the most suitable regions for its establishment. Thirty-one occurrence points, elevation, and the bioclimatic variables bio1, bio3, bio4, bio7, bio8, bio12, bio14, bio15, and bio17 were used to infer the models. Hundred percent of the occurrence points coincided with the Cfa climate distributed in Acrisol (60.9%), Andosol (17.4%), Cambisol (13%), Fluvisol (4.35%), and Gleysol (4.35%) soil. The Maxent model reported the highest AUC value (0.93), while the GLM obtained the best TSS value (0.84); the SVM model reported the largest suitability area ≥ 0.5 with 100,394.4 km2. Temperature-related variables were the major contributors to the models and the ones explaining the distribution limits of S. edule in Japan. The coastal eastern prefectures of Kantō, Chūbu, Kinki, Chūgoku, Kyūshū, and Shikoku regions showed a suitability ≥ 0.5.