Nanoparticles can be used to tune the properties of polyelectrolyte brushes, and polyelectrolyte brushes can be used to control the interaction between nanoparticles and substrates. In the present paper, we investigate the polyelectrolyte brushes immersed in a nanoparticle solution within the analytical strong-stretching theoretical framework. The theoretical model does not take the excluded volume interaction between any two components into account. When there is no nanoparticle loaded, the polyelectrolyte brush is assumed to be an osmotic brush. Local electroneutral approximation is assumed to be still valid after the nanoparticles have been loaded. The loaded nanoparticles are not big enough to deform the grafted polyelectrolyte chains laterally. Analytical formulae for density profiles of each component and brush thickness are derived. The loaded nanoparticles always compress the polyelectrolyte brush. By analyzing the limiting case, a scaling-type diagram for behaviors of the nanoparticle-loading polyelectrolyte brush is constructed. Two characteristic nanoparticle controlling regimes are shown. When the charge of the nanoparticle, Z, is not very large, charged nanoparticles penetrate into the brush and the brush thickness is scaled by $H \sim (Z\varPhi)^{-1/3}$ , where $\varPhi$ is the nanoparticle volume fraction. When the nanoparticle charge Zis large enough, nanoparticles are mainly distributed outside the brush and the brush thickness is scaled by $H \sim (Z\varPhi)^{-1}$ . In the former case, the Coulombic repulsion between the grafted polyelectrolyte chains is screened by the counterions and the nanoparticles, and the brush behavior is determined by the balance between the chain elasticity and the osmotic pressure of the counterions and the nanoparticles. In the latter case, the electrostatic screening is executed by the counterions, and the chain elasticity is balanced by the osmotic pressure of the counterions. The two regimes are divided into subregimes which are dominated respectively by electrostatic or non-electrostatic interaction. The effects of size polydispersity of the nanoparticles are also investigated. It is found that the behaviors of the grafted polyelectrolyte chains are mainly determined by the ratio between the first two moments of the nanoparticle size distribution function. The polyelectrolyte brush is compressed more by the polydispere nanoparticles than by the monodisperse ones. Possible improvement in the present theory is discussed in the conclusion section.