In this work, we study the neutrino mixing sum rules arising from discrete symmetries and the class of Littlest Seesaw (LS) neutrino models. These symmetry-based approaches all offer predictions for the cosine of the leptonic CP phase cosδ in terms of the mixing angles, θ13, θ12, θ23, while the LS models also predict the sine of the leptonic CP phase sinδ, as well as making other predictions. In particular, we study the solar neutrino mixing sum rules, arising from charged lepton corrections to tri-bimaximal (TB), bimaximal (BM), golden ratio (GR) and hexagonal (HEX) neutrino mixing, and the atmospheric neutrino mixing sum rules, arising from preserving one of the columns of these types of mixing—for example, the first or second column of the TB mixing matrix (TM1 or TM2)—and we confront them with an up-to-date global fit of the neutrino oscillation data. We show that some mixing sum rules, such as an atmospheric neutrino mixing sum rule arising from a version of neutrino golden ratio mixing (GRa1), are already excluded at 3σ, and we determine the remaining models allowed by the data. We also consider the more predictive LS models (which obey the TM1 sum rules and offer further predictions) based on constrained sequential dominance CSD(n) with n≈3. We compare for the first time the three cases n=2.5, n=3 and n=1+6≈3.45, which are favored by theoretical models, using a new type of analysis to accurately predict the observables θ12, θ23 and δ. We study all the above approaches, solar and atmospheric mixing sum rules and LS models, together so that they may be compared and to give an up-to-date analysis of the predictions of all of these possibilities, when confronted with the most recent global fits.