Biochemical oxygen demand (BOD) is one of the most important parameters used for water quality assessment. Alternative methods are essential for accurately prediction of this parameter because the traditional method in predicting the BOD is time-consuming and it is inaccurate due to inconstancies in microbial multiplicity. In this study, the applicability of four hybrid neuro-fuzzy (ANFIS) methods, ANFIS with genetic algorithm (GA), ANFIS with particle swarm optimization (PSO), ANFIS with sine cosine algorithm (SCA), and ANFIS with marine predators algorithm (MPA), was investigated in predicting BOD using distinct input combinations such as potential of hydrogen (pH), dissolved oxygen (DO), electrical conductivity (EC), water temperature (WT), suspended solids (SS), chemical oxygen demand (COD), total nitrogen (TN), and total phosphorus (T-P) acquired from two river stations, Gongreung and Gyeongan, South Korea. The applicability of multi-variate adaptive regression spline (MARS) in determination of the best input combination was examined. The ANFIS-MPA was found to be the best model with the lowest root mean square error and mean absolute error and the highest determination coefficient. It improved the root mean square error of ANFIS-PSO, ANFIS-GA, and ANFIS-SCA models by 13.8%, 12.1%, and 6.3% for Gongreung Station and by 33%, 25%, and 6.3% for Gyeongan Station in the test stage, respectively.