The mesh generation in the standard BEM is still not trivial as one may imagine, especially for high-dimensional moving boundary problems. To overcome this difficulty, the boundary-type RBF collocation methods have been proposed and endured a fast development in the recent decade thanks to being integration-free, spectral convergence, easy-to-use, and inherently truly meshless. First, this chapter introduces the basic concepts of the method of fundamental solutions (MFS). Then a few recent boundary-type RBF collocation schemes are presented to tackle the issue of the fictitious boundary in the MFS, such as boundary knot method (BKM), regularized meshless method, and singular boundary method. Following this, an improved multiple reciprocity method (MRM), the recursive composite MRM (RC-MRM), is introduced to establish a boundary-only discretization of nonhomogeneous problems. Finally, numerical demonstrations show the convergence rate and stability of these boundary-type RBF collocation methods for several benchmark examples.Keywords Meshless Á Integration-free Á Collocation Á Fundamental solutions Á Singularity Á Method of fundamental solutions Á Boundary knot method Á Regularized meshless method Á Singular boundary method Á Boundary particle method During the past two decades we have witnessed a research boom on the boundarytype meshless techniques since the construction of a mesh in the standard BEM is not trivial. Among the typical techniques are the boundary node method, the local boundary integral equation method, the boundary cloud method, the boundary point method, the boundary point interpolation method (BPIM), and the method of fundamental solutions (MFS). The essence of all these techniques, excluding the MFS and BPIM, is basically a combination of the moving least square (MLS) technique with various boundary element schemes, whereas the MFS is a boundary-type RBF collocation scheme.Such MLS-based methods involve singular integration and are mathematically complicated, and their low-order approximations also depress computational efficiency. The numerical integration requires a background mesh, and thus the