I think I now understand why this is sometimes said (and why it is inaccurate).
Our polarization in its full form is written as:
$$\frac{2 \text{$\Omega $}_p (\text{$\gamma $}_{12}-i (\text{$\Delta $}_c-\text{$\Delta $p}))}{4 (\text{$\gamma $}_{13}+i \text{$\Delta $p}) (-i \text{$\gamma $}_{12}-\text{$\Delta $c}+\text{$\Delta $p})-i \text{$\Omega $c}^2}$$
In EIT the control field is typically very strong. But we can consider the situation when it's weak and taylor expand $\Omega_c$:
$$\frac{\text{$\Omega $p}}{2 (\text{$\Delta $p}-i \text{$\gamma $13})}+\frac{i \text{$\Omega $c}^2 \text{$\Omega $p}}{8 (\text{$\Delta $p}-i \text{$\gamma $13})^2 (\text{$\gamma $12}-i \text{$\Delta $c}+i \text{$\Delta $p})}+O\left(\text{$\Omega $c}^3\right)$$
(I've added some simple mathematica code that does this if anyone wants it at the end.)
So now we can see that the polarizability can be broken down into a linear sum of different terms. The first term is:
$$ \frac{\text{$\Omega $p}}{2 (\text{$\Delta $p}-i \text{$\gamma $13})} $$
Here we can see that this term can only represent transitions between $|1\rangle$ and $|3\rangle$ (since no $\Omega_c$ is involved). (the atom is in a Lambda configuration where the probe $E_p$ causes transitions from $|1\rangle$ and $|3\rangle$ and the control causes transitions between $|2\rangle$ and $|3\rangle$. $|1\rangle$ and $|2\rangle$ are ground states and $|3\rangle$ is the excited state.)
So when the control field is off, we can think of the strength of the polarization as being due to the transitions from $|1\rangle$ and $|3\rangle$.
The second term has the form:
$$\frac{i \text{$\Omega $c}^2 \text{$\Omega $p}}{8 (\text{$\Delta $p}-i \text{$\gamma $13})^2 (\text{$\gamma $12}-i \text{$\Delta $c}+i \text{$\Delta $p})}$$
Now if we think of each unit of $\Omega_c$ or $\Omega_p$ as being a single photon that causes the atom to transition to another state, we can think of $\text{$\Omega $}_p\text{$\Omega $}_c^2 = \text{$\Omega $}_p\text{$\Omega $}_c \text{$\Omega $}_c$ as the situation where the atom goes from $|1\rangle$ to $|3\rangle$ via $\Omega_p$ then $|3\rangle$ to $|2\rangle$ via $\Omega_c$ and finally $|2\rangle$ back to $|3\rangle$ via the second $\Omega_c$.
Now intuitively, our atom has two pathways to transition to the state $|3\rangle$, one via $|1\rangle \rightarrow |3\rangle$ and one via $|1\rangle \rightarrow |3\rangle \rightarrow |2\rangle \rightarrow |3\rangle$. Now we have a quantum system that can travel two different pathways, and therefore there can be some quantum interference between the two possible paths, like the double-slit experiment or a quantum Mach-Zhender interferometer.
Now from our expansion we know the probability amplitudes. The first path will have an amplitude $\propto \frac{1}{2 (\text{$\Delta $p}-i \text{$\gamma $13})}$ while the second path will have an amplitude $\propto \frac{i }{8 (\text{$\Delta $p}-i \text{$\gamma $13})^2 (\text{$\gamma $12}-i \text{$\Delta $c}+i \text{$\Delta $p})}$. I admit though that it gets a bit confusing at this point, because we are referring to $\rho_{13}$,the quantum superposition strength of $|1\rangle$ and $|3\rangle$, and how that value changes depending on quantum interference.
In any case, this is what I believe to be what is meant when people referring to EIT being due to quantum interference between different transitions. Obviously, EIT with its strong control field cannot be Taylor expanded like this, so this particular path interference only works in this very weak limit. But, intuitively we can still technically think of interference by taking the infinite limit of the sum of all of the contributing terms.
Mathematica code for anyone who wants it:
Series[(2*(\[Gamma]12 -
I (\[CapitalDelta]c - \[CapitalDelta]p)) \[CapitalOmega]p)/(
4*(\[Gamma]13 +
I \[CapitalDelta]p) (-I \[Gamma]12 - \[CapitalDelta]c + \
\[CapitalDelta]p) - I \[CapitalOmega]c^2), {\[CapitalOmega]c, 0, 2}]