Using the fact that the neutrino mixing matrix U = (UeU nu)-U-dagger, where U-e and U-v, result from the diagonalisation of the charged lepton and neutrino mass matrices, we analyse the sum rules which the Dirac phase delta present in U satisfies when U-nu has a form dictated by, or associated with, discrete symmetries and U-e has a "minimal" form (in terms of angles and phases it contains) that can provide the requisite corrections to U-nu, so that reactor, atmospheric and solar neutrino mixing angles theta(13), theta(23) and theta(12) have values compatible with the current data. The following symmetry forms are considered: i) tri-bimaximal (IBM), ii) bimaximal (BM) (or corresponding to the conservation of the lepton charge L' = L-e - L-mu - L-tau (LC)), iii) golden ratio type A (GRA), iv) golden ratio type B (GRB), and v) hexagonal (HG). We investigate the predictions for delta in the cases of TBM, BM (LC), GRA, GRB and HG forms using the exact and the leading order sum rules for cos delta proposed in the literature, taking into account also the uncertainties in the measured values of sin(2) theta(12), sin(2) theta(23) and sin(2) theta(13). This allows us, in particular, to assess the accuracy of the predictions for cos delta based on the leading order sum rules and its dependence on the values of the indicated neutrino mixing parameters when the latter are varied in their respective 3 sigma experimentally allowed ranges. (C) 2015 The Authors. Published by Elsevier B.V.