A neutral structure in the DD[over ¯]^{*} system around the DD[over ¯]^{*} mass threshold is observed with a statistical significance greater than 10σ in the processes e^{+}e^{-}→D^{+}D^{*-}π^{0}+c.c. and e^{+}e^{-}→D^{0}D[over ¯]^{*0}π^{0}+c.c. at sqrt[s]=4.226 and 4.257 GeV in the BESIII experiment. The structure is denoted as Z_{c}(3885)^{0}. Assuming the presence of a resonance, its pole mass and width are determined to be [3885.7_{-5.7}^{+4.3}(stat)±8.4(syst)] MeV/c^{2} and [35_{-12}^{+11}(stat)±15(syst)] MeV, respectively. The Born cross sections are measured to be σ[e^{+}e^{-}→Z_{c}(3885)^{0}π^{0},Z_{c}(3885)^{0}→DD[over ¯]^{*}]=[77±13(stat)±17(syst)] pb at 4.226 GeV and [47±9(stat)±10(syst)] pb at 4.257 GeV. The ratio of decay rates B[Z_{c}(3885)^{0}→D^{+}D^{*-}+c.c.]/B[Z_{c}(3885)^{0}→D^{0}D[over ¯]^{*0}+c.c.] is determined to be 0.96±0.18(stat)±0.12(syst), consistent with no isospin violation in the process, Z_{c}(3885)^{0}→DD[over ¯]^{*}.