The integrality gap is a useful indicator of how well an IP can be approximated. It might be better to think of it in an informal, intuitive way. A high integrality gap implies that certain methods won't work. Certain primal/dual methods, for example, depend on a small integrality gap. For the standard primal Vertex Cover LP, the dual LP asks for a maximum matching. In this case, we can do the following:
- find an optimum fractional solution $\boldsymbol{y}$ to the dual LP (a maximum fractional matching)
- multiply the solution $\boldsymbol{y}$ by a factor of 2 (double all edge weights)
- convert this to a feasible integral $\boldsymbol{x}$ for the primal LP (each edge gives half of its weight from the $2\boldsymbol{y}$ vector to each of its endpoints in the $\boldsymbol{x}$ vector, then each $x_i$ is replaced with $\min(\lfloor x_i\rfloor, 1)$).
In this case this simple strategy works and we end up with a feasible integral solution to the primal LP whose weight is no more than twice the weight of a feasible solution for the dual LP. Since the weight of a feasible solution for the dual LP is a lower bound for OPT, this is a 2-approximation algorithm.
Now, where does the integrality gap come in? The IG is 2 in this case, but that alone doesn't imply that the algorithm will work. Rather, it suggests that it might work. And if the IG was more than 2, it would guarantee that the simple strategy would not always work. At the very least we would have to multiply the dual solution by the IG. So the integrality gap sometimes tells us what won't work. The integrality gap can also indicate what kind of approximation factor we can hope for. A small integrality gap suggests that investigating rounding strategies, etc., might be a worthwhile approach.
For a more interesting example, consider the Hitting Set problem and the powerful technique of approximating the problem using $\varepsilon$-nets (Brönnimann & Goodrich, 1995). Many problems can be formulated as instances of Hitting Set, and a strategy that has been successful for many problems is to do this, then just find a good net finder, i.e., an algorithm to construct small $\varepsilon$-nets, and crank everything through the B&G meta-algorithm. So people (myself included) try to find net finders for restricted instances of Hitting Set that, for any $\varepsilon$, can build an $\varepsilon$-net of size $f(1/\varepsilon)$, where the function $f$ should be as small as possible. Having $f(1/\varepsilon) = \mathcal{O}(1/\varepsilon)$ is a typical goal; this would give a $\mathcal{O}(1)$-approximation.
As it turns out, the best possible function $f$ is bounded by the integrality gap of a certain LP for Hitting Set (Even, Rawitz, Shahar, 2005). Specifically, the optimum integral and fractional solutions satisfy $\mathrm{OPT}_I \leq f(\mathrm{OPT}_f)$. For unrestricted instances of Hitting Set the integrality gap is $\Theta(\log(m))$, but when formulating another problem as Hitting Set, the IG can be lower. In this example the authors show how to find $\varepsilon$-nets of size $\mathcal{O}((1/\varepsilon) \log \log (1/\varepsilon))$ for the restricted instances of Hitting Set that correspond to the problem of hitting axis-parallel boxes. In this way they improve upon the best known approximation factor for that problem. It's an open problem whether or not this can be improved. If, for these restricted Hitting Set instances, the IG for the Hitting Set LP is $\Theta(\log \log m)$, it would be impossible to design net finder guaranteeing $\varepsilon$-nets of size $o((1/\varepsilon) \log \log (1/\varepsilon))$, since doing so would imply the existence of an algorithm that guarantees integral hitting sets of size $o(\mathrm{OPT}_f \log\log \mathrm{OPT}_f)$, but since $\mathrm{OPT}_f\leq m$ this would imply a smaller integrality gap. So if the integrality gap is large, proving it could prevent people from wasting their time looking for good net finders.