Models with range-free frustrated Ising spin interaction and additional Hubbard interaction are treated exactly by means of the discrete time slicing method of Grassmann-field theory. Critical and tricritical points, spin and charge correlations, and the fermion propagator are derived as a function of temperature and chemical potential , of the Hubbard coupling U and the spin-glass energy J. U is allowed to be either repulsive (U Ͼ0) or attractive (UϽ0). Cuts through the multidimensional phase diagram are obtained. Analytical and numerical evaluations take important replica symmetry-breaking ͑RSB͒ effects into account. Results for the ordered phase are given at least in a one-step approximation ͑1RSB͒, and for Tϭ0 we report the first two-, three-, and four-step calculations ͑4RSB͒ for fermionic spin glasses. The use of exact relations and invariances under RSB together with 2RSB calculations for all fillings and 4RSB solutions for half filling allow to model exact solutions by interpolation. For Tϭ0, our numerical results provide strong evidence that the exact spin-glass pseudogap obeys (E)ϭc 1 ͉EϪE F ͉ for energies close to the Fermi level with c 1 Ϸ0.13. Rapid convergence of Ј(E F ) under increasing orders of RSB is observed and Љ(E) is evaluated to estimate subleading powers. Over a wide range of the pseudogap and after a small transient regime (E) regains a linear shape with a larger slope and a small S-like perturbation. The leading term resembles the Efros-Shklovskii Coulomb pseudogap of two-dimensional localized disordered fermionic systems. Beyond half filling we obtain a Ϫ1ϳ(ϪU) 2 , уU dependence of the fermion filling factor . We find a half filling transition between a phase for UϾ, where the Fermi level lies inside the Hubbard gap, into a phase where (ϾU) is located at the center of the upper spin-glass ͑SG͒ pseudogap ͑SG gap͒. For ϾU the Hubbard gap combines with the lower one of two SG gaps ͑phase I͒, while for ϽU it joins the sole SG gap which exists in this half-filling regime ͑phase II͒. Shoulders of the combined gaps are shaped by RSB due to spin-glass order. We predict scaling behavior at the half-filling transition which becomes continuous due to ϱRSB. Implications of the half-filling transition between the deeper insulating phase II and phase I for the eventual delocalization by additional hopping processes in itinerant model extensions are discussed. Possible metalinsulator transition scenarios are described. Generalizations to random Hubbard coupling and alloy models as well as frustrated magnetic interactions with ferro-or antiferromagnetic components are discussed.