  The candidate PeVatron MGRO J1908+06, which shows a hard spectrum beyond 100
TeV, is one of the most peculiar $\gamma$-ray sources in the Galactic plane.
Its complex morphology and some possible counterparts spatially related with
the VHE emission region, preclude to distinguish between a hadronic or leptonic
nature of the $\gamma$-ray emission. In this paper we illustrate a new
multiwavelength analysis of MGRO J1908+06, with the aim to shed light on its
nature and the origin of its ultra high-energy emission. We performed an
analysis of the $^{12}$CO and $^{13}$CO molecular line emission demonstrating
the presence of dense molecular clouds spatially correlated with the source
region. We also analyzed 12-years of Fermi-LAT data between 10 GeV and 1 TeV
finding a counterpart with a hard spectrum ($\Gamma \sim 1.6$). Our reanalysis
of XMM-Newton data allowed us to put a more stringent constraint on the X-ray
flux from this source. We demonstrate that a single accelerator cannot explain
the whole set of multiwavelength data, regardless of whether it accelerates
protons or electrons, but a 2-zone model is needed to explain the emission from
MGRO J1908+06. The VHE emission seems most likely the superposition of a TeV
PWN powered by PSR J1907+0602, in the southern part, and of the interaction
between the supernova remnant G40.5-0.5 and the molecular clouds towards the
northern region.

High Energy Astrophysical Phenomena