Let us consider scalar electrodynamics (Klein-Gordon-Maxwell electrodynamics) with the Lagrangian:
\begin{eqnarray}\label{eq:pr6}
-\frac{1}{4}F^{\mu\nu}F_{\mu\nu}+\frac{1}{2}(\psi^*_{,\mu}-ieA_\mu\psi^*)(\psi^{,\mu}+ieA^\mu\psi)-\frac{1}{2}m^2\psi^*\psi
\end{eqnarray}
and the equations of motion
\begin{equation}\label{eq:pr7}
(\partial^\mu+ieA^\mu)(\partial_\mu+ieA_\mu)\psi+m^2\psi=0,
\end{equation}
\begin{equation}\label{eq:pr8}
\Box A_\mu-A^\nu_{,\nu\mu}=j_\mu,
\end{equation}
\begin{equation}\label{eq:pr9}
j_\mu=ie(\psi^*\psi_{,\mu}-\psi^*_{,\mu}\psi)-2e^2 A_\mu\psi^*\psi.
\end{equation}
The complex charged matter field $\psi$ can be made real by a gauge transform (at least locally), and the equations of motion in the relevant gauge (unitary gauge) for the transformed four-potential of electromagnetic field $B^{\mu}$ and real matter field $\varphi$ are as follows ($E.~Schr\ddot{o}dinger, {Nature}$, 169:538, 1952):
\begin{equation}\label{eq:pr10}
\Box\varphi-(e^2 B^\mu B_\mu-m^2)\varphi=0,
\end{equation}
\begin{equation}\label{eq:pr11}
\Box B_\mu-B^\nu_{,\nu\mu}=j_\mu,
\end{equation}
\begin{equation}\label{eq:pr12}
j_\mu=-2e^2 B_\mu\varphi^2.
\end{equation}
Schrödinger made the following comment: "That the wave function ... can be made real by a change of gauge is but a truism, though it contradicts the widespread belief about 'charged' fields requiring complex representation."
These equations of motion can be obtained from the Lagrangian of work (T. Takabayasi(1953), Progr. Theor. Phys.,9,187):
\begin{equation}\label{eq:pr12a}
-\frac{1}{4}F^{\mu\nu}F_{\mu\nu}+\frac{1}{2}e^2 B_\mu B^\mu \phi^2+\frac{1}{2}(\varphi_{,\mu}\varphi^{,\mu}-m^2\varphi^2).
\end{equation}
Actually, it coincides with the previous Lagrangian up to the replacement of the complex scalar field by a real one. Any solution of the equations of motion for the first Lagrangian has a physically equivalent solution of the equations of motion for the second Lagrangian.
I did not consider quantization here.
EDIT(02/10/2018):
I would like to add that the above system of interacting real scalar field and electromagnetic field has some amazing properties. One can see from the equations of motion that the real scalar field can be algebraically eliminated, and it turns out that the resulting equations for the electromagnetic field describe its independent evolution (my article in European Physical Journal C http://link.springer.com/content/pdf/10.1140%2Fepjc%2Fs10052-013-2371-4.pdf and references there).
Apparently, it is also possible to introduce a Lorentz-invariant Lagrangian with higher derivatives that does not include the matter field, but is largely equivalent to the Takabayasi Lagrangian (my article https://arxiv.org/abs/1006.2578 ; the significance of some special cases, e.g., $\varphi=0$ and $B^\mu B_\mu=0$ (see below) is unclear (different particles?)). To this end, the Takabayasi Lagrangian can be expressed in terms of $\Phi=\varphi^2$, rather than $\varphi$, using, e.g., the following:
\begin{equation}\label{eq:pr16a}
\varphi_{,\mu}\varphi^{,\mu}=\frac{1}{4}\frac{\Phi_{,\mu}\Phi^{,\mu}}{\Phi},
\end{equation}
and then $\Phi$ can be replaced by the following expression obtained from the equations of motion:
\begin{equation}\label{eq:pr16b}
\Phi=-\frac{1}{2e^2} \frac{B^\mu(\Box B_\mu-B^\nu_{,\nu\mu})}{B^\mu B_\mu}.
\end{equation}
Thus, a Lagrangian including only electromagnetic field describes pretty much the same physics as scalar electrodynamics.
