The effect of bootstrap current on the beta limit of MHD equilibria is studied systematically byan iterative calculation of MHD equilibrium and the consistent bootstrap current in high betaheliotron plasmas. The LHD machine is treated as a standardconfiguration heliotron with an L = 2 planar axis. The effects of vacuum magneticconfigurations, pressure profiles and the vertical field control method are studied. Theequilibrium beta limit with consistent bootstrap current is quite sensitive to the magnetic axislocation for finite beta, compared with the currentless cases. For a vacuum configuration withthe magnetic axis shifted inwards in the torus, even in the high beta regimes, the bootstrap currentflows to increase the rotational transform, leading to an increase in the equilibrium beta limit. Onthe contrary, for a vacuum configuration with the magnetic axis shifted outwards in the torus, even inthe low beta regimes, the bootstrap current flows so as to reduce the rotational transform;therefore, there is an acceleration of the Shafranov shift increase as beta increases, leading to adecrease in the equilibrium beta limit. The pressure profiles and vertical field controlmethods influence the equilibrium beta limit through the location of the magnetic axis for finite beta.These characteristics are independent of both device parameters, such as magnetic fieldstrength, and device size in the low collisional regime.