In the presence of diffusion, stability of three-layer Hele-Shaw flows which models enhanced oil recovery processes by polymer flooding is studied for the case of variable viscosity in the middle layer. This leads to the coupling of the momentum equation and the species advection-diffusion equation the hydrodynamic stability study of which is presented in this paper. Linear stability analysis of a potentially unstable three-layer rectilinear Hele-Shaw flow is used to examine the effects of species diffusion on the stability of the flow. Using a weak formulation of the disturbance equations, upper bounds on the growth rate of individual disturbances and on the maximal growth rate over all possible disturbances are found. Analytically, it is shown that a short-wave disturbance if unstable can be stabilized by mild diffusion of species, where as an unstable long-wave disturbance can always be stabilized by strong diffusion of species. Thus, an otherwise unstable three-layer Hele-Shaw flow can be completely stabilized by a large enough diffusion, i.e., by increasing enough the magnitude of the species diffusion coefficient. The magnitude of this diffusion coefficient required to completely stabilize the flow will depend on the magnitude of interfacial viscosity jumps and the viscosity gradient of the basic viscous profile of the middle layer.