  We present Subaru/Suprime-Cam deep V and I imaging of seven fields in the
outer regions of M33. Our aim is to search for stellar structures corresponding
to extended HI clouds found in a recent 21-cm survey of the galaxy. Three
fields probe a large HI complex to the southeastern (SE) side of the galaxy. An
additional three fields cover the northwestern (NW) side of the galaxy along
the HI warp. A final target field was chosen further north, at a projected
distance of approximately 25 kpc, to study part of the large stellar plume
recently discovered around M33. We analyse the stellar population at R > 10 kpc
by means of V, I colour magnitude diagrams reaching the red clump. Evolved
stellar populations are found in all fields out to 120' (~ 30 kpc), while a
diffuse population of young stars (~ 200 Myr) is detected out to a
galactocentric radius of 15 kpc. The mean metallicity in the southern fields
remains approximately constant at [M/H] = -0.7 beyond the edge of the optical
disc, from 40' out to 80'. Along the northern fields probing the outer \hi
disc, we also find a metallicity of [M/H] = -0.7 between 35' and 70' from the
centre, which decreases to [M/H] = -1.0 at larger angular radii out to 120'. In
the northernmost field, outside the disc extent, the stellar population of the
large stellar feature possibly related to a M33-M31 interaction is on average
more metal-poor ([M/H] = -1.3) and older (> 6 Gyr). An exponential disc with a
large scale-length (~ 7 kpc) fits well the average distribution of stars
detected in both the SE and NW regions from a galactocentric distance of 11 kpc
out to 30 kpc. The stellar distribution at large radii is disturbed and,
although there is no clear correlation between the stellar substructures and
the location of the HI clouds, this gives evidence for tidal interaction or
accretion events.
