If the question is is there a real (deeper) reason and connection between the Ising model (in 2D) and free fermion (in 1D) the answer is yes. I don't know how deep. I believe this was shown for the first time neatly in H. Lieb's solution of the 2D classical Ising model.
The steps are the following
1) Write the partition function of the classical 2D Ising model (on a, say, $N\times N$ square lattice) as
$$
Z := \sum_{\sigma} e^{-\beta E[\sigma]} = \mathrm{Tr} (T^N),
$$
where $\sigma$ is a classical configuration and $E[\sigma]$ is its energy. The above equation defines the transfer matrix $T$.
The matrix $T$ turns out to be a product of Pauli operators defined on a line (i.e., one dimension less than 2D). In particular one obtains something like $T=T_1 T_2$
$$
T_1 \ \sim \ \exp{\left ( K_{\parallel} \sum_i \sigma_i^z \sigma_{i+1}^z \right ) }\\
T_2 \ \sim \ \exp{\left ( K^{\ast} \sum_i \sigma_i^x \right ) }
$$
where $K_{||}, K^{\ast}$ are functions of the original parameters of the 2D Ising model ($J_{\parallel}, J_\perp$).
3) There is a mapping (named after Jordan-Wigner) that maps a system of $N$ spins to a system of $N$ fermions. This is not too surprising giving that both Hilbert spaces are of dimension $2^N$. So the next step is to write out $Z$ in terms of fermion operators $a_i, a_i^\dagger$.
4) It turns out that the resulting expression is integrable. What this mean is somehow involved but we can say it amounts to the fact that, essentially, it can be written as
$$
Z = \mathrm{Tr} e^{\Gamma[A_1]} \dots e^{\Gamma[A_n]}
$$
where $\Gamma[A] :=\sum_{i,j} a^\dagger_i A_{i,j} a_j$.
5) Integrability is now expressed by the following identity
$$
\mathrm{Tr} \left( e^{\Gamma[A_1]} \dots e^{\Gamma[A_n]} \right) = \det \left(1 +e^{A_1} \dots e^{A_1} \right)
$$
namely a trace over a $2^N$ dimensional space is reduced to a trace over an $N$ dimensional one (the one-particle space).
It should be pointed out that not every expression of $N$ Pauli operators is mapped into free fermions. This lucky coincidence is what is responsible for the integrability of problem.
EDIT
I just checked and it turns out that the original paper is by T. D. SCHULTZ, D. C. MATTIS, and E. H. LIEB and can be found here. Not surprisingly its title is Two-Dimensional Ising Model as a Soluble Problem of Many Fermions.