See Section 5.1.2, starting on page 108. See, for instance, Example 3, p. 112. This example reveals that the MOL approach with an RK4 solver gives a competitive performance for smooth kernels.
The method of lines is used for sloving PDEs. Most often it refers to the construction or analysis of numerical methods for partial differential equations that proceeds by first discretizing the spatial derivatives only and leaving the time variable continuous. This leads to a system of ordinary differential equations to which a numerical method for initial value ordinary equations can be applied
There are many references you may refere to:
Causley, M., Christlieb, A., Ong, B., & Van Groningen, L. (2014). Method of lines transpose: An implicit solution to the wave equation. Mathematics of Computation.
Dereli, Y., & Schaback, R. (2013). The meshless kernel-based method of lines for solving the equal width equation. Applied Mathematics and Computation, 219(10), 5224-5232.
Northrop, P. W., Ramachandran, P. A., Schiesser, W. E., & Subramanian, V. R. (2013). A robust false transient method of lines for elliptic partial differential equations. Chemical Engineering Science, 90, 32-39.
Furzeland, R. M., Verwer, J. G., & Zegeling, P. A. (1990). A numerical study of three moving-grid methods for one-dimensional partial differential equations which are based on the method of lines. Journal of Computational Physics, 89(2), 349-388.
Reddy, S. C., & Trefethen, L. N. (1992). Stability of the method of lines. Numerische Mathematik, 62(1), 235-267.
Rektorys, K. (1982). The method of discretization in time and partial differential equations. Equadiff 5, 293-296.