Four state-of-the-art metaheuristic algorithms including the genetic algorithm (GA), particle swarm optimization (PSO), differential evolutionary (DE), and ant colony optimization (ACO) are applied to an adaptive neuro-fuzzy inference system (ANFIS) for spatial prediction of landslide susceptibility in Qazvin Province (Iran). To this end, the landslide inventory map, composed of 199 identified landslides, is divided into training and testing landslides with a 70:30 ratio. To create the spatial database, thirteen landslide conditioning factors are considered within the geographic information system (GIS). Notably, the spatial interaction between the landslides and mentioned conditioning factors is analyzed by means of frequency ratio (FR) theory. After the optimization process, it was shown that the DE-based model reaches the best response more quickly than other ensembles. The landslide susceptibility maps were developed, and the accuracy of the models was evaluated by a ranking system, based on the calculated area under the receiving operating characteristic curve (AUROC), mean absolute error, and mean square error (MSE) accuracy indices. According to the results, the GA-ANFIS with a total ranking score (TRS) = 24 presented the most accurate prediction, followed by PSO-ANFIS (TRS = 17), DE-ANFIS (TRS = 13), and ACO-ANFIS (TRS = 6). Due to the excellent results of this research, the developed landslide susceptibility maps can be applied for future planning and decision making of the related area.