Well, that's anyway not a very good explanation. First of all, in relativistic quantum theory, a one-particle description doesn't make sense (except for free particles), and thus one should immediately use quantized fields rather than classical fields.
Then, the idea is to get relavistically covariant field equations (for the field operators). Thus the most convenient starting point is a Poincare invariant action with a Lagrange density that depends only on the fields and the first space-time derivatives. For the free field there should be no higher powers of field operators and their derivatives than 2nd. Thus, the most simple possibility is
[tex]S[\phi]=\int \mathrm{d}^4 x [(\partial_{\mu} \phi^\dagger)(\partial^{\mu} \phi) - m^2 \phi^{\dagger} \phi ].[/tex]
Hamilton's principle leads to the Klein-Gordan equation,
[tex](\Box+m^2) \phi=(\Box+m^2) \phi^{\dagger}=0.[/tex]
Now, if you want to interpret the field operator to describe, e.g., the electromagnetics of charged scalar bosons, you need to couple this to the electromagnetic field in a gauge invariant way since the electromagnetic field is described by a massless vector field, and if you do not describe it as a gauge field already the representation theory of the Poincare group tells you said all hell breaks loose. So better describe it as a gauge field.
The standard description is to substitute
[tex]\partial_{\mu} \rightarrow D_{\mu}=\partial_{\mu}+\mathrm{i} q A_{\mu}.[/tex]
Then the variation of the corresponding action and setting [tex]A_{\mu}=0[/tex] leads to the conserved electromagnetic current for a free Klein-Gordon field,
[tex]j_{\mu} = \mathrm{i} q (\phi^{\dagger} \partial_{\mu} \phi - (\partial \mu \phi^{\dagger}) \phi).[/tex]
Imho that's the most physically motivated explanation for using this expression as the four-current density of a Klein-Gordon field, and its time component is (the operator representing) the charge density of the corresponding Bose particles.