We study the effect of deformation on the two-neutrino double beta decay for ground state to ground state transitions in all the nuclei whose half-lives have been measured. Our theoretical framework is a deformed QRPA based in Woods-Saxon or Hartree-Fock mean fields. We are able to reproduce at the same time the main characteristics of the two single beta branches, as well as the double beta matrix elements. We find a suppression of the double beta matrix element with respect to the spherical case when the parent and daughter nuclei have different deformations.It has been observed that the matrix elements M 2ν in all the measured 2ν ββ cases for the 0 + → 0 + transition are quenched with respect to those predicted from pure quasiparticle transitions. The physical mechanisms responsible for this reduction have been a subject of great interest. In this work, we study the effect of the parent and daughter deformations on the M 2ν matrix elements.We use a deformaed quasiparticle random phase approximation (pnQRPA) to describe simultaneously the energy distributions of the single β Gamow-Teller (GT) strength and the matrix element M 2ν of all the double beta emitters whose half-lives have been measured: 48 Ca, 76 Ge, 82 Se, 96 Zr, 100 Mo, 116 Cd, 128 Te, 130 Te, 136 Xe, and 150 Nd. The formalism includes a deformed quasiparticle basis and residual spin-isospin separable interactions in both particle-hole (ph) and particle-particle (pp) channels. We consider two different deformed mean fields: a phenomenological Woods-Saxon (WS) potential and a selfconsistent Hartree-Fock (HF) with Skyrme interactions [1]. While the deformation in the HF calculation is obtained selfconsistently, in the WS case it is an input parameter taken to reproduce the experimental quadrupole moments.We describe the 2ν ββ process as two successive GT transitions via intermediate 1 + states,1