We construct explicit cohomogeneity 2 metrics of G2 holonomy, which are foliated by twistor spaces. The twistor spaces are S2 bundles over four-dimensional Bianchi IX Einstein metrics with self-dual (or anti-self-dual) Weyl tensor. Generically the 4-metric is of triaxial Bianchi IX type, with SU(2) isometry. We derive the first-order differential equations for the metric coefficients, and obtain the corresponding superpotential governing the equations of motion, in the general triaxial Bianchi IX case. In general our metrics have singularities, which are of orbifold or cosmic-string type. For the special case of biaxial Bianchi IX metrics, we give a complete analysis of their local and global properties, and the singularities. In the triaxi...