Neural network adaptive iterative learning control (ILC) is developed in this article to treat strict-feedback nonlinear systems with unknown state delays and input saturation. These delays are treated by constructing the Lyapunov-Krasovskii (L-K) functions for each subsystem. A command filter is employed to avoid the derivative explosion caused by continuous differentiation of the virtual controller. Corresponding auxiliary systems are designed and integrated into the backstepping procedure to compensate input saturation and the unimplemented part of the filter. Hyperbolic tangent functions and radial basis function neural networks (RBF NNs) are employed to treat singularity and related unknown terms, respectively. The convergence of the resultant strict-feedback systems is ensured in the framework of composite energy function (CEF). Finally, a simulation example is adopted to substantiate the validity of the proposed algorithm.