Orthogonal time frequency space (OTFS) can provide better error performance than orthogonal frequency division modulation in the high-speed scenario. However, the two-dimensional convolution between the information-bearing symbols and the channel response in the delay-Doppler domain induces high-complexity detection, which hinders the implementation of the OTFS. In this paper, expectation propagation (EP) is introduced in the OTFS system to improve the reliability compared to the message passing approaches. The low-complexity EP detection with a log-linear order is proposed by observing the sparsity and the block quasi-banded matrix structure of the time-domain matrix. The exploitation of the lower-upper factorization for the banded matrix, the forward or the back substitution algorithm for the lower or the upper triangular matrices reduces the complexity of the large-scale matrix inversion involved in the EP detection from a cubic order to a linear order. The error performance of the proposed scheme is further analyzed by utilizing the state evolution and discussed under the different frame sizes. Moreover, the proposed low-complexity detector is then explored together with the likelihood decoder in an iterative manner to further improve the reliability of the proposed receiver. Finally, the better error performance and the lower computational complexity of the proposed receiver are validated by the numerical results.