Abstract: We present a lattice QCD calculation of the up, down, strange and charm quark masses performed using the gauge configurations produced by the European Twisted Mass Collaboration with N-f = 2 + 1 + 1 dynamical quarks, which include in the sea, besides two light mass degenerate quarks, also the strange and charm quarks with masses close to their physical values. The simulations are based on a unitary setup for the two light quarks and on a mixed action approach for the strange and charm quarks. The analysis uses data at three values of the lattice spacing and pion masses in the range 210-450 MeV, allowing for accurate continuum limit and controlled chiral extrapolation. The quark mass renormalization is carried out non-perturbatively using the RI'-MOM method. The results for the quark masses converted to the (MS) over bar scheme are: m(ud) (2 GeV) = 3.70(17) MeV, m(s)(2 GeV) = 99.6(4.3) MeV and m(c)(m(c)) = 1.348(46) GeV. We obtain also the quark mass ratios m(s)/m(ud) = 26.66(32) and m(c)/m(s) = 11.62(16). By studying the mass splitting between the neutral and charged kaons and using available lattice results for the electromagnetic contributions, we evaluate m(u)/m(d) = 0.470(56), leading to m(u) = 2.36(24) MeV and m(d) = 5.03(26) MeV.