Although the least squares meshfree collocation formulation (LSMC) is very popular among various meshfree collocation schemes, a precise elucidation on the temporal stability of LSMC for transient analysis is still absent due to its formulation complexity. In this talk, we present a thorough theoretical investigation for the temporal stability of LSMC particularly referring to the linear transient heat conduction problems. As a matter of fact, it is disclosed that the semi-discrete LSMC for transient heat conduction analysis essentially corresponds to a non-symmetric discrete formulation of the underlying generalized eigenvalue problem, which is totally different from the symmetric scenario for steady and free vibration analysis. This non-symmetric feature of the generalized eigenvalue problem for LSMC poses severe difficulty to perform the modal reduction via the classical eigenmode orthogonality. This obstacle is removed herein by directly employing the eigenmode independence to decouple the system of discrete equations for LSMC. The generic decoupled characteristic equation along with the generalized trapezoidal rule (GTR) is then used for the temporal stability assessment, and a criterion is accordingly established to estimate the stable time step. Regarding the temporal stability, it is found that when the eigenvalues are real and positive, the stability requirements of GTR for LSMC are almost identical to those derived from the symmetric formulation, except that the eigenvalues now also rely on the time step. In the case that the eigenvalues are complex, the situation becomes more involved and can be divided into two circumstances. If the real parts of the eigenvalues are all non-negative, the stable time step is bounded from the above by both real and imaginary parts of the eigenvalues. On the other hand, when the real parts of certain eigenvalues are negative, a reversed stability condition can be possibly obtained since the stable time step is bounded from the bottom on this specific occasion. Meanwhile, it is shown that LSMC with GTR is first order temporally accurate except the mid-point rule that is second order accurate. In addition, it is also testified that LSMC could effectively alleviate the occurrence of negative-real-part eigenvalues compared with the direct meshfree collocation formulation, which is beneficial for the temporal stability in transient heat conduction analysis.