Nuclei interacting with electrons in dense plasmas acquire electronic bound states, modify continuum states, generate resonances and hopping electron states, and generate short-range ionic order. The mean ionization state (MIS), i.e, the mean charge Z of an "average ion" in such plasmas is a valuable concept: pseudo-potentials, pair-distribution functions, equations of state, transport properties, energy-relaxation rates, opacity, radiative processes, etc., can all be formulated using the MIS of the plasma more concisely than with an "all-electron" description. However, the MIS does not have a unique definition and is used and defined differently in different statistical models of plasmas. Here, using the MIS formulations of several average-atom models based on density functional theory, we compare numerical results for Be, Al, and Cu plasmas for conditions inclusive of incomplete atomic ionization and partial electron degeneracy. By contrasting modern orbital-based models with orbital-free Thomas-Fermi models, we quantify effects of shell structure, continuum resonances, the role of exchange and correlation, and the effects of different choices of the fundamental cell and boundary conditions. Finally, the role of the MIS in plasma applications is illustrated in the context of X-ray Thomson scattering in warm dense matter.