Role of vector interaction and axial anomaly in the PNJL modeling of the QCD phase diagram